Stereoselectivity, Different Oxidation States, and Multiple Spin States in the Cyclopropanation of Olefins Catalyzed by Fe–Porphyrin Complexes

Iron Heme Enzyme-Catalyzed Cyclopropanations with Diazirines as Carbene Precursors: Computational Explorations of Diazirine Activation and Cyclopropanation Mechanism

The mechanism of cyclopropanations with diazirines as air-stable and user-friendly alternatives to commonly employed diazo compounds within iron heme enzyme-catalyzed carbene transfer reactions has been studied by means of density functional theory (DFT) calculations of model systems, quantum mechanics/molecular mechanics (QM/MM) calculations, and molecular dynamics (MD) simulations of the iron carbene and the cyclopropanation transition state in the enzyme active site. The reaction is initiated by a direct diazirine-diazo isomerization occurring in the active site of the enzyme. In contrast, an isomerization mechanism proceeding via the formation of a free carbene intermediate in lieu of a direct, one-step isomerization process was observed for model systems. Subsequent reaction with benzyl acrylate takes place through stepwise C-C bond formation via a diradical intermediate, delivering the cyclopropane product. The origin of the observed diastereo- and enantioselectivity in the enzyme was investigated through MD simulations, which indicate a preferred formation of the cis-cyclopropane by steric control.

Does an Enol Pathway Preclude High Stereoselectivity in Iron-Catalyzed Indole C-H Functionalization via Carbene Insertion?

C-H functionalization of indoles via Fe carbenoids presents an attractive strategy to obtain biologically important structural motifs. However, obtaining good stereoselectivity with Fe has been a significant challenge. It is unclear whether the low selectivity is due to a radical pathway or an ionic mechanism involving metal-free species. We therefore present a density functional theory (DFT) study of indole alkylation with diazoacetates catalyzed by Fe(ClO₄)TMEDA/spirobisoxazoline and myoglobin. We explore three mechanistic pathways: nucleophilic, radical, and oxocarbenium routes. The nucleophilic pathway is the most feasible with the formation of an enol species that tautomerizes to furnish the alkylated indole. While this mechanism is routinely proposed, the stereochemical model has been conspicuously absent until now. We show that the conventionally invoked enol pathway is not responsible for the low enantiomeric excess. The enol intermediate can stay coordinated to the catalyst via different binding sites placing the enol in proximity to the chiral environment and affecting the stereoselective proton transfer. Both the binding strength and the chiral environment are crucial for obtaining high selectivity. Our study provides the much needed insights for the modest-low selectivities of Fe systems and could help in expediting the discovery of an efficient catalytic system. These mechanistic underpinnings could also be applicable to other metal (Rh, Pd, Cu, etc.)-catalyzed X-H insertion reactions.

On the Sensitivity to Density-Functional Approximations for CO Binding Energies of Single-Atom Catalysts in Nitrogen-Doped Graphene

Density functional theory (DFT) methods are the working horse in screening new catalytic materials. They are widely used to predict trends in binding energies, which are then used to compare the activity of different materials. The binding strength of CO is an important descriptor to the CO₂ reduction catalytic activity of the single transition metal atoms embedded on nitrogen-doped graphene (TM/NG). In this work, however, we show that CO binding strengths in different TM/NG has very different sensitivity to DFT methods. Specifically, Fe/NG CO binding energy changes dramatically with the percentage of exact exchange in the functional; Co/NG does less so, while Ni/NG nearly has no change. Such varying behaviors is a direct result of different local spin configurations, similar to the performance of DFT methods for metal porphyrin complexes. Therefore, caution should be exercised when using DFT binding energies for quantitative predictions in TM/NG single atom catalysis.

Porphyrin N-Pincer Pd(II)-Complexes in Water: A Base-Free and Nature-Inspired Protocol for the Oxidative Self-Coupling of Potassium Aryltrifluoroborates in Open-Air

