Kinetic Isotope Effect Probes the Reactive Spin State, As Well As the Geometric Feature and Constitution of the Transition State during H-Abstraction by Heme Compound II Complexes

Spin Crossover and Its Application in Organometallic Catalysis: Concepts and Recent Progress

Spin crossover is one of the most important properties of open-shell metal complexes. In organometallic catalytic reactions, catalysts can alter reaction kinetics by spin crossover, i. e., accelerating or hindering the reaction progression, as well as altering reaction pathways, modulating the reaction selectivity or promoting new reactions. This personal account outlines the introduction and development of important concepts such as "two-state reactivity" involving spin crossover, and proposes a new concept, "spin-responsive catalysis" to summarize the catalytic processes in which spin effects are present. Finally, the electronic mechanism of spin crossover accelerating the reaction and the role of spin crossover in changing the reaction path and regulating the reaction selectivity are introduced by taking two recent typical iron-catalyzed reactions recently reported by our group as examples.

C-H Bond Activation by Iridium(III) and Iridium(IV) Oxo Complexes

Oxidation of an iridium(III) oxo precursor enabled the structural, spectroscopic, and quantum-chemical characterization of the first well-defined iridium(IV) oxo complex. Side-by-side examination of the proton-coupled electron transfer thermochemistry revealed similar driving forces for the isostructural oxo complexes in two redox states due to compensating contributions from H+ and e- transfer. However, C-H activation of dihydroanthracene revealed significant hydrogen tunneling for the distinctly more basic iridium(III) oxo complex. Our findings complement the growing body of data that relate tunneling to ground state properties as predictors for the selectivity of C-H bond activation.

Two-State Hydrogen Atom Transfer Reactivity of Unsymmetric [Cu2(O)(NO)]2+ Complexes

We report the temperature-dependent spin switching of dicopper oxo nitrosyl [Cu2(O)(NO)]2+ complexes and their influence on hydrogen atom transfer (HAT) reactivity. Electron paramagnetic resonance (EPR) and Evans method analysis suggest that [Cu2(O)(NO)]2+ complexes transition from the S = 1/2 to the S = 3/2 state around ca. 202 K. At low temperatures (198 K) where S = 3/2 dominates, a strong correlation between the rate of HAT (kHAT) and the population of the S = 1/2 state was identified (R2 = 0.988), suggesting that the HAT by [Cu2(O)(NO)]2+ complexes proceeds by the S = 1/2 isomer. Installation of functional groups that introduce an unsymmetric secondary coordination environment accelerates the HAT rates through perturbation of the spin equilibria. Given the often unsymmetric coordination sphere of bimetallic active sites in natural proteins, we anticipate that similar strategies could be employed by metalloenzymes to control HAT reactions.

Insights into the electronic structure and mechanism of norcarane hydroxylation by OxoMn(V) porphyrin complexes: A density functional theory study

Norcarane hydroxylation by neutral [PorMn(V)O-L] (L═OH^- , F^- ) and cationic [PorMn(V)O-L]⁺ (L═H₂ O, imidazole) oxoMn(V) porphyrin complex models has been investigated by density functional theory calculations to better understand the reaction mechanism and electronic structure. We found that the energy barriers of norcarane hydroxylation by cationic oxoMn(V) porphyrin complexes are lower than those by neutral oxoMn(V) porphyrin complexes. This indicates that cationic oxoMn(V) porphyrin complexes enhance norcarane hydroxylation compared with neutral oxoMn(V) porphyrin complexes. According to electronic structure analysis, in the C─H activation step, electron transfer occurs through initial interaction between the σ_CH and rich-oxygen π(Mn═O) orbitals to form real donor orbitals, followed by transfer to the acceptor π*(Mn═O) orbitals. Moreover, single electron shifts from norcarane (CH) to Mn atom during C─H activation. The positive charge of the cationic complex stabilizes the acceptor orbital more than the donor orbital, reducing the energy gap between these orbitals, thus lowering the reaction barrier.

