Ligand-Mediated Regioselective Rhodium-Catalyzed Benzotriazole-Allene Coupling: Mechanistic Exploration and Quantum Chemical Analysis.

The ligand-controlled rhodium-catalyzed regioselective coupling of 1,2,3-benzotriazoles and allenes was investigated by DFT calculations. Because allylation can occur at either the N1 or N2 position of the 1,2,3-benzotriazole, the complete Gibbs free energy profiles for both pathways were computed. A kinetic preference emerged for the experimentally observed N1 allylation with the JoSPOphos ligand, whereas N2 allylation was favored with DPEphos. Analysis of the regiodetermining oxidative addition step by using the activation strain model in conjunction with a matching energy decomposition analysis has revealed that the unprecedented N2 reaction regioselectivity is dictated by the strength of the electrostatic interactions between the 1,2,3-benzotriazole and the rhodium catalyst. The nature of the electrostatic interaction was rationalized by analysis of the electrostatic potential maps and Hirshfeld charges: a stabilizing electrostatic interaction was found between the key atoms involved in the oxidative addition for the N2 pathway, analogous interactions are weaker in the N1 case.

Controlling the regioselectivity of fundamental organic transformations is an outstanding challenge in synthesis.1 The N‐selective functionalization of 1,2,3‐benzotriazoles is a particular challenge due to the dynamic equilibrium between two tautomeric forms of the benzotriazole (Scheme 1 a).2 Despite of the greater stability of tautomer 1 a, alkylation of 1,2,3‐benzotriazole has typically been observed to be moderately selective at N1/N3 position and extremely unselective at N2.

Scheme 1

(a) Tautomeric equilibrium of 1,2,3‐benzotriazole; (b) Synthetic protocol of regioselective Rh catalyzed addition of 1,2,3‐benzotriazoles to allenes.

Major advances have been forged for N1‐selective methods.3 N2‐selective modification is underdeveloped, primarily because of difficulties in shifting the equilibrium toward 1 b tautomer. As a consequence, the structural space around valuable N2‐substituted 1,2,3‐benzotriazoles that often occur in pharmaceutical agents, agrochemicals and advanced materials is limited.3, 4

The rare examples of N2 procedures reported during the last decade mostly rely on structural modification of the coupling partners.3 These methods are usually based on sterically bulky groups on the n class="Chemical">benzotriazole or alkylating agent that block access to the terminal nitrogens N1/N3, thereby making the N2 position the preferred site of nucleophilic attack. Despite being practically successful, these strategies are restricted in functional‐group tolerance and require prefunctionalized reactants. Transition‐metal catalysis offers an alternative approach that circumvents these previously mentioned issues as it involves the direct activation of X−H bonds for functionalizing molecules. In 2014, we proposed an unprecedented rhodium (Rh) catalyzed allylation of 1,2,3‐benzotriazoles. This methodology was attractive in the sense that it gave high N1/N2 regioselectivity and had a high tolerance toward a broad range of 1,2,3‐benzotriazole and allene derivatives.5 Importantly, the N1 or N2 regioselectivity could be tuned by judicious choice of the ligand on the Rh center. Among a series of diphosphine ligands used for the initial screening, we found that JoSPOphos (L1) gave the preference of N1 product, whereas DPEphos (L2) gave exclusive formation of the desired N2 substituted benzotriazole (Scheme 1 b). Intrigued by this fact, we performed additional experiments5 to probe the reaction mechanism and proposed the catalytic cycle illustrated in Scheme 2. The reaction involves three steps and begins with oxidative addition of benzotriazole tautomers 1 a/b to the RhI complex. The generated RhIII−H intermediates A and B proceed through hydrometalation and transform into σ‐allyl complexes C and D. Intermediates C or D then undergo reductive elimination with the contaminant release of N‐allylated products 4 and 3, respectively.

Scheme 2

Proposed reaction mechanism for Rh catalyzed coupling of 1,2,3‐benzotriazoles with allenes.

Despite the plausibility of the reaction mechanism outlined in the Scheme 2, the factors determining the unusual N2 regioselectivity remain unclear. This fact has motivated us to study the reaction mechanism using density functional theory (DFT) calculations. To provide general concepts for the rational design of other complementary approaches for the regioselective synthesis of valuable n class="Chemical">N‐substituted synthons, the following questions must be addressed: 1) What is the reaction mechanism? 2) Which step of the catalytic cycle is regiodetermining? 3) What are the physical factors behind the regioselectivity? To answer these questions, we have performed a systematic in‐depth theoretical investigation based on the accounts of the previously reported experimental data.5

