Mechanism of the Molybdenum-Mediated ... - ACS Publications

Jun 28, 2018 - time.1−5 Traditionally, most of the interest surrounding this kind of oxygen .... molybdenum-catalyzed OAT reaction between Ar-NO2 an...
1 downloads 0 Views 1MB 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 Cite This: ACS Omega 2018, 3, 7019−7026

Mechanism of the Molybdenum-Mediated Cadogan Reaction Marta Castiñeira Reis,† Marta Marín-Luna,† Carlos Silva López,*,† and Olalla Nieto Faza*,‡ †

Departamento de Quı ́mica Orgánica, Universidade de Vigo, Campus Lagoas-Marcosende, 36310 Vigo, Spain Departamento de Quı ́mica Orgánica, Universidade de Vigo, Campus As Lagoas, 32004 Orense, Spain



Downloaded via 80.82.77.83 on June 30, 2018 at 23:53:52 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

S Supporting Information *

ABSTRACT: Oxygen atom transfer reactions are receiving increasing attention because they bring about paramount transformations in the current biomass processing industry. Significant efforts have therefore been made lately in the development of efficient and scalable methods to deoxygenate organic compounds. One recent alternative involves the modification of the Cadogan reaction in which a Mo(VI) core catalyzes the reduction of o-nitrostyrene derivatives to indoles in the presence of PPh3. We have used density functional theory calculations to perform a comprehensive mechanistic study on this transformation, in which we find two clearly defined stages: an associative path from the nitro to the nitroso compound, characterized by the reduction of the catalyst in the first step, and a peculiar mechanism involving oxazaphosphiridine and nitrene intermediates leading to an indole product, where the metal catalyst does not participate.



the initial step of the mechanism36 diverging from previous studies on similar reactions.14,23 In 2007, Sanz et. al. reported a promising new synthetic method for the preparation of complex organic molecules bearing an indole,37 a scaffold with a strong relevance to medicinal and pharmaceutical chemistry.38−40 This procedure is becoming very popular to prepare indole derivatives41−43 and it involves the modification of the classical Cadogan reaction in which a MoO2Cl2(dmf)2 complex catalyzes the reduction of o-nitrostyrene derivatives in the presence of PPh3 under thermal conditions (see Figure 1).44,45 Sanz et. al.

INTRODUCTION Heteroatom transfer reactions (particularly those involving oxygen) have captivated the scientific community for a long time.1−5 Traditionally, most of the interest surrounding this kind of oxygen transfer has pivoted around the development of reactions oxygenating hydrocarbons (epoxydations, dihydroxylations, etc.), with the main goal of increasing the functionalization of simple petrol-derived building blocks.6−10 The interest in the generation of renewable feedstocks for the chemical industry has prompted attention to the reverse procedure, deoxygenation of oxygen-rich biomass sub(products), to produce highly reduced building blocks, analogous to those obtained from petrol sources.11−17 In the development of oxygen transfer reactions, nature has provided us with a model: metalloenzymes.18−26 Notable key examples of these enzymes are the high-valent manganese species implicated in the photosynthetic water splitting,27,28 or cytochrome P-450 as the responsible agent for the metabolism of steroids.29 Among these metalloenzymes, molybdenum is found in different types of oxotransferases sharing a common active center that consists of one Mo(VI) atom with at least one oxo ligand.30,31 Synthetic molybdenum systems mimicking the chemical behavior of these oxotransferases have been reported, the vast majority involving penta- and hexacoordinated molybdenum complexes featuring sulphurated bidentate ligands.25,32,33 There is limited information on systems in which the coordination sphere is reduced to four ligands. In all these systems, the [MoO2]2+ active core has been only taken into account considering its oxidizing power and more scarcely as a Lewis acid.34,35 On the basis of this, we have recently published a theoretical study on the deoxygenation of dimethylsulfoxide catalyzed by MoO2Cl2 that showcases the unusual behavior of the molybdenum core as a Lewis acid in © 2018 American Chemical Society

Figure 1. Conversion of a general o-nitrostyrene scaffold into its indole derivative via a molybdenum-catalyzed reduction.

proposed that molybdenum only partially catalyzes the overall process exerting its catalytic activity during the reduction of the nitro to nitroso group. The subsequent reductive cyclization process furnishing the indole is considered to evolve without the participation of the catalyst. In light of the potential of this reaction in the synthesis of bio-active compounds41,43 and its mechanistic uncertainties,46−50 we decided to study the mechanism of the Momediated deoxygenation of ortho-nitrostyrene in the presence Received: June 7, 2018 Accepted: June 13, 2018 Published: June 28, 2018 7019

DOI: 10.1021/acsomega.8b01278 ACS Omega 2018, 3, 7019−7026

Article

ACS Omega of a phosphine as a reducing agent.37,51 The possible initial roles of [MoO2]2+ as a Lewis acid and as an oxidizing agent have been explored, resulting in five alternative reaction paths. The subsequent transformation of the transient nitroso derivative into the final indole has also been analyzed and several mechanistic questions on this second part of the mechanism have been addressed. In our recent work, we have found that similar chemistry operates through multisurface pathways, involving spin-crossing events along the reaction mechanisms.14,36 We have, therefore, taken particular account of the spin multiplicity of the stationary points along the studied reaction mechanisms and when changes in spin were present, the minimum energy crossing points (MECP) on the seam formed between the crossing potential energy surfaces were computed to obtain a qualitative picture of the probability of such forbidden process to take place (for a quantitative picture spin−orbit coupled density functional theory should be considered).52