Ligand Architecture Perturbation Influences the Reactivity of Nonheme Iron(V)-Oxo Tetraamido Macrocyclic Ligand Complexes: A Combined Experimental and Theoretical Study

Iron(V)-oxo complexes bearing negatively charged tetraamido macrocyclic ligands (TAMLs) have provided excellent opportunities to investigate the chemical properties and the mechanisms of oxidation reactions of mononuclear nonheme iron(V)-oxo intermediates. Herein, we report the differences in chemical properties and reactivities of two iron(V)-oxo TAML complexes differing by modification on the "Head" part of the TAML framework; one has a phenyl group at the "Head" part (1), whereas the other has four methyl groups replacing the phenyl ring (2). The reactivities of 1 and 2 in both C-H bond activation reactions, such as hydrogen atom transfer (HAT) of 1,4-cyclohexadiene, and oxygen atom transfer (OAT) reactions, such as the oxidation of thioanisole and its derivatives, were compared experimentally. Under identical reaction conditions, 1 showed much greater reactivity than 2, such as a 10²-fold decrease in HAT and a 10⁵-fold decrease in OAT by replacing the phenyl group (i.e., 1) with four methyl groups (i.e., 2). Then, density functional theory calculations were performed to rationalize the reactivity differences between 1 and 2. Computations reproduced the experimental findings well and revealed that the replacement of the phenyl group in 1 with four methyl groups in 2 not only increased the steric hindrance but also enlarged the energy gap between the electron-donating orbital and the electron-accepting orbital. These two factors, steric hindrance and the orbital energy gap, resulted in differences in the reduction potentials of 1 and 2 and their reactivities in oxidation reactions.

EPR spectroscopy elucidates the electronic structure of [FeV(O)(TAML)] complexes

The complete hyperfine tensor of ¹⁷O of the Fe^V-oxo moeity was probed by ENDOR spectroscopy. The EPR spectroscopic results reported here provide a conclusive experimental basis for elucidating the electronic structure of the Fe^V-oxo complex.

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.

DFT Mechanistic Account for the Site Selectivity of Electron-Rich C(sp3)-H Bond in the Manganese-Catalyzed Aminations

DFT study suggests that the C-H cleavage involved in the C(sp³)-H amination catalyzed by manganese perchlorophthalocyanine complex [Mn^III(ClPc)]⁺ proceeds via hydride transfer (HYT), instead of hydrogen atom transfer (HAT), thus elucidating why the reaction favors aminating the electron-rich C-H bond, rather than that with smaller bond dissociation energy preferred by HAT. Detailed analyses indicate the HYT actually occurs via concerted electron and hydrogen atom transfers and the redox-active ClPc ligand enables the HYT.

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.

Computational Exploration of Chiral Iron Porphyrin-Catalyzed Asymmetric Hydroxylation of Ethylbenzene Where Stereoselectivity Arises from π-π Stacking Interaction

The mechanism and origins of stereoselectivity of chiral iron porphyrin-catalyzed asymmetric hydroxylation of ethylbenzene were explored with density functional theory. The hydrogen atom abstraction is the rate- and stereoselectivity-determining step. In good agreement with experimental results, the formation of the (R)-1-phenylethanol product is found to be the most favorable pathway. The transition state of hydrogen atom abstraction which leads to the (S)-1-phenylethanol product is unfavorable by 1.7 kcal/mol compared to the corresponding transition state which leads to the (R)-1-phenylethanol product. Enantioselectivity arises from an attractive π-π stacking interaction between the phenyl group of ethylbenzene substrate and the naphthyl group of the porphyrin ligand.

Conformational turn triggers regio-selectivity in the bioactivation of thiophene-contained compounds mediated by cytochrome P450

