Copper(I)-Catalyzed 1,3-Dipolar Cycloaddition of Ketonitrones to

Apr 10, 2017 - The metal-involving process is much more asynchronous and profitable from both kinetic and thermodynamic viewpoints than the hypothetic...
0 downloads 10 Views 3MB Size
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.

Article http://pubs.acs.org/journal/acsodf

Copper(I)-Catalyzed 1,3-Dipolar Cycloaddition of Ketonitrones to Dialkylcyanamides: A Step toward Sustainable Generation of 2,3Dihydro-1,2,4-oxadiazoles Anna A. Melekhova, Andrey S. Smirnov, Alexander S. Novikov, Taras L. Panikorovskii, Nadezhda A. Bokach,* and Vadim Yu. Kukushkin* Saint Petersburg State University, 7/9 Universitetskaya Nab., 199034 Saint Petersburg, Russian Federation S Supporting Information *

ABSTRACT: CuI-catalyzed cycloaddition (CA) of the ketonitrones, Ph2CN+(R′)O− (R′ = Me, CH2Ph), to the disubstituted cyanamides, NCNR2 (R = Me2, Et2, (CH2)4, (CH2)5, (CH2)4O, C9H10, (CH2Ph)2, Ph(Me)), gives the corresponding 5-amino-substituted 2,3-dihydro-1,2,4-oxadiazoles (15 examples) in good to moderate yields. The reaction proceeds under mild conditions (CH2Cl2, RT or 45 °C) and requires 10 mol % of [Cu(NCMe)4](BF4) as the catalyst. The somewhat reduced yields are due to the individual properties of 2,3-dihydro-1,2,4-oxadiazoles, which easily undergo ring opening via NO bond splitting. Results of density functional theory calculations reveal that the CA of ketonitrones to CuIbound cyanamides is a concerted process, and the copper-catalyzed reaction is controlled by the predominant contribution of the HOMOdipole−LUMOdipolarophile interaction (group I by Sustmann’s classification). The metal-involving process is much more asynchronous and profitable from both kinetic and thermodynamic viewpoints than the hypothetical metal-free reaction.



INTRODUCTION Although 2,3-dihydro-1,2,4-oxadiazoles1 (DHODs; Scheme 1) are relevant to 1,2,4-oxadiazoles with broad application of the latter in material and medicinal chemistry,2 it is still an almost unexplored class of heterocyclic systems because synthetic routes to DHODs are poorly elaborated. In the context of known DHOD properties, platinum(II) species featuring ligated DHODs are of biological importance, exhibiting antitumor properties.3 In addition, NR2 substituents in 5amino-substituted DHODs could be potentially considered as a useful tool for design of bioisosteres (e.g., carboxylic acid, amide, and esters) based on this type of heterocycles.4 Hence, novel synthetic approaches to DHODs warrant investigation. The first group of reactions leading to DHODs includes metal-free cycloaddition (CA) of nitrones to nitriles (Scheme 1, A), and this reaction proceeds under harsh conditions.1b,c,5 Although this atom-economic method is perhaps the simplest route to obtain DHODs, it is restricted exclusively to electrondeficient nitriles bearing strong acceptor substituents, such as CCl3, and it also utilizes quite reactive nitrones. The reaction of nitriles employing oxaziridines (B), which serve as a hidden and reactive form of nitrones, allows high-yield syntheses of DHODs, but even these reactions involve only rather electron-deficient nitriles, in particular, those with R = Ar.6 CA between oxazol-5-ones and nitrosobenzene gives 3carboxylic-substituted DHODs (R3 = CO2H) (C).7 Another group of reactions leading to DHODs utilizes metal centers to promote CA of nitrones to nitriles (Scheme 1, right). © 2017 American Chemical Society

Metal-mediated processes gradually lose their initial popularity because metals serve as a potential toxic waste source. However, in many instances, metal-involving syntheses (particularly metal-catalyzed) are still preferred over relevant metal-free routes because metal centers could strongly activate reactants and/or substantially reduce the number of steps leading to target compounds. These metal-involving methods were applied for generation of DHODs, and it appeared that some metals act as efficient activators of RCN dipolarophiles in the CA (D).8 The preparation of metal-free DHODs via routes (D)−(E) requires an additional step (E) of ligand liberation.1d,8a−i,9 This method (D−E), despite its generality, employs rather expensive platinum and palladium species. Recently, we have demonstrated (F) that disubstituted cyanamides, NCNRR′, being coordinated to a ZnII center, can be employed in the CA with acyclic N-alkyl ketonitrones Ph2CN+(O−)R′ (R′ = Me, CH2Ph), and this reaction proceeds easily achieving the respective metal-free heterocycles.10 We succeeded in replacing platinum and palladium in the CA with more favorable zinc(II), but the developed method employs stoichiometric rather than catalytic amounts of Zn(OTf)2. Although CA with stoichiometric amounts of cheap zinc is an obvious development in the generation of DHODs, our goal Received: February 6, 2017 Accepted: March 24, 2017 Published: April 10, 2017 1380

DOI: 10.1021/acsomega.7b00130 ACS Omega 2017, 2, 1380−1391

ACS Omega

Article

Taking into account our general interest in metal-mediated reactions of substrates bearing a CN triple bond (for our reviews, see refs 9c and 12) and, in particular, in their metalcatalyzed organic transformations (for recent works, see ref 13), we focused our efforts on the search for a practically useful catalytic system for generation of DHODs. After many unsuccessful attempts, we found that copper(I) exhibits catalytic properties in the CAs of ketonitrones to push−pull nitriles, such as cyanamides, giving 5-amino-substituted DHODs, and all our results disclosing the first catalytic system for generation of DHOD’s are consistently disclosed in sections that follow.

Scheme 1. Metal-Free (Left) and Metal-Mediated (Right) Synthetic Approaches to 2,3-Dihydro-1,2,4-oxadiazoles



RESULTS AND DISCUSSION Copper(I)-Catalyzed CA of Ketonitrones to Cyanamides. Optimization of the Catalytic System and Reaction Conditions. Initially, dimethylcyanamide and N(diphenylmethylene)methanamine oxide, Ph2CN+(Me)O−, were addressed as model substrates for optimization of the reaction conditions. Our variations included choice of copper species, load of catalyst, solvents, time, and reaction temperatures. The results are summarized in Table S1 (see the Supporting Information). We tested copper(I and II) species, namely, [Cu(NCMe)4](BF4) and Cu(OTf)2 (Table S1, entry 15), which are soluble in the used organic solvents. The complex [Cu(NCMe)4](BF4) was found to be the best choice, and no reaction occurred without copper species (Table S1, entries 1, 2, and 14). Variation of load of catalyst (Table S1, entries 3−6) revealed that the highest yield was obtained when 10 mol % of [Cu(NCMe)4](BF4) was employed. The lowering of the yield of 1 upon decreasing the load of catalyst is probably due to decreased concentration of activated substrates. More unusual is the yield reduction upon increasing the load of catalyst above 10 mol %. This phenomenon can be rationalized by the occurrence of some side reactions at a higher catalyst concentration, for instance, CuI-mediated reduction of the ketonitrone. The reaction proceeds in the temperature range from 20 to 80 °C (Table S1, entries 7−11 and 17, in CH2Cl2 or toluene), and increasing the temperature leads to an enhancement of the reaction rate. Effects of concentration (Table S1, entry 12) and microwave irradiation (Table S1, entries 13 and 14) were also studied, but they only slightly affected the reaction rate. The yield was decreased when dimethylcyanamide was taken in 10-fold excess compared with the nitrone (Table S1, entry 16), and the optimal molar ratio between these substrates is 1:1. We believe that a large excess of NCNMe2 could lead to transformation of the copper catalyst to its less active form. Solvent variation was studied for the system Ph2CN+(Me)O−:NCNMe2 in the presence of 10 mol % of the copper(I) catalyst (Table S1, entries 4 and 18−22). The yields of 1 are based on NMR integration with 1,4dimethoxybenzene taken as an internal standard. Only a small increase in the yield of 1 (from 60 to 68%) upon increasing the reaction time from 2 to 24 h (Table S1, entries 4 and 8−11) is probably explained by the consumption of the catalyst in the reaction. In summary, the results of these preliminary tests demonstrate that the optimal conditions for the model reaction include [Cu(NCMe)4](BF4) (10 mol %), molar ratio between the reactants 1:1, temperature 45 °C and application of CH2Cl2 as the solvent. The order of addition of the reactants is important, and the highest yield (68%; 45 °C, 24 h) of 1 was achieved when NCNMe2 was added to a solution of

was to find a catalytic system for the CA and synthesis of these heterocycles. At first glance, this task may seem simple, but it proves surprisingly challenging as the dominant part of useful catalytic systems utilizes Pearson’s “hard” metal centers. These centers, in turn, preferably coordinate with the hard oxygen center of nitrones, thus blocking these dipoles toward the CA, and also promote nitrone hydrolytic and/or deoxygenation decomposition (Scheme 2).8d,11 Scheme 2. Metal-Mediated Reactions of Nitrones

1381

DOI: 10.1021/acsomega.7b00130 ACS Omega 2017, 2, 1380−1391

ACS Omega

Article