Metalloporphyrins (and porphyrins) are well known as pigments of life in nature, since representatives of this group include chlorophylls (Mg-porphyrins) and heme (Fe-porphyrins). Hence, the construction of chemistry based on these substances can be based on the imitation of biological systems. Inspired by nature, in this article we present the preparation of five different porphyrin, meso-tetraphenylporphyrin (TPP), meso-tetra(p-anisyl)porphyrin (TpAP), tetrasodium meso-tetra(p-sulfonatophenyl)porphyrin (TSTpSPP), meso-tetra(m-hydroxyphenyl)porphyrin (TmHPP), and meso-tetra(m-carboxyphenyl)porphyrin (TmCPP) as well as their N-pincer Pd(II)-complexes such as Pd(II)-meso-tetraphenylporphyrin (PdTPP), Pd(II)-meso-tetra(p-anisyl)porphyrin (PdTpAP), Pd(II)-tetrasodium meso-tetra(p-sulfonatophenyl)porphyrin (PdTSTpSPP), Pd(II)-meso-tetra(m-hydroxyphenyl)porphyrin (PdTmHPP), and Pd(II)-meso-tetra(m-carboxyphenyl)porphyrin (PdTmCPP). These porphyrin N-pincer Pd(II)-complexes were studied and found to be effective in the base-free self-coupling reactions of potassium aryltrifluoroborates (PATFBs) in water at ambient conditions. The catalysts and the products (symmetrical biaryls) were characterized using their spectral data. The high yields of the biaryls, the bio-mimicking conditions, good substrate feasibility, evading the use of base, easy preparation and handling of catalysts, and the application of aqueous media, all make this protocol very attractive from a sustainability and cost-effective standpoint.

Combined Experimental and Computational Study of Ruthenium N-Hydroxyphthalimidoyl Carbenes in Alkene Cyclopropanation Reactions

A combined experimental-computational approach has been used to study the cyclopropanation reaction of N-hydroxyphthalimide diazoacetate (NHPI-DA) with various olefins, catalyzed by a ruthenium-phenyloxazoline (Ru-Pheox) complex. Kinetic studies show that the better selectivity of the employed redox-active NHPI diazoacetate is a result of a much slower dimerization reaction compared to aliphatic diazoacetates. Density functional theory calculations reveal that several reactions can take place with similar energy barriers, namely, dimerization of the NHPI diazoacetate, cyclopropanation (inner-sphere and outer-sphere), and a previously unrecognized migratory insertion of the carbene into the phenyloxazoline ligand. The calculations show that the migratory insertion reaction yields an unconsidered ruthenium complex that is catalytically competent for both the dimerization and cyclopropanation, and its relevance is assessed experimentally. The stereoselectivity of the reaction is argued to stem from an intricate balance between the various mechanistic scenarios.

A computational approach to understand the role of metals and axial ligands in artificial heme enzyme catalyzed C-H insertion

Engineered heme enzymes such as myoglobin and cytochrome P450s metalloproteins are gaining widespread importance due to their efficiency in catalyzing non-natural reactions. In a recent strategy, the naturally occurring Fe metal in the heme unit was replaced with non-native metals such as Ir, Rh, Co, Cu, etc., and axial ligands to generate artificial metalloenzymes. Determining the best metal-ligand for a chemical transformation is not a trivial task. Here we demonstrate how computational approaches can be used in deciding the best metal-ligand combination which would be highly beneficial in designing new enzymes as well as small molecule catalysts. We have used Density Functional Theory (DFT) to shed light on the enhanced reactivity of an Ir system with varying axial ligands. We look at the insertion of a carbene group generated from diazo precursors via N2 extrusion into a C-H bond. For both Ir(Me) and Fe systems, the first step, i.e., N2 extrusion is the rate determining step. Strikingly, neither the better ligand overlap with 5d orbitals on Ir nor the electrophilicity on the carbene centre play a significant role. A comparison of Fe and Ir systems reveals that a lower distortion in the Ir(Me)-porphyrin on moving from the reactant to the transition state renders it catalytically more active. We notice that for both metal porphyrins, the free energy barriers are affected by axial ligand substitution. Further, for Fe porphyrin, the axial ligand also changes the preferred spin state. We show that for the carbene insertion into the C-H bond, Fe porphyrin systems undergo a stepwise HAT (hydrogen atom transfer) instead of a concerted hydride transfer process. Importantly, we find that the substitution of the axial Me ligand on Ir to imidazole or chloride, or without an axial substitution changes the rate determining step of the reaction. Therefore, an optimum ligand that can balance the barriers for both steps of the catalytic cycle is essential. We subsequently used the QM cluster approach to delineate the protein environment's role and mutations in improving the catalytic activity of the Ir(Me) system.