and OPMe3. Our calculations show that the triplet state in this complex (2) is 7.3 kcal/mol more stable than the singlet, suggesting that a switch between the singlet and triplet potential energy surfaces may occur along the reaction path during the evolution from 1 to 2. In this sense, we have located a MECP near 2 (see Figure 2) analogous to that described for the Mo(VI)-mediated reduction of sulfoxides with phosphine.36 The initial reduction of the catalyst therefore involves a barrier of 23.6 kcal/mol and a late post-transitional spin flip process. The ulterior release of the OPMe3 ligand trough TS2−3 (ΔG‡ = 7.7 kcal/mol) results in the formation of a trigonal complex, 3 (5.7 kcal/mol), which features an extra coordination vacancy prompted to be filled by nitrostyrene with the formation of complex 4 (through a barrier of 10.8 kcal/mol). This energetically uphill path leading to 4, further evolves toward the more stable chelate 5 via a facile coordination of the dangling oxygen atom of the nitro group to the metal center (see Figure 2). The difference in energy between the triplet and singlet in 5 is here reduced to just 0.13 kcal/mol from 2.27 in 4, indicating that a crossing between the singlet and triplet surfaces may lie nearby. Then, an oxygen atom transfer reaction from the nitro group to the molybdenum center takes place through an N−O bond cleavage (TS5−6) involving an activation barrier of 19.3 kcal/ mol. The ground state of this transition state is singlet, confirming that the spin switch has occurred before the transition state. A low-energy crossing point could actually be located between 5 and TS5−6 (see Figure 2). In the last step of this catalytic cycle, the departure of the reduced nitroso ligand allows the recovery of the catalyst through TS6−1 in an exergonic process (−24.7 kcal/mol). Using the energetic span model,61 we found that the rate determining stationary points of this catalytic cycle are 2 (as the turnover determining intermediate) and TS3−4 (as the turnover determining transition state) with an energy span of 28.0 kcal/mol. This value is in fair agreement with the thermal conditions required in this reaction (see Figure 1). The uphill energy values found for the tricoordinated Mo(IV) species 3 and the region around it, led us to reconsider the reactivity of complex 2. This complex, which still displays a coordination vacancy, could behave as a Lewis acid, thus, opening the door to alternative mechanisms, as depicted in Figure 3. The Lewis acid reactivity of polyoxomolybdenum(VI) compounds in redox processes is not usually taken into account, although there is some literature on this topic.34−36,62,63 In this context, we found it relevant to also explore the possibility of the catalyst behaving as a Lewis acid before (complex 1 in Figure 2) and after its reduction to Mo(IV) has taken place (complex 2 in Figure 2). Figure 3, therefore, contains the complex network of reaction pathways available for both complexes when they react as Lewis acids. The top section (blue) highlights the chemistry in the Mo(VI) state, whereas the bottom section (yellow) contains the chemistry in the Mo(IV) state. A manifold of steps connecting both regions via PMe3 addition or OPMe3 release has also been considered. The alternative pathways starting from Mo(VI) complex 1 (see Figure 3, blue region) can be classified in terms of how late the reduction of molybdenum occurs along the path. Hence, how late they merge with the Mo(IV) path (yellow section in Figure 3). The first step, common to all these paths, involves the direct coordination of nitrostyrene to the metal center, yielding 11

RESULTS AND DISCUSSION

Different mechanisms can be envisaged for the molybdenumcatalyzed monodeoxygenation of the nitro group in the presence of PMe3. We first followed the classical approach for the oxygen atom transfer (OAT) reaction followed in nature,23,53 that is: the reduction of the catalyst as the first step, followed by the generation of a coordination vacancy on the metal center and a subsequent coordination of the substrate. The reduction of the molybdenum catalyst is usually proposed to take place through the nucleophilic attack of PMe3 on the π*Mo−O orbital in the Mo(VI) catalyst 1, leading to the formation of Mo(IV) complex 2 (see Figure 2).36,54−57 The P−O distance found in 2 (dP−O = 1.56 Å) is very similar to that found in isolated phosphine oxide compounds (dP−O = 1.52 Å),58−60 indicating a weak interaction between the metal

Figure 2. Enzymatic-type path: proposed mechanism for the molybdenum-catalyzed OAT reaction between Ar-NO2 and PMe3, considering the reduction of the catalyst as the first step. Red color indicates structures with triplet ground states, whereas black color is used for structures in the singlet state. Computations were performed at the UB3PW91/Def2-SVPD (polarizable continuum model (PCM), toluene) theoretical level. Gibbs free energies are reported in kcal/mol (1 atm and 298 K), relative to 1 + PMe3 + Ar-NO2. 7020

DOI: 10.1021/acsomega.8b01278 ACS Omega 2018, 3, 7019−7026

Article

ACS Omega

Figure 3. Lewis acid-type path: available reaction pathways for the molybdenum-catalyzed OAT reaction between Ar-NO2 and PMe3 considering the Lewis acid character of the catalyst. Mo(VI) pathways in blue (top) and Mo(IV) paths in yellow (bottom) with the connecting manifold between the two alternatives as vertical arrows on a blank background. Entry/exit points to the traditional mechanism illustrated in Figure 2 and are indicated as boldface structures (red color for triplets and black color for singlets as in Figure 2). Computations were performed at the UB3PW91/ Def2-SVPD (PCM, toluene) theoretical level. Gibbs free energies in kcal/mol (1 atm and 298 K) relative to 1 + PMe3 + Ar-NO2 are shown.