reaction mixture after 2 days. This observation additionally demonstrates the greater reactivity of cyanamides toward CA as compared to that of conventional nitriles; this trend was highlighted previously in our review.1a Aldonitrone, 4MeC 6H 4 CHN +(Me)O −, did not react with NCNMe 2 under the reaction conditions (10 mol % [Cu(NCMe)4](BF4), CH2Cl2, 45 °C, 2 days), which additionally confirmed the previously established greater reactivity of C,C-diaryl ketonitrones than that of C-aryl aldonitrones.8b The higher reactivity of ketonitrones is related to the different electronic effect of phenyl groups in aldo- and ketonitrones. In aldonitrones, the aromatic substituent at the C atom is involved in the conjugation and acts as an electron-acceptor substituent, decreasing the reactivity of the 1,3-dipoles. In C,C-diaryl ketonitrones, both aromatic substituents are out of the CNO plane and are not involved in the conjugation and therefore, could not act as electron acceptors. 2,3-Dihydro-1,2,4-oxadiazoles are unstable in the presence of copper(I), and heating of 1 with [Cu(NCMe)4](BF4) (50 mol %; nitromethane, 70 °C, 2 h) leads to almost complete decomposition of the heterocycle and formation of a mixture of benzophenone (isolated yield 25%) and 3-(diphenylmethylene)-1,1-dimethylurea (isolated yield 56%) (Scheme 4); these

[Cu(NCMe) 4 ](BF 4 ) followed by addition of Ph 2 C N+(Me)O−. If the ketonitrone was added to a solution of the catalyst followed by addition of dimethylcyanamide, the yield of the target heterocycle was substantially lower (ca. 10%). This behavior could be explained taking into account the competition of the ketonitrone and NCNMe2 for the copper(I) coordination site. We studied the reaction of ketonitrones with CuI in a separate experiment and found that their interaction leads to oxidation of copper(I) and coordination of the nitrone O atom followed by decomposition of the complexes thus formed (see later). Hence, a combination of coordinated dimethylcyanamide and free dipole leads to the CA;1a,14 therefore, the order of reactant addition indicated above is preferable. Our theoretical calculations (see later) agree with the experimental data and they support group I by Sustmann’s classification15 of CA. Substrate Variation. As the next step, we extended the scope of the substrate to various disubstituted cyanamides and two aryl ketonitrones (Scheme 3), and the obtained results are summarized in Table 1. Scheme 3. Copper(I)-Catalyzed CA of Ketonitrones to Cyanamides

Scheme 4. Copper(I)-Mediated Decomposition of the Heterocycle

For most combinations of cyanamides and ketonitrones, the corresponding 2,3-dihydro-1,2,4-oxadiazoles were obtained and isolated in 28−73% yields. In all reaction mixtures, unreacted ketonitrone was identified, and this reflects an incomplete conversion of substrates. The target product was not isolated for the NCNR2/Ph2CN+(CH2Ph)O− (R = Me, Et) systems, whereas signals corresponding to 4 were determined by HRESI+-MS monitoring of the reaction mixtures. Product 4 was not obtained, probably due to steric restriction of the NEt2 group that prevents cyanamide−nitrone CA. When the reaction with conventional nitriles RCN (R = Me, Ph, 4-Cl-C6H4, 4-CF3C6H4) was attempted under the same conditions (Ph2CN+(Me)O−, 10 mol % [Cu(NCMe)4](BF4), CH2Cl2, 45 °C), no CA products were detected in the

species are probably derived from the known10 2,3-dihydo1,2,4-oxadiazole ring opening. Heterocycle 1 remains intact under the same conditions in the absence of the copper(I) complex. The observed metal-mediated ring opening of 2,3-dihydro1,2,4-oxadiazoles explains the moderate isolated yields of these heterocycles. Characterization of 6 and 8−16. Heterocycles 1−3, 5, and 7 were identified by comparison of their 1H NMR and high-resolution ESI+-MS spectra with those known from the

Table 1. Compound Numbering Scheme and Isolated Yields (%) of Heterocycles R′ of nitrone/reaction conditions R2 of NCNR2 Me2 Et2 (CH2)5 (CH2)4O (CH2)4 tetrahydroisoquinolin-2-yl (C9H10) (CH2Ph)2 Ph(Me)

Me 1 3 5 7 9 11 13 15

(68) (30) (60) (73) (57) (45) (68) (45)

reaction conditions

CH2Ph

2 h, 45 °C 2 h, 45 °C 2 h, 45 °C 2 h, 45 °C 2 h, 45 °C reaction time 24 h, RT 2 h, 45 °C 2 h, 45 °C

2 (30) 4 tracesa 6 (50) 8 (55) 10 (34) 12 (28) 14 (53) 16 (50)

reaction conditions 2 h, 45 °C 2 h, 45 °C 2 h, 45 °C 2 h, 45 °C reaction time 24 h, RT 24 h, 45 °C 48 h, 45 °C

a We were unable to isolate the cycloadduct when the reaction was carried out under standard conditions (2 h, 45 °C); upon prolonged treatment (24 h, 45 °C) under harsh conditions (MW irradiation, 100 °C, 30 min), we only observed peaks corresponding to the cycloadduct in HRESI+-MS of the reaction mixture.

1382

DOI: 10.1021/acsomega.7b00130 ACS Omega 2017, 2, 1380−1391

ACS Omega

Article

Figure 1. Molecular structure of 12 with the atomic numbering scheme.

Figure 2. Molecular structure of [Cu{ON(R′′)CPh2}4](BF4)2 (17) with the atomic numbering scheme. Thermal ellipsoids are drawn at the 50% probability level. Selected bond lengths (Å) and angles (deg): Cu1O1 1.9333(9), Cu1O1′ 1.9286(9), O1N1 1.3439(14), O1′N1′ 1.3509(14), N1C2 1.2987(17), N1′C2′ 1.3021(16), O1CuO1′ 90.54(4), Cu1O1N1 118.10(7), O1N1C2 123.00(11).

literature.10 New heterocycles 6, 8−16 were characterized by elemental analyses (C, H, N), HRESI+-MS, 1H and 13C{1H} NMR, Fourier transform infrared (FTIR), and also by X-ray diffraction (for 9, 12, and 16). Heterocycle 4 was not isolated from the reaction mixture due to low yield and was identified by HRESI+-MS in the reaction mixtures (m/z: 386.2246 [M + H]+, calcd 386.2227). Compounds 6 and 8−16 give satisfactory

C, H, and N elemental analyses for the proposed formulas. The HRESI+-MS of these species exhibit peaks corresponding to [M + H]+. The FTIR spectra of all complexes display the ν(CN) bands in the range 1655−1660 cm−1, which is specific for relevant DHODs.8c,16 The protons of the N-methyl group of 9, 11, 13, and 15 appeared as a singlet in the interval 2.33−2.49 ppm; the methylene protons (NCH2Ph) of 6, 8, 10, 12, 14, and 1383

DOI: 10.1021/acsomega.7b00130 ACS Omega 2017, 2, 1380−1391

ACS Omega

Article

Scheme 5. Proposed Mechanisms for the CA

that were studied by X-ray diffraction (17; Figure 2; for description of the X-ray structure, see the SI). Moreover, we also found that the nitrones are unstable in the presence of a catalytic amount of the copper(I) complex. The heating of Ph2CN+(O−)Me with [Cu(NCMe)4](BF4) (10 mol %) under optimized conditions (45 °C, 2 h) leads to a partial hydrolytic decomposition of the dipole giving benzophenone Ph2CO (preparative yield 22%). Thus, upon treatment of ketonitrones with CuI we observed that the reactant interplay leads to oxidation of copper(I) and a broad spectrum of products. In the context of the copper(I)catalyzed CA, these synthetic experiments give an idea that the generation of 2,3-dihydro-1,2,4-oxadiazoles is likely to proceed via ligation of the nitriles to the CuI center, followed by the reaction of the metal-bound dipolarophiles with the nitrones. This assumption is fully supported by theoretical calculations whose results are presented below. Theoretical Study of the Cu I -Catalyzed CA of Ketonitrones to Cyanamides. To understand the mechanism of the copper-catalyzed CA of ketonitrones to cyanamides we have carried out quantum-chemical calculations at the density functional theory (DFT) level of theory of model CA of Me2 CN +(Me)O− to uncomplexed and copper-bound NCNMe2. We have successfully used this approach in studies of similar metal-assisted and metal-free 1,3-dipolar CA reactions of nitrones to isocyanides20 and nitriles,21 and metal-assisted coupling of oximes and nitriles.22 We proposed three types of mechanisms for CA of Me2C N+(Me)O− to NCNMe2, two of them are stepwise (Scheme 5 A,C) and the third is concerted (B). CA is initiated by the endoergonic formation of the orientation complex OC (+4.2 and +2.3 kcal/mol for metal-free and CuI-catalyzed reactions, respectively). Stepwise pathways involve the formation of acyclic zwitterionic intermediates INTa or INTb (via TSs1a or TSs1b), which then undergo the cyclization (via TSs2a or TSs2b, respectively) giving the cyclic product P, whereas the