Accuracy of DLPNO-CCSD(T): Effect of Basis Set and System Size

The DLPNO-CCSD(T) method is designed to study large molecular systems at significantly reduced cost relative to its canonical counterpart. However, the error in this approach is also size-extensive and relies on cancellation of errors for the calculation of relative energies. This work provides a direct comparison of canonical CCSD(T) and TightPNO DLPNO-CCSD(T) calculations of reaction energies and barriers of a broad range of chemical reactions. The dataset includes acidities, anion binding affinities, enolization, Diels-Alder, nucleophilic substitution, and atom transfer reactions and complements existing theoretical datasets in terms of system size as well as new reaction types (e.g., anion binding affinities and chlorine atom transfer reactions). The performance of DLPNO-CCSD(T) was further examined with respect to systematic variation of basis set and system size and amounts of nonbonded interaction present in the system. The errors in the DLPNO-CCSD(T) were found to be relatively insensitive to the choice of basis set for small systems but increase monotonically with system size. Additionally, calculations of barriers appear to be more challenging than reaction energies with errors exceeding 5 kJ mol^-1 for many Diels-Alder reactions. Further tests on three realistic organic reactions reveal the impact of the DLPNO approximation in calculating absolute and relative barriers that are important for predictions such as stereoselectivity.

Theoretical insight into the opposite redox activity of iron complexes toward the ring opening polymerization of lactide and epoxide

The origin of opposite reactivity in the ring-opening polymerization of lactide (LA) and cyclohexene oxide (CHO) catalyzed by redox-switchable bis(imino)pyridine iron complexes has been computationally elucidated.

Computational Studies on the Mechanism and Origin of the Different Regioselectivities of Manganese Porphyrin-Catalyzed C-H Bond Hydroxylation and Amidation of Equilenin Acetate

The manganese porphyrin-catalyzed C-H bond hydroxylation and amidation of equilenin acetate developed by Breslow and his co-worker have been investigated with density functional theory (DFT) calculations. The hydroxylation of C(sp²)-H bond of equilenin acetate leading to the 6-hydroxylated product is more favorable than the hydroxylation of C(sp³)-H bond of equilenin acetate, leading to the 11β-hydroxylation product. The computational results suggest that the C(sp²)-H bond hydroxylation of equilenin acetate undergoes an oxygen-atom-transfer mechanism, which is more favorable than the C(sp³)-H bond hydroxylation undergoing the hydrogen-atom-abstraction/oxygen-rebound (HAA/OR) mechanism by 1.6 kcal/mol. That is why, the 6-hydroxylated product is the major product and the 11β-hydroxylated product is the minor product. In contrast, the 11β-amidated product is the only observed product in manganese porphyrin-catalyzed amidation reaction. The benzylic amidation undergoes a hydrogen-atom-abstraction/nitrogen-rebound (HAA/NR) mechanism, in which hydrogen atom abstraction is followed by nitrogen rebound, leading to the 11β-amidated product. The benzylic C(sp³)-H bond amidation at the C-11 position is more favorable than aromatic amidation at the C-6 position by 4.9 kcal/mol. Therefore, the DFT computational results are consistent with the experiments that manganese porphyrin-catalyzed C-H bond hydroxylation and amidation of equilenin acetate have different regioselectivities.

Bridging and axial carbene binding modes in cobalt corrole complexes: effect on carbene transfer

Catalytically relevant intermediates in carbene transfer reactions from a diazo precursor were investigated using cobalt corrole complexes. Two divergent mechanisms are proposed depending on the oxidation state of the cobalt center. Mechanistically driven factors for the usage of cobalt corroles in carbene transfer reactions are discussed.

An In-Depth Computational Study of Alkene Cyclopropanation Catalyzed by Fe(porphyrin)(OCH3) Complexes. The Environmental Effects on the Energy Barriers