(11.0 kcal/mol). The long Mo−ONO bond distance, 2.54 Å, and the invariance of the bond length between the oxygen and nitrogen atoms dO−NO = 1.23 Å with respect to that of the free nitro-system suggest a weak coordination between this ligand and the [MoO2]2+ core. This may be the result of the strong trans effect caused by the oxo ligand.64a Starting from intermediate 11, we found three different reduction paths: 1. Path A: this path avoids altogether merging with the Mo(IV) mechanism. The reduction step occurs via a direct attack of PMe3 on an oxygen of the nitro group (instead of attacking an oxo ligand on the metal center), which leads to intermediate 13 (Figure 3). This intermediate can then evolve toward 6 through the release of OPMe3 with an almost barrierless transition state (0.6 kcal/mol). It could also be proposed that complex 13 releases the organic ligand, thus, connecting with the energetically more demanding noncatalyzed process (see the Supporting Information (SI) for more details). Once 6 has been formed, it evolves through a low-lying transition state (ΔG‡ = 4.4 kcal/mol) to the aryl nitroso compound recovering the active catalyst, 1. 2. Path B: the nitro ligand chelates the metal core through TS11−12 with a very low barrier (0.6 kcal/mol) yielding 12. This complex is subsequently reduced to 8 in the presence of PMe3, thus, entering the Mo(IV) pathway. This reduction step is however energetically very demanding (32.0 kcal/mol). 3. Path C: PMe3 attacks the π* MoO bond trans to the nitro ligand in 11 through TS11−7 (ΔG‡ = 21.2 kcal/

mol) and enters the Mo(IV) pathway (Figure 3, in yellow) with the formation of complex 7. The significant decrease of the Mo−ONO distance (from 2.54 Å in 11 to 2.00 Å in 7) and the increase of the O−NO distance (1.27 Å) with respect to the free ligand are associated with the impact of the softer trans effect created by the phosphine oxide. The selectivity found in TS11−7 can be rationalized on the grounds of the steric hindrance imparted by the nitro ligand and the presumed lower apicophilicity of the oxo ligand.65,66 (See the SI for further discussion on the stereoselectivity of the described processes.) The energy spans in these three mechanisms are 32.2, 32.0, and 30.9 kcal/mol for paths A, B, and C, respectively. These values are too high for these mechanisms to be competitive with the enzymatic-type mechanism provided in Figure 2 (28.0 kcal/mol). If we consider that the action of molybdenum as a Lewis acid takes place after reduction of complex 1 to 2, the first step would be the complexation of a nitrostyrene molecule to 2, leading to the formation of complex 7 (−8.1 kcal/mol, see Figure 3). The ground state for complex 2 is a triplet, whereas 7 shows near degeneracy between the singlet and triplet states (the latter lying only 0.6 kcal/mol above the former). A spin flip area should therefore be located near 7. Indeed, a minimum energy crossing point between the two surfaces, which qualitatively determines the activation barrier, can be located so close to 7 that it is isoenergetic with this minimum at the 0.1 kcal/mol precision level commonly employed in 7021

DOI: 10.1021/acsomega.8b01278 ACS Omega 2018, 3, 7019−7026

Article

ACS Omega density functional theory calculations (and also in most experimental kinetic measurements). Once complex 7 is obtained, the nitro ligand chelates easily, through a barrier of 3.5 kcal/mol, leading to octahedral complex 8 (−15.4 kcal/ mol). At this point, intermediate 8 can release a molecule of OPMe3 to furnish complex 5 (thus, connecting this path to that shown in Figure 2), but it can also suffer from (2 + 2) cycloreversion through TS8−9 by which the initial nitro group loses an oxygen atom (ΔG‡ = 23.8 kcal/mol) leading to 9. This last alternative is not only kinetically but also thermodynamically preferred, because 9 is a rather stable intermediate compared with 5 (−27.8 vs −5.1 kcal/mol, for 9 and 5, respectively, see Figure 3). Complex 9 can then evolve either through cleaving OPMe3, which yields 6 (−21.9 kcal/ mol) or through cleavage of the nitrosoaromatic ligand, yielding 10 (−31.3 kcal/mol), the latter being clearly favored. Finally, the departure of the OPMe3 ligand in 10 through TS10−1 recovers the initial catalyst which can start a new cycle. This last reaction pathway involving the initial reduction of the complex and proceeding via 9 has an energy span of 23.8 kcal/mol, and hence, is more favorable than either of the alternatives considered above, including the classical mechanism illustrated in Figure 2. These results suggest that the classical pathway, involving a dissociative mechanism and occurring via severely deficient species in terms of coordination number, should be replaced by an associative mechanism, first involving a reduction step and then the addition of nitrostyrene in a Lewis acid/base step. Once the deoxygenation of ortho-nitrostyrene is completed, the reduction and concomitant cyclization of the nitrosostyrene derivative to furnish the final indole is expected, nevertheless, the mechanism is far from being well-known. For this second step, we have discarded the participation of the metal as a catalyst taking into account previous experimental results.67 Cadogan described how 2-nitrosobiphenyl and other derivatives can be reduced to the corresponding carbazoles in the presence of phosphines or phosphites very rapidly even under mild conditions (a few minutes at 0−5 °C). This is indicative of low activation barriers in the absence of the catalyst. Our own exploratory calculations agreed with these results.b Some studies in the literature suggest a tandem cyclization/reduction process in which the indole is transformed first into an N-oxo-indole derivative and ultimately into a hydroxyindole.41,68,69 However, the high reaction barriers computed for the reduction of the former and the absence of a complete description for the reduction of the latter led us to discard these alternatives (see the SI for further discussion of these pathways) and to consider the approach described for the Cadogan−Sundberg reaction.70 In this approach the proposed order of the steps is inverted with respect to the tandem cyclization/reduction and a deoxygenation reaction is followed by a cyclization step (Figure 4). We propose that the reducing agent, PMe3, initially attacks the nitrogen atom of 14, through a barrier of just 13.1 kcal/mol, leading to 19, which can further evolve toward oxazaphosphiridine 20. The N−P bond in 20 is subsequently cleaved, allowing the formation of a nitrene intermediate, 22, via the loss of OPMe3. All our efforts to locate TS20−21 were unfruitful, although the chemistry associated seems relatively simple (cleaving a single N−P bond). To explore the potential energy surface in greater detail in this region, we decided to apply the nudged elastic band (NEB) method. This method provides the lowest energy

Figure 4. Proposed mechanisms for the transformation of nitroso compound 14 into indole 24. Computations were performed at the UB3PW91/Def2-SVPD theoretical level. Gibbs free energies in kcal/ mol (1 atm and 298 K) relative to 14 + PMe3 are shown.

