Next Article in Journal
A Review on the Morphology, Cultivation, Identification, Phytochemistry, and Pharmacology of Kitagawia praeruptora (Dunn) Pimenov
Next Article in Special Issue
Synthesis of Trifluoromethylated Pyrimido[1,2-b]indazole Derivatives through the Cyclocondensation of 3-Aminoindazoles with Ketoester and Their Functionalization via Suzuki-Miyaura Cross-Coupling and SNAr Reactions
Previous Article in Journal
Construction of Ketoenamine-Based Covalent Organic Frameworks with Electron-Rich Sites for Efficient and Rapid Removal of Iodine from Solution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Question of Zwitterionic Intermediates in the [3+2] Cycloaddition Reactions between Aryl Azides and Ethyl Propiolate

by
Ewa Dresler
1,
Przemysław Woliński
2,
Aneta Wróblewska
3 and
Radomir Jasiński
2,*
1
Łukasiewicz Research Network—Institute of Heavy Organic Synthesis “Blachownia”, Energetyków 9, 47-225 Kędzierzyn-Koźle, Poland
2
Institute of Organic Chemistry and Technology, Cracow University of Technology, Warszawska 24, 31-155 Cracow, Poland
3
Department of Organic Chemistry, Faculty of Chemistry, University of Lódź, Tamka 12, 91-403 Łódź, Poland
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(24), 8152; https://doi.org/10.3390/molecules28248152
Submission received: 28 November 2023 / Revised: 13 December 2023 / Accepted: 16 December 2023 / Published: 18 December 2023

Abstract

:
The molecular mechanism of the [3+2] cycloaddition reactions between aryl azides and ethyl propiolate was evaluated in the framework of the Molecular Electron Density Theory. It was found that independently of the nature of the substituent within the azide molecule, the cycloaddition process is realized via a polar but single-step mechanism. All attempts of localization as postulated earlier by Abu-Orabi and coworkers’ zwitterionic intermediates were not successful. At the same time, the formation of zwitterions with an “extended” conformation is possible on parallel reaction paths. The ELF analysis shows that the studied cycloaddition reaction leading to the 1,4-triazole proceeds by a two-stage one-step mechanism. It also revealed that both zwitterions are created by the donation of the nitrogen atom’s nonbonding electron densities to carbon atoms of ethyl propiolate.

1. Introduction

The chemistry of heterocyclic compounds is one of the most dynamically developing areas of science [1,2,3,4,5]. Among a number of compounds of this type, structures from the triazole group enjoy undying interest [6,7,8]. In practice, the triazole heterocyclic ring is a key molecular segment of many compounds that exhibit different types of important activities such as anti-tubercular [9,10], antiplasmodial, antimalarial [11], anticancer [12], antifungal [13,14,15], anticonvulsant [16], antiviral, antineoplastic, antihypertensive, antianxiety, and other [17].
The most universal methodology of the preparation of five-membered heterocycles is [3+2] cycloaddition reactions (32CA reactions) with the participation of different types of three-atom components (TACSs) [18,19,20]. The preparation of the 1,2,3-triazole molecular segment is also possible in this way due to 32CA between azides and alkynes.
In general, the 32CA processes between organic azides and alkynes are realized under relatively difficult conditions and with lower regioselectivity [21] than similar processes with the participation of other nitrogen-containing TACs such as analogs of diazomethane [22], nitrile N-oxides [23,24], nitrones [25,26] nitronates [27], or azomethine ylides [28]. For example, the benzyl azide reacts with the 3-phenoxypropyne at 92 °C under solvent-free conditions yielding a mixture of both possible regioisomeric cycloadducts: 3-benzyl-4-phenoxymethyl-1,2,3-triazole and 3-benzyl-5-phenoxymethyl-1,2,3-triazole (Scheme 1) [29].
Relatively easier are analogous 32CAs with the participation of cyclooctyne analogs [30,31]. This is, however, an incidental example, determined by the specific, strained nature of the alkyne. Similar, atypical reactivity is also observed in the case of azide-fullerene C20 cycloadditions (Scheme 2) [32].
Full regioselective synthesis of 1,2,3-triazoles via the 32CA reaction and under milder conditions is possible via the catalytic protocol, including the presence of cooper or rhutenium catalysts. For example, benzyl azide reacts with simple alkynes at 60 °C. The reaction is realized with full regioselectivity (Scheme 3) [33].
Some years ago, Abu-Orabi and coworkers [34] described an interesting case of the 32CA reactions between aryl azides (1ae) and ethyl propiolate (2) (Scheme 4). Independently of the absence of a catalyst, all analyzed reactions were realized under relatively mild conditions and (most importantly) with full regioselectivity. The authors explained this phenomenon by assuming a stepwise mechanism via a zwitterionic intermediate. The presence of the zwitterionic intermediates is, however, not confirmed in any way.
It should be underlined that, at this moment, the actual state of knowledge shed new light on the earlier interpretation of the 32CA reaction mechanism as one-step “concerted” independently of the nature of starting molecules. In particular, at this moment, several different mechanisms can be assumed within the exploration of the 32CA mechanisms; in particular, (i) polar mechanisms (the one-step mechanism, two-stage one-step mechanism, and stepwise zwitterionic mechanism) [35,36]) or (ii) non-polar mechanism (the one-step mechanism, two-stage one-step mechanism asynchronous mechanism, and stepwise biradical mechanism) [37,38]. Which of them is implemented in practice is determined by the interplay of two independent factors—the nature of electronic interactions [39] and the power of steric interactions [40]. Next, Domingo generally undermined the term “pericyclic” regarding all 32CA reactions based on the recent progress of the exploration of reaction profiles via BET and ELF techniques [41,42]. Therefore, the mechanistic aspects of the title processes evidently require deeper exploration.

2. Results and Discussion