The coordination of Rh to L1′ was initially investigated to determine the active form of the catalyst. Figure 1 illustrates the two possible binding modes of the JoSPOphos ligand (L1′) to the Rh center.6 The large difference in computed Gibbs free energies (ΔΔG=17.4 kcal mol−1) results from a favorable P‐type of binding in Rh‐L1′. Having determined the structure of the catalyst, we then investigated the reaction mechanism.

Figure 1

Coordination modes of L1′ to the Rh. To reduce computational cost the tert‐butyl group of the ligand (L1) was replaced with methyl (L1′). Relative energies [kcal mol−1] were computed at SMD(DCE)‐B3LYP‐D3/6‐31G(d)/LANL2DZ.

All relevant mechanistic pathways (for details, see the Supporting Information, Figure S1) were computationally7 explored and the Gibbs free energy profiles of the most favorable catalytic process were constructed (Scheme 3 a). DFT calculations revealed a kinetic preference for N1 allylation (black color pathway) over N2, which is in line with experimental results of N1‐substituted benzotriazole (3) being the major product when JoSPOphos is employed.

Scheme 3

(a) Gibbs free energy profile for coupling of the 1,2,3‐benzotriazoles with allenes catalyzed by Rh‐L1′ at 80 °C and a standard state of 1 mol L−1; (b) Regioselectivity determining transition states; (c) Rate‐determining transition states. Bond lengths in Å, Gibbs free energies in kcal mol−1. To reduce computational cost the cyclohexyl group of allene (2) was replaced by the methyl group (2′) and tert‐butyl group of the ligand (L1) was replaced with methyl (L1′). All data computed at SMD(DCE)‐B3LYP‐D3/6‐311++G(d,p)/SDD//SMD(DCE)‐B3LYP‐D3/6‐31G(d)/LANL2DZ.

The reaction is initiated by coordination of the substrate 1 a/b to Rh‐L1′ to form the pre‐reactant complexes n class="Gene">INT1 and INT2. Notably, the hydroxy group on the Rh‐L1′ catalyst acts as a directing group for the incoming benzotriazole and results in an OH⋅⋅⋅N1/N2 hydrogen‐bonding interaction. Either the N1−H (1 a tautomer) or N2−H bond (1 b tautomer) can undergo oxidative addition depending on the tautomeric form of the benzotriazole. The oxidative addition involves the classical 3‐membered TSs (TS1‐3, TS2‐4) with relative free energies of 20.3 and 19.2 kcal mol−1. Computed structures of TS1‐3 and TS2‐4 are illustrated in Scheme 3 b and reveal a late, product‐like character for TS1‐3 with a nearly formed Rh−N bond at 2.36 Å and an earlier TS structure for TS2‐4, in which the Rh−N bond is still relatively long 3.04 Å. Because formation of INT3 and INT4 is downhill with respect to INT1, INT2 by 15.4 kcal mol−1 in both cases, the oxidative addition step is predicted to be irreversible. To provide a coordination site necessary for the subsequent hydrometalation of the allene, the Cl ligand migrates to an axial position, whereas the allene binds to the Rh via an external C=C double bond in the equatorial plane. This ligand transformation raises energies to −1.0 and 0.6 kcal mol−1 to form INT5 and INT6, respectively. The catalytic cycle continues with the hydrometalation and passes over transition states TS5‐7, TS6‐8 with Gibbs free activation barriers of 8.9 and 4.9 kcal mol−1, respectively. The resulting intermediates INT7, INT8 are asymmetrical π‐type complexes located at −11.9 and −15.3 kcal mol−1, which further isomerize to more stable σ‐allyl intermediates8 INT9, INT10 with relative free energies of −21.8 and −23.7 kcal mol−1. The reductive elimination traversing transition states TS9‐11, TS10‐12 and affords the product complexes INT11, INT12. It is worth pointing out that the calculated activation barrier of reductive elimination for the N1‐selective pathway (Scheme 3 a, in black) is significantly lower (ΔΔG ≠ =5.4 kcal mol−1) compared to the N2 route. The optimized transition structures (Scheme 3 c) reveal a stabilizing hydrogen‐bonding interaction in TS10‐12 that is absent in TS9‐11, which may be the origin of the lower barrier for the former. Release of products 3′ and 4′ allows the Rh catalyst to reenter the catalytic cycle and bind a new substrate. The Gibbs free energy profile provided in Scheme 3 a reveals that the oxidative addition is the regioselectivity determining step, while the reductive elimination is the rate‐determining step. Our DFT calculations support the experimental observation5 of N1‐allylated regioisomer (3′) being a preferred product as a TS2‐4 associated with N1‐selective path is more stable than a competitive TS1‐3 for 1.1 kcal mol−1.