16 appear as a singlet located between 3.37 and 3.64 ppm. The 13 C{1H} NMR spectra display resonances in the interval 94.7− 96.0 ppm corresponding to the quaternary C3 atoms of the heterocycles and signals at 158.0−160.6 ppm attributed to the C5 atom from the CN moiety. In the 13C{1H} NMR spectra, the peaks from the CPh2 moiety emerge as one broad singlet due to a dynamic process. Heterocycles 9, 12, and 16 were characterized by X-ray crystallography (Figures 1, S1, and S2, see the Supporting Information). Bond distances and angles for these structures are close to those in the previously reported structures of DHODs (Table S2).10 Generation of Nitrone Complexes [Cu{ON(R′′)CPh2}4](BF4)2 (17). If the reaction of copper(I) centers with nitriles are well known,17 those of nitrones are poorly studied. Copper(II) nitrone complexes documented in the literature feature bidentately coordinated ligands with the nitrone moiety18 and the monodentate 2,5,5-trimethyl-1-pyrroline-N-oxide.18c For a deeper understanding of the copper(I)-catalyzed CA, we studied the reaction between the nitrone dipoles Ph2C N+(O−)R″ (R″ = Me, CH2Ph) and [Cu(NCMe)4](BF4). Four equivalents of any one of the nitrones were added to [Cu(NCMe)4](BF4) in CH2Cl2 and the reaction mixture was stirred at RT. HRESI+-MS monitoring of the reaction mixture after 24 h allows the identification of peaks corresponding to the fragment ions, namely, [MCl(nitrone) 2 ] + , [MCl(nitrone)3]+, and [M(nitrone)2]+. IR monitoring of the reaction mixture allowed the observation of the CN stretching bands of the (nitrone)CuI complexes at ca. 1660 cm−1 (R″ = Me, Ph), whereas the uncomplexed nitrones Ph2CN+(O−)R″ (R = Me, Ph) exhibit a medium-to-weak band at ca. 1530−1520 cm−1 of the CCAr and/or CN stretches.19 The resulting mixture was diluted with hexane and concentrated in vacuo at RT giving a dark greenish oily residue that, on the next day, solidified furnishing a few black needle-like crystals (R″ = Me) 1384

DOI: 10.1021/acsomega.7b00130 ACS Omega 2017, 2, 1380−1391

ACS Omega

Article

concerted CA proceeds via a cyclic transition state (TSc). For both metal-free and metal-involving processes no minima for structures corresponding to any acyclic intermediates were located upon a detailed search for the potential energy surface. Despite numerous attempts, in all cases, optimization led to the formation of initial species, which are separated from each other, or to product P. Thus, stepwise pathways have been excluded from consideration. On the other hand, the TSc structures, corresponding to concerted CA, were successfully found for both metal-free and metal-involving reactions (TScfree and TScCu, respectively). The cyclic nature of these transition states (TSs) has been confirmed by the topological analysis of the electron density distribution within the formalism of the quantum theory of atoms in molecules (QTAIM) method23 (we successfully used this approach in studies of noncovalent interactions and properties of coordination bonds in various transition metal complexes20,24). The Poincare−Hopf relation in both cases is satisfied, thus all critical points have been found. The contour line diagrams of the Laplacian distribution ∇2ρ(r), bond paths, and selected zero-flux surfaces are shown in Figure S48. The ring critical point (3, +1) for the 2,3-dihydro-1,2,4-oxadiazole cycle and five suitable bond critical points (3, −1; BCPs) were determined. For TScCu, the electron density value ρ(r) at the O···Ccyanamide BCP is significantly higher than that at the Cnitrone···Ncyanamide BCP, but for TScfree, these parameters at both BCPs are comparable (Table S4). The energy density Hb at the O··· Ccyanamide and Cnitrone···Ncyanamide BCPs in TScfree and at O··· Ccyanamide BCP in TScCu is clearly negative demonstrating some covalent contribution in these interactions, whereas Hb at the Cnitrone···Ncyanamide BCP in TScCu is virtually zero and hence, no covalent contribution in this contact was detected. These observations reveal that the concerted metal-involving CA reaction is significantly more asynchronous as compared to the metal-free process, which is consistent with the general trend: coordination of a Lewis acid to any reactant leads to substantial decreasing of the reaction synchronicity [CA of nitrones to the CN bond: refs 1a, 20, 21, 24d, 25; CA of nitrones to the C C bond: ref 26]. Indeed, appropriate S y parameters (quantitative measure of the synchronicity of concerted CAs proposed by Moyano et al.27) are 0.68 for metal-involving CA versus 0.93 for metal-free CA. Inspection of the calculated activation and reaction energies (Figure 3, Table S6) indicates, first, that the activation barrier for the Cu-catalyzed CA is significantly lower than that for the metal-free coupling (by 5.4 kcal/mol in terms of Gibbs free energies in solution). Second, metal-involving process is more preferably from a thermodynamic viewpoint. An interaction between the frontier MOs of the reactants is the driving force of CA reactions. Our calculations indicate that the CA reaction of Me2CN+(Me)O− to free NCNMe2 belongs to group II processes in the Sustman classification (energy gaps between HOMOdipole−LUMOdipolarophile and HOMOdipolarophile−LUMOdipole are 0.26 and 0.25 au, respectively, see Figure 4). The coordination of NCNMe2 to the copper(I) metal center results in a decrease of both HOMOπ(CN) and LUMOπ*(CN) energies and leads to significant contraction of the HOMOdipole−LUMOdipolarophile energy gap and expansion of the HOMOdipolarophile−LUMOdipole energy gap. Thus, the simple qualitative MO consideration suggests that (i) the coordination of NCNMe2 should accelerate the nitrone CA and (ii) the

Figure 3. Energy profiles for metal-free (blue) and metal-involving (red) processes.

metal-involving CA belongs to normal electron demand processes (group I by Sustmann’s classification). Literature data confirm our observations and, indeed, CAs of nitrones often belong to group II reactions,28 but it is easily possible to change the type of reaction by the variation of substituents in the reactants or its coordination to the metal center. Thus, CA of nitrones to electron-deficient dipolarophiles (namely, acrylonitrile,29 acrolein,30 methyl propiolate,31 methyl acrylate,32 nitroalkenes,33 fluoroalkenes, and fluoroalkynes34) as well as metal-involving processes35 is determined by the interaction of HOMOdipole−LUMOdipolarophile (group I, normal electron demand reactions), whereas usage of dipolarophiles with electron-donating substituents (e.g., methyl vinyl ether)26a leads to switching of the type of reaction to inverse electron demand (HOMOdipolarophile−LUMOdipole, group III by Sustmann’s classification).



CONCLUDING REMARKS We found a novel synthetic approach to 2,3-dihydro-1,2,4oxadiazoles that includes CuI-catalyzed CA of ketonitrones to disubstituted cyanamides and gives 5-amino-substituted ring systems. The reaction proceeds under mild conditions (RT or 45 °C) and requires 10 mol % of [Cu(NCMe)4](BF4) as the catalyst. The application of the developed method furnishes the heterocycles in good to moderate yields. The somewhat reduced yields are due to the individual properties of 2,3dihydro-1,2,4-oxadiazoles, which easily undergo ring opening via the NO bond splitting. All previous approaches to 2,3-dihydro-1,2,4-oxadiazoles bearing donor substituents at the fifth position of the ring were based on the usage of stoichiometric amounts of activating metal centers (ZnII, PtII, PdII, or PtIV) and these methods included a two-step procedure followed by (for platinum and palladium) conventional recycling of these expensive metals. The discovery of a simple catalytic system for generation of 2,3dihydro-1,2,4-oxadiazoles is certainly a step toward sustainable synthesis of this yet unexplored class of heterocycles. Inspection of our experimental data indicates that the combination of coordinated cyanamide dipolarophile and uncomplexed nitrone dipole, not vice versa, leads to the CA. Theoretical calculations are agreeable with this conclusion. Results of the DFT calculations reveal that the CA of ketonitrones to CuI-bound cyanamides is a concerted process 1385

DOI: 10.1021/acsomega.7b00130 ACS Omega 2017, 2, 1380−1391

ACS Omega

Article

Figure 4. Energies of the interacting frontier MOs of Me2CN+(Me)O− and uncomplexed or Cu-bound NCNMe2.

the CrysAlisPro41 program complex using spherical harmonics, implemented in the SCALE3 ABSPACK scaling algorithm. Supporting crystallographic data for this article have been deposited at the Cambridge Crystallographic Data Centre (CCDC 1456546, 1456547, 1470262, 1468157) and can be obtained free of charge via www.ccdc.cam.ac.uk/data_request/ cif. Computational Details. The full geometry optimization of all structures and TSs has been carried out at the DFT level of theory using the M06 functional42 with the help of the Gaussian-0943 program package. No symmetry operations have been applied. The calculations were carried out using the multielectron fit fully relativistic energy-consistent pseudopotential MDF10 of the Stuttgart/Cologne group that described 10 core electrons and the appropriate contracted basis set for the copper atom44 and the 6-31G(d) basis set for other atoms. The Hessian matrix was calculated analytically for the optimized structures to prove the location of correct minima (no imaginary frequencies) or saddle points (only one imaginary frequency) and to estimate the thermodynamic parameters, the latter being calculated at 25 °C. The nature of all TSs was studied by the analysis of vectors associated with the imaginary frequency and by calculations of the intrinsic reaction coordinates using the Gonzalez−Schlegel method.45 The total energies corrected for solvent effects (Es) were estimated by single-point calculations on the basis of equilibrium gas-phase geometries at the same level of theory using the SMD continuum solvation model by Truhlar and coworkers46 with CH2Cl2 as the solvent. The entropic term in CH2Cl2 solution (Ss) was calculated according to the procedure described by Wertz47 and Cooper and Ziegler48 using eqs 1−4

and the copper-catalyzed reaction is controlled by the predominant contribution of the HOMOdipole−LUMOdipolarophile interaction (group I by Sustmann’s classification). The metalinvolving process is much more asynchronous and favorable from both kinetic and thermodynamic viewpoints than the hypothetical metal-free reaction.