The title reaction can be theoretically realized according to two competitive regioisomeric channels leading to the 3-aryl-5-carboethoxy-1,2,3-pyrazoles (3ae) or 3-aryl-4-carboethoxy-1,2,3-pyrazoles (4ae), respectively (Scheme 5). Within both reaction channels, however, two alternative mechanisms are possible: One-step mechanisms (paths A and B) or stepwise zwitterionic mechanisms (paths C and D).
In the initial step, we decided to analyze the nature of the intermolecular interactions between respective pairs of reagents (Table 1, Figure 1 and Figure 2). It was found that ethyl propiolate is characterized by global electrophilicity equal to 1.48 eV. Therefore, according to the unique Domingo electrophilicity scale [43,44], this component should be treated as a moderate electrophile. The analysis of the distribution of local electrophilicities subsequently shows that the most activated reaction center is located on the terminal carbon atom of acetylenic moiety (ω1 = 0.67 eV). In comparison, the local electrophilicity at the second carbon atom of the acetylenic moiety is equal to 0.04 eV. Therefore, in light of the Molecular Electron Density Theory [45], the course of the addition reaction with the participation of ethyl propiolate should be determined by the attack of nucleophilic agents on the terminal position of the ethyne moiety. On the other hand, the considered organic azides are characterized by global nucleophilicities in the range of 2.12–3.45 eV. Relatively stronger nucleophilicity properties are observed in the case of EDG-substituted azides, whereas weaker nucleophilicities are assigned to EWG-substituted azides (Table 1). It is interesting that local nucleophilicities on both important reaction centers within the NNN moiety are not significantly different and exist in the range of 0.5–0.7 eV. This suggests that in light of electronic effects, both regioisomeric cycloaddition paths should be kinetically allowed.
Within the next research step, we explored theoretically possible channels of reactions between aryl azides and ethyl propiolate, starting from the model addition 1c+2. The results of the wb97xd/6-311+G(d)(PCM) calculations show that from a qualitative point of view, enthalpy profiles of reactions leading to 3-phenyl-5-carboethoxy-1,2,3-pyrazoles (A) or 3-phenyl-4-carboethoxy-1,2,3-pyrazoles (B) are very similar. In both cases, between the valley of the individual reagents and the valley of the respective cycloadduct, two critical points were detected. These are connected by the existence of the pre-reaction molecular complex (MC) and transition state (TS), respectively.
Interactions between addend molecules initially lead to the formation of the pre-reaction molecular complex MC (MCA and MCB for paths A and B, respectively). This transformation of the reaction system is accompanied by reductions in enthalpy of 3.2 and 5.0 kcal/mol for MCA and MCB, respectively (Table 2, Figure 3). At the same time, however, a substantial reduction in the entropy of the reaction system is observed. Consequently, Gibbs free energies of the formation of MCs exhibit positive values. This excludes the possibility of the existence of MCs as thermodynamically stable intermediates. In the framework of MCs, substructures of addends adopt specific orientations, which stimulate stabilization via the coulombic interactions (Figure 4). In general, the relative orientation of the reaction centers at this stage determines the further observed regioselectivity. Therefore, localized intermediates can be interpreted as orientation complexes [46]. Within MCs, any new bonds are not formed (Table 3). Key interatomic distances of N3-C4 and C5-N1 exist beyond the area of the range typical for new C-N bonds within transition states [47,48,49]. These types of MCs were experimentally observed earlier regarding 32CAs of ozone with ethene and ethyne [50,51]. It should be underlined that both MCs do not exhibit the nature of the charge-transfer complexes [52,53]. This was confirmed by analysis of global electron density transfer (GEDT [54]) values (GEDT = 0.00e).
The further transformation of MCs on both considered reaction paths leads to TS (TSA and TSB on reaction paths A and B, respectively). Kinetic factors clearly favored the formation of the TSA structure, which is connected to the 3-phenyl-5-carboethoxy-1,2,3-triazole 3c. This adduct was detected experimentally in the post-reaction mixture. It is interesting that the competitive reaction path B should be treated as disfavored, but not forbidden from a kinetic point of view, because the activation barrier in this case is only 1.1 kcal/mol lower than in the case of path A. However, one should wonder whether the authors of the work [34] carefully analyzed the post-reaction mass. The published procedure only mentions the evaporation of the solvent and recrystallization of the post-reaction mass. No detailed composition analyses were performed using, e.g., the HPLC technique [55,56,57]. We believe that the second regioisomer is likely formed there as a minority, with a yield of several percent. This is very likely because the isolated product is produced with a yield of 80–90% and the vast majority of 32CA reactions occur with a total yield close to 95–99%. This hypothesis correlates well with the above-described analysis of the local reactivity of title addends.
Both optimized TSs are characterized by imaginary frequencies (−484.72 and −524.70 cm−1, respectively) and are structurally similar. In particular, within both TSs, key interatomic distances (N3-C4 and C5-N1) are substantially reduced in comparison to the respective MCs (Table 4). However, the kinetically favored TSA transition state is slightly less unsymmetrical. It should be underlined that both considered TSs exhibit evidently polar natures. This was confirmed by an examination of GEDT values (0.31 and 0.37e for TSA and TSB, respectively). The IRC analysis confirms without any doubt that optimized TSs are connected directly to the valleys of the respective MCs and respective products (Figure 3). All attempts at the localization of hypothetical zwitterionic intermediates were not successful. Therefore, postulated reaction channels C and D cannot be realized in practice. At the same time, however, we detected the possibility of the formation of other zwitterionic structures as a consequence of the interactions between azide 1c and alkyne 2 (Scheme 6). These are not postulated intermediates I1C and I2D, but other, different adducts I3E and I4F (Figure 5 and Figure 6) characterized by an “extended” conformation [58,59]. This type of conformation excludes the possibility of the cyclization of zwitterions to cycloadducts. Their conversion to triazoles is realized via dissociation to individual reagents, and, in the next step, via one-step cycloaddition according to the A or B path. Similar dissociations are observed during the rotation around the NNCC bond.
The formation of zwitterions I3E and I4F are realized via transition states TSE and TSF, respectively. Within the considered TSs, only one new single bond is formed—N3-C4 in the case of TSE and N1-C5 in the case of TSF. It should be underlined, however, that the formation of intermediates I3E and I4F is not favored from a thermodynamic point of view, because the Gibbs free energies of the formation of I3E and I4F are substantially positive. Next, from the kinetic point of view, both these transformations should be considered forbidden due to the very high activation barriers in comparison to reaction paths A and B.
In a similar way, we examined similar reactions with the participation of azides substituted by different types of EDG or EWG groups. In all cases, we detected analogous molecular mechanisms as in the case of the addition 1c+2. Therefore, the proposed mechanism can be assumed to be general for a defined group of 32CA reactions.
In the last step, we decided to characterize the mechanism of 32CA of phenylazide 1c with ethyl propiolate 2 as a representative model with the bonding evolution theory (BET) [60] study along the more favorable reaction path. In Table 4 and Scheme 7, detailed BET data of the critical points of the reaction are presented. The choice requirement is the change in the electronic structure, such as the creation or disappearance of a basin compared to the previous point.
The BET study of 716 points along the reaction path revealed twelve phases following intrinsic reaction coordinates. The most important observations can be summed up as follows:
(1)
In phases III, disynaptic basins V(N1,N2) and V′(N1,N2) go through topological changes leading to the creation of a new monosynaptic basin V(N2). First, in P1, disynaptic basin V′(N1,N2) disappears, transferring its electron density into V(N1,N2) and integrating 4.20e. Next, in P2, monosynaptic basin V(N2) representing nonbonding electron density is created with a population of 0.62e originating from V(N1,N2). The energy increases by 9.4 kcal/mol.
(2)
In phases IIIVI, events leading to the creation of the first new single bond can be observed. During phase III, the integration of disynaptic basins V(N1,N2) and V(N1,N2) decreases by 0.81e and 0.42e, respectively, increasing the population of monosynaptic basin V(N2). In P3, a pseudoradical center is created on C5 represented by the V(C5) monosynaptic basin, integrating 0.14e. An increase of 9.9 kcal/mol in energy can be observed, and the GEDT is −0.08e. Next, a short-phase V starts with the creation of a pseudoradical center on C4 integrating 0.05e; at the same time, an increase in the population in basins V(C5) and V(N2) can be seen. In P5, monosynaptic basin V(C4) disappears, with its electron density being transferred to the monosynaptic basin V(N3) and integrating 3.48e.
(3)
Phase VII, d(N1-C5) = 2.104 Å and d(N3-C4) = 1.878 Å, starts with the creation of the first new C-N bond by donation of the nonbonding electron density of N3 to C4. The disynaptic basin V(N3,C4) is created with an initial population of 1.45e, while the V(N3) monosynaptic basin’s integration decreases to 2.06e, as shown by structures P6′ and P6 in Figure 6.
(4)
In P7, a new monosynaptic basin V′(N1) integrating 0.38e is created with its population originating from V(N1) now integrating 3.45e. Meanwhile, a transfer of population from disynaptic basins V(C4,C5) and V′(C4,C5) of the C4-C5 double bond to V(N3,C4) and V(C5) can be observed. The GEDT is 0.14e, and the energy decreases by 15.3 kcal/mol.
(5)
Phase IX starts with the creation of a monosynaptic basin V′(N3), integrating 0.50e, and at the same time the population of the disynaptic basin V(N3,C4) decreases to 1.33e. Next, at P9, basin V′(N3) disappears transferring its population back to V(N3,C4), integrating now 2.02e. Additionally, an increase in population of the monosynaptic basin V′(N1) of 0.30e originating from V(N1) can be seen. The energy decreases by 20.2 kcal/mol and the GEDT increases to 0.24e.
(6)
The second new C-N bond is created at P10, d(N1-C5) = 1.720 Å and d(N3-C4) = 1.561 Å, by shearing of the pseudoradical center C5 and nonbonding electron density of N1. The monosynaptic basin V(C5) and V′(N1), both integrating 0.75e, merge to create the new disynaptic basin V(N1,C5) with a population of 1.50e, as shown by structures P10′ and P10 in Figure 6.
(7)
Phase XII starts with the disappearance of V′(C4,C5) with the transfer of its population to V(C4,C5), now integrating 3.67e. A decrease in energy of 27.3 kcal/mol and a GEDT of 0.29e can be observed.
(8)
The 32CA of phenylazide 1c with alkyne 2 takes place via a two-stage one-step mechanism of the attack of nucleophilic N3 of the azide on the most electrophilic carbon C4 of the alkyne.
In the same way, the process of N3-C4 zwitterion creation has been evaluated (Table 5, Scheme 8). Analysis of hte 264 points along the IRC revealed four phases, of which the most important aspects are:
(1)
In P1, a monosynaptic basin V(C5) integrating 0.24e is created with its population originating from two disynaptic basins V(C4,C5) and V′(C4,C5), as shown by P1′ and P1 in Figure 7. The energy increases by 17.0 kcal/mol and the GEDT is −0.03e.
(2)
At P2, the most significant topological change can be observed, and a new bond is created by sharing the N3 nonbonding electron density with C4. The monosynaptic basin V(N3) is depopulated and a new disynaptic basin V(N3,C4) is created, integrating 1.68e, as shown by structures P2′ and P2 in Figure 7. The GEDT increases to −0.31e, and another significant rise in energy of 17.7 kcal/mol can be seen.
(3)
The last phase IV starts with the creation of monosynaptic basin V′(N3) integrating 0.69e. The energy decreases by 0.9 kcal/mol, while the GEDT increases to −0.51e. The energy of the final structure of the zwitterion I3E is 0.5 kcal/mol lower and the GEDT reaches −0.53e.
Scheme 8. Simplified representation of the molecular mechanism of the reaction of phenyl azide 1c and ethyl propiolate 2 by Lewis-like structures based on the topological analysis of ELF along the reaction path leading to I3E.
Scheme 8. Simplified representation of the molecular mechanism of the reaction of phenyl azide 1c and ethyl propiolate 2 by Lewis-like structures based on the topological analysis of ELF along the reaction path leading to I3E.
Molecules 28 08152 sch008
Figure 7. ELF valence basins attractors of structures of the most important topological changes, creation of negative charge on C5 and creation of N3-C4 single bond, during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to I3E. Structures with apostrophes represent the first structure before the critical point.
Figure 7. ELF valence basins attractors of structures of the most important topological changes, creation of negative charge on C5 and creation of N3-C4 single bond, during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to I3E. Structures with apostrophes represent the first structure before the critical point.
Molecules 28 08152 g007
Table 5. ELF valence basin populations of the IRC points, MCEI3E, defining the four different phases characterizing the reaction of the phenyl azide 1c and ethyl propiolate 2. The stationary points 1c, 2, MCE, TSE, and I3E are also included. Distances are given in angstroms, Å, electron populations in the average number of electrons, [e], relative energies in kcal·mol−1, and GEDT values in an average number of electrons, [e].
Table 5. ELF valence basin populations of the IRC points, MCEI3E, defining the four different phases characterizing the reaction of the phenyl azide 1c and ethyl propiolate 2. The stationary points 1c, 2, MCE, TSE, and I3E are also included. Distances are given in angstroms, Å, electron populations in the average number of electrons, [e], relative energies in kcal·mol−1, and GEDT values in an average number of electrons, [e].
Structures1c2MCE P1 P2TSEP3 I3E
Phases I II III IV
d1(N3-C4) 3.460 2.474 1.7981.7101.532 1.518
GEDT 0.00 −0.03 −0.31−0.38−0.51 −0.53
dE −6.2 10.8 28.528.827.6 27.1
V(N1)3.76 3.77 3.80 3.673.643.59 3.59
V(N1.N2)2.34 2.32 2.50 2.672.652.37 2.32
V′(N1.N2)1.80 1.81 1.67 1.701.752.02 2.06
V(N2.N3)2.50 2.51 2.47 2.282.252.27 2.28
V(N3)3.37 3.37 3.28 1.701.541.25 1.10
V(N3.C7)1.83 1.84 1.88 1.811.801.80 1.79
V(C4.C5) 2.562.56 2.49 2.062.001.83 1.81
V′(C4.C5) 2.732.72 2.59 2.102.041.88 1.85
V(C5.C6) 2.402.40 2.41 2.312.322.32 2.29
V′(N3) 0.69 0.84
V(N3.C4) 1.681.901.59 1.61
V(C4) 0.24 1.351.471.81 1.90
The analysis of the mechanism of formation of less favorable zwitterion N1-C5 (Table 6, Scheme 9) revealed ten phases during which the most important topological changes are:
(1)
The first topological change takes place in P1 where disynaptic basin V′(C4,C5) merges with V(C4,C5), now integrating 5.19e, and an increase in energy of 27.8 kcal/mol can be observed. Next, phase III starts with the creation of a monosynaptic basin V(C4), with an initial population of 0.26e, which represents the developing negative charge of the zwitterion, as shown in Figure 8 by P2′ and P2. The energy increases by 8.0 kcal/mol and the GEDT is −0.06e.
(2)
With the start of phase IV, another mayor jump in energy can be seen; it rises by 17.4 kcal/mol. The new bond between C5 and N1 is created by shearing the nonbonding electron density of N1 with C5. The disynaptic basin V(N1,C5) is created with the initial population of 1.00e coming from the monosynaptic basin V(N1), as shown in Figure 8 by P3′ and P3. The GEDT increases to −0.26e.
(3)
In phases V-X, further reorganization of the electron density can be observed. The disynaptic basin V(C4,C5) and monosynaptic basins V(N1) and V(N3) gradually depopulate while the disynaptic basins V(N1,C5), V(N2,N3), and V′(N2,N3) and monosynaptic basin V(C4) increase in population. The final structure of the zwitterion I4F possesses 5.1 kcal/mol lower energy than the transition state.
Scheme 9. Simplified representation of the molecular mechanism of the reaction of phenyl azide 1c and ethyl propiolate 2 by Lewis-like structures based on the topological analysis of the ELF along the reaction path leading to I4F.
Scheme 9. Simplified representation of the molecular mechanism of the reaction of phenyl azide 1c and ethyl propiolate 2 by Lewis-like structures based on the topological analysis of the ELF along the reaction path leading to I4F.
Molecules 28 08152 sch009
Figure 8. ELF valence basins attractors of structures of the most important topological changes, creation of negative charge on C4, and creation of N1-C5 single bond during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to I4F. Structures with apostrophes represent the first structure before the critical point.
Figure 8. ELF valence basins attractors of structures of the most important topological changes, creation of negative charge on C4, and creation of N1-C5 single bond during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to I4F. Structures with apostrophes represent the first structure before the critical point.
Molecules 28 08152 g008
Table 6. ELF valence basin populations of the IRC points, MCFI4F, defining the ten different phases characterizing the reaction of the phenyl azide 1c and ethyl propiolate 2. The stationary points 1c, 2, MCF, TSF, and I4F are also included. Distances are given in angstroms, Å, electron populations in the average number of electrons, [e], relative energies in kcal·mol-1, and GEDT values in an average number of electrons, [e].
Table 6. ELF valence basin populations of the IRC points, MCFI4F, defining the ten different phases characterizing the reaction of the phenyl azide 1c and ethyl propiolate 2. The stationary points 1c, 2, MCF, TSF, and I4F are also included. Distances are given in angstroms, Å, electron populations in the average number of electrons, [e], relative energies in kcal·mol-1, and GEDT values in an average number of electrons, [e].
Structures1c2MCF P1 P2 P3 P4 P5TSFP6 P7 P8 P9 I4F
Phases I II III IV V VI VII VIII IX X
d1(N1-C5) 3.740 2.189 2.075 1.817 1.778 1.7221.7001.558 1.492 1.446 1.443 1.428
GEDT 0.00 −0.03 −0.06 −0.26 −0.30 −0.36−0.38−0.49 −0.51 −0.50 −0.50 −0.48
dE −6.6 21.2 29.2 46.6 47.7 48.448.547.1 46.1 44.8 44.5 43.4
V(N1)3.76 3.74 3.74 3.74 2.89 2.78 2.682.652.53 2.49 2.46 2.45 2.41
V(N1.N2)1.80 1.81 1.90 1.90 1.79 1.77 1.741.731.70 3.71 2.90 3.79 3.84
V′(N1.N2)2.34 2.37 2.24 2.23 2.10 2.08 2.052.042.01
V(N2.N3)2.50 2.47 2.54 2.57 2.80 2.84 2.902.921.53 1.50 1.49 1.51 2.85
V(N3)3.37 3.39 3.34 3.29 3.07 3.03 2.982.952.83 2.79 2.78 2.78 2.82
V(N3.C7)1.83 1.83 1.84 1.86 1.89 1.90 1.911.911.94 1.95 1.97 1.97 1.95
V(C4.C5) 2.732.78 5.19 4.97 4.07 1.95 3.903.863.63 3.53 3.42 3.41 3.37
V′(C4.C5) 2.562.56 2.06
V(C5.C6) 2.402.41 2.47 2.49 2.46 2.43 2.412.402.37 2.36 2.36 2.36 2.36
V(N1.C5) 1.00 1.17 1.351.401.68 1.80 1.92 1.93 2.01
V(N2) 0.90
V′(N2.N3) 1.47 1.49 1.41 1.39 1.39
V(C4) 0.26 1.62 1.73 1.871.912.11 2.16 2.19 2.19 2.20