pathway connecting two minima and it thus offers an approximated value for the activation energy associated with this bond cleavage.71 With this method (see the SI), we estimated the reaction barrier to be about 11 kcal/mol and found why a transition state could not be located with standard algorithms.c In this step the P−N bond cleavage is strongly coupled with a twist of the aromatic ring, as in an fouetté en tournant pirouette, whose ends rotated 180° with respect to its initial position (Figure 5). Then 21 releases a molecule of

Figure 5. Strong coupled motion in the P−N bond cleavage of 20 (See the Web enhanced object for motion picture). Ballerina picture reproduced with permission from GFYCAT.com (unknown artist).

OPMe3 providing a nitrene intermediate which can undergo cyclization. For the transformation of 14 into 22, we also explored the direct attack of the reducing agent on the oxygen atom on the nitroso group but this resulted in a noncompetitive path (see the SI, Figure S2). Cadogan proposed that the cyclization step from nitrene 22 could occur on the most stable triplet surface.67,70,72,73 We confirmed that the triplet state is indeed the ground state of the nitrene, lying 9.3 kcal/mol lower than the singlet state, but found that the barrier for cyclization in this spin state is incompatible with the low thermal requirements of this reduction/cyclization tandem (0−5 °C).67 If, on the contrary, 7022

DOI: 10.1021/acsomega.8b01278 ACS Omega 2018, 3, 7019−7026

Article

ACS Omega

have checked that the reaction barrier for the rate determining step is not severely affected by this simplification: the reaction barrier obtained for TS1−2 with PPh3 is 34.11 kcal/mol and the reaction energy −7.83 kcal/mol at the UB3PW91/Def2-SVPP level in toluene. Although the catalyst used in the experimental conditions is a hexacoordinated species (MoO2Cl2(dmf)2, dmf = dimethylformamide), we have used the tetracoordinated MoO2Cl2 complex as the starting point of all our catalytic cycles. In previous work, with similar Mo(VI)-catalyzed oxygen transfer reactions, we analyzed the reactivity of MoO2Cl2(dmf)2, MoO2Cl2(dmf), and MoO2Cl2 and found that the latter is the most active catalyst.36 To simplify the reaction schemes, we have chosen a color code: the structures and data colored in red are those for which the most stable spin state is the triplet state and the data corresponds to this same state. On the other hand, those structures and data colored in black are those for which the singlet state is the most stable state and the data corresponds to this same state.

the reaction is considered to proceed adiabatically, nitrene 22 would evolve through cyclization via transition state TS22−23 in the singlet state, with an associated activation barrier of only 5.8 kcal/mol. This small activation energy is not only compatible with the experimental conditions, but also with the nitrene being such a fleeting intermediate so as to prevent it from relaxing electronically to the triplet state due to the spin forbidden nature of this process.



CONCLUSIONS In summary, we have disclosed here a detailed theoretical study on the reductive cyclization of ortho-nitrostyrene to furnish indole catalyzed by MoO2Cl2 and with the participation of a trialkyl-phosphine molecule as a reductant. Two series of mechanistic alternatives have been considered for the monodeoxygenation of the nitro to nitroso group mainly diverging on the role of the metal complex in the initial steps of the mechanism: as a redox agent or as Lewis acid. Our computations predicted that the most favorable mechanistic alternative is that in which the molybdenum(VI) catalyst acts as a redox agent rather than a Lewis acid in the first step and it immediately switches to a Lewis acid role for the second step. The catalyst is thus acting in a manner similar to that proposed for molybdoenzymes, but we have found that an associative mechanism is more favorable than the commonly proposed dissociative path. The subsequent conversion of the nitroso substrate to the final indole system occurs preferentially following a reduction/cyclization tandem process. In this second step the reductant is involved in a nucleophilic attack on the NO bond leading to an unusual oxazaphosphiridine intermediate. These results, in agreement with the previous hypothesis based on the experimental observations, reveal that the participation of a PMe3 molecule is a requirement in both sections of the mechanism (the nitro reduction and the reduction/cyclization tandem process) but molybdenum seems to be only necessary for the reduction of the initial nitro compound.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsomega.8b01278. Cartesian coordinates, electronic energies and the number of imaginary frequencies for all the structures reported in this work, additional schemes and discussion of alternative mechanisms that have been considered (PDF) Strong coupled motion in the P−N bond cleavage of 20 (MPG)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected]. Phone: +34 986813268 (C.S.L.). *E-mail: [email protected] (O.N.F.).



COMPUTATIONAL METHODS Geometries of all the stationary points were fully optimized by using the UB3PW91 functional74,75 with the Def2-SVPD basis set.76,77 The SDD(28, MWB) effective core potential associated with this basis set was used to describe the inner electrons of the Mo atom.78 The effect of the solvent (toluene) was modeled using the polarizable continuum model (PCM)79 with the default parameters implemented in the Gaussian 09 package.80 The Orca program81 was used, at the same computational level, to obtain the minimum energy crossing points between triplet and singlet surfaces when a change in multiplicity is needed along the reaction path. In these cases, the conductor-like screening model82 was applied to take into account the solvation effect. To check the accuracy of this approximation, we optimized some of the minima and checked their relative energies with both programs resulting in similar values, as expected. Harmonic analysis was used to establish the nature of all optimized structures as either minima or transition structures. For all stationary points, the stability of the wave function was also confirmed. To reduce the computational cost associated with both the extra electrons and a larger conformational space, we are using the smaller trimethylphosphine as the reducing agent, instead of the actual triphenylphosphine used in the experiments. We

ORCID