In the present work, we performed Density Functional Theory calculations to explore the bioactivation mechanism of thiophene-containing molecules mediated by P450s. For this purpose, relatively large size compounds, 2,5-diaminothiophene derivatives were selected particularly for this investigation. Here we found the successive regio-selectivity triggered by conformational turn played a significant role in the occurrence of bioactivation. 2,5-Diaminothiophene was oxidized to a 2,5-diimine thiophene-reactive intermediate by Compound I (Cpd I) through successive activations of two N-H bonds (H3-N11 and H1-N6). This reaction exhibited three special characteristics: (1) self-controlled regio-selectivity during the oxidation process. There was a large scale of conformational turn in the abstraction of the first H atom which triggers the selection of the second H for abstraction. (2) Proton-shuttle mechanism. In high spin (HS) state, proton-shuttle mechanism was observed for the abstraction of the second H atom. (3) Spin-selective manner. In protein environment, the energy barrier in HS state was much lower than that in low spin state. The novel proposed bioactivation mechanism of 2,5-diaminothiophene compounds can help us in rational design of thiophene-contained drugs avoiding the occurrence of bioactivation.

Quantum Mechanics/Molecular Mechanics Studies on the Relative Reactivities of Compound I and II in Cytochrome P450 Enzymes

The cytochromes P450 are drug metabolizing enzymes in the body that typically react with substrates through a monoxygenation reaction. During the catalytic cycle two reduction and protonation steps generate a high-valent iron (IV)-oxo heme cation radical species called Compound I. However, with sufficient reduction equivalents present, the catalytic cycle should be able to continue to the reduced species of Compound I, called Compound II, rather than a reaction of Compound I with substrate. In particular, since electron transfer is usually on faster timescales than atom transfer, we considered this process feasible and decided to investigate the reaction computationally. In this work we present a computational study using density functional theory methods on active site model complexes alongside quantum mechanics/molecular mechanics calculations on full enzyme structures of cytochrome P450 enzymes. Specifically, we focus on the relative reactivity of Compound I and II with a model substrate for O⁻H bond activation. We show that generally the barrier heights for hydrogen atom abstraction are higher in energy for Compound II than Compound I for O⁻H bond activation. Nevertheless, for the activation of such bonds, Compound II should still be an active oxidant under enzymatic conditions. As such, our computational modelling predicts that under high-reduction environments the cytochromes P450 can react with substrates via Compound II but the rates will be much slower.

Stereoselective Metabolism of Omeprazole by Cytochrome P450 2C19 and 3A4: Mechanistic Insights from DFT Study

The efficacy of S-omeprazole as a proton pump inhibitor compared with that of its enantiomer R-omeprazole is studied using density functional theoretical calculations. The pharmacokinetic studies suggest that the efficacy of S-omeprazole presumably depends on metabolic pathway and excretion from the human body. The density functional theory calculations at SMD_water-B3LYP-D3/6-311+G(d,p)/LANL2DZ//B3LYP/6-31G(d)/LANL2DZ with triradicaloid model active species, [Por^•+Fe^IV(SH)O], of CYP2C19 enzyme with high-spin quartet and low-spin doublet states demonstrate C-H bond activation mechanism through a two-state rebound process for the hydroxylation of R-omeprazole and S-omeprazole. The calculated activation free energy barriers for the hydrogen abstraction are 15.7 and 17.5 kcal/mol for R-omeprazole and S-omeprazole, respectively. The hydroxylation of R-omeprazole and S-omeprazole is thermodynamically favored; however, the hydroxylated intermediate of S-omeprazole further disintegrates to metabolite 5- O-desmethylomeprazole with a higher kinetic barrier. We have examined the sulfoxidation of S-omeprazole to omeprazole sulfone metabolite by CYP3A4, and the observed activation free energy barrier is 9.9 kcal/mol. The computational results reveal that CYP2C19 exclusively metabolizes R-omeprazole to hydroxyomeprazole, which is hydrophilic and can easily excrete, whereas CYP3A4 metabolizes S-omeprazole to lipophilic sulfone; hence, the excretion of this metabolite would be relatively slower from the body. The spin density analysis and molecular orbital analysis performed using biorthogonalization calculations indicate that R-omeprazole favors high-spin pathway for metabolism process whereas S-omeprazole prefers the low-spin pathway.