EXPERIMENTAL SECTION Materials and Instrumentation. The dialkylcyanamides, NCNR2 (R = Me, Et, 1/2(CH2)5, 1/2(CH2)4O, 1/2(CH2)4; Aldrich), and solvents were obtained from a commercial source and used as received. The dialkylcyanamides NCN(CH2Ph)2, NCN(Me)Ph, and NCNC9H10,36 the copper(I) complex [Cu(NCMe)4](BF4),37 and the nitrones38 were synthesized in accord with the published recipes. The HRESI mass spectra were obtained on a Bruker micrOTOF spectrometer equipped with an electrospray ionization source, and MeOH was employed as the solvent. The instrument was operated in positive ion mode using an m/z range of 50−3000. The capillary voltage of the ion source was set at −4500 V (ESI+ MS) and the capillary exit was set at ±(70−150) V. In the isotopic pattern, the most intensive peak is reported. Infrared spectra were recorded using a Bruker FTIR TENSOR 27 instrument in KBr pellets. 1H and 13C{1H} NMR spectra were measured using a Bruker Avance III 400/100 MHz spectrometer at ambient temperature. X-ray Structure Determinations. Crystals of 9, 12, 16, and 17 were measured on a Agilent Technologies SuperNova diffractometer at a temperature of 100 K using monochromated Cu Kα radiation. All structures were solved by direct methods by means of the SHELX program39 incorporated into the OLEX2 program package.40 For crystallographic data and refinement parameters, see Table S3. The carbon-bound H atoms were placed in calculated positions and were included in the refinement in the “riding” model approximation, with Uiso(H) set to 1.5Ueq(C) and CH 0.98 Å for CH3 groups, with Uiso(H) set to 1.2Ueq(C) and CH 0.99 Å for CH2 groups, and with Uiso(H) set to 1.2Ueq(C) and CH 0.95 Å for CH groups. Empirical absorption correction was applied in

ΔS1 = R ln V s m,liq /Vm,gas

(1)

ΔS2 = R ln V o m/V s m,liq

(2)

α= 1386

S o,s liq − (S o,s gas + R ln V s m,liq /Vm,gas) (S o,s gas + R ln V s m,liq /Vm,gas)

(3)

DOI: 10.1021/acsomega.7b00130 ACS Omega 2017, 2, 1380−1391

ACS Omega

Article

from further fractions eluted by EtOAc/hexane 1:2 → pure EtOAc (5−15% recovered). Previously reported1 heterocycles 1−3, 5, and 7 were identified by comparison of their 1H NMR, HRESI+-MS, and TLCs with those of the known species. Full characterization of novel 2,3-dihydro-1,2,4-oxadiazoles is given below. 2-Benzyl-3,3-diphenyl-5-(piperidin-1-yl)-2,3-dihydro-1,2,4oxadiazole (6). Yield 50%. Mp 128−132 °C. Anal. Calcd for C26H27N3O (%): C, 78.56; H, 6.85; N, 10.57. Found (%): C, 78.74; H, 6.90; N 10.61. HRESI+-MS, m/z: 398.2212 ([M + H]+, calcd 398.2227). IR spectrum, selected bands, cm−1: ν(CH), ν(NCN). 1H NMR in CDCl3 δ: 1.62 (m, 6H, piperidine), 3.35 (m, 4H, piperidine), 3.48 (s, 2H, CH2, benzyl), 7.26−7.38 (m, 11H, Ph), 7.69 (br, 4H, Ph). 13C{1H} NMR in CDCl3, δ: 24.3 (CH2CH2CH2), 25.4 (CH2CH2CH2), 47.5 (CH2NCH2), 58.6 (CH2, benzyl), 95.34 (C3), 127.1, 128 br, 129.0, 138.4 (Ph), 160.1 (C5N). 4-(2-Benzyl-3,3-diphenyl-2,3-dihydro-1,2,4-oxadiazol-5yl)morpholine (8). Yield 55%. Mp 137−140 °C. Anal. Calcd for C25H25N3O2 (%): C, 75.16; H, 6.31; N, 10.52. Found (%): C, 75.14; H, 6.32, N 10.48. HRESI+-MS, m/z: 420.2033 ([M + H]+, calcd 400.2020). IR spectrum in KBr, selected bands, cm−1: ν(CH), ν(NCN). 1H NMR in CDCl3 δ: 3.37 (m, 4H, morpholine), 3.47 (s, 2H, benzyl), 3.65 (m, 4H, morpholine), 7.24−7.36 (m, 11H, Ph), 7.65 (br, 4H, Ph). 1H NMR in DMSO-d6 δ: 3.26 (m, 4H, morpholine), 3.37 (br s, 2H, benzyl), 3.57 (m, 4H, morpholine), 7.25 (m, 11H, Ph), 7.56, 7.58 (two br, s, 4H, Ph). 13C{1H} NMR in DMSO-d6, δ: 46.1 (CH2NCH2), 57.7 (CH2), 65.3 (CH2OCH2), 94.7 (C3), 127.1, 127.3 br, 128.0, 128.2, 128.9, 137.6 (Ph), 159.4 (C5 N). 2-Methyl-3,3-diphenyl-5-(pyrrolidin-1-yl)-2,3-dihydro1,2,4-oxadiazole (9). Yield 57%. Mp 148−150 °C. Anal. Calcd for C19H21N3O (%): C, 74.24; H, 6.89; N, 13.67. Found (%): C 74.77; H 6.87; N 13.68. HRESI+-MS, m/z: 308.1760 ([M + H]+, calcd 308.1758). IR spectrum in KBr, selected bands, cm−1: 3060, 2980, 2870 m, ν(CH), 1660 s, ν(NCN). 1 H NMR in CDCl3 δ: 1.87 (m, 4H, pyrrolidine), 2.33 (s, 3H, Me), 3.42 (m, br, 4H, pyrrolidine), 7.20−7.31 (m, 6H, Ph), 7.61 (br, 4H, Ph). 13C{1H} NMR in CDCl3, 25.6 (CH2), 42.4 (NMe), 47.6 (NCH2), 95.8 (C3), 127.2, 128.0 (br; Ph), 158.0 (C5N). 2-Benzyl-3,3-diphenyl-5-(pyrrolidin-1-yl)-2,3-dihydro1,2,4-oxadiazole (10). Yield 34%. Mp 166−168 °C. Anal. Calcd for C25H25N3O (%): C, 78.30; H, 6.57; N, 10.96. Found (%): C, 78.90; H 6.25; N 11.46. HRESI+-MS, m/z: 384.2067 ([M + H]+, calcd 384.2071). IR spectrum in KBr, selected bands, cm−1: 3055, 3025, 2975, 2870 m, ν(CH), 1660 s, ν(NCN). 1 H NMR in CDCl 3 δ: 1.86 (m, 4H, pyrrolidine), 3.36 (m, 4H, pyrrolidine), 3.47 (s, 2H, CH2, benzyl), 7.22−7.36 (m, 6H, Ph), 7.67 (br, 4H, Ph,). 13C{1H} NMR in CDCl3, δ: 25.6 (CH2 pyrrolidine), 47.6 (NCH2, pyrrolidine), 58.6 (NCH2, benzyl), 95.7 (C3), 127.1, 127.4 (br), 128.3, 129.0, 138.4 (Ph), 158.3 (C5N). 5-(3,4-Dihydroisoquinolin-2(1H)-yl)-2-methyl-3,3-diphenyl-2,3-dihydro-1,2,4-oxadiazole (11). Yield 45%. Mp 145−148 °C (dec). Anal. Calcd for C24H23N3O (%): C, 78.02; H, 6.27; N, 11.37. Found (%): C, 78.53; H, 6.30; N, 11.59. HRESI+-MS, m/z: 370.1902 ([M + H]+, calcd 370.1914). IR spectrum in KBr, selected bands, cm−1: 3060, 3000, 2915, 2855 m, ν(C H), 1655 s, ν(NCN). 1H NMR in CDCl3, δ: 2.35 (s, 3H, Me), 2.89 (br, 2H, NCH2CH2), 3.69 (br, 2H, NCH2CH2), 4.66 (s, 2H, CH2). 7.12−7.33 (m, Ph, 10H), 7.62 (br, Ph, 4H).

Ss = Sg + ΔSsol = Sg + [ΔS1 + α(Sg + ΔS1) + ΔS2] = Sg + [(− 11.80 cal/mol K) − 0.21(Sg − 11.80 cal/mol K) + 5.45 cal/mol K]

(4)

where Sg is the gas-phase entropy of the solute, ΔSsol is the solvation entropy, S°,sliq, S°,sgas, and Vsm,liq are the standard entropies and molar volume of the solvent in liquid or gas phases (173.84 and 270.28 J/mol K and 64.15 mL/mol, respectively, for CH2Cl2), Vm,gas is the molar volume of ideal gas at 25 °C (24 450 mL/mol), V°m is the molar volume of the solution corresponding to the standard conditions (1000 mL/ mol). The enthalpies and Gibbs free energies in solution (Hs and Gs) were estimated by using expressions 5 and 6 Hs = Es − Eg + Hg

(5)

Gs = Hs − TSs

(6)

where Es, Eg, and Hg are the total energies in solution and in the gas phase and gas-phase enthalpy, respectively. A topological analysis of the electron density distribution with the help of the QTAIM formalism developed by Bader23 has been performed using the Multiwfn program (version 3.3.4).49 The Wiberg bond indices were computed using the natural bond orbital partitioning scheme.50 The synchronicity of concerted CAs (Sy)27 was calculated using eqs 7−9 n