To understand why simply changing the JoSPOphos ligand (L1) with DPEphos (L2) completely switches the regioselectivity from n class="Chemical">N1 to N2, we again sampled various mechanistic possibilities (see the Supporting Information, Figures S2, S3) and outlined the lowest‐energy reaction profiles corresponded to N1 and N2 pathways (Scheme 4 a). In contrast to L1′ (Scheme 3 a), now the N2 reaction channel (red curve, Scheme 4 a) is kinetically favored, which is in line with the experimental regioselective allylation at N2 atom of the benzotriazole.

Scheme 4

(a) Gibbs free energy profile for coupling of the 1,2,3‐benzotriazoles with allenes catalyzed by Rh‐L2 at 80 °C and a standard state of 1 mol L−1; (b) Regioselectivity determining transition states. Bond lengths in Å, Gibbs free energies in kcal mol−1. To reduce computational cost the cyclohexyl group of allene (2) was replaced by the methyl group (2′). All data computed at SMD(DCE)‐B3LYP‐D3/6‐311++G(d,p)/SDD//SMD(DCE)‐B3LYP‐D3/6‐31G(d)/LANL2DZ.

The first elementary step of the catalytic cycle (Scheme 4 a), namely, oxidative addition takes place from INT13 or INT14, in which the substrates 1 a or 1 b are pre‐coordinated to the catalyst via a weak hydrogen‐bond interaction between N1/2−H⋅⋅⋅Cl. From the substrate‐coordinated complexes INT13 and INT14, the N1/N2−H activation is taking place via TS13‐15/TS14‐16 located at 20.9 and 22.8 kcal mol−1. Optimized structures of TS13‐15, TS14‐16 are provided in Scheme 4 b and are late with respect to the N−H bond breaking, displaying N−H distances of approximately 2.0 Å. The RhIII hydride intermediates (INT15, INT16) furnished after the oxidative addition adopt a pentacoordinate structure bearing a hydride in the apical and benzotriazole fragment in the equatorial positions. The formation of these species is exergonic relative to the pre‐reactive complexes (INT13, INT14) by 1.9 and 5.5 kcal mol−1, respectively, suggesting oxidative addition to be irreversible. Binding of the incoming allene 2′ during the second stage of the catalytic cycle is uphill in terms of Gibbs free energy by 9.9 (INT17) and 10.3 (INT18) kcal mol−1. During this process, the coordination geometry at the Rh center changes from square pyramidal to octahedral, in which Cl and H ligands occupy apical positions. The reaction continues with a fast hydrometalation process via TS17‐19, TS18‐20 with activation free energies of 6.0 and 6.9 kcal mol−1 to give σ‐allyl complexes (INT19, INT20). Finally, the hydrometalation adducts undergo reductive elimination with barriers of 13.1 and 13.4 kcal mol−1 (TS19‐21 and TS20‐22), and deliver the product complexes INT21 and INT22. The catalytic cycle completes upon dissociation of the N‐allylated products 3′ or 4′ from the coordination sphere of the metal with the regeneration of the catalyst. Considering the full Gibbs free energy profile illustrated in Scheme 4 a, the key step determining the regioselectivity of the allylation is the oxidative addition step. Our DFT results show TS13‐15 to be more stable than TS14‐16 by 1.9 kcal mol−1. This ΔΔG ≠ between the regiodetermining TSs corresponds to a ratio of N1/N2 4:96 with predominance of N2 allylated regioisomer 4′, which is in perfect agreement with the experimentally isolated ratio of N1/N2 6:94.

To understand why DPEphos (L2) leads to N2 regioselectivity, we quantitatively analyzed the regiodetermining oxidative addition step catalyzed by n class="Chemical">Rh‐L2 using the activation strain model ASM (Figure 2).9 Within the ASM, the relative energy in solution ΔE solution along the reaction energy profile is separated into the energy of the solute ΔE solute (reaction system in vacuum with the solution phase geometry) and the solvation energy ΔE solvation[Eq. (1)].

Figure 2

(a) Activation strain analyses and (b) energy decomposition analyses for the oxidative addition of 1,2,3‐benzotriazole 1 a (N2‐selective pathway) or 1 b (N1‐selective pathway) to Rh‐L2 computed at ZORA‐B3LYP‐D3/TZ2P level. Dots represent TSs. Electrostatic potential maps along with the Hirshfeld charges (m a.u.) at N1, N2, N3 atoms of the 1,2,3‐benzotriazole fragment (1a/1b) in both equilibrium (1a/1b) and transition state (TS13‐15/TS14‐16) geometries correspond to (c) N2‐ and (d) N1‐selective pathways. ESP maps are plotted on the total electron density of the 1,2,3‐benzotriazole fragment from the ZORA‐B3LYP‐D3/TZ2P calculations by using a consistent surface potential range of −0.03 a.u. to 0.30 a.u. and isovalue 0.045.