3. Computational Details

All calculations reported in this paper were performed using the “Ares” infrastructure in the “Cyfronet” computational center in Cracow. Hybrid functional wb97xd with the 6-311+G(d) basis set included in the GAUSSIAN 09 package [61] was used. Previously, we found that the wb97xd/6-311+G(d) calculations illustrate the molecular mechanisms of several different 32CA reactions well [25,49,62,63].
Global electronic properties of reactants were estimated according to the equations recommended by Parr and Domingo [64,65,66]. According to Domingo’s recommendation, for this purpose, the B3LYP/6-31G(d) level of theory was used. In particular, the electronic chemical potentials (μ) and chemical hardness (η) were evaluated in terms of one-electron energies of FMO (EHOMO and ELUMO) using Equations (1) and (2):
μ ≈ (EHOMO + ELUMO)/2
η ≈ ELUMO − EHOMO
Next, the values of μ and η were then used for the calculation of global electrophilicity (ω) according to the Formula (3) [66], and the global nucleophilicity (N) [67] can be expressed in terms of Equation (4)
ω = μ2/2η
N = EHOMO − ELUMO (rerracyanoethene)
The local electrophilicity (ωk) [68] condensed to atom k was calculated by projecting the index ω onto any reaction center k in the molecule using the Parr function P+k:
ωk = P+k·ω
The local nucleophilicity (Nk) [68] condensed to atom k was calculated using global nucleophilicity N and Paar function Pk according to the formula:
Nk = Pk·N
For structure optimization of the reactants, intermediates, and products, the Berny algorithm was applied. First-order saddle points were localized using the QST2 or QST3 procedure. Stationary points were checked by vibrational frequency analyses to determine whether they constituted minima or maxima on the potential energy surface (PES). All transition structures showed a single imaginary frequency (νi), whereas reactants, products, and pre-reaction complexes had none. The intrinsic reaction coordinate (IRC) path was traced in order to check the energy profiles connecting each transition structure to the two associated minima of the proposed mechanism. The calculations were carried out for the simulated presence of ethanol as the reaction medium (the PCM model [69] was used). All calculations were performed for 298 K and 1 atm pressure. The absolute entropies of critical structures (S298) were estimated from the complete vibrational analysis. Enthalpies (H298) were corrected to Gibbs free energies (G298) using the calculated entropies.
Global electron density transfer (GEDT) [54] was calculated according to Formula (7):
GEDT = −ΣqA
where qA is the net charge and the sum is taken for all the atoms of nitroalkene.
Indexes of σ-bonds development (l) were calculated according to Formula (8) [36]:
l X Y = 1 r X Y TS r X Y P r X Y P
where rTSX−Y is the distance between the reaction centers X and Y in the transition structure and rPX−Y is the same distance in the corresponding product.
The kinetic parameters and essential properties of critical structures are presented in Table 2 and Table 3. Consistent with the previous convention [70,71], in this paper, the letters MC and TS are designated to the pre-reaction complex and TS, respectively. They are distinguished by appending the letters AD depending on the reaction pathway.
The Electron Localization Function (ELF) [72] analysis was performed using the TopMod 09 package [73] at the standard cubical grid of step size of 0.1 Bohr. To visualize the molecular geometries and ELF basin attractors, the GaussView program [74] was used.