Sy = 1−(2n − 2)−1 ·∑ i=1

|δBi − δBav | δBav

(7)

where n is the number of bonds directly involved in the reaction (n = 5 for 1,3-dipolar CAs) and δBi is the relative variation of a given Wiberg bond index Bi at the TS relative to reactants (R) and products (P), and it is calculated as

δBi =

BiTS − BiR BiP − BiR

(8)

If the δBi value is negative, it is assumed to be zero. The average value of δBi (δBav) is defined as n

δBav = n−1∑ δBi i=1

(9)

Sy is 0 for stepwise CAs and is 1 for the fully concerted reactions. The Cartesian atomic coordinates of the calculated equilibrium structures are presented in the Table S7. Synthetic Work. CuI-Catalyzed Synthesis of 2,3-Dihydro1,2,4-oxadiazoles. A solution of any one of NCNR2 (0.473 mmol) in CH2Cl2 (2 mL) was added to solid [Cu(NCMe)4](BF4) (15 mg, 0.047 mmol), and the formed mixture was stirred for 5 min for homogenization, whereupon Ph2C N+(O−)R′ (0.473 mmol; R = Me, 100 mg; CH2Ph 0.135 mg) was added, and the formed reaction mixture was stirred under conditions given in Table S1. The reaction mixture was subjected to chromatographic separation on silica gel. Heterocycles 1−16 were isolated from the first fraction (the fractions were identified by TLC on Merck 60 F254 plates under UV light; hexane/EtOAc 10:1, v/v); unreacted nitrone was isolated 1387

DOI: 10.1021/acsomega.7b00130 ACS Omega 2017, 2, 1380−1391

ACS Omega



C{1H} NMR in CDCl3, δ: 28.7 (NCH2CH2), 42.4 (Me), 43.6 (NCH2CH2), 47.8 (NCH2), 95.7 (C3), 126.4, 126.5, 126.6, 127.4 (br), 128.0, 128.9, 133.1, 134.3 (Ph), 159.4 (C5N). 2-Benzyl-5-(3,4-dihydroisoquinolin-2(1H)-yl)-3,3-diphenyl2,3-dihydro-1,2,4-oxadiazole (12). Yield 28%. Mp 122−125 °C. Anal. Calcd for C30H27N3O (%): C, 80.87; H, 6.11; N, 9.43. Found (%): C, 80.41; H, 6.30; N, 9.86. HRESI+-MS, m/z: 446.2227 ([M + H]+, calcd 446.2227). IR spectrum in KBr, selected bands, cm−1: 3080, 3060, 3025, 2925, 2895, 2845 ν(CH), 1660 s ν(NCN). 1H NMR in CDCl3, δ: 2.88 (br, 2H, NCH2CH2), 3.48 (br, 2H, CH2, benzyl), 3.62 (br, 2H, NCH2CH2), 4.61 (s, 2H, NCH2), 7.10−7.36 (m, Ph, 15H), 7.71 (br, Ph, 4H). 13C{1H} NMR in CDCl3, δ: 28.6 (NCH2CH2), 43.6 (NCH2CH2), 47.8 (NCH2), 58.7 (CH2, benzyl) 95.6 (C3), 126.4, 126.4, 126.7, 127.2, 127.5 (br), 128.0, 128.3, 128.9, 129.0, 133.1, 134.3, 138.2 (Ph), 159.6 (C5N). N,N-Dibenzyl-2-methyl-3,3-diphenyl-2,3-dihydro-1,2,4-oxadiazol-5-amine (13). Yield 68%. Mp 94−96 °C. Anal. Calcd for C29H27N3O (%): C, 80.34; H, 6.28; N, 9.69. Found (%): C, 80.14; H, 6.48; N, 9.81. HRESI+-MS, m/z: 434.2245 ([M + H]+, calcd 434.2227). IR spectrum in KBr, selected bands, cm−1: 3085, 3060, 3028, 2910, 2860 ν(CH), 1655 ν(N CN). 1H NMR in CDCl3, δ: 2.49 (Me, 3H), 4.49 (br, 2H, CH2), 4.60 (br, 2H, CH2), 7.26−7.33 (m, Ph, 4H), 7.32−7.44 (m, Ph, 12H), 7.78 (br, Ph, 4H). 13C{1H} NMR in CDCl3, δ: 42.1 (NMe), 51.2 (CH2, benzyl), 96.0 (C3), 127.9, 128.0 br, 128.4, 128.6, 137.21 (Ph), 160.3 (C5N). N,N,2-Tribenzyl-3,3-diphenyl-2,3-dihydro-1,2,4-oxadiazol5-amine (14). Yield 53%. Mp 119−122 °C. Anal. Calcd for C35H31N3O (%): C, 82.48; H, 6.13; N, 8.25. Found (%): C, 82.41, H 6.22, N 8.33. HRESI+-MS, m/z: 509.2523 ([M + H]+, calcd 510.2540). IR spectrum in KBr, selected bands, cm−1: 3085, 3060, 3030, 2920, 2870 sh ν(CH), 1655 s ν(NC N). 1H NMR in CDCl3, δ: 3.64 (br, 2H, N5CH2), 4.45 (br, 4H, CH2), 7.16−7.18 (m, 4H, Ph), 7.27−7.36 (m, 13H, Ph), 7.39− 7.43 (m, 4H, Ph). 13C{1H} NMR in CDCl3, δ: 51.1 (CH2, benzyl), 58.4 (N5CH2), 95.8 (C3), 127.0, 127.5, 127.7 br, 128.0, 128.2, 128.6, 128.7 (Ph), 160.6 (C5N). N,2-Dimethyl-N,3,3-triphenyl-2,3-dihydro-1,2,4-oxadiazol5-amine (15). Yield 45%. Mp 131−133 °C. Anal. Calcd for C22H21N3O (%): C, 76.94; H, 6.16; N, 12.24. Found (%): C, 77.40; H, 6.32; N, 12.61. HRESI+-MS, m/z: 344.1759 ([M + H]+, calcd 344.1758). IR spectrum in KBr, selected bands, cm−1: 3060, 2920, 2850 ν(CH), 1660 ν(NCN). 1H NMR in CDCl3, δ: 2.37 (s, 3H, Me), 3.44 (s, 3H, Me), 7.09− 7.32 (m, Ph, 11H), 7.61 (br, Ph, 4H). 13C{1H} NMR in CDCl3, δ: 39.7 (Me), 42.2 (Me), 95.8 (C3), 123.9, 125.1, 127.4 br, 128.9, 143.7 (Ph), 158.7 (C5N). 2-Benzyl-N-methyl-N,3,3-triphenyl-2,3-dihydro-1,2,4-oxadiazol-5-amine (16). Yield 50%. Mp 150−152 °C. Anal. Calcd for C28H25N3O (%): C, 80.16; H, 6.01; N, 10.02. Found (%): C, 80.43; H, 6.25; N, 10.41. HRESI+-MS, m/z: 420.2072 ([M + H]+ calcd 420.2071). IR spectrum in KBr, selected bands, cm−1: 3060, 3025, 2950, 2890 ν(CH), 1655 ν(NCN). 1 H NMR in CDCl3, δ: 3.41 (s, 3H, Me), 3.55 (s, br, 2H, CH2), 7.07−7.35 (m, 16H, Ph), 7.70 (br, 4H, Ph). 13C{1H} NMR in CDCl3, δ: 39.2 (Me), 58.5 (CH2), 95.7 (C3), 123.4, 124.9, 127.1, 127.6, 128.1, 128.2, 128.7 br, 128.8, 137.9, 143.5 (Ph), 158.9 (C5N). 13

Article

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsomega.7b00130. Experimental procedures, spectral and crystallographic data, computational details, and Cartesian coordinates for all optimized structures (PDF) Crystallographic data for 9, 12, 16, and Ph2CN(O)Me (CIF)(CIF)(CIF)(CIF)(CIF)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected] (N.A.B.). *E-mail: [email protected] (V.Y.K.). ORCID

Anna A. Melekhova: 0000-0002-1066-2177 Alexander S. Novikov: 0000-0001-9913-5324 Taras L. Panikorovskii: 0000-0002-2323-1413 Nadezhda A. Bokach: 0000-0001-8692-9627 Vadim Yu. Kukushkin: 0000-0002-2253-085X Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors are grateful to the Russian Foundation for Basic Research (grants 14-03-00080-a and 16-33-60063) and RAS Program 1.14P for support of their studies. A.S.S. is much obliged to the Scientific Council of the President of the Russian Federation (grant MK-3228.2015.3). A.S.N. is much obliged to the Government of Saint Petersburg for the Grant for young scientists in 2016. Physicochemical studies were performed at the Center for Magnetic Resonance, Center for X-ray Diffraction Studies, Center for Chemical Analysis and Materials Research, and the Center for Thermogravimetric and Calorimetric Research (all belong to Saint Petersburg State University).



REFERENCES