Carlos Silva López: 0000-0003-4955-9844 Olalla Nieto Faza: 0000-0001-8754-1341 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors thank the Centro de Supercomputación de Galicia (CESGA) for the allocation of computational resources. M.C.R. is thankful to the Spanish Ministerio de Ministerio de Educación Cultura y Deporte for her FPU fellowship. M.M.-L. thanks the Xunta de Galicia for a post-doctoral research contract (ED481B 2016/166-0). This work has been ́ Industria y funded by the Ministerio de Economi a, Competitividad (grant CTQ2016-75023-C2-2P) and the Xunta de Galicia (grant EM2014/040).



ADDITIONAL NOTES Despite being aware of the basic character of PMe3, we have omitted the study of those paths in which it is directly coordinated to the metal center due to the high barriers reported elsewhere.36 As a consequence, only those routes in a

7023

DOI: 10.1021/acsomega.8b01278 ACS Omega 2018, 3, 7019−7026

Article

ACS Omega

(18) Young, C. G.; Laughlin, L. J.; Colmanet, S.; Scrofani, S. D. B. Oxygen Atom Transfer, Sulfur Atom Transfer, and Correlated Electron-Nucleophile Transfer Reactions of Oxo- and Thiomolybdenum(IV) Complexes: Synthesis of Oxothiomolybdenum(VI) and (Hydroxo)oxomolybdenum(V) Species. Inorg. Chem. 1996, 35, 5368−5377. (19) Most, K.; Koepke, S.; Dall’Antonia, F.; Moesch-Zanetti, N. C. The First Molybdenum Dioxo Compounds With η2 -Pyrazolate Ligands: Crystal Structure and Oxo Transfer Properties. Chem. Commun. 2002, 1676−1677. (20) Eagle, A. A.; Gable, R. W.; Thomas, S.; Sproules, S. A.; Young, C. G. Sulfur Atom Transfer Reactions of Tungsten(Vi) and Tungsten(IV) Chalcogenide Complexes. Polyhedron 2004, 23, 385− 394. (21) Guretzki, S.; Papenbrock, J. Characterization of the Sulfurtransferase Family from Oryza Sativa ,L. Plant Physiol. Biochem. 2011, 49, 1064−1070. (22) Farwell, C. C.; McIntosh, J. A.; Hyster, T. K.; Wang, Z. J.; Arnold, F. H. Enantioselective Imidation of Sulfides Via EnzymeCatalyzed Intermolecular Nitrogen-Atom Transfer. J. Am. Chem. Soc. 2014, 136, 8766−8771. (23) Li, J.; Andrejić, M.; Mata, R. A.; Ryde, U. A Computational Comparison of Oxygen Atom Transfer Catalyzed by Dimethyl Sulfoxide Reductase with Mo and W. Eur. J. Inorg. Chem. 2015, 2015, 3580−3589. (24) Farwell, C. C.; Zhang, R. K.; McIntosh, J. A.; Hyster, T. K.; Arnold, F. H. Enantioselective Enzyme-Catalyzed Aziridination Enabled by Active-Site Evolution of a Cytochrome P450. ACS Cent. Sci. 2015, 1, 89−93. (25) Elrod, L. T.; Kim, E. Lewis Acid Assisted Nitrate Reduction with Biomimetic Molybdenum Oxotransferase Complex. Inorg. Chem. 2018, 57, 2594−2602. (26) Heinze, K. Bioinspired functional analogs of the active site of molybdenum enzymes: Intermediates and mechanisms. Coord. Chem. Rev. 2015, 300, 121−141. (27) Yachandra, V. K.; Sauer, K.; Klein, M. P. Manganese Cluster in Photosynthesis: Where Plants Oxidize Water to Dioxygen. Chem. Rev. 1996, 96, 2927−2950. (28) Najafpour, M. M.; Ghobadi, M. Z.; Haghighi, B.; Tomo, T.; Shen, J.-R.; Allakhverdiev, S. I. Comparison of Nano-Sized Mn Oxides With the Mn Cluster of Photosystem II as Catalysts for Water Oxidation. Biochim. Biophys. Acta, Bioenerg. 2015, 1847, 294−306. (29) Lorbek, G.; Lewinska, M.; Rozman, D. Cytochrome P450s in the Synthesis of Cholesterol and Bile Acids - from Mouse Models to Human Diseases. FEBS J. 2012, 279, 1516−1533. (30) Hille, R. Molybdenum Enzymes. Essays Biochem. 1999, 34, 125−137. (31) Hille, R. The Molybdenum Oxotransferases and Related Enzymes. Dalton Trans. 2013, 42, 3029−3042. (32) Colorado-Peralta, R.; Sánchez-Vázquez, M.; HernándezAhuactzi, I. F.; Sánchez-Ruiz, S. A.; Contreras, R.; Flores-Parra, A.; Castillo-Blum, S. E. Structural study of molybdenum(VI) complexes containing bidentate ligands: Synthesis, characterization and DFT calculations. Polyhedron 2012, 48, 72−79. (33) Sugimoto, H.; Sato, M.; Asano, K.; Suzuki, T.; Mieda, K.; Ogura, T.; Matsumoto, T.; Giles, L. J.; Pokhrel, A.; Kirk, M. L.; Itoh, S. A Model for the Active-Site Formation Process in DMSO Reductase Family Molybdenum Enzymes Involving Oxido-Alcoholato and Oxido-Thiolato Molybdenum(VI) Core Structures. Inorg. Chem. 2016, 55, 1542−1550. (34) de Noronha, R. G.; Costa, P. J.; Romão, C. C.; Calhorda, M. J.; Fernandes, A. C. MoO2Cl2 as a Novel Catalyst for C-P Bond Formation and for Hydrophosphonylation of Aldehydes. Organometallics 2009, 28, 6206−6212. (35) de Noronha, R. G.; Romão, C. C.; Fernandes, A. C. MoO2Cl2 as a Novel Catalyst for the Synthesis of α-Aminophosphonates. Catal. Commun. 2011, 12, 337−340.