The intrinsic energy of the solute ΔE solute is split further into two terms: ΔE solute,strain (energy required for deformation of fragments from their starting geometries to the geometries they obtain over the course of the reaction) and ΔE solute,int (energy of interaction between deformed fragments) [Eq. (2)].

Next, the interaction energy ΔE solute,int is decomposed into three terms using a canonical energy decomposition analysis (EDA) [Eq. (3)]:

in which ΔV solute,elstat corresponds to the electrostatic interaction between unperturbed charges of the deformed reactants, ΔE solute,Pauli accounts for the repulsion between occupied orbitals, ΔE solute,oi is responsible for interaction between occupied and unoccupied orbitals and polarization. The ASM and EDA terms were projected on the N1/N2−H bond stretch, as this geometrical parameter is critically involved in the reaction and undergoes a well‐defined change9 over the course of the oxidative addition step.

Figure 2 a shows that the interaction energies (ΔE solute,int) control the reaction regiochemistry, whereas the strain energy required for deformation of catalyst (Rh‐L2) and substrate (one of the two 1,2,3‐benzotriazole tautomers, 1 a/b) during oxidative addition (ΔE solute,strain) is nearly identical for both N1‐ (black) and N2‐selective (red) pathways. Similar ΔE solute,strain curves is a consequence of breaking the same N−H bond at either N1 or N2 during the oxidative addition. Because the ΔE solute,int is decisive in determining the trend in ΔE solute and thus the observed reactivity trends, it was further analyzed by using the EDA, and the results are plotted in Figure 2 b. Differences in the orbital‐interaction curves (ΔE solute,oi) are minimal. Thus, it becomes evident that the predominance of the N2‐selective pathway can be attributed to the more stabilizing electrostatic interaction (ΔV solute,elstat), which effectively overrules the less destabilizing Pauli repulsion (ΔE solute,Pauli) preference for the N1 pathway.

Next, to understand the trend in ΔV solute,elstat, we examined the electrostatic potential maps (ESP) and Hirshfeld charges10 of Rh and N1, N2, N3 atoms of 1,2,3‐benzotriazole fragments (1 a/b). Figure 2 c and d illustrate the ESP and charge analysis for the 1,2,3‐benzotriazole tautomers at their equilibrium and TSs geometries. Our analysis has identified a causal relationship between the charge of the N atom of the benzotriazole involved in the regiodetermining oxidative addition step and the degree of electrostatic stabilization: a more electronegative N atom leads to a stronger, more stabilizing electrostatic interaction with the electropositive Rh metal center. Thus, the more reactive tautomer 1 a has an N1 atom positioned between a carbon atom and nitrogen N2. This results in a relatively electronegative N1 (−36 m a.u.), which leads to a strongly stabilizing electrostatic interaction with Rh (20 m a.u.) metal center (Figure 2 c). In contrast, the less reactive tautomer 1 b has a partially positive (44 m a.u.) N2 atom (Figure 2 d) as a result of being positioned between two nitrogens (N1 and N3 atoms). This leads to less stabilizing electrostatic interaction between the N2 atom and the electropositive Rh. The electrostatic interactions are more stabilizing for the reaction of N2 compared with N1 over the entire reaction coordinate (Figure 2 b), thus this charge and molecular electrostatic potential (MEP) analysis can be performed at any point.

In summary, we have computationally analyzed the reaction mechanism of a newly developed coupling of 1,2,3‐benzotriazoles with n class="Chemical">allenes catalyzed by Rh‐diphosphine ligand (JoSPOphos/DPEphos) complexes. Our DFT calculations reveal that regardless of the ligand being used, the reaction proceeds via a three‐step catalytic cycle, comprised of first an oxidative addition, then hydrometalation, and finally a reductive elimination. Our activation strain and energy decomposition analyses on the regiodetermining oxidative addition step catalyzed by Rh‐DPEphos show that the preference for the N2‐selective reaction channel arises from a previously unrecognized electronic mechanism, namely, a more stabilizing electrostatic interaction between the electron‐enriched region at the N1 atom of tautomer 1 a with the positively charged Rh during N1−Rh bond formation in the rate‐determining N1−H oxidative addition step. We envisage that the newly identified electrostatic interactions can be used to further tune the 1,2,3‐benzotriazole and catalyst interaction through the introduction of various functionalities, which may allow new reactions with tailored N2 regioselectivity trend.

Conflict of interest

The authors declare no conflict of interest.

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

Supplementary

Click here for additional data file.