(1) (a) Bokach, N. A.; Kuznetsov, M. L.; Kukushkin, V. Y. 1,3Dipolar cycloaddition of nitrone-type dipoles to uncomplexed and metal-bound substrates bearing the CN bond. Coord. Chem. Rev. 2011, 255, 2946−2967. (b) Hermkens, P. H. H.; Van Maarseveen, J. H.; Kruse, C. G.; Scheeren, H. W. 1,3-Dipolar cycloaddition of nitrones with nitriles. Scope and mechanistic study. Tetrahedron 1988, 44, 6491−6504. (c) Plate, R.; Hermkens, P. H. H.; Smits, J. M. M.; Nivard, R. J. F.; Ottenheijm, H. C. J. Employment of nitriles in the stereoselective cycloaddition to nitrones. J. Org. Chem. 1987, 52, 1047−1051. (d) Wagner, G.; Garland, T. Synthesis of 5trichloromethyl-Δ4-1,2,4-oxadiazolines and their rearrangement into formamidine derivatives. Tetrahedron Lett. 2008, 49, 3596−3599. (2) (a) Pace, A.; Pierro, P. The new era of 1,2,4-oxadiazoles. Org. Biomol. Chem. 2009, 7, 4337−4348. (b) Burns, A. R.; Kerr, J. H.; Kerr, W. J.; Passmore, J.; Paterson, L. C.; Watson, A. J. B. Tuned methods for conjugate addition to a vinyl oxadiazole; synthesis of pharmaceutically important motifs. Org. Biomol. Chem. 2010, 8, 2777−2783. (3) Coley, H. M.; Sarju, J.; Wagner, G. Synthesis and Characterization of Platinum(II) Oxadiazoline Complexes and their In Vitro Antitumor Activity in Platinum-Sensitive and -Resistant Cancer Cell Lines. J. Med. Chem. 2008, 51, 135−141. (4) Meanwell, N. A. Synopsis of Some Recent Tactical Application of Bioisosteres in Drug Design. J. Med. Chem. 2011, 54, 2529−2591. (5) Yu, Y.; Fujita, H.; Ohno, M.; Eguchi, S. 1,3-Dipolar cycloaddition of nitrile functions with nitrones under high-pressure conditions − a 1388

DOI: 10.1021/acsomega.7b00130 ACS Omega 2017, 2, 1380−1391

ACS Omega

Article

new and direct synthesis of 2,3-dihydro-1,2,4-oxadiazole derivatives. Synthesis 1995, 498−500. (6) (a) Troisi, L.; Caccamese, S.; Pilati, T.; Videtta, V.; Rosato, F. 3 + 2-Cycloaddition Reaction between Chiral Oxaziridines and Nitriles: Enantioselective Synthesis of 2,3-Dihydro-1,2,4-oxadiazoles. Synthesis 2010, 3211−3216. (b) Troisi, L.; Ronzini, L.; Rosato, F.; Videtta, V. 3 + 2 Cycloaddition of Oxaziridines with Nitriles: Synthesis of 2,3Dihydro-1,2,4-Oxadiazoles. Synlett 2009, 1806−1808. (c) Videtta, V.; Perrone, S.; Rosato, F.; Alifano, P.; Tredici, S. M.; Troisi, L. Condensed Oxaziridine-Mediated 3 + 2 Cycloaddition: Synthesis of Polyhetero-bicyclo Compounds. Synlett 2010, 2781−2783. (7) (a) Pavez, H.; Marquez, A.; Navarrete, P.; Tzichinovsky, S.; Rodriguez, H. Regiospecific syntheses of delta-4-1,2,4-oxadiazolin-3carboxylic acids. Tetrahedron 1987, 43, 2223−2228. (b) Rodriguez, H.; Pavez, H.; Marquez, A.; Navarrete, P. Reaction between delta-2oxazolin-5-ones and nitrosobenzene − formation of 1,2,4-oxadiazolines. Tetrahedron 1983, 39, 23−27. (8) (a) Kritchenkov, A. S.; Bokach, N. A.; Haukka, M.; Kukushkin, V. Y. Unexpectedly efficient activation of push-pull nitriles by a PtII center toward dipolar cycloaddition of Z-nitrones. Dalton Trans. 2011, 40, 4175−4182. (b) Kritchenkov, A. S.; Bokach, N. A.; Kuznetsov, M. L.; Dolgushin, F. M.; Tung, T. Q.; Molchanov, A. P.; Kukushkin, V. Y. Facile and Reversible 1,3-Dipolar Cycloaddition of Aryl Ketonitrones to Platinum(II)-Bound Nitriles: Synthetic, Structural, and Theoretical Studies. Organometallics 2012, 31, 687−699. (c) Kritchenkov, A. S.; Bokach, N. A.; Starova, G. L.; Kukushkin, V. Y. A Palladium(II) Center Activates Nitrile Ligands toward 1,3-Dipolar Cycloaddition of Nitrones Substantially More than the Corresponding Platinum(II) Center. Inorg. Chem. 2012, 51, 11971−11979. (d) Bokach, N. A.; Krokhin, A. A.; Nazarov, A. A.; Kukushkin, V. Y.; Haukka, M.; Frausto da Silva, J. J. R.; Pombeiro, A. J. L. Interplay between nitrones and (nitrile)PdII complexes: Cycloaddition vs. Complexation followed by cyclopalladation and deoxygenation reactions. Eur. J. Inorg. Chem. 2005, 3042−3048. (e) Wagner, G.; Haukka, M. Stereoselective [2 + 3] cycloaddition of nitrones to platinum-bound organonitriles. First enantioselective synthesis of Δ4-1,2,4-oxadiazolines. J. Chem. Soc., Dalton Trans. 2001, 2690−2697. (f) Wagner, G.; Haukka, M.; Fraústo da Silva, J. J.; Pombeiro, A. J. L.; Kukushkin, V. Y. [2 + 3] cycloaddition of nitrones to platinum-bound organonitriles: effect of metal oxidation state and of nitrile substituent. Inorg. Chem. 2001, 40, 264−271. (g) Wagner, G.; Pombeiro, A. J. L.; Kukushkin, V. Y. Platinum(IV)-Assisted [2 + 3] Cycloaddition of Nitrones to Coordinated Organonitriles. Synthesis of Δ4-1,2,4-Oxadiazolines. J. Am. Chem. Soc. 2000, 122, 3106−3111. (h) Desai, B.; Danks, T. N.; Wagner, G. Cycloaddition of nitrones to free and coordinated (E)cinnamonitrile: effect of metal coordination and microwave irradiation on the selectivity of the reaction. Dalton Trans. 2003, 2544−2549. (i) Desai, B.; Danks, T. N.; Wagner, G. Ligand discrimination in the reaction of nitrones with [PtCl2(PhCN)2]. Selective formation of mono-oxadiazoline and mixed bis-oxadiazoline complexes under thermal and microwave conditions. Dalton Trans. 2004, 166−171. (j) Coley, H. M.; Sarju, J.; Wagner, G. Synthesis and Characterization of Platinum(II) Oxadiazoline Complexes and their In Vitro Antitumor Activity in Platinum-Sensitive and -Resistant Cancer Cell Lines. J. Med. Chem. 2008, 51, 135−141. (k) Charmier, M. A. J.; Haukka, M.; Pombeiro, A. J. L. Unprecedented single-pot synthesis of nitrilederived ketoimino platinum(II) complexes by ring opening of Δ41,2,4-oxadiazolines. Dalton Trans. 2004, 2741−2745. (9) (a) Sarju, J.; Arbour, J.; Sayer, J.; Rohrmoser, B.; Scherer, W.; Wagner, G. Synthesis and characterization of mixed ligand Pt(II) and Pt(IV) oxadiazoline complexes. Dalton Trans. 2008, 5302−5312. (b) Kritchenkov, A. S.; Gurzhiy, V. V.; Bokach, N. A.; Kalibabchuk, V. A. Dichloridobis[3-(4-chlorophenyl)-2,N,N-trimethyl-2,3-dihydro1,2,4-oxadiazole-5-amine-κN4]platinum(II)-4-chlorobenzaldehyde (1/ 1). Acta Crystallogr., Sect. E: Struct. Rep. Online 2013, 69, m446−m447. (c) Kritchenkov, A. S.; Lavnevich, L. V.; Starova, G. L.; Bokach, N. A.; Kalibabchuk, V. A. Dichloridobis[3-(4-methoxyphenyl)-2-methyl-5(piperidin-1-yl)-2,3-dihydro-1,2,4-oxadiazole-κN4]platinum(II). Acta Crystallogr., Sect. E: Struct. Rep. Online 2013, 69, m435−m436.