which the nitro compound is in the inner coordination sphere of the metal, in the first step, have been explored. b We performed preliminary calculations considering the metalcatalyzed reduction of nitrosostyrene to indole and found that the catalyst had deleterious effects on the kinetics of this process (produced higher activation barriers). c The energy obtained from the application of NEB is electronic energy.



REFERENCES

(1) Woo, L. K. Intermetal Oxygen, Sulfur, Selenium, and Nitrogen Atom Transfer Reactions. Chem. Rev. 1993, 93, 1125−1136. (2) Johnson, M. J. A.; Lee, P. M.; Odom, A. L.; Davis, W. M.; Cummins, C. C. Atom-bridged Intermediates in P- and P-Atom Transfer Reactions. Angew. Chem., Int. Ed. Engl. 1997, 36, 87−91. (3) Nemykin, V. N.; Laskin, J.; Basu, P. Isolation, Characterization of an Intermediate in an Oxygen Atom-Transfer Reaction, and the Determination of the Bond Dissociation Energy. J. Am. Chem. Soc. 2004, 126, 8604−8605. (4) Ding, M.; Rouzieres, M.; Losovyj, Y.; Pink, M.; Clerac, R.; Smith, J. M. Partial Nitrogen Atom Transfer: A New Synthetic Tool to Design Single-Molecule Magnets. Inorg. Chem. 2015, 54, 9075−9080. (5) Paul, T.; Rodehutskors, P. M.; Schmidt, J.; Burzlaff, N. Oxygen Atom Transfer Catalysis with Homogenous and Polymer-Supported N,N- and N,N,O-Heteroscorpionate Dioxidomolybdenum(VI) Complexes. Eur. J. Inorg. Chem. 2016, 2016, 2595−2602. (6) Ishii, Y.; Sakaguchi, S.; Iwahama, T. Innovation of Hydrocarbon Oxidation with Molecular Oxygen and Related Reactions. Adv. Synth. Catal. 2001, 343, 393−427. ́ (7) Shulpin, G. B. Metal-Catalyzed Hydrocarbon Oxygenations in Solutions: the Dramatic Role of Additives: a Review. J. Mol. Catal. A: Chem. 2002, 189, 39−66. (8) Qi, B.; Lou, L.-L.; Yu, K.; Bian, W.; Liu, S. Selective Epoxidation of Alkenes With Hydrogen Peroxide over Efficient and Recyclable Manganese Oxides. Catal. Commun. 2011, 15, 52−55. (9) Pariyar, A.; Bose, S.; Biswas, A. N.; Das, P.; Bandyopadhyay, P. Catalytic Hydrocarbon Oxidation by Iron Complex of 5,10,15Tris(Difluorophenyl)Corrole Via Activation of Hydroperoxides. Catal. Commun. 2013, 32, 23−27. (10) Qi, Y.; Luan, Y.; Yu, J.; Peng, X.; Wang, G. Nanoscaled Copper Metal Organic Framework (MOF) Based on Carboxylate Ligands as an Efficient Heterogeneous Catalyst for Aerobic Epoxidation of Olefins and Oxidation of Benzylic and Allylic Alcohols. Chem. - Eur. J. 2015, 21, 1589−1597. (11) Na, J.-G.; Yi, B. E.; Kim, J. N.; Yi, K. B.; Park, S.-Y.; Park, J.-H.; Kim, J.-N.; Ko, C. H. Hydrocarbon Production from Decarboxylation of Fatty Acid Without Hydrogen. Catal. Today 2010, 156, 44−48. (12) Harms, R. G.; Herrmann, W. A.; Kühn, F. E. Organorhenium dioxides as oxygen transfer systems: Synthesis, reactivity, and applications. Coord. Chem. Rev. 2015, 296, 1−23. (13) Lupp, D.; Christensen, N. J.; Dethlefsen, J. R.; Fristrup, P. DFT Study of the Molybdenum-Catalyzed Deoxydehydration of Vicinal Diols. Chem. - Eur. J. 2015, 21, 3435−3442. (14) de Vicente Poutás, L. C.; Castiñeira Reis, M.; Sanz, R.; López, C. S.; Faza, O. N. A Radical Mechanism for the Vanadium-Catalyzed Deoxydehydration of Glycols. Inorg. Chem. 2016, 55, 11372−11382. (15) Lup, A. N. K.; Abnisa, F.; Daud, W. M. A. W.; Aroua, M. K. A Review on Reaction Mechanisms of Metal-Catalyzed Deoxygenation Process in Bio-Oil Model Compounds. Appl. Catal. A 2017, 541, 87− 106. (16) Sandbrink, L.; Beckerle, K.; Meiners, I.; Liffmann, R.; Rahimi, K.; Okuda, J.; Palkovits, R. Supported Molybdenum Catalysts for the Deoxydehydration of 1,4-Anhydroerythritol into 2,5-Dihydrofuran. ChemSusChem 2017, 10, 1375−1379. (17) Gossett, J.; Srivastava, R. Rhenium-Catalyzed Deoxydehydration of Renewable Biomass Using Sacrificial Alcohol as Reductant. Tetrahedron Lett. 2017, 58, 3760−3763. 7024

DOI: 10.1021/acsomega.8b01278 ACS Omega 2018, 3, 7019−7026

Article

ACS Omega