Kinetic Isotope Effect Determination Probes the Spin of the Transition State, Its Stereochemistry, and Its Ligand Sphere in Hydrogen Abstraction Reactions of Oxoiron(IV) Complexes

This Account outlines interplay of theory and experiment in the quest to identify the reactive-spin-state in chemical reactions that possess a few spin-dependent routes. Metalloenzymes and synthetic models have forged in recent decades an area of increasing appeal, in which oxometal species bring about functionalization of hydrocarbons under mild conditions and via intriguing mechanisms that provide a glimpse of Nature's designs to harness these reactions. Prominent among these are oxoiron(IV) complexes, which are potent H-abstractors. One of the key properties of oxoirons is the presence of close-lying spin-states, which can mediate H-abstractions. As such, these complexes form a fascinating chapter of spin-state chemistry, in which chemical reactivity involves spin-state interchange, so-called two-state reactivity (TSR) and multistate reactivity (MSR). TSR and MSR pose mechanistic challenges. How can one determine the structure of the reactive transition state (TS) and its spin state for these mechanisms? Calculations can do it for us, but the challenge is to find experimental probes. There are, however, no clear kinetic signatures for the reactive-spin-state in such reactions. This is the paucity that our group has been trying to fill for sometime. Hence, it is timely to demonstrate how theory joins experiment in realizing this quest. This Account uses a set of the H-abstraction reactions of 24 synthetic oxoiron(IV) complexes and 11 hydrocarbons, together undergoing H-abstraction reactions with TSR/MSR options, which provide experimentally determined kinetic isotope effect (KIE_exp) data. For this set, we demonstrate that comparing KIE_exp results with calculated tunneling-augmented KIE (KIE_TC) data leads to a clear identification of the reactive spin-state during H-abstraction reactions. In addition, generating KIE_exp data for a reaction of interest, and comparing these to KIE_TC values, provides the mechanistic chemist with a powerful capability to identify the reactive-TS in terms of not only its spin state but also its geometry and ligand-sphere constitution. Since tunneling "cuts through" barriers, it serves as a chemical selectivity factor. Thus, we show that in a family of oxoirons reacting with one hydrocarbon, the tunneling efficiency increases as the ligands become better electron donors. This generates counterintuitive-reactivity patterns, like antielectrophilic reactivity, and induces spin-state reactivity reversals because of differing steric demands of the corresponding ^2S+1TS species, etc. Finally, for the same series, the Account reaches intuitive understanding of tunneling trends. It is shown that the increase of ligand's donicity results in electrostatic narrowing of the barrier, while the decrease of donicity and increase of bond-order asymmetry in the TS (inter alia due to Bell-Evans-Polanyi effects) broadens the barrier. Predictions are made that usage of powerful electron-donating ligands may train H-abstractors to activate the strongest C-H bond in a molecule. The concepts developed here are likely to be applicable to other oxometals in the d- and f-blocks.

Privileged Role of Thiolate as the Axial Ligand in Hydrogen Atom Transfer Reactions by Oxoiron(IV) Complexes in Shaping the Potential Energy Surface and Inducing Significant H-Atom Tunneling

An H/D kinetic isotope effect (KIE) of 80 is found at -20 °C for the oxidation of 9,10-dihydroanthracene by [Fe^IV(O)(TMCS)]⁺, a complex supported by the tetramethylcyclam (TMC) macrocycle with a tethered thiolate. This KIE value approaches that previously predicted by DFT calculations. Other [Fe^IV(O)(TMC)(anion)] complexes exhibit values of 20, suggesting that the thiolate ligand of [Fe^IV(O)(TMCS)]⁺ plays a unique role in facilitating tunneling. Calculations show that tunneling is most enhanced (a) when the bond asymmetry between C-H bond breaking and O-H bond formation in the transition state is minimized, and (b) when the electrostatic interactions in the O---H---C moiety are maximal. These two factors-which peak for the best electron donor, the thiolate ligand-afford a slim and narrow barrier through which the H-atom can tunnel most effectively.