(10) Smirnov, A. S.; Yandanova, E. S.; Bokach, N. A.; Starova, G. L.; Gurzhiy, V. V.; Avdontceva, M. S.; Zolotarev, A. A.; Kukushkin, V. Y. Zinc(II)-mediated generation of 5-amino substituted 2,3-dihydro1,2,4-oxadiazoles and their further Zn-II-catalyzed and O-2-involving transformations. New J. Chem. 2015, 39, 9330−9344. (11) Wagner, G. Reactions of nitrones with transition metal nitrile complexes: cycloaddition, ligand substitution, or hydrolysis. Inorg. Chim. Acta 2004, 357, 1320−1324. (12) (a) Bolotin, D. S.; Bokach, N. A.; Kukushkin, V. Y. Coordination chemistry and metal-involving reactions of amidoximes: Relevance to the chemistry of oximes and oxime ligands. Coord. Chem. Rev. 2016, 313, 62−93. (b) Boyarskiy, V. P.; Bokach, N. A.; Luzyanin, K. V.; Kukushkin, V. Y. Metal-Mediated and Metal-Catalyzed Reactions of Isocyanides. Chem. Rev. 2015, 115, 2698−2779. (13) (a) Rassadin, V. A.; Boyarskiy, V. P.; Kukushkin, V. Y. Facile Gold-Catalyzed Heterocyclization of Terminal Alkynes and Cyanamides Leading to Substituted 2-Amino-1,3-Oxazoles. Org. Lett. 2015, 17, 3502−3505. (b) Timofeeva, S. A.; Kinzhalov, M. A.; Valishina, E. A.; Luzyanin, K. V.; Boyarskiy, V. P.; Buslaeva, T. M.; Haukka, M.; Kukushkin, V. Y. Application of palladium complexes bearing acyclic amino(hydrazido)carbene ligands as catalysts for copper-free Sonogashira cross-coupling. J. Catal. 2015, 329, 449−456. (c) Rocha, B. G. M.; Valishina, E. A.; Chay, R. S.; Guedes da Silva, M. F. C.; Buslaeva, T. M.; Pombeiro, A. J. L.; Kukushkin, V. Y.; Luzyanin, K. V. ADCmetal complexes as effective catalysts for hydrosilylation of alkynes. J. Catal. 2014, 309, 79−86. (d) Kinzhalov, M. A.; Luzyanin, K. V.; Boyarskiy, V. P.; Haukka, M.; Kukushkin, V. Y. ADC-Based Palladium Catalysts for Aqueous Suzuki−Miyaura Cross-Coupling Exhibit Greater Activity than the Most Advantageous Catalytic Systems. Organometallics 2013, 32, 5212−5223. (14) Bokach, N. A. Cycloaddition of nitrones to metal-activated nitriles and isocyanides. Russ. Chem. Rev. 2010, 79, 89−100. (15) (a) Sustmann, R. Simple model for substituent effects in cycloaddition reactions. II. Diels-Alder reaction. Tetrahedron Lett. 1971, 2721−2724. (b) Sustmann, R. Simple model for substituent effects in cycloaddition reactions. I. 1,3-Dipolar cycloadditions. Tetrahedron Lett. 1971, 2717−2720. (16) Eberson, L.; McCullough, J. J.; Hartshorn, C. M.; Hartshorn, M. P. 1,3-Dipolar cycloaddition of α-phenyl-N-tert-butylnitrone (PBN) to dichloro- and dibromo-malononitrile, chlorotricyanomethane and tetracyanomethane. Structure of products and kinetics of their formation. J. Chem. Soc., Perkin Trans. 2 1998, 41−48. (17) (a) Zhang, Y.; Sun, W.; Freund, C.; Santos, A. M.; Herdtweck, E.; Mink, J.; Kuhn, F. E. Cationic copper(I) and silver(I) nitrile complexes with fluorinated weakly coordinating anions: Metal-nitrile bond strength and its influence on the catalytic performance. Inorg. Chim. Acta 2006, 359, 4723−4729. (b) Nelson, I. V.; Iwamoto, R. T.; Kleinberg, J. The formation of copper(1) chloride by the action of copper(1) ion on carbon tetrachloride in 2-butanol. J. Am. Chem. Soc. 1964, 86, 364. (c) Farha, F.; Iwamoto, R. T. The Preparation and Infrared Examination of the 2-, 3-, and 4-Cyanopyridine Complexes of Copper(I), Silver(I), and Gold (I) Perchlorates. lnorg. Chem. 1965, 4, 844−848. (d) Kubota, M.; Johnston, D. L. Rotational Conformers of Glutaronitrile in Complexes with Metal Ions. J. Am. Chem. Soc. 1966, 88, 2451−2455. (e) Kubota, M.; Johnston, D. L. Nitrile complexes of copper(I) perchlorate. J. Inorg. Nucl. Chem. 1967, 29, 769−773. (f) Meerwein, H.; Hederich, V.; Wunderlich, K. Studies with anhydro silver fluoroborate and copper(I) fluoborate. Arch. Pharm. 1958, 291, 541−554. (g) Bergerhoff, G. Darstellung von Kupfer(I)- und Gold(I)Verbindungen in Acetonitril. Z. Anorg. Allg. Chem. 1964, 327, 139− 142. (h) Hemmerich, P.; Sigwart, C. Cu(CH3CN)2+, ein Mittel zum Studium homogener Reaktionen des einwertigen Kupfers in wässriger Lösung. Experientia 1963, 19, 488−489. (i) Hathaway, B. J.; Holoh, D. G.; Postlethwaite, J. D. The preparation and properties of some tetrakis(methylcyanide)copper(I) complexes. J. Chem. Soc. 1961, 3215−3218. (18) (a) Sivasubramanian, S.; Manisankar, P.; Palaniandavar, M.; Arumugam, N. Donor properties of the nitrone function in copper(II) complexes of some 2-hydroxy-1-naphthylnitrones. Transition Met. 1389

DOI: 10.1021/acsomega.7b00130 ACS Omega 2017, 2, 1380−1391

ACS Omega

Article

Chem. 1982, 7, 346−349. (b) Thirumalaikumar, M.; Sivakolunthu, S.; Ponnuswamy, A.; Sivasubramanian, S. Synthesis and characterization of cobalt(II), nickel(II) and copper(II) bis-chelates. Indian J. Chem. 1999, 38A, 720−722. (c) Villamena, F. A.; Crist, R. D. Metal-nitrone complexes: spin trapping and solution characterization. J. Chem. Soc., Dalton Trans. 1998, 4055−4062. (d) Petkova, E. G.; Domasevitch, K. V.; Gorichko, M. V.; Zub, V. Y.; Lampeka, R. D. New Coordination Compounds Derived from Nitrone Ligands: Copper(II) Complexes with 8-Hydroxyquinoline-2 carbaldehyde- and Pyridine-2-carbaldehyde-N-methylnitrones. Z. Naturforsch., B: J. Chem. Sci. 2001, 56, 1264−1270. (19) (a) Baliah, V.; Nair, V. C. Infrared spectra of some nitrones, azoxybenzenes, O-methyl oximes and sydnones. J. Indian Chem. Soc. 1989, 66, 42−44. (b) Smith, P. A. S.; Robertson, J. Some Factors Affecting the Site of Alkylation of Oxime Salts. J. Am. Chem. Soc. 1962, 84, 1197. (c) Hamer, J.; Macaluso, A. Nitrones. Chem. Rev. 1964, 64, 473−495. (20) Novikov, A. S.; Kuznetsov, M. L.; Pombeiro, A. J. L. Theory of the Formation and Decomposition of N-Heterocyclic Aminooxycarbenes through Metal-Assisted [2 + 3]-Dipolar Cycloaddition/ Retro-Cycloaddition. Chem. − Eur. J. 2013, 19, 2874−2888. (21) Kuznetsov, M. L.; Kukushkin, V. Y.; Pombeiro, A. J. L. Comparative Theoretical Study of 1,3-Dipolar Cycloadditions of AllylAnion Type Dipoles to Free and Pt-Bound Nitriles. J. Org. Chem. 2010, 75, 1474−1490. (22) Kuznetsov, M. L.; Bokach, N. A.; Kukushkin, V. Y.; Pakkanen, T.; Wagner, G.; Pombeiro, A. J. L. Metal-assisted coupling of oximes and nitriles: a synthetic, structural and theoretical study. J. Chem. Soc., Dalton Trans. 2000, 4683−4693. (23) Bader, R. F. W. A quantum theory of molecular structure and its applications. Chem. Rev. 1991, 91, 893−928. (24) (a) Ding, X.; Tuikka, M. J.; Hirva, P.; Kukushkin, V. Y.; Novikov, A. S.; Haukka, M. Fine-tuning halogen bonding properties of diiodine through halogen−halogen charge transfer − extended [Ru(2,2′-bipyridine)(CO)2X2]·I2 systems (X = Cl, Br, I). CrystEngComm 2016, 18, 1987−1995. (b) Ivanov, D. M.; Kirina, Y. V.; Novikov, A. S.; Starova, G. L.; Kukushkin, V. Y. Efficient π-stacking with benzene provides 2D assembly of trans-[PtCl2(p-CF3C6H4CN)2]. J. Mol. Struct. 2016, 1104, 19−23. (c) Ivanov, D. M.; Novikov, A. S.; Ananyev, I. V.; Kirina, Y. V.; Kukushkin, V. Y. Halogen bonding between metal centers and halocarbons. Chem. Commun. 2016, 52, 5565−5568. (d) Novikov, A. S.; Kuznetsov, M. L. Nucleophilicity of Oximes Based upon Addition to a Nitrilium closo-Decaborate Cluster. Inorg. Chim. Acta 2012, 380, 78−89. (e) Serebryanskaya, T. V.; Novikov, A. S.; Gushchin, P. V.; Haukka, M.; Asfin, R. E.; Tolstoy, P. M.; Kukushkin, V. Y. Identification and H(D)-bond energies of C H(D)···Cl interactions in chloride−haloalkane clusters: a combined Xray crystallographic, spectroscopic, and theoretical study. Phys. Chem. Chem. Phys. 2016, 18, 14104−14112. (f) Mikhaylov, V. N.; Sorokoumov, V. N.; Korvinson, K. A.; Novikov, A. S.; Balova, I. A. Synthesis and Simple Immobilization of Palladium(II) Acyclic Diaminocarbene Complexes on Polystyrene Support as Efficient Catalysts for Sonogashira and Suzuki−Miyaura Cross-Coupling. Organometallics 2016, 35, 1684−1697. (g) Mikherdov, A. S.; Kinzhalov, M. A.; Novikov, A. S.; Boyarskiy, V. P.; Boyarskaya, I. A.; Dar’in, D. V.; Starova, G. L.; Kukushkin, V. Y. Difference in Energy between Two Distinct Types of Chalcogen Bonds Drives Regioisomerization of Binuclear (Diaminocarbene)PdII Complexes. J. Am. Chem. Soc. 2016, 138, 14129−14137. (h) Bolotin, D. S.; Burianova, V. K.; Novikov, A. S.; Demakova, M. Y.; Pretorius, C.; Mokolokolo, P. P.; Roodt, A.; Bokach, N. A.; Suslonov, V. V.; Zhdanov, A. P.; Zhizhin, K. Y.; Kuznetsov, N. T.; Kukushkin, V. Y. Nucleophilicity of Oximes Based upon Addition to a Nitrilium closoDecaborate Cluster. Organometallics 2016, 35, 3612−3623. (i) Melekhova, A. A.; Novikov, A. S.; Bokach, N. A.; Avdonceva, M. S.; Kukushkin, V. Y. Characterization of Cu-ligand bonds in trispyrazolylmethane isocyanide copper(I) complexes based upon combined X-ray diffraction and theoretical study. Inorg. Chim. Acta 2016, 450, 140−145.