Molybdenum(VI) to Phosphorus(III). Chem. - Eur. J. 2005, 11, 3255−3267. (54) Ž mirić, A.; Zaric, S. D. Theoretical Study on Intermediate in Oxygen Transfer Reaction in Molybdoenzyme Model System. Inorg. Chem. Commun. 2002, 5, 446−448. (55) Pietsch, M. A.; Hall, M. B. Theoretical Studies on Models for the Oxo-Transfer Reaction of Dioxomolybdenum Enzymes. Inorg. Chem. 1996, 35, 1273−1278. (56) Antony, S.; Bayse, C. A. Theoretical Studies of Models of the Active Site of the Tungstoenzyme Acetylene Hydratase. Organometallics 2009, 28, 4938−4944. (57) Fey, N.; Orpen, A. G.; Harvey, J. N. Building Ligand Knowledge Bases for Organometallic Chemistry: Computational Description of Phosphorus(III)-Donor Ligands and the MetalPhosphorus Bond. Coord. Chem. Rev. 2009, 253, 704−722. (58) Basu, P.; Nemykin, V. N.; Sengar, R. S. Substituent Effect on Oxygen Atom Transfer Reactivity from Oxomolybdenum Centers: Synthesis, Structure, Electrochemistry, and Mechanism. Inorg. Chem. 2009, 48, 6303−6313. (59) Smith, P. D.; Millar, A. J.; Young, C. G.; Ghosh, A.; Basu, P. Detection, Isolation, and Characterization of Intermediates in Oxygen Atom Transfer Reactions in Molybdoenzyme Model Systems. J. Am. Chem. Soc. 2000, 122, 9298−9299. (60) Nemykin, V. N.; Laskin, J.; Basu, P. Isolation, Characterization of an Intermediate in an Oxygen Atom-Transfer Reaction, and the Determination of the Bond Dissociation Energy. J. Am. Chem. Soc. 2004, 126, 8604−8605. (61) Kozuch, S.; Shaik, S. How to Conceptualize Catalytic Cycles? The Energetic Span Model. Acc. Chem. Res. 2011, 44, 101−110. (62) Guggilapu, S. D.; Prajapti, S. K.; Nagarsenkar, A.; Lalita, G.; Naidu Vegi, G. M.; Babu, B. N. MoO2Cl2 Catalyzed Efficient Synthesis of Functionalized 3,4-Dihydropyrimidin-2(1H)-ones/thiones and Polyhydroquinolines: Recyclability, Fluorescence and Biological Studies. New J. Chem. 2016, 40, 838−843. (63) Jeyakumar, K.; Chand, D. K. Application of Molybdenum(VI) Dichloride Dioxide (MoO2Cl2) in Organic Transformations. J. Chem. Sci. 2009, 121, 111−123. (64) Hartley, F. R. The cis- and trans-Effects of Ligands. Chem. Soc. Rev. 1973, 2, 163−179. (65) Westheimer, F. H. Pseudo-Rotation in the Hydrolysis of Phosphate Esters. Acc. Chem. Res. 1968, 1, 70−78. (66) López Silva, C.; Faza Nieto, O.; de Lera, A. R.; York, D. M. Pseudorotation Barriers of Biological Oxyphosphoranes: A Challenge for Simulations of Ribozyme Catalysis. Chem. - Eur. J. 2005, 11, 2081−2093. (67) Bunyan, P. J.; Cadogan, J. I. G. 7. The Reactivity of Organophosphorus Compounds. Part XIV. Deoxygenation of Aromatic C-Nitroso-Compounds by Triethyl Phosphite and Triphenylphosphine: a New Cyclisation Reaction. J. Chem. Soc. 1963, 42− 49. (68) Davies, I. W.; Guner, V. A.; Houk, K. N. Theoretical Evidence for Oxygenated Intermediates in the Reductive Cyclization of Nitrobenzenes. Org. Lett. 2004, 6, 743−746. (69) Penoni, A.; Palmisano, G.; Zhao, Y.-L.; Houk, K. N.; Volkman, J.; Nicholas, K. M. On the Mechanism of Nitrosoarene Alkyne Cycloaddition. J. Am. Chem. Soc. 2009, 131, 653−661. (70) Cadogan, J. I. G.; Todd, M. J. Reduction of Nitro- and NitrosoCompounds by Tervalent Phosphorus Reagents. Part IV. Mechanistic Aspects of the Reduction of 2,4,6-Trimethyl-2′-Nitrobiphenyl, 2Nitrobiphenyl, and Nitrobenzene. J. Chem. Soc. C 1969, 2808−2813. (71) Sheppard, D.; Terrell, R.; Henkelman, G. Optimization Methods for Finding Minimum Energy Paths. J. Chem. Phys. 2008, 128, 134106−134110. (72) Smolinsky, G.; Wasserman, E.; Yager, W. A. The E.P.R. of Ground State Triplet Nitrenes. J. Am. Chem. Soc. 1962, 84, 3220− 3221. (73) McConaghy, J. S.; Lwowski, W. Singlet and Triplet Nitrenes. I. Carbethoxynitrene Generated by .Alpha. Elimination. J. Am. Chem. Soc. 1967, 89, 2357−2364.