4. Conclusions

The wb97xd/6-311+G(d) computational study evidently underlies the hypothesis of the stepwise mechanism via a zwitterionic intermediate of the 32CA reactions between aryl azides and ethyl propiolate. Our investigations show that all title reactions are realized via a polar one-step reaction mechanism without the intervention of zwitterionic intermediates. It should be underlined, however, that the new single bond within the formed triazole ring is always formed faster at the activated electrophilically terminal carbon atom on the acetylenic moiety. Some zwitterionic intermediates can be theoretically formed in competitive reaction paths. These transformations are, however, not favored from thermodynamic and kinetic points of view. Additionally, our quantumchemical calculations suggest that within the reaction course, a minor amount of 3-aryl-4-carboethoxy-1,2,3-triazole not isolated and not identified by Abu-Orabi and co-workers can be formed. This suggests a further, deeper experimental reexamination of title reactions should be performed. This will be the subject of our future research. The ELF analysis revealed that the 32CA reaction takes place via a non-concerted two-stage one-step mechanism during which the first single bond is created by donation of the nonbonding electron density of N3 to C4, and the second single bond is created by shearing of the pseudoradical center C5 with the nonbonding electron density of N1. Meanwhile, new bonds in the theoretically possible zwitterions are created by the donation of the nonbonding electron density of the appropriate nitrogen atom to the carbon atom.

Author Contributions

Methodology, R.J.; formal analysis, E.D. and A.W.; investigation, R.J. and P.W.; data curation, E.D. and A.W.; writing—original draft, R.J. and P.W.; supervision, R.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

