Carla
Gomes
a,
Mariana
Costa
a,
Susana M. M.
Lopes
a,
Bernardo Albuquerque
Nogueira
a,
Rui
Fausto
ab,
José A.
Paixão
c,
Teresa M. V. D.
Pinho e Melo
a,
Luísa M. D. R. S.
Martins
d and
Marta
Pineiro
*a
aCQC-IMS and Department of Chemistry, University of Coimbra, 3004-535, Coimbra, Portugal. E-mail: mpineiro@qui.uc.pt
bFaculty of Sciences and Letters, Department of Physics, Istanbul Kultur University, Ataköy Campus, Bakirköy 34156, Istanbul, Turkey
cCFisUC and Department of Physics, University of Coimbra, 3004-516 Coimbra, Portugal
dCentro de Química Estrutural, Institute of Molecular Sciences, Departamento de Engenharia Química, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal
First published on 15th December 2023
The development of sustainable processes requires the integration of the principles of green chemistry. In this work, we develop the solvent-free synthesis of bis- and tris-(pyrazol-1-yl)methanes from pyrazole and halooximes via 1,4-addition to conjugated α-halogenated nitrosoalkenes generated in situ under ball milling. The pro-ligands were complexed with copper, using a solvent-free procedure with equimolar quantities of reactants, to give catalysts copper(II) C-scorpionates with the oxime substituent at the C-center. The new catalysts were used to efficiently catalyze azide–alkyne cycloaddition reactions under solvent-free mechanochemical conditions. Through the overall process, we combine atom economy, the use of safer solvents and auxiliaries, design for energy efficiency and reduction of derivatives and catalysis (the 1st, 5th, 6th, 8th, and 9th principles of Green Chemistry). The sustainability metrics show that mechanochemistry provides a greener and more sustainable approach for the preparation of the catalysts and for the copper-catalyzed azide–alkyne cycloaddition (CuAAC) reactions. The development of novel C-scorpionates and their metal complexes, and their subsequent application in CuAAC reactions under mechanochemical conditions, presents an exciting avenue for the future of sustainable click chemistry.
The copper-catalyzed azide–alkyne cycloaddition (CuAAC) reaction yielding triazoles is a widely used metal catalyzed reaction and a truly “click reaction”. The unique features of the CuAAC reaction, including its biorthogonality, high regioselectivity and its characteristics in line with the requirements of green chemistry,2–4 make it a popular choice for applications in diverse fields, such as medicinal and polymer chemistries, as well as materials and nanomaterials science. This was publicly recognized with the award of Nobel Prize in Chemistry to Morten Meldal, Carolyn R. Bertozzi and K. Barry Sharpless in 2022.5–8 Researchers in this area have focused on developing more sustainable conditions by seeking greener solvents and copper precursors that can generate copper(I) efficiently, while decreasing or eliminating the need for reducing agents and bases.
Mechanochemical click CuAAC reactions enable the preparation of triazoles in good yields and with high regioselectivity under solvent-free conditions, using diverse catalytic systems. Applying Cu(OAc)2 as the copper source, sodium ascorbate as the reducing agent, and SiO2 as the grinding auxiliary in a planetary ball mill, triazoles were obtained in good yields, with almost no side products and in a short reaction time (∼10 min).9 Other studies carried out using a mixer milling system showed moderate yields of quinolines containing the phenyl-1,2,3-triazole unit when Cu(OAc)2 was used as a catalyst without any reducing agent. Using the same type of milling, the reaction yield increased when copper(I) was used as a catalyst in the presence of a tertiary amine, and when milling vessels and balls of brass (an alloy of copper and zinc) were used as catalysts, with N,N-di-isopropylethylamine (DIPEA), together with a small amount of acetic acid as auxiliaries.10 Furthermore, copper powder was also selected as a catalyst for the reaction between benzyl azides and alkynes under planetary ball milling.11 In copper vials as reaction containers, the catalyzed reaction between aryl and alkyl alkynes, benzyl bromine, and sodium azide afforded target triazoles, although longer reaction times were required.12 This study also allowed concluding on the safety of manipulating sodium azide under mechanochemical conditions (high-energy milling conditions), highlighting that azide reactions under mechanochemical conditions are indeed much safer than when conducted through traditional solution-based chemistry in a glass apparatus.12 Propargyl glycoside derivatives were also shown to be efficiently coupled with azides,13 by using CuSO4/sodium ascorbate and tBuOH/H2O, while applying similar experimental conditions, the preparation of cancer cell targeting polyethylene glycol-functionalized nanographene oxide was reported,14 pointing out the capability of mechanochemistry to solve the problem of promoting the reaction between reactants with different solubilities.
N-ligands protect copper(I) from oxidation, by modulating the redox potential of the Cu(II)/Cu(I) pair, and they influence the solvation energy of copper ions and also act as a base to promote the initial formation of copper(I) acetylide, thus increasing the CuAAC reaction rate. Monodentate N-ligands (amines) and polydentate N-ligands, such as 1,1,4,7,7-pentamethyldiethylenetriamine (PMDETA), significantly accelerate the CuAAC reaction.5
Tris(pyrazol-1-yl)methanes, also known as C-scorpionates, are among the most versatile classes of N-heterocyclic ligands,15–20 however the application of their copper complexes as catalysts for the CuAAC reaction is scarce. To the best of our knowledge, only two examples of such use of the copper complexes of C-scorpionates have been reported hitherto. Cano and co-workers used copper(I) C-scorpionates to efficiently promote the formation of N-sulfonyl-1,2,3-triazoles from sulfonylazides and alkynes in chloroform, at 40 °C for 24 h,21 while Mahmoud and co-workers used copper(I) C-scorpionates and copper(II) C-scorpionates as catalysts for the three-component CuAAC reaction in aqueous medium, under microwave irradiation, at 125 °C for 30 min.22
Having in mind the development of genuine sustainable processes, in the present investigation we explored the utilization of mechanochemistry for the synthesis of a novel class of C-scorpionates and of their metal complexes, and the subsequent application of the new Cu(II) complexes in click chemistry, in particular in the multicomponent CuAAC reaction, under mechanochemical conditions, to afford triazoles.
Scheme 1 Strategies for the synthesis of tris(pyrazol-1-yl)methane featuring an oxime group at the C-center. |
To obtain the desired compound, 3a, as depicted in Scheme 1, equivalent amounts of oxime 2, synthesized using a previously described methodology,23 pyrazole (1a) and sodium carbonate were mixed in a stainless steel jar under ball milling conditions. Specifically, high-speed vibrations at 20 Hz for 30 min with two stainless steel balls (7 mm of diameter) were utilized. Using this approach, 33% yield of the desired complex was obtained (entry 1, Table 1). By increasing the milling frequency to 25 Hz, the reaction yield reached 52% (entry 2, Table 1). However, increasing the milling frequency (up to 30 Hz), the number of milling balls, or extended reaction times did not improve the reaction yield (entries 3 to 5, Table 1). Interestingly, the order in which reactants were added to the milling jar was found to have a notable influence on the reaction outcome. It was observed that sodium carbonate should be added on top of pyrazole. Adding it on top of the oxime resulted in the immediate onset of the elimination reaction, with the formation of HCl, which led to the depletion of the reactant and subsequent dropping of the reaction yield to only 12% (entry 6, Table 1).
Mechanochemistry conditions | ||||
---|---|---|---|---|
Entry | Frequency (Hz) | Number of balls | Time (min) | Yielda (%) |
a Isolated yields after column chromatography [eluent: AcOEt/hexane (2:1)]. b Order of addition of the reactants: oxime, pyrazole and then sodium carbonate. c Order of addition of the reactants: oxime, sodium carbonate and then pyrazole. | ||||
1 | 20 | 2 | 30 | 33b |
2 | 25 | 2 | 30 | 52b |
3 | 30 | 2 | 30 | 51b |
4 | 25 | 2 | 60 | 52b |
5 | 25 | 3 | 30 | 49b |
6 | 25 | 2 | 30 | 12c |
The optimization of reaction conditions under mechanical action allowed to obtain the desired compound in 52% yield, after ball milling at 25 Hz during 30 min (see Scheme 1). The elimination of dichloromethane (use restricted in accordance with the European directive24) as well as the reduction of the reaction time from 16 h to 30 min is in line with the concerns for the development of more sustainable processes. According to previous studies, in solution, the reaction mechanism is highly dependent on both the basicity of the N-heterocycle and the nature of the oxime.25 In the case of the reaction of pyrazole and 2,2,2-trichloroacetaldehyde oxime to afford tris(pyrazol-1-yl)methane, the proposed mechanism consists of three consecutive conjugated addition reactions (Michael addition) of pyrazole to the in situ generated nitrosoalkenes.
The X-ray analysis of the crystals of 3a obtained after the recrystallization of the reaction product synthesized under mechanical action in ethanol (Fig. 1) clearly established that the (E)-oxime was obtained exclusively. The 1H-NMR spectrum of the reaction product (ESI,† Fig. S7) shows only one peak at 11.98 ppm, i.e., at the typical region of the proton of an oxime group with E configuration.26 The E configuration was not evident in the bidimensional NOESY spectrum (ESI,† Fig. S8).
The performed single-crystal X-ray crystallography study unambiguously confirmed that the E configuration is adopted by 3a in the crystalline state, as depicted in the ORTEP graph shown in Fig. 1 (left). All H-atoms could be located using the difference electron density maps. The compound crystallizes in the centrosymmetric monoclinic space group P21/n with cell parameters a = 10.2989(3)Å, b = 9.5354(3) Å, c = 12.3745(4) Å, α = 90°, β = 91.395(2)°, and γ = 90°. Valence angles and bond lengths are unexceptional, with average values Nsp2–Nsp2 [1.355(2) Å], Csp2–Nsp2 [1.32(5) Å], Nsp2–Csp3 [1.453(4) Å], Csp3–Csp2 [1.5159(2) Å], Csp3–Csp3 [1.37(2) Å], Nsp2–O [1.3906(14) Å]. The acetaldehyde oxime group is nearly planar, as shown by the torsion angle O1–N1–C2–C1 [178.99(10)°]. There is a small but significant deviation from the ideal tetrahedral geometry around C1, with valence angles N2–C1–N6 [107.92(9)°], N4–C1–N6 [108.91(9)°], N4–C1–N2 [107.85(9)°], N2–C1–C2 [108.91(9)°] N4–C1–C2 [110.38(9)°] and N6–C1–C2 [112.72(10)°]. The three pyrazole rings are planar with the rms deviation of the atoms from the least squares planes not exceeding 0.006 Å. The rings adopt a propeller-like conformation as shown by angles ∠(N2–N3–C3–C4–C5/N4–N5–C6–C7–C8) [88.06(5)°], ∠(N2–N3–C3–C4–C5/N6–N7–C9–C10–11) [81.39(6)°], ∠(N4–N5–C6–C7–C8/N6–N7–C9–C10–C11) [79.28(6)°] and depicted in Fig. 1 (center). The torsion angle N1–C2–C1–N2 [166.71(11)°] shows that the aldoxime group is rotated by 13° from the ideal bisecting position arising from a rotation around the single C1–C2 bond likely caused by a combination of steric effects and the involvement of O1 in intermolecular hydrogen bonding.
The main synthon occurring in the crystal structure is a pair of molecules related by an inversion center, linked by a bifurcated hydrogen bond, as shown in Fig. 1 (right). The donated hydroxyl proton is shared towards the lone pairs of two bare N atoms of the symmetry-related molecule, one belonging to the nearest pyrazole rings, and the other being the N atom of the oxime group, with geometry O1–H1⋯N1i [i = 1 − x, 1 − y, 1 − z; D–H 0.961(19); H⋯A: 2.190(19) Å; D⋯A: 3.0044(14) Å; ∠D–H–A: 141.7(15)°]; O1–H1⋯N7i [i = 1 − x, 1 − y, 1 − z; D–H 0.961(19); H⋯A: 2.162(18) Å; D⋯A: 2.9274(15) Å; ∠D–H–A: 135.7(14)°]. These pairs of molecules are further interconnected in a 3D network by weaker interactions of the C–H⋯N type involving the remaining two bare N atoms of the pyrazole rings. The packing is compact and there are no sizeable voids in the crystal structure accessible to solvent molecules. For the XRD see Tables S10–S22 in the ESI.†
The infrared (IR) and Raman spectra of the synthesized ligand 3a are presented in Fig. 2, together with the corresponding calculated spectra (B3LYP/6-311++G(d,p), scaled as described in the Experimental section – computational details). The theoretical spectra were calculated after the optimization of the dimeric structure present in the crystal of the compound, using the crystallographic data as initial geometry.
As shown in Fig. 2, the calculated and experimental spectra exhibited a remarkable resemblance, with the major spectral features observed experimentally being well reproduced by the calculations (for assignments, see the ESI,† Table S1). This fact should be highlighted, because the calculations were undertaken considering only the isolated dimeric structural unit of the compound, allowing to conclude that the inter-dimer interactions in the crystal are weak enough not to disturb extensively the intra-dimer vibrational potential. The general good reproducibility of the experimental spectra of the studied model ligand 3a by the calculations also suggested the possibility to use successfully vibrational data obtained theoretically to help to undertake the general structural characterization of the complexes of the oxime ligand investigated in the present study, as demonstrated in Section 2.2.
Besides the chosen model ligand 3a, other oxime derivatives of tris(pyrazolyl)methanes and bis(pyrazol-1-yl)methanes were synthesized via mechanochemistry in the present investigation (see Scheme 2, where the optimized conditions for the synthesis are also provided). The reproducibility of the synthetic methods using mechanochemistry is high. For example, the synthesis of ligand 3a yielded 52, 53 and 50% in three independent experiments. It is interesting to point out that the yields of the mechanochemical syntheses performed using oximes from aldehydes were significantly higher than those obtained using oximes derived from ketones, thus following the same trend observed for the corresponding reactions in solution.23
Scheme 2 Synthesis of tris(pyrazol-1-yl)methanes and bis(pyrazol-1-yl)methanes under mechanical action and the corresponding reaction yields. The precursor oximes, 1,1-dichloropropano-2-one oxime (4a) and 1-(1-p-chlorophenyl-1H-1,2,3-triazol-4-yl)-2,2-dibromoethanone oxime (4b), were prepared using a methodology described previously27 and through α-bromination with Br2/HOAc of 1-(1-p-chlorophenyl-1H-1,2,3-triazol-4-yl)ethanone28 followed by oximation with hydroxylamine hydrochloride in ethanol (yield: 81%), respectively (see structural characterization in the ESI,† Fig. S1–S6). |
Scheme 3 Mechanochemical synthesis of copper complexes 6a and 6b of (E)-2,2,2-tris(1H-pyrazol-1-yl)acetaldehyde oxime (3a). |
For comparison, the two model copper complexes were also synthesized in solution, following the general procedure described.29 The melting points of the complexes prepared under ball milling and in solution were found to be similar (151.7–152.6 °C and 151.9–152.6 °C, respectively, in the case of 6a, and 118.9–119.5 °C and 118.5–119.0 °C, respectively, for 6b). The high-resolution mass spectrum (HRMS) of 6a shows a peak at m/z = 379.0446, corresponding to the loss of an acetate group, and the elemental analysis data match the molecular formula C15H17CuN7O5 + H2O. For copper complex 6b, the high-resolution mass spectrum shows a peak at m/z = 390.977, which matches the molecular formula C11H11CuCl2N7O, and the elemental analysis correlates with the formula C11H11CuN7OCl2 + H2O.
Owing to the impossibility to obtain 1H-NMR spectra due to the paramagnetic nature of the copper(II) complexes, we endeavoured to acquire high-quality crystals of the synthesized complexes for the determination of their single crystal X-ray structure. Regrettably, our efforts proved unsuccessful. Consequently, we resorted to IR and Raman spectroscopies to investigate the structure of the copper complexes. The analysis of the experimental spectroscopic data was supported by the theoretical prediction of the vibrational spectra of the complexes (see the Experimental section for details).
The structures of the isolated molecule of the ligand and of complexes 6a and 6b were first investigated using quantum chemical calculations performed at the DFT(B3LYP)/6-311++G(d,p) level. The optimized structure of the ligand (Fig. 3, left) closely resembles that observed experimentally in crystalline phase (Fig. 1), with the aldoxime group being almost planar (O1–N1–C2–C1 torsion angle = 179.9°). There is also a deviation from the tetrahedral geometry around the C1 atom predicted by the calculation that is comparable to that found in the crystal, with the calculated N2–C1–N6, N4–C1–N6, N4–C1–N2, N2–C1–C2, N4–C1–C2 and N6–C1–C2 angles presenting values in the range 107.7–114.1°. Additionally, also in agreement with the crystallographic data, the calculations predict the plane of the oxime group slightly deviated (by ∼9.0°, as defined by the N1–C2–C1–N2 torsion angle) from the plane that bisects perpendicularly the angle N4–C1–N6.
Fig. 3 Optimized DFT(B3LYP)/6-311++G(d,p) calculated structures for pro-ligand 3 (left) and copper complexes 6a (centre) and 6b (right). |
In the optimized structures of complexes 6a and 6b (Fig. 3), the calculated bond angles and dihedral angles of the oxime ligand exhibit values similar to those found for the free ligand. The N1–C2–C1–N2 torsion, which represents the deviation from the plane that bisects perpendicularly the angle N4–C1–N6, is equal to 10.8° in both complexes 6a and 6b. On the other hand, the angles formed between the pyrazole rings and the oxime moiety plane differ in the two complexes and in the free ligand (see Fig. 3). In the free ligand, the orientation of the rings (as defined by the dihedral angles C2–C1–N2–N3, C2–C1–N4–N5 and C2–C1–N6–N7) confers the structure a propeller-like arrangement, with dihedral angles equal to 40.5, 42.2 and 63.8°, respectively. In the predicted structure of the 6a complex, the corresponding values of these dihedral angles are 70.2°, 174.1° and 164.4°, while in that of complex 6b, these dihedral angles are 72.2°, 172.1° and 172.2°, respectively. As shown in Fig. 3, this means that in complexes 6a and 6b, two of the pyrazole rings (rings 2 and 3) are oriented in such a way that their bare nitrogen atom (N5 and N7) faces the copper atom, while ring 1 is oriented resembling the correspondent ring in the free ligand structure.
The theoretical IR and Raman spectra of complexes 6a and 6b were obtained for the optimized structures and were compared with the corresponding calculated spectrum of pro-ligand 3a (Fig. 4). Both calculated Raman and IR spectra of the complexes present noticeable differences compared to the computed spectra of the ligand, some of which are highlighted in Fig. 4 by labelling the bands with the corresponding maximum intensity wavenumber. The tentative assignment of the major features of the experimental and calculated spectra of 6a and 6b is provided in the ESI,† Tables S2 and S3, respectively.
Fig. 4 Comparison between the IR (left) and Raman (right) scaled calculated spectra of the pro-ligand 3a and the complexes 6a and 6b, obtained at the B3LYP/6-311++G(d,p) level of theory. |
As expected, the calculated Raman spectra of the complexes present several intense additional bands, compared to the spectrum of the C-scorpionate pro-ligand. In the case of 6a, these bands are mainly due to the presence of the complex of the acetate (OAc) ligands, such as those predicted at 300 cm−1 (corresponding to the Cu–OAc stretching mode), 473 cm−1 (C–C bending mode of OAc) and 1596/1608 cm−1 (CO stretching of OAc). These bands can then be expected to be possible marker bands for this complex. On the other hand, the calculated Raman spectrum of complex 6b presents its highest intensity band at 315 cm−1, which is ascribable to the Cu–Cl stretching vibration and, thus, also appears as a possible marker band for this complex.
The calculated IR spectra of the isolated C-scorpionate pro-ligand and those of the two complexes are in general more similar to each other than the corresponding Raman spectra (excluding the CO stretching modes of OAc in 6a, in the 1595–1610 cm−1 region). However, in the range of 1350–1405 cm−1, it is also possible to observe some noticeable spectral differences in the IR spectra. While in the ligand spectrum, only one band is present (1383 cm−1, OH bending vibration), in the spectrum of 6a, there are bands at 1365 (OAc CH3 rocking mode), 1373 (OAc C–C stretching vibration) and at 1399 cm−1 (OH bending vibration), and in the spectrum of complex 6b, bands are observed at 1387 and at 1399 cm−1, which are assigned to the CH (aromatic) and OH bending modes, respectively.
Fig. 5 compares the experimental spectra (both IR and Raman) of ligand 3a and of complexes 6a and 6b, obtained under ball milling (the same comparison, with the spectra of the complexes obtained in solution is presented in Fig. S19, in the ESI†). The above highlighted marker bands of the 6a and 6b complexes, identified based on their calculated spectra, are, as expected, clearly visible in their experimental spectra (compare Fig. 4, 5 and Fig. S19 in the ESI†). In the Raman experimental spectrum of complex 6a, the Cu–OAc stretching mode appears at 365 cm−1, and the C–C bending vibration and CO stretching modes of the OAc ligands give rise to bands at 478 cm−1 and 1572 cm−1 (single broad band), respectively, in good correspondence with the theoretically predicted spectrum. In the case of 6b, the Cu–Cl stretching marker band vibration is observed in the experimental spectrum at 290 and 329 cm−1 (calculated: 315 cm−1), respectively. The appearance of the two bands due to the Cu–Cl stretching mode in the experimental spectra (for both the spectra of the complex obtained using mechanochemical and conventional procedures, see Fig. 5 and Fig. S19, ESI†) might indicate that the two chlorine atoms have different positions relative to the metal atom in the complex, with one being more tightly bonded than the other. In the case of the IR spectra, the above-mentioned identified theoretical marker bands for the ligand and the complexes in the 1350–1405 cm−1 spectral range are observed in the experimental spectra of the ligand at 1384 cm−1, at 1365 and 1388 cm−1 for complex 6a (the latter being a broad band comprising both the OAc C–C stretching vibration and the OH bending mode), and at 1384 and 1407 cm−1 for complex 6b, showing very good correspondence with the predicted frequency values.
Fig. 5 Comparison between the IR (left) and Raman (right) experimental spectra of pro-ligand 3a and complexes 6a and 6b, obtained under ball milling. |
The resemblance of the experimental spectra obtained for the complexes obtained under ball milling and from the solution (compare Fig. 5 and Fig. S19 in the ESI†) is highlighted. The few noticeable spectral differences might be attributed to the different crystalline/amorphous phases (or degree of crystallinity) resulting from the two complexation procedures.
Following the success in preparing the C-scorpionate complexes 6a and 6b, a set of other copper complexes of triazole-oxime ligands was prepared, using the above described mechanochemistry methodology, in high yields (Chart 1; characterization provided in the ESI†).
Scheme 4 Three-component AAC reaction of 3-butynol (10), benzyl bromide (11) and sodium azide (12) used as the model reaction. |
Ball milling was applied for all the catalyzed reactions, using equimolar amounts of the three reactants in stainless steel jars with two stainless steel balls (7 mm) and a frequency vibration of 30 Hz. Ball milling of the three reactants under the above described conditions for 3 h did not result in any reaction product, indicating that the metal in stainless steel jars and balls cannot catalyze the reaction (entry 1, Table 3). Nonetheless, the use of 1.5 mol% (relative to alkyne) of copper acetate monohydrate as the catalyst led to the reaction product in 37% yield, after isolation and purification through column chromatography (entry 2, Table 3). This result highlights the crucial role of copper as a catalyst, and the mechanochemical action leads to the formation of catalytic Cu(I) species due to the reduction of Cu(OAc)2 in the presence of sodium azide (12), benzyl bromide (11) and the alkyne (10). Moreover, the reduction of Cu(II) was promoted under ball milling conditions through the Glaser reaction30–33 without requiring reducing agents, additional bases or water.
To determine the influence of basicity of the oxime group [pKa (aldoxime) = 28.5, Bordwell pKa Table34] in the catalytic cycle, the reaction was performed adding copper acetate and ligand 3a to the azide/alkyne/benzyl mixture. After milling for 90 min under the above conditions, 1,2,3-triazole 13 was obtained in 33% yield (entry 3, Table 2), proving that the basicity of the oxime does not influence the reaction.
Entry | Catalyst | Reaction conditions | Yield (%) |
---|---|---|---|
1 | — | Ball milling 30 Hz, 3 h | 0 |
2 | Cu(OAc)2·H2O | Ball milling 30 Hz, 90 min | 37 |
3 | Cu(OAc)2·H2O + 3a | Ball milling 30 Hz, 90 min | 33 |
4 | 6a | Ball milling 30 Hz, 90 min | 78 |
5 | 6a | H2O/MeOH MW (125 °C, 30 min) | 64 |
6 | 6a | LAG (H2O/MeOH) Ball milling 30 Hz, 90 min | 53 |
7 | 6b | Ball milling 30 Hz, 90 min | 49 |
8 | 7 | Ball milling 30 Hz, 90 min | 59 |
9 | 8 | Ball milling 30 Hz, 90 min | 68 |
10 | 9 | Ball milling 30 Hz, 90 min | 55 |
When 1.5 mol% (relative to alkyne) of copper C-scorpionate 6a was used as a catalyst, the triazole was obtained in 78% yield (entry 4, Table 2), pointing out the capability of the copper C-scorpionate as a catalyst and its ability to increase the reaction rate.
E-Factor | EcoScale | Yield (%) | |
---|---|---|---|
Tris-pyrazole ligand 3a | |||
Mechanochemistry (this work) | 5.7 | 51.5 | 53 |
Pinho e Melo method23 | 68.4 | 46.5 | 87 |
Copper C-Scorpionate 6a | |||
Mechanochemistry (this work) | 0.4 | 72.5 | 79 |
Pettinari method29 | 320.8 | 34.5 | 43 |
CuAAC reaction for the synthesis of triazole 13 | |||
Mechanochemistry (this work) | 0.5 | 46 | 78 |
CuAAC reaction for the synthesis of triazole 1322 | 29.0 | 35.5 | 79 |
Martins and co-workers have successfully used copper(I) scorpionates and copper(II) scorpionates as catalysts for the three-component CuAAC reaction under microwave irradiation.22 In that study, it was observed that under similar reaction conditions (solvent: H2O/MeOH (1:1), MW, 125 °C, 30 min) but using 6a as the catalyst, triazole 2-(1-benzyl-1H-1,2,3-triazol-4-yl)ethanol was obtained in 64% yield (entry 5, Table 2). Knowing the effect that oxidizable alcohols and water could have in the reduction of Cu(II), a H2O/MeOH (1:1) mixture was chosen to perform the reaction under liquid-assisted mechanochemistry conditions. Using 15 μL of H2O/MeOH (1:1) and milling for 90 min at 30 Hz, the desired triazole 13 was obtained in 53% yield (entry 6, Table 2), indicating that neither water nor methanol has influence on the catalytic cycle, reinforcing the oxidative homocoupling of terminal alkynes as the main mechanism to form the active Cu(I) catalytic species.
Although acetate was previously described as a co-ligand with a good compromise between basicity and non-nucleophilicity, favoring metalation of the terminal alkyne without disfavoring the formation of the dinuclear intermediate and being very efficient for the protodemetallation step,35 the C-scorpionate complex 6b was tested as a catalyst of the reaction model (Scheme 4). The substitution of acetate by chloride resulted in a significant decrease in reaction yield, dropping to 49%. The methyl group at the C-4 of pyrazole also decreased the efficiency of the catalyst (entry 8, Table 2). The copper complex of the unsubstituted bis-pyrazole oxime ligand was found to produce good reaction yields (entry 9, Table 2), while the presence of the triazole unit at the C-scorpionate, complex 9, decreases the yield of the catalyzed reaction (entry 10, Table 2). The results show the influence of the structure of the ligand in the complex catalytic activity, where the increase of the size of the pro-ligand leads to a decrease in the activity, while the presence of unsubstituted pyrazole units leads to higher yields of triazole.
The CuAAC reaction under ball milling using the best reaction conditions (entry 4, Table 2) was further studied by carrying out the formal cycloaddition using aliphatic and aromatic alkynes, leading to triazole derivatives 13–18 in good yields, (Chart 2, see Fig. S20–S25 in the ESI† for structural characterization data).
Atom economy is determined using the chemical equation and is independent of experimental reaction conditions. While the ligand synthesis shows modest atom economy due to the use of 3 equivalents of sodium carbonate, the complexation reaction and the CuAAC reaction exhibit almost perfect and high atom economy, respectively (Fig. 6).
The values obtained for the E-factor and the EcoScale are presented in Table 3. The data used to estimate both metrics are presented in the ESI† (Tables S4–S7). For comparison, previously described methods for the synthesis of the same compounds or structurally similar compounds (in the case of the new compounds) were chosen.
The E-factor is the ratio of the mass of waste per mass of product: E-factor = (total mass of waste/mass of the product). Therefore, the lower the value, the more sustainable the process can be considered. The EcoScale is a post-synthetic tool that evaluates other aspects of the synthetic process, namely, toxicity of the reactants, energy input (time and temperature of the reaction) and procedure safety (see Table S4 in the ESI†). All the procedures to be analyzed start with 100 points, and the points are subtracted as each of the evaluated parameters deviates from the ideal of sustainability. Therefore, the higher the value, the more sustainable the process will be.
The synthesis of the pro-ligands, and therefore of the catalyst, requires the prior preparation of trichloroacetaldehyde oxime. The latter was prepared under conventional conditions following a known procedure. The E-factor of this reaction is 6, a good ratio of waste/product, and the EcoScale is 47.5, (ESI,† Table S5) a moderate value, far from the 100 points ideals for this metric, mostly due to safety issues of the reactants used.
The synthesis of tris-pyrazole ligands under free-solvent ball milling conditions (this work) yields the desired compound in lower yield than the synthesis in solution, using dichloromethane as a solvent as described by Pinho e Melo et al.23 However, the E-factor for the mechanochemical transformation (E-Factor = 6) is ten times lower than that obtained for the synthesis in solution (E-Factor = 68), which reflects the influence of the mass of solvent in the reaction waste. The almost equal EcoScale values achieved for both methods (51.5 and 46.5 for mechanochemical and solution method, respectively) show that the decrease of the reaction yield under ball milling is compensated from the point of view of sustainability by the elimination of the chlorinated solvent (ESI,† Table S6).
The preparation of the copper(II) C-scorpionate complexes under mechanochemical conditions is highly effective, and the use of equimolar amounts of the copper salt and the absence of a solvent led to an almost ideal E-factor value (E-Factor = 0.4) reflecting the low amount of waste obtained from the reaction and far from that obtained when the method in solution previously described is used to synthesize a similar copper complex (E-Factor = 320.8).29 A high value of EcoScale was also obtained for the complex synthesis under ball milling (72.5), close to the ideal 100 points, which are unachievable mostly due to the inherent toxicity of the copper acetate salt (ESI,† Table S7).
The three-component AAC reaction using copper(II) C-scorpionate complexes as catalysts gave similar yields under mechanochemical action and under microwave irradiation, however, the synthesis of the triazole under ball milling gave an impressive E-Factor value of 0.5, an improved performance when compared to the reaction carried out under microwave irradiation (29.0), in which the amount of solvent (MeOH/H2O) was already greatly reduced. The increase of sustainability using the mechanochemical conditions is also evidenced in the recovery of the EcoScale values from 35.5 to 46 points, reflecting the elimination of the solvent and the formation of sodium bromide as only secondary product (ESI,† Table S8).
Overall, the reported results show that the development of synthetic methods under mechanochemical conditions can reduce the use of reactants in excess, eliminate solvents, and maintain the efficiency of the synthetic process. This reduction in waste generation is reflected in the obtained sustainability metrics.
Overall, the development of synthetic methods using mechanochemistry in combination with catalysis offers advantages such as reduction of the use of excess reactants, the elimination of solvents and an increase in the efficiency of the catalytic process. These benefits are reflected in the sustainability metrics obtained for the developed processes.
All commercially acquired reagents were high-grade chemicals and were used without any additional purification: pyrazole (288-13-1, Fluorochem, 99%), 4-methylpyrazole (7554-65-6, Fluorochem, 97%), sodium carbonate anhydrous (497-19-8, Merck, 99%), copper(II) acetate monohydrate (142-71-2, J. T. Baker), sodium azide (26628-22-8, Acros Organics, 99%), benzyl bromide (100-39-0, Alfa Aesar, 99%), 3-butyn-1-ol (927-74-2, Alfa Aesar, 98%), ethynyltrimethylsilane (1066-54-2, Merck, 98%), 1,1-diphenyl-2-propyn-1-ol (3923-52-2, Sigma-Aldrich, 99%), phenylacetylene (536-74-3, Alfa Aesar, 98+%), 4-ethynyltoluene (766-97-2, Sigma-Aldrich, 97%), 3-butyn-2-ol (2028-63-9, Fluorochem, 97%).
Footnote |
† Electronic supplementary information (ESI) available. CCDC 2096452. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d3nj05017c |
This journal is © The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2024 |