(36) Castiñeira Reis, M.; Marín Luna, M.; Silva López, C.; Nieto Faza, O. [MoO2]2+-Mediated Oxygen Atom Transfer via an Unusual Lewis Acid Mechanism. Inorg. Chem. 2017, 56, 10570−10575. (37) Sanz, R.; Escribano, J.; Pedrosa, M. R.; Aguado, R.; Arnáiz, F. J. Dioxomolybdenum(Vi)-Catalyzed Reductive Cyclization of Nitroaromatics. Synthesis of Carbazoles and Indoles. Adv. Synth. Catal. 2007, 349, 713−718. (38) Juwarker, H.; Suk, J.-m.; Jeong, K.-S. In Anion Recognition in Supramolecular Chemistry; Gale, P. A., Dehaen, W., Eds.; Springer: Berlin, Heidelberg, 2010; pp 177−204. (39) Prudhomme, M. Biological Targets of Antitumor Indolocarbazoles Bearing a Sugar Moiety. Curr. Med. Chem.: Anti-Cancer Agents 2004, 4, 509−521. (40) Stolc, S. Indole Derivatives as Neuroprotectants. Life Sci. 1999, 65, 1943−1950. (41) Formenti, D.; Ferretti, F.; Ragaini, F. Synthesis of NHeterocycles by Reductive Cyclization of Nitro Compounds using Formate Esters as Carbon Monoxide Surrogates. ChemCatChem 2018, 10, 148−152. (42) Yamamoto, Y.; Ohkubo, E.; Shibuya, M. Synthesis of 3-Aryl-2(trifluoromethyl)indoles via Copper-Catalyzed Hydroarylation and Subsequent Cadogan Cyclization. Adv. Synth. Catal. 2017, 359, 1747−1751. (43) Tong, S.; Xu, Z.; Mamboury, M.; Wang, Q.; Zhu, J. Aqueous Titanium Trichloride Promoted Reductive Cyclization of o-Nitrostyrenes to Indoles: Development and Application to the Synthesis of Rizatriptan and Aspidospermidine. Angew. Chem., Int. Ed. 2015, 54, 11809−11812. (44) Cadogan, J. I. G.; Cameron-Wood, M.; Mackie, R. K.; Searle, R. J. G. 896. The Reactivity of Organophosphorus Compounds. Part Xix. Reduction of Nitro-Compounds by Triethyl Phosphite: a Convenient New Route to Carbazoles, Indoles, Indazoles, Triazoles, and Related Compounds. J. Chem. Soc. 1965, 4831−4837. (45) Cadogan, J. I. G. Organophosphorus Reagents in Organic Synthesis; Academic Press, 1979; p 608. (46) Nguyen, S. T.; Kwasny, S. M.; Ding, X.; Williams, J. D.; Peet, N. P.; Bowlin, T. L.; Opperman, T. J. Synthesis and Antifungal Evaluation of Head-to-Head and Head-to-Tail Bisamidine Compounds. Bioorg. Med. Chem. 2015, 23, 5789−5798. (47) Nguyen, S. T.; Williams, J. D.; Butler, M. M.; Ding, X.; Mills, D. M.; Tashjian, T. F.; Panchal, R. G.; Weir, S. K.; Moon, C.; Kim, H.O.; Marsden, J. A.; Peet, N. P.; Bowlin, T. L. Synthesis and Antibacterial Evaluation of New, Unsymmetrical Triaryl Bisamidine Compounds. Bioorg. Med. Chem. Lett. 2014, 24, 3366−3372. (48) Ansari, N. H.; Taylor, M. C.; Söderberg, B. C. Syntheses of Three Naturally Occurring Polybrominated 3,3-bi-1H-Indoles. Tetrahedron Lett. 2017, 58, 1053−1056. (49) Ansari, N. H.; Söderberg, B. C. Short syntheses of the indole alkaloids alocasin A, scalaridine A, and hyrtinadine A-B. Tetrahedron 2016, 72, 4214−4221. (50) Williams, J. D.; Nguyen, S. T.; Gu, S.; Ding, X.; Butler, M. M.; Tashjian, T. F.; Opperman, T. J.; Panchal, R. G.; Bavari, S.; Peet, N. P.; Moir, D. T.; Bowlin, T. L. Potent and Broad-Spectrum Antibacterial Activity of Indole-Based Bisamidine Antibiotics: Synthesis and SAR of Novel Analogs of MBX 1066 and MBX 1090. Bioorg. Med. Chem. 2013, 21, 7790−7806. (51) Freeman, A. W.; Urvoy, M.; Criswell, M. E. Triphenylphosphine-Mediated Reductive Cyclization of 2-Nitrobiphenyls: A Practical and Convenient Synthesis of Carbazoles. J. Org. Chem. 2005, 70, 5014−5019. (52) Alberto, G. C.; Leonardo, B.; Francesco, T.; Harvey, J. N.; Paola, B. Spin Forbidden Reactions: Adiabatic Transition States Using Spin-Orbit Coupled Density Functional Theory. Chem. - Eur. J. 2017, 24, 5006−5015. (53) Millar, A. J.; Doonan, C. J.; Smith, P. D.; Nemykin, V. N.; Basu, P.; Young, C. G. Oxygen Atom Transfer in Models for Molybdenum Enzymes: Isolation and Structural, Spectroscopic, and Computational Studies of Intermediates in Oxygen Atom Transfer from 7025

DOI: 10.1021/acsomega.8b01278 ACS Omega 2018, 3, 7019−7026

Article

ACS Omega (74) Becke, A. D. Density-Functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648−5652. (75) Perdew, J. P. Electronic Structure of Solids; Akademie-Verlag: Berlin, 1991; pp 11−20. (76) Rappoport, D.; Furche, F. Property-Optimized Gaussian Basis Sets for Molecular Response Calculations. J. Chem. Phys. 2010, 133, No. 134105. (77) Weigend, F.; Ahlrichs, R. Balanced Basis Sets of Split Valence, Triple Zeta Valence and Quadruple Zeta Valence Quality for H to Rn: Design and Assessment of Accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297−3305. (78) Andrae, D.; Häußermann, U.; Dolg, M.; Stoll, H.; Preuß, H. Energy-Adjusted Ab Initio Pseudopotentials for the Second and Third Row Transition Elements. Theor. Chim. Acta 1990, 77, 123−141. (79) Tomasi, J.; Mennucci, B.; Cammi, R. Quantum Mechanical Continuum Solvation Models. Chem. Rev. 2005, 105, 2999−3093. (80) Frisch, M. J.;et al. et al.Gaussian 09, revision C.1; Gaussian Inc.: Wallingford, CT, 2009. (81) Neese, F. The ORCA Program System. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2012, 2, 73−78. (82) Sinnecker, S.; Rajendran, A.; Klamt, A.; Diedenhofen, M.; Neese, F. Calculation of Solvent Shifts on Electronic G-Tensors With the Conductor-Like Screening Model (COSMO) and its SelfConsistent Generalization to Real Solvents (Direct COSMO-RS). J. Phys. Chem. A 2006, 110, 2235−2245.

7026

DOI: 10.1021/acsomega.8b01278 ACS Omega 2018, 3, 7019−7026