We gratefully acknowledge the Polish high-performance computing infrastructure PLGrid (HPC Centers: ACK Cyfronet AGH) for providing computer facilities and support within computational grant no. PLG/2023/016659.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Katritzky, A.R. Introduction: Heterocycles. Chem. Rev. 2004, 104, 2125–2126. [Google Scholar] [CrossRef]
  2. Kerru, N.; Gummidi, L.; Maddila, S.; Gangu, K.K.; Jonnalagadda, S.B. A Review on Recent Advances in Nitrogen-Containing Molecules and Their Biological Applications. Molecules 2020, 25, 1909. [Google Scholar] [CrossRef] [PubMed]
  3. Saliyeva, L.M.; Dyachenko, I.V.; Danyliuk, I.Y.; Vovk, M.V. Di-, tetra-, and perhydropyrrolo[1,2-a]imidazoles: The Methods of Synthesis and Some Aspects of Application. Chem. Heterocycl. Compd. 2022, 58, 661–680. [Google Scholar] [CrossRef]
  4. Hyjek, K.; Jodłowski, P. Metal-organic frameworks for efficient drug adsorption and delivery. Sci. Rad. 2023, 2, 118–189. [Google Scholar] [CrossRef]
  5. Komkov, A.V.; Sukhanova, A.А.; Menchikov, L.G.; Zavarzin, I.V. o-Aminopyrimidine Aldehydes and Ketones: Synthesis and use as Precursors to Fused Pyrimidines. Chem. Heterocycl. Compd. 2022, 58, 547–566. [Google Scholar] [CrossRef]
  6. Danyliuk, I.Y.; Vovk, M.V. Tetrahydroazepines with an annulated five-membered heteroaromatic ring. Chem. Heterocycl. Compd. 2022, 58, 567–581. [Google Scholar] [CrossRef]
  7. Obernikhina, N.V.; Kachaeva, M.V.; Kachkovsky, O.D.; Brovarets, V.S. In silico Study of Conjugated Nitrogen Heterocycles Affinity in their Biological Complexes. Chem. Heterocycl. Compd. 2022, 58, 412–420. [Google Scholar] [CrossRef]
  8. Khomenko, D.M.; Doroshchuk, R.O.; Ohorodnik, Y.M.; Ivanova, H.V.; Zakharchenko, B.V.; Raspertova, I.V.; Vaschenko, O.V.; Dobrydnev, A.V.; Grygorenko, O.O.; Lampeka, R.D. Expanding the chemical space of 3(5)-functionalized 1,2,4-triazoles. Chem. Heterocycl. Compd. 2022, 58, 116–128. [Google Scholar] [CrossRef]
  9. Zhang, S.; Xu, Z.; Gao, C.; Ren, Q.-C.; Chang, L.; Lv, Z.-S.; Feng, L.-S. Triazole derivatives and their anti-tubercular activity. Eur. J. Med. Chem. 2017, 138, 501–513. [Google Scholar] [CrossRef]
  10. Keri, R.S.; Patil, S.A.; Budagumpi, S.; Nagaraja, B.M. Triazole: A Promising Antitubercular Agent. Chem. Biol. Drug. Des. 2015, 86, 410–423. [Google Scholar] [CrossRef]
  11. Chu, X.M.; Wang, C.; Wang, W.L.; Liang, L.L.; Liu, W.; Gong, K.K.; Sun, K.L. Triazole derivatives and their antiplasmodial and antimalarial activities. Eur. J. Med. Chem. 2019, 15, 206–223. [Google Scholar] [CrossRef] [PubMed]
  12. Xia, Y.; Liu, Y.; Wan, J.; Wang, M.; Rocchi, P.; Qu, F.; Iovanna, J.L.; Peng, L. Novel Triazole Ribonucleoside Down-Regulates Heat Shock Protein 27 and Induces Potent Anticancer Activity on Drug-Resistant Pancreatic Cancer. J. Med. Chem. 2009, 52, 6083–6096. [Google Scholar] [CrossRef] [PubMed]
  13. Groll, A.H.; Townsend, R.; Desai, A.; Azie, N.; Jones, M.; Engelhardt, M.; Schmitt-Hoffman, A.; Brüggemann, R.J.M. Drug-drug interactions between triazole antifungal agents used to treat invasive aspergillosis and immunosuppressants metabolized by cytochrome P450 3A4. Transpl. Infect. Dis. 2017, 19, e12751. [Google Scholar] [CrossRef] [PubMed]
  14. Lass-Flörl, C. Triazole Antifungal Agents in Invasive Fungal Infections. Drugs 2011, 71, 2405–2419. [Google Scholar] [CrossRef] [PubMed]
  15. Groll, A.H.; Townsend, R.; Desai, A.; Azie, N.; Jones, M.; Engelhardt, M.; Schmitt-Hoffman, A.; Brüggemann, R.J.M. The Enzymatic Basis of Drug-Drug Interactions with Systemic Triazole Antifungals. Clin. Pharmacokinet. 2008, 47, 779–792. [Google Scholar]
  16. Ayati, A.; Emami, S.; Foroumadi, A. The importance of triazole scaffold in the development of anticonvulsant agents. Eur. J. Med. Chem. 2016, 109, 380–392. [Google Scholar] [CrossRef]
  17. Kharb, R.; Sharma, P.C.; Yar, M.S. Pharmacological significance of triazole scaffold. J. Enzyme Inhib. Med. Chem. 2011, 26, 1–21. [Google Scholar] [CrossRef]
  18. Kras, J.; Sadowski, M.; Zawadzińska, K.; Nagatsky, R.; Woliński, P.; Kula, K.; Łapczuk, A. Thermal [3+2] cycloaddition reactions as most universal way for the effective preparation of five-membered nitrogen containing heterocycles. Sci. Rad. 2023, 2, 247–267. [Google Scholar] [CrossRef]
  19. Tiwari, G.; Khanna, A.; Mishra, V.K.; Sagar, R. Recent developments on microwave-assisted organic synthesis of nitrogen- and oxygen-containing preferred heterocyclic scaffolds. RSC Adv. 2023, 13, 32858–32892. [Google Scholar] [CrossRef]
  20. Appukkuttan, P.; Mehta, V.P.; Van der Eycken, E.V. Microwave-assisted cycloaddition reactions. Chem. Soc. Rev. 2009, 39, 1467–1477. [Google Scholar] [CrossRef]
  21. Domingo, L.R.; Ríos-Gutiérrez, M.; Pérez, P. Why is phenyl azide so unreactive in [3+2] cycloaddition reactions? Demystifying Sustmann’s paradigmatic parabola. Org. Chem. Front. 2023, 10, 5579–5591. [Google Scholar] [CrossRef]
  22. Gomonov, K.A.; Pilipenko, I.A. Formation of Five- and Six-membered Oxygen-containing Heterocycles on the Basis of 1-halo-1-nitroalkenes. Chem. Heterocycl. Compd. 2023, 59, 1–4. [Google Scholar] [CrossRef]
  23. Milišiūnaitė, V.; Plytninkienė, E.; Bakšienė, R.; Bieliauskas, A.; Krikštolaitytė, S.; Račkauskienė, G.; Arbačiauskienė, E.; Šačkus, A. Convenient Synthesis of Pyrazolo[4′,3′:5,6]pyrano[4,3-c][1,2]oxazoles via Intramolecular Nitrile Oxide Cycloaddition. Molecules 2021, 26, 5604. [Google Scholar] [CrossRef] [PubMed]
  24. Zawadzińska, K.; Ríos-Gutiérrez, M.; Kula, K.; Woliński, P.; Mirosław, B.; Krawczyk, T.; Jasiński, R. The Participation of 3,3,3-Trichloro-1-nitroprop-1-ene in the [3+2] Cycloaddition Reaction with Selected Nitrile N-Oxides in the Light of the Experimental and MEDT Quantum Chemical Study. Molecules 2021, 26, 6774. [Google Scholar] [CrossRef]
  25. Sadowski, S.; Mudyna, A.; Knap, K.; Demchuk, O.M.; Łapczuk, A. Synthesis of (Z)-N-aryl-C-(pyrid-3-yl)-nitrones. Sci. Rad. 2023, 2, 319–324. [Google Scholar] [CrossRef]
  26. Jasiński, E.; Żmigrodzka, M.; Dresler, E.; Kula, K. A Full Regioselective and Stereoselective Synthesis of 4-Nitroisoxazolidines via Stepwise [3+2] Cycloaddition Reactions between (Z)-C-(9-Anthryl)-N-arylnitrones and (E)-3,3,3-Trichloro-1-nitroprop-1-ene: Comprehensive Experimental and Theoretical Study. J. Heterocycl. Chem. 2017, 54, 3314–3320. [Google Scholar] [CrossRef]
  27. Woliński, P.; Kącka-Zych, A.; Dziuk, B.; Ejsmont, K.; Łapczuk-Krygier, A.; Dresler, E. The structural aspects of the transformation of 3-nitroisoxazoline-2-oxide to 1-aza-2,8-dioxabicyclo[3.3.0]octane derivatives: Experimental and MEDT theoretical study. J. Mol. Struct. 2019, 1192, 27–34. [Google Scholar] [CrossRef]
  28. Żmigrodzka, M.; Sadowski, M.; Kras, J.; Desler, E.; Demchuk, O.M.; Kula, K. Polar [3+2] cycloaddition between N-methyl azomethine ylide and trans-3,3,3-trichloro-1-nitroprop-1-ene. Sci. Rad. 2022, 1, 26–35. [Google Scholar] [CrossRef]
  29. Rostovtsev, V.V.; Green, L.G.; Fokin, V.V.; Sharpless, K.B. A stepwise huisgen cycloaddition process: Copper(I)-catalyzed regioselective “ligation” of azides and terminal alkynes. Angew. Chem. Int. Ed. Engl. 2002, 41, 2596–2599. [Google Scholar] [CrossRef]
  30. Levandowski, B.J.; Gamache, R.F.; Murphy, J.M.; Houk, K.N. Readily Accessible Ambiphilic Cyclopentadienes for Bioorthogonal Labeling. J. Am. Chem. Soc. 2018, 140, 6426–6431. [Google Scholar] [CrossRef]
  31. Yoshida, S.; Kuribara, T.; Ito, H.; Meguro, T.; Nishiyama, Y.; Karaki, F.; Hatakeyama, Y.; Koike, Y.; Kii, I.; Hosoya, T. A facile preparation of functional cycloalkynes via an azide-to-cycloalkyne switching approach. Chem. Commun. 2019, 19, 3556–3559. [Google Scholar] [CrossRef] [PubMed]
  32. Siadati, S.A.; Rezazadeh, S. The extraordinary gravity of three atom 4π-components and 1,3-dienes to C20-nXn fullerenes; a new gate to the future of Nano technology. Sci. Rad. 2022, 1, 46–68. [Google Scholar] [CrossRef]
  33. Boren, B.C.; Narayan, S.; Rasmussen, L.K.; Zhang, L.; Zhao, H.; Lin, Z.; Jia, G.; Fokin, V.V. Ruthenium-Catalyzed Azide−Alkyne Cycloaddition: Scope and Mechanism. J. Am. Chem. Soc. 2008, 130, 8923–8930. [Google Scholar] [CrossRef] [PubMed]
  34. Abu-Orabi, S.T.; Atfah, M.A.; Jibril, I.; Mari’i, F.M.; Ali, A.A.-S. Dipolar cycloaddition reactions of organic azides with some acetylenic compounds. J. Heterocycl. Chem. 1989, 26, 1461–1468. [Google Scholar] [CrossRef]
  35. Kula, K.; Kącka-Zych, A.; Łapczuk-Krygier, A.; Wzorek, Z.; Nowak, A.K.; Jasiński, R. Experimental and Theoretical Mechanistic Study on the Thermal Decomposition of 3,3-diphenyl-4-(trichloromethyl)-5-nitropyrazoline. Molecules 2021, 26, 1364. [Google Scholar] [CrossRef] [PubMed]
  36. Jasiński, R. A stepwise, zwitterionic mechanism for the 1,3-dipolar cycloaddition between (Z)-C-4-methoxyphenyl-N-phenylnitrone and gem-chloronitroethene catalysed by 1-butyl-3-methylimidazolium ionic liquid cations. Tetrahedron Lett. 2015, 56, 532–535. [Google Scholar] [CrossRef]
  37. Jasiński, R. A new insight on the molecular mechanism of the reaction between (Z)-C,N-diphenylnitrone and 1,2-bismethylene-3,3,4,4,5,5-hexamethylcyclopentane. J. Mol. Graph. Model. 2020, 94, 107461. [Google Scholar] [CrossRef]
  38. Domingo, L.R.; Saéz, J.A.; Zaragozá, R.J.; Arnó, M. Understanding the Participation of Quadricyclane as Nucleophile in Polar [2σ + 2σ + 2π] Cycloadditions toward Electrophilic π Molecules. J. Org. Chem. 2008, 73, 8791–8799. [Google Scholar] [CrossRef]
  39. Domingo, L.R.; Ríos-Gutiérrez, M. A Useful Classification of Organic Reactions Based on the Flux of the Electron Density. Sci. Rad. 2023, 2, 1–24. [Google Scholar] [CrossRef]
  40. Mondal, A.; Mohammad-Salim, H.A.; Acharjee, N. Unveiling substituent effects in [3+2] cycloaddition reactions of benzonitrile N-oxide and benzylideneanilines from the molecular electron density theory perspective. Sci. Rad. 2023, 2, 75–92. [Google Scholar] [CrossRef]
  41. Ríos-Gutiérrez, M.; Domingo, L.R. Unravelling the Mysteries of the [3+2] Cycloaddition Reactions. Eur. J. Org. Chem. 2019, 2019, 267–282. [Google Scholar] [CrossRef]
  42. Domingo, L.R.; Ríos-Gutiérrez, M.; Pérez, P. Unveiling the Chemistry of Higher-Order Cycloaddition Reactions within the Molecular Electron Density Theory. Chemistry 2022, 4, 735–752. [Google Scholar] [CrossRef]
  43. Domingo, L.R.; Aurell, M.J.; Pérez, P.; Contreras, R. Quantitative characterization of the global electrophilicity power of common diene/dienophile pairs in Diels–Alder reactions. Tetrahedron 2002, 58, 4417–4423. [Google Scholar] [CrossRef]
  44. Domingo, L.R.; Ríos-Gutiérrez, M.; Pérez, P. Applications of the Conceptual Density Functional Theory Indices to Organic Chemistry Reactivity. Molecules 2016, 21, 748. [Google Scholar] [CrossRef]
  45. Domingo, L.R. Molecular Electron Density Theory: A Modern View of Reactivity in Organic Chemistry. Molecules 2016, 21, 1319. [Google Scholar] [CrossRef] [PubMed]
  46. Tolstikov, G.A.; Shults, E.E.; Malikova, T.S.; Spirikhin, L.V. The Importance of Preliminary Orientation in [4+2]-Cycloadditions of Dienes and Dienophiles with Complex Structures. Mendeleev Commun. 1994, 4, 60–62. [Google Scholar] [CrossRef]
  47. Aitouna, A.O.; Barhoumi, A.; Zeroual, A. A Mechanism Study and an Investigation of the Reason for the Stereoselectivity in the [4+2] Cycloaddition Reaction between Cyclopentadiene and Gem-substituted Ethylene Electrophiles. Sci. Rad. 2023, 2, 217–228. [Google Scholar] [CrossRef]
  48. Mondal, A.; Acharjee, N. Unveiling the exclusive stereo and site selectivity in [3+2] cycloaddition reactions of a tricyclic strained alkene with nitrile oxides from the molecular electron density theory perspective. Chem. Heterocycl. Comp. 2023, 59, 145–154. [Google Scholar] [CrossRef]
  49. Kula, K.; Sadowski, M. Regio- and stereoselectivity of [3+2] cycloaddition reactions between (Z)-1-(anthracen-9-yl)-N-methyl nitrone and analogs of trans-β-nitrostyrene on the basis of MEDT computational study. Chem. Heterocycl. Comp. 2023, 59, 138–144. [Google Scholar] [CrossRef]
  50. Gillies, C.W.; Gillies, J.Z.; Suenram, D.R.; Lovas, F.J.; Kraka, E.; Cremer, D. Van der Waals complexes in 1,3-dipolar cycloaddition reactions: Ozone-ethylene. J. Am. Chem. Soc. 1991, 113, 2412–2421. [Google Scholar] [CrossRef]
  51. Gillies, J.Z.; Gillies, C.W.; Lovas, F.J.; Matsumura, K.; Suenram, R.D.; Kraka, E.; Cremer, D. Van der Waals complexes of chemically reactive gases: Ozone-acetylene. J. Am. Chem. Soc. 1991, 113, 6408–6415. [Google Scholar] [CrossRef]
  52. Lupinski, J.H. The chargé transfer complex between β-carotene and iodine. II. Characterization of complex. J. Phys. Chem. 1963, 67, 2725–2728. [Google Scholar] [CrossRef]
  53. Bahnick, D.A.; Person, W.B. Raman Intensity Study of Charge-Transfer Complexes of ICN. J. Chem. Phys. 1968, 48, 5637–5645. [Google Scholar] [CrossRef]
  54. Domingo, L.R. A new C–C bond formation model based on the quantum chemical topology of electron density. RSC Adv. 2014, 4, 32415–32428. [Google Scholar] [CrossRef]
  55. Woliński, P.; Kącka-Zych, A.; Wróblewska, A.; Wielgus, E.; Dolot, R.; Jasiński, R. Fully Selective Synthesis of Spirocyclic-1,2-oxazine N-Oxides via Non-Catalysed Hetero Diels-Alder Reactions with the Participation of Cyanofunctionalysed Conjugated Nitroalkenes. Molecules 2023, 28, 4586. [Google Scholar] [CrossRef] [PubMed]
  56. Kula, K.; Łapczuk, A.; Sadowski, M.; Kras, J.; Zawadzińska, K.; Demchuk, O.M.; Gaurav, G.K.; Wróblewska, A.; Jasiński, R. On the Question of the Formation of Nitro-Functionalized 2,4-Pyrazole Analogs on the Basis of Nitrylimine Molecular Systems and 3,3,3-Trichloro-1-Nitroprop-1-Ene. Molecules 2022, 27, 8409. [Google Scholar] [CrossRef] [PubMed]
  57. Woliński, P.; Kącka-Zych, A.; Demchuk, O.M.; Łapczuk-Krygier, A.; Mirosław, B.; Jasiński, R. Clean and molecularly programmable protocol for preparation of bis-heterobiarylic systems via a domino pseudocyclic reaction as a valuable alternative for TM-catalyzed cross-couplings. J. Clean. Prod. 2020, 275, 122086. [Google Scholar] [CrossRef]
  58. Jasiński, R. β-Trifluoromethylated nitroethenes in Diels-Alder reaction with cyclopentadiene: A DFT computational study. J. Fluor. Chem. 2018, 206, 1–7. [Google Scholar] [CrossRef]
  59. Jasiński, R. Searching for zwitterionic intermediates in Hetero Diels–Alder reactions between methyl α,p-dinitrocinnamate and vinyl-alkyl ethers. Comput. Theor. Chem. 2014, 1046, 93–98. [Google Scholar] [CrossRef]
  60. Krokidis, K.; Noury, S.; Silvi, B. Characterization of Elementary Chemical Processes by Catastrophe Theory. J. Phys. Chem. A. 1997, 101, 7277–7282. [Google Scholar] [CrossRef]
  61. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  62. Kącka-Zych, A. The Molecular Mechanism of the Formation of Four-Membered Cyclic Nitronates and Their Retro (3 + 2) Cycloaddition: A DFT Mechanistic Study. Molecules 2021, 26, 4786. [Google Scholar] [CrossRef] [PubMed]
  63. Fryźlewicz, A.; Kącka-Zych, A.; Demchuk, O.M.; Mirosław, B.; Woliński, P.; Jasiński, R. Green synthesis of nitrocyclopropane-type precursors of inhibitors for the maturation of fruits and vegetables via domino reactions of diazoalkanes with 2-nitroprop-1-ene. J. Clean. Prod. 2021, 292, 126079. [Google Scholar] [CrossRef]
  64. Pérez, P.; Domingo, L.R.; Aizman, A.; Contreras, R. The electrophilicity index in organic chemistry. Theor. Aspects Chem. React. 2008, 19, 139–201. [Google Scholar]
  65. Perez, P.; Domingo, L.R.; Aurell, M.J.; Contreras, R. Quantitative characterization of the global electrophilicity patternof some reagents involved in 1,3-dipolar cycloaddition reactions. Tetrahedron 2003, 59, 3117–3125. [Google Scholar] [CrossRef]
  66. Parr, R.G.; Szentpály, L.V.; Liu, S. Electrophilicity Index. J. Am. Chem. Soc. 1999, 121, 1922–1924. [Google Scholar] [CrossRef]
  67. Domingo, L.R.; Chamorro, E.; Pérez, P. Understanding the Reactivity of Captodative Ethylenes in Polar Cycloaddition Reactions. A Theoretical Study. J. Org. Chem. 2008, 73, 4615–4624. [Google Scholar] [CrossRef] [PubMed]
  68. Domingo, L.R.; Pérez, P.; Sáez, J.A. Understanding the local reactivity in polar organic reactions through electrophilic and nucleophilic Parr functions. RSC Adv. 2013, 3, 1486–1494. [Google Scholar] [CrossRef]
  69. Cossi, M.; Rega, N.; Scalmani, G.; Barone, V. Energies, Structures, and Electronic Properties of Molecules in Solution with the C-PCM Solvation Model. J. Comp. Chem. 2003, 24, 669–681. [Google Scholar] [CrossRef]
  70. Kącka-Zych, A. Understanding the uniqueness of the stepwise [4 + 1] cycloaddition reaction between conjugated nitroalkenes and electrophilic carbene systems with a molecular electron density theory perspective. Int. J. Quant. Chem. 2021, 121, e26440. [Google Scholar] [CrossRef]
  71. Jasiński, R. Stepwise, zwitterionic course of hetero-Diels–Alder reaction between 1,2,4-triazine molecular systems and 2-cyclopropylidene-1,3-dimethylimidazoline. Chem. Heterocycl. Comp. 2022, 58, 260–262. [Google Scholar] [CrossRef]
  72. Becke, A.D.; Edgecombe, K.E. A simple measure of electron localization in atomic and molecular systems. J. Chem. Phys. 1990, 92, 5397–5403. [Google Scholar] [CrossRef]
  73. Noury, S.; Krokidis, X.; Fuster, F.; Silvi, B. Computational Tools for the Electron Localization Function Topological Analysis. Comput. Chem. 1999, 23, 597–604. [Google Scholar] [CrossRef]
  74. Dennington, R.; Keith, T.A.; Millam, J.M. GaussView Version 6; Semichem Inc.: Shawnee Mission, KS, USA, 2016. [Google Scholar]