Iron porphyrin methoxy complexes, of the general formula [Fe(porphyrin)(OCH₃)], are able to catalyze the reaction of diazo compounds with alkenes to give cyclopropane products with very high efficiency and selectivity. The overall mechanism of these reactions was thoroughly investigated with the aid of a computational approach based on density functional theory calculations. The energy profile for the processes catalyzed by the oxidized [Fe^III(Por)(OCH₃)] (Por = porphine) as well as the reduced [Fe^II(Por)(OCH₃)]^- forms of the iron porphyrin was determined. The main reaction step is the same in both of the cases, that is, the one leading to the terminal-carbene intermediate [Fe(Por)(OCH₃)(CHCO₂Et)] with simultaneous dinitrogen loss; however, the reduced species performs much better than the oxidized one. Contrarily to the iron(III) profile in which the carbene intermediate is directly obtained from the starting reactant complex, the favored iron(II) process is more intricate. The initially formed reactant adduct between [Fe^II(Por)(OCH₃)]^- and ethyl diazoacetate (EDA) is converted into a closer reactant adduct, which is in turn converted into the terminal iron porphyrin carbene [Fe(Por)(OCH₃)(CHCO₂Et)]^-. The two corresponding transition states are almost isoenergetic, thus raising the question of whether the rate-determining step corresponds to dinitrogen loss or to the previous structural and electronic rearrangement. The ethylene addition to the terminal carbene is a downhill process, which, on the open-shell singlet surface, presents a defined but probably short-living diradicaloid intermediate, though other spin-state surfaces do not show this intermediate allowing a direct access to the cyclopropane product. For the crucial stationary points, the more complex catalyst [Fe(2)(OCH₃)], in which a sterically hindered chiral bulk is mounted onto the porphyrin, was investigated. The corresponding computational data disclose the very significant effect of the porphyrin skeleton on the reaction energy profile. Though the geometrical features around the reactive core of the system remain unchanged, the energy barriers become much lower, thus revealing the profound effects that can be exerted by the three-dimensional organic scaffold surrounding the reaction site.

Spin-inversion mechanisms in O2 binding to a model heme complex revisited by density function theory calculations

Spin-inversion mechanisms in O₂ binding to a model heme complex, consisting of Fe(II)-porphyrin and imidazole, were investigated using density-functional theory calculations. First, we applied the recently proposed mixed-spin Hamiltonian method to locate spin-inversion structures between different total spin multiplicities. Nine spin-inversion structures were successfully optimized for the singlet-triplet, singlet-quintet, triplet-quintet, and quintet-septet spin-inversion processes. We found that the singlet-triplet spin-inversion points are located around the potential energy surface region at short Fe-O distances, whereas the singlet-quintet and quintet-septet spin-inversion points are located at longer Fe-O distances. This suggests that both narrow and broad crossing models play roles in O₂ binding to the Fe-porphyrin complex. To further understand spin-inversion mechanisms, we performed on-the-fly Born-Oppenheimer molecular dynamics calculations. The reaction coordinates, which are correlated to the spin-inversion dynamics between different spin multiplicities, are also discussed.

Mechanism and stereoselectivity of benzylic C-H hydroxylation by Ru-porphyrin: a computational study

The mechanism and origin of the stereoselectivity of asymmetric benzylic C-H hydroxylation by Ru-porphyrin were elucidated with density functional theory calculations. The reaction proceeds via a hydrogen-atom abstraction/oxygen-rebound pathway, wherein a high-valent ruthenium-oxo species abstracts a hydrogen atom from ethylbenzene to generate a radical pair intermediate, followed by the oxygen-rebound process to form 1-phenylethanol. The hydrogen-atom abstraction step is the rate- and stereoselectivity-determining step. Based on the mechanistic model, the computed stereoselectivity is in agreement with the experimental observations. Analysis of the distortion/interaction model suggests that stereoselectivity is determined by both the distortion energy of the ethylbenzene and the interaction energy between the ethylbenzene and the chiral Ru-porphyrin. The steric repulsion between the phenyl group of ethylbenzene and the bulky substituent of Ru-porphyrin is the leading cause of chiral induction.

Iron catalysts with N-ligands for carbene transfer of diazo reagents

This review provides an overview of the catalytic activity of iron complexes of nitrogen ligands in driving carbene transfers towards CC, C–H and X–H bonds. The reactivity of diazo reagents is discussed as well as the proposed reaction mechanisms.

Iron- and cobalt-catalyzed C(sp3)–H bond functionalization reactions and their application in organic synthesis

This review highlights the developments in iron and cobalt catalyzed C(sp³)–H bond functionalization reactions with emphasis on their applications in organic synthesis, i.e. natural products and pharmaceuticals synthesis and/or modification.