(25) (a) Menezes, F. M. C.; Kuznetsov, M. L.; Pombeiro, A. J. L. Isocyanide Complexes with Platinum and Palladium and Their Reactivity toward Cycloadditions with Nitrones to Form Aminooxycarbenes: A Theoretical Study. Organometallics 2009, 28, 6593− 6602. (b) Kuznetsov, M. L.; Kukushkin, V. Y.; Dement’ev, A. I.; Pombeiro, A. J. L. 1,3-Dipolar Cycloaddition of Nitrones to Free and Pt-Bound Nitriles. A Theoretical Study of the Activation Effect, Reactivity, and Mechanism. J. Phys. Chem. A 2003, 107, 6108−6120. (26) (a) Domingo, L. R. Theoretical Study of 1,3-Dipolar Cycloaddition Reactions with Inverse Electron Demand − A DFT Study of the Lewis Acid Catalyst and Solvent Effects in the Reaction of Nitrones with Vinyl Ethers. Eur. J. Org. Chem. 2000, 2000, 2265−2272. (b) Tanaka, J.; Kanemasa, S. Ab initio study of Lewis acid catalyzed nitrone cycloaddition to electron deficient alkenes. Does a Lewis acid catalyst change the reaction mechanism? Tetrahedron 2001, 57, 899− 905. (c) Gothelf, K. V.; Hazell, R. G.; Jørgensen, K. A. Control of Diastereo- and Enantioselectivity in Metal-Catalyzed 1,3-Dipolar Cycloaddition Reactions of Nitrones with Alkenes. Experimental and Theoretical Investigations. J. Org. Chem. 1996, 61, 346−355. (27) (a) Moyano, A.; Pericas, M. A.; Valenti, E. A theoretical study on the mechanism of the thermal and the acid-catalyzed decarboxylation of 2-oxetanones (β-lactones). J. Org. Chem. 1989, 54, 573−582. (b) Lecea, B.; Arrieta, A.; Roa, G.; Ugalde, J. M.; Cossio, F. P. Catalytic and Solvent Effects on the Cycloaddition Reaction between Ketenes and Carbonyl Compounds To Form 2-Oxetanones. J. Am. Chem. Soc. 1994, 116, 9613−9619. (c) Morao, I.; Lecea, B.; Cossío, F. P. In-Plane Aromaticity in 1,3-Dipolar Cycloadditions. J. Org. Chem. 1997, 62, 7033−7036. (d) Cossío, F. P.; Morao, I.; Jiao, H.; Schleyer, P. V. R. InPlane Aromaticity in 1,3-Dipolar Cycloadditions. Solvent Effects, Selectivity, and Nucleus-Independent Chemical Shifts. J. Am. Chem. Soc. 1999, 121, 6737−6746. (28) Kuznetsov, M. L. Theoretical studies on [3 + 2]-cycloaddition reactions. Russ. Chem. Rev. 2006, 75, 935−960. (29) Carda, M.; Portolés, R.; Murga, J.; Uriel, S.; Marco, J. A.; Domingo, L. R.; Zaragozá, R. J.; Röper, H. Stereoselective 1,3-Dipolar Cycloadditions of a Chiral Nitrone Derived from Erythrulose. An Experimental and DFT Theoretical Study. J. Org. Chem. 2000, 65, 7000−7009. (30) Silva, M. A.; Goodman, J. M. Nitrone cyclisations: the development of a semi-quantitative model from ab initio calculations. Tetrahedron 2002, 58, 3667−3671. (31) Yeung, M. L.; Li, W. K.; Liu, H. J.; Wang, Y.; Chan, K. S. Kinetic and Theoretical Studies of the [3 + 2] Cycloadditions of Alkynyl Fischer Carbene Complexes with N-Alkyl Nitrones. J. Org. Chem. 1998, 63, 7670−7673. (32) Merino, P.; Anoro, S.; Merchan, F.; Tomas, T. 1,3-Dipolar Cycloaddition between Hetaryl Nitrones and Methyl Acrylate: Theoretical Study and Application to the Synthesis of Functionalized Pyrrolidines. Heterocycles 2000, 53, 861−875. (33) Baranski, A. Substituent effect on the [2 + 3] cycloaddition of (E)-β-Nitrostyrene with (Z)-C,N-Diarylnitrones. J. Phys. Org. Chem. 2002, 15, 78−82. (34) (a) Naji, N.; Marakchi, K.; Kabbaj, O. K.; Komiha, N.; Joffre, J.; Chraïbi, M.; Soufiaoui, M. Approche théorique de la cycloaddition dipolaire-1,3 sur des dipolarophiles fluorés. J. Fluorine Chem. 1999, 94, 127−133. (b) Marakchi, K.; Kabbaj, O. K.; Komiha, N. Etude DFT du mécanisme des réactions de cycloaddition dipolaire-1,3 de la C,Ndiphénylnitrone avec des dipolarophiles fluorés de type éthylénique et acétylénique. J. Fluorine Chem. 2002, 114, 81−89. (35) Bokach, N. A.; Kuznetsov, M. L.; Kukushkin, V. Y. 1,3-Dipolar cycloaddition of nitrone-type dipoles to uncomplexed and metalbound substrates bearing the CN triple bond. Coord. Chem. Rev. 2011, 255, 2946−2967. (36) Stolley, R. M.; Maczka, M. T.; Louie, J. Nickel-Catalyzed [2 + 2 + 2] Cycloaddition of Diynes and Cyanamides. Eur. J. Org. Chem. 2011, 3815−3824. (37) Kubas, G. J. AQ8: Please provide a DOI number for ref 37 or indicate if one doesn’t 1390

DOI: 10.1021/acsomega.7b00130 ACS Omega 2017, 2, 1380−1391

ACS Omega

Article

exist.Tetrakis(Acetonitrile)Copper(I) Hexafluorophosphate. Inorg. Synth. 1979, 19, 90−92. (38) (a) Exner, O. A new synthesis of N-methyl ketoximes. Collect. Czech. Chem. Commun. 1951, 16, 258−267. (b) Ovens, A. C. Inhibitors of procollagen C-proteinase. Dissertation, University of Manchester, UK, 1999; 2012:1449371, Chem. Abs. 157:519778. (39) Sheldrick, G. M. A short history of SHELX. Acta Crystallogr., Sect. A: Found. Adv. 2008, 64, 112−122. (40) Dolomanov, O. V.; Bourhis, L. J.; Gildea, R. J.; Howard, J. A. K.; Puschmann, H. OLEX2: A complete structure solution, refinement and analysis program. J. Appl. Crystallogr. 2009, 42, 339−341. (41) CrysAlisPro, 1.171.36.20; Agilent Technologies. (42) Zhao, Y.; Truhlar, D. G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215−241. (43) 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.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Keith, T.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, O.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09, revision C.01; Gaussian, Inc.: Wallingford, CT, 2010. (44) Figgen, D.; Rauhut, G.; Dolg, M.; Stoll, H. Energy-consistent pseudopotentials for group 11 and 12 atoms: adjustment to multiconfiguration Dirac−Hartree−Fock data. Chem. Phys. 2005, 311, 227− 244. (45) (a) Gonzalez, C.; Schlegel, H. B. Improved algorithms for reaction path following: Higher-order implicit algorithms. J. Chem. Phys. 1991, 5853−5860. (b) Gonzalez, C.; Schlegel, H. B. An improved algorithm for reaction path following. J. Chem. Phys. 1989, 90, 2154−2161. (c) Gonzalez, C.; Schlegel, H. B. Reaction path following in mass-weighted internal coordinates. J. Phys. Chem. 1990, 94, 5523−5527. (46) Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. Universal Solvation Model Based on Solute Electron Density and on a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions. J. Phys. Chem. B 2009, 113, 6378−6396. (47) Wertz, D. H. Relationship between the gas-phase entropies of molecules and their entropies of solvation in water and 1-octanol. J. Am. Chem. Soc. 1980, 102, 5316−5322. (48) Cooper, J.; Ziegler, T. A Density Functional Study of SN2 Substitution at Square-Planar Platinum(II) Complexes. Inorg. Chem. 2002, 41, 6614−6622. (49) Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580−592. (50) Glendening, E. D.; Landis, C. R.; Weinhold, F. Natural bond orbital methods. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2012, 2, 1− 42.

1391

DOI: 10.1021/acsomega.7b00130 ACS Omega 2017, 2, 1380−1391