Scheme 1. The regioselectivity of the 32CA reactions between benzyl azide and phenoxyethyne.
Scheme 1. The regioselectivity of the 32CA reactions between benzyl azide and phenoxyethyne.
Molecules 28 08152 sch001
Scheme 2. The 32CA reactions between azides and strained alkynes.
Scheme 2. The 32CA reactions between azides and strained alkynes.
Molecules 28 08152 sch002
Scheme 3. The rhutenium-catalyzed 32CA reactions between benzyl azide and simple alkynes.
Scheme 3. The rhutenium-catalyzed 32CA reactions between benzyl azide and simple alkynes.
Molecules 28 08152 sch003
Scheme 4. The regioselectivity and reaction conditions of the 32CAs between aryl azides 1ae and ethyl propiolate 2.
Scheme 4. The regioselectivity and reaction conditions of the 32CAs between aryl azides 1ae and ethyl propiolate 2.
Molecules 28 08152 sch004
Scheme 5. Enthalpy profiles for the reactions between cyclopentadiene 1 with isopropyl 3-nitroprop-2-enate 2 in the DCM environment in light of the wb97xd/6-311+G(d) (PCM) calculations Theoretical possible ways of the 32CA reactions between aryl azides 1ae and ethyl propiolate 2.
Scheme 5. Enthalpy profiles for the reactions between cyclopentadiene 1 with isopropyl 3-nitroprop-2-enate 2 in the DCM environment in light of the wb97xd/6-311+G(d) (PCM) calculations Theoretical possible ways of the 32CA reactions between aryl azides 1ae and ethyl propiolate 2.
Molecules 28 08152 sch005
Figure 1. The distribution of local electrophilicities along the azide 1c molecule.
Figure 1. The distribution of local electrophilicities along the azide 1c molecule.
Molecules 28 08152 g001
Figure 2. The distribution of local nucleophilicities along the ethyl propiolate 2 molecule.
Figure 2. The distribution of local nucleophilicities along the ethyl propiolate 2 molecule.
Molecules 28 08152 g002
Figure 3. Gibbs free energy profiles for all explored paths of the reaction between phenyl azide 1c and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.
Figure 3. Gibbs free energy profiles for all explored paths of the reaction between phenyl azide 1c and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.
Molecules 28 08152 g003
Figure 4. Views of critical structures of the 32CA reaction between phenyl azide 1c and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.
Figure 4. Views of critical structures of the 32CA reaction between phenyl azide 1c and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.
Molecules 28 08152 g004
Scheme 6. Detected ways of reactions between aryl azides 1ae and ethyl propiolate 2.
Scheme 6. Detected ways of reactions between aryl azides 1ae and ethyl propiolate 2.
Molecules 28 08152 sch006
Figure 5. Views of critical structures of the creation of zwitterions via reaction between phenyl azide 1c and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.
Figure 5. Views of critical structures of the creation of zwitterions via reaction between phenyl azide 1c and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.
Molecules 28 08152 g005
Figure 6. ELF valence basins attractors of structures of the most important topological changes, creation of N3-C4 and N1-C5 single bonds, during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to 3c. Structures with apostrophes represent the first structure before the critical point.
Figure 6. ELF valence basins attractors of structures of the most important topological changes, creation of N3-C4 and N1-C5 single bonds, during the reaction of phenyl azide 1c and ethyl propiolate 2 along the reaction path leading to 3c. Structures with apostrophes represent the first structure before the critical point.
Molecules 28 08152 g006
Scheme 7. Simplified representation of the molecular mechanism of 32CA of phenyl azide 1c and ethyl propiolate 2 by Lewis-like structures based on the topological analysis of the ELF along the reaction path leading to 3c.
Scheme 7. Simplified representation of the molecular mechanism of 32CA of phenyl azide 1c and ethyl propiolate 2 by Lewis-like structures based on the topological analysis of the ELF along the reaction path leading to 3c.
Molecules 28 08152 sch007
Table 1. Global and local nucleophilicities at azide molecules 1ac (in eV).
Table 1. Global and local nucleophilicities at azide molecules 1ac (in eV).
Molecules 28 08152 i001
NoRNN1N3
1aMeO3.450.650.47
1bMe3.110.710.57
1cH2.920.720.64
1dCl2.800.600.51
1eNO22.120.530.52
Table 2. Energetic characterization of key transitions within 32CA reaction between aryl azides 1ae and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations (ΔH, ΔG are in kcal/mol; ΔS are in cal/molK).
Table 2. Energetic characterization of key transitions within 32CA reaction between aryl azides 1ae and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations (ΔH, ΔG are in kcal/mol; ΔS are in cal/molK).
ReactionPathTransitionΔHΔSΔG
1a+2A1a+2→MCA−3.3−25.54.3
1a+2→TSA17.4−40.729.6
1a+2→3a−76.6−48.6−62.1
B1a+2→MCB−3.7−39.38.0
1a+2→TSB17.1−44.330.3
1a+2→4a−74.4−45.4−60.9
E1a+2→MCE−5.9−37.35.2
1a+2→TSE28.3−35.338.8
1a+2→I3E27.2−39.038.8
F1a+2→MCF−5.7−35.34.9
1a+2→TSF46.8−41.759.3
1a+2→I4F42.6−43.655.6
1b+2A1b+2→MCA−5.7−34.14.4
1b+2→TSA17.8−42.430.4
1b+2→3b−76.1−49.2−61.5
B1b+2→MCB−6.0−39.15.7
1b+2→TSB17.0−49.231.7
1b+2→4b−74.3−53.5−58.4
E1b+2→MCE−6.0−39.15.7
1b+2→TSE28.3−43.941.3
1b+2→I3E28.1−36.839.0
F1b+2→MCF−5.7−34.14.4
1b+2→TSF47.6−39.759.5
1b+2→I4F43.8−41.656.2
1c+2A1c+2→MCA−3.2−29.05.4
1c+2→TSA18.0−41.030.2
1c+2→3c−75.7−48.6−61.3
B1c+2→MCB−5.0−32.64.7
1c+2→TSB18.0−44.831.3
1c+2→4c−73.2−45.1−59.7
E1c+2→MCE−4.9−34.45.3
1c+2→TSE29.4−37.540.5
1c+2→I3E28.7−35.939.4
F1c+2→MCF−5.3−34.75.0
1c+2→TSF48.3−40.560.4
1c+2→I4F44.6−42.857.4
1d+2A1d+2→MCA−3.4−30.45.7
1d+2→TSA18.2−41.730.6
1d+2→3d−75.3−48.7−60.8
B1d+2→MCB−5.9−36.65.0
1d+2→TSB18.1−44.131.2
1d+2→4d−73.1−47.1−59.0
E1d+2→MCE−5.9−36.34.9
1d+2→TSE29.8−37.441.0
1d+2→I3E29.4−38.040.7
F1d+2→MCF−5.6−34.54.7
1d+2→TSF48.8−41.361.1
1d+2→I4F45.3−41.757.7
1e+2A1e+2→MCA−4.2−35.46.3
1e+2→TSA18.7−38.930.2
1e+2→3e−73.5−48.7−58.9
B1e+2→MCB−5.2−35.05.2
1e+2→TSB18.9−42.631.6
1e+2→4e−71.0−46.7−57.1
E1e+2→MCE−5.6−36.85.4
1e+2→TSE31.9−38.943.4
1e+2→I3E32.0−33.942.1
F1e+2→MCF−5.3−34.85.0
1e+2→TSF51.0−40.963.2
1e+2→I4F48.0−40.760.2
Table 3. Key parameters of critical structures of the 32CA reaction between aryl azides 1ae and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.
Table 3. Key parameters of critical structures of the 32CA reaction between aryl azides 1ae and ethyl propiolate 2 in the ethanolic solution according to the wb97xd/6-311+G(d)(PCM) calculations.
ReactionPathStructureInteratomic Distances r [Å]GEDT
N1-N2N2-N3N3-C4C4-C5C5-N1[e]
1a+2A1a1.1271.226
2 1.200
MCA1.1271.2263.4471.1993.403
TSA1.1491.2542.0311.2322.2790.30
3a1.2901.3491.3431.3751.359
BMCB1.1251.2273.6551.2003.798
TSB1.1601.2512.1911.2312.0750.25
4a1.3021.3321.3611.3741.350
EMCE1.1271.2263.5231.1994.030
TSE1.1191.2461.7201.2604.5610.75
I3E1.1141.2571.5151.3034.5480.89
FMCF1.1261.2273.9071.1993.893
TSF1.1561.1954.4591.2841.7070.55
I4F1.1771.1924.6611.3551.4280.70
1b+2A1b1.1261.227
MCA1.1261.2273.4611.1993.416
TSA1.1491.2532.0331.2322.2720.31
3b1.2891.3491.3431.3751.359
MCB1.1261.2273.6491.1993.534
BTSB1.1591.2502.1931.2312.0660.26
4b1.3021.3321.3611.3741.349
MCE1.1261.2273.4911.1994.047
ETSE1.1191.2461.7171.2614.5470.74
I3E1.1141.2561.5161.3034.5500.90
MCF1.1251.2283.9491.1994.052
FTSF1.1541.1954.5021.2851.7020.55
I4F1.1751.1924.6631.3541.4280.68
1c+2A1c1.1251.228
MCA1.1251.2283.4301.1993.412
TSA1.1491.2542.0381.2322.2640.31
3c1.2891.3491.3441.3741.360
MCB1.1251.2283.6841.2003.326
HBTSB1.1591.2512.1951.2312.0610.27
4c1.3021.3321.3611.3741.349
MCE1.1251.2283.4591.1994.056
ETSE1.1181.2471.7101.2634.5250.76
I3E1.1141.2561.5181.3034.5370.87
MCF1.1251.2293.8861.1993.740
FTSF1.1531.1964.5381.2861.7000.55
I4F1.1731.1934.6651.3541.4280.68
1d+2A1d1.1241.229
MCA1.1251.2293.4451.1993.426
TSA1.1491.2552.0461.2322.2530.29
3d1.2881.3501.3441.3741.360
MCB1.1241.2293.5981.1993.568
BTSB1.1591.2522.1971.2312.0600.24
4d1.3011.3331.3611.3741.350
MCE1.1241.2293.4331.1994.036
ETSE1.1181.2481.7061.2644.5170.75
I3E1.1141.2571.5211.3024.5260.89
MCF1.1241.2303.9491.1994.079
FTSF1.1521.1964.5481.2871.6950.55
I4F1.1731.1944.6631.3551.4260.68
1e+2A1e1.1221.234
MCA1.1221.2343.4191.2003.464
TSA1.1491.2582.0731.2312.2170.27
3e1.2861.3521.3461.3721.361
MCB1.1221.2343.7411.2003.416
BTSB1.1571.2562.2091.2312.0410.22
4e1.2991.3351.3621.3731.351
MCE1.1211.2353.3881.1994.099
ETSE1.1161.2521.6901.2684.4540.75
I3E1.1131.2611.5301.3034.3390.87
MCF1.1211.2353.9781.1994.218
FTSF1.1491.1994.6151.2891.6810.55
I4F1.1691.1974.6651.3551.4240.65
Table 4. ELF valence basin populations of the IRC points, MC 1c3c, defining the twelve different phases characterizing the reaction of the phenyl azide 1c and ethyl propiolate 2. The stationary points 1c, 2, MCA, TSA, and 3c are also included. Distances are given in angstroms, Å, electron populations in an average number of electrons, [e], relative energies in kcal·mol−1, and GEDT values in an average number of electrons, [e].
Table 4. ELF valence basin populations of the IRC points, MC 1c3c, defining the twelve different phases characterizing the reaction of the phenyl azide 1c and ethyl propiolate 2. The stationary points 1c, 2, MCA, TSA, and 3c are also included. Distances are given in angstroms, Å, electron populations in an average number of electrons, [e], relative energies in kcal·mol−1, and GEDT values in an average number of electrons, [e].
Structures1c2MCA P1 P2 P3TSAP4 P5 P6 P7 P8 P9 P10 P11 3c
Phases I II III IV V VI VII VIII IX X XI XII
d2(N1-C5) 3.431 3.024 2.693 2.3352.2632.179 2.108 2.104 1.931 1.910 1.743 1.720 1.477 1.360
d1(N3-C4) 3.137 2.717 2.433 2.1052.0371.953 1.883 1.878 1.718 1.701 1.575 1.561 1.435 1.344
GEDT 0.00 −0.03 −0.06 −0.08−0.07−0.03 0.01 0.02 0.14 0.15 0.24 0.24 0.29 0.25
dE −4.5 0.7 6.9 16.817.516.3 13.2 12.9 −2.4 −4.8 −25.0 −28.6 −55.9 −73.1
V(N1)3.76 3.77 3.76 3.74 3.703.703.72 3.75 3.75 3.45 3.42 3.20 3.19 3.00 2.89
V(N1.N2)2.34 2.32 4.20 3.77 2.962.872.79 2.72 2.72 2.56 2.54 2.41 2.39 2.18 2.04
V′(N1.N2)1.80 1.82
V(N2.N3)2.50 2.51 2.46 2.31 1.891.831.77 1.74 1.73 1.70 1.69 1.70 1.70 1.72 1.71
V(N3)3.37 3.37 3.31 3.26 3.253.273.30 3.48 2.06 1.89 1.87 1.73 1.72 1.43 0.59
V(N3.C7)1.83 1.84 1.87 1.9 1.911.891.88 1.86 1.85 1.83 1.82 1.82 1.82 1.86 2.02
V(C4.C5) 2.562.52 2.48 2.38 2.432.392.31 2.20 2.19 2.04 2.02 1.93 1.92 3.67 3.36
V′(C4.C5) 2.732.75 2.78 2.87 2.562.422.32 2.21 2.21 2.05 2.04 1.98 1.97
V(C5.C6) 2.42.40 2.43 2.46 2.512.452.42 2.4 2.4 2.35 2.35 2.31 2.3 2.25 2.38
V(N1.C5) 1.50 1.93 2.23
V(N3.C4) 1.45 1.78 1.33 2.02 2.04 2.29 3.00
V′(N1) 0.38 0.43 0.73
V(N2) 0.62 1.882.062.22 2.35 2.35 2.59 2.62 2.77 2.79 3.00 3.16
V′(N3) 0.50
V(C4) 0.05
V(C5) 0.140.360.49 0.58 0.59 0.71 0.72 0.75
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dresler, E.; Woliński, P.; Wróblewska, A.; Jasiński, R. On the Question of Zwitterionic Intermediates in the [3+2] Cycloaddition Reactions between Aryl Azides and Ethyl Propiolate. Molecules 2023, 28, 8152. https://doi.org/10.3390/molecules28248152

AMA Style

Dresler E, Woliński P, Wróblewska A, Jasiński R. On the Question of Zwitterionic Intermediates in the [3+2] Cycloaddition Reactions between Aryl Azides and Ethyl Propiolate. Molecules. 2023; 28(24):8152. https://doi.org/10.3390/molecules28248152

Chicago/Turabian Style

Dresler, Ewa, Przemysław Woliński, Aneta Wróblewska, and Radomir Jasiński. 2023. "On the Question of Zwitterionic Intermediates in the [3+2] Cycloaddition Reactions between Aryl Azides and Ethyl Propiolate" Molecules 28, no. 24: 8152. https://doi.org/10.3390/molecules28248152

APA Style

Dresler, E., Woliński, P., Wróblewska, A., & Jasiński, R. (2023). On the Question of Zwitterionic Intermediates in the [3+2] Cycloaddition Reactions between Aryl Azides and Ethyl Propiolate. Molecules, 28(24), 8152. https://doi.org/10.3390/molecules28248152

Article Metrics

Back to TopTop