Keto-Enol Tautomerization Triggers an Electrophilic Aldehyde Deformylation Reaction by a Nonheme Manganese(III)-Peroxo Complex

What is the Origin of the Regioselective C3-Hydroxylation of L-Arg by the Nonheme Iron Enzyme Capreomycin C?

The nonheme iron dioxygenase capreomycin C (CmnC) hydroxylates a free L-arginine amino acid regio- and stereospecifically at the C3-position as part of the capreomycin antibiotics biosynthesis. Little is known on its structure, catalytic cycle and substrate specificity and, therefore, a comprehensive computational study was performed. A large QM cluster model of CmnC was created of 297 atoms and the mechanisms for C3-H, C4-H and C5-H hydroxylation and C3-C4 desaturation were investigated. All low-energy pathways correspond to radical reaction mechanisms with an initial hydrogen atom abstraction followed by OH rebound to form alcohol product complexes. The work is compared to alternative L-Arg hydroxylating nonheme iron dioxygenases and the differences in active site polarity are compared. We show that a tight hydrogen bonding network in the substrate binding pocket positions the substrate in an ideal orientation for C3-H activation, whereby the polar groups in the substrate binding pocket induce an electric field effect that guides the selectivity.

Tutorial Review on the Set-Up and Running of Quantum Mechanical Cluster Models for Enzymatic Reaction Mechanisms

Enzymes turnover substrates into products with amazing efficiency and selectivity and as such have great potential for use in biotechnology and pharmaceutical applications. However, details of their catalytic cycles and the origins surrounding the regio- and chemoselectivity of enzymatic reaction processes remain unknown, which makes the engineering of enzymes and their use in biotechnology challenging. Computational modelling can assist experimental work in the field and establish the factors that influence the reaction rates and the product distributions. A popular approach in modelling is the use of quantum mechanical cluster models of enzymes that take the first- and second coordination sphere of the enzyme active site into consideration. These QM cluster models are widely applied but often the results obtained are dependent on model choice and model selection. Herein, we show that QM cluster models can give highly accurate results that reproduce experimental product distributions and free energies of activation within several kcal mol-1, regarded that large cluster models with >300 atoms are used that include key hydrogen bonding interactions and charged residues. In this tutorial review, we give general guidelines on the set-up and applications of the QM cluster method and discuss its accuracy and reproducibility. Finally, several representative QM cluster model examples on metal-containing enzymes are presented, which highlight the strength of the approach.

Dioxygenase Chemistry in Nucleophilic Aldehyde Deformylations Utilizing Dicopper O2-Derived Peroxide Complexes

The chemistry of copper-dioxygen complexes is relevant to copper enzymes in biology as well as in (ligand)Cu-O2 (or Cu2-O2) species utilized in oxidative transformations. For overall energy considerations, as applicable in chemical synthesis, it is beneficial to have an appropriate atom economy; both O-atoms of O2(g) are transferred to the product(s). However, examples of such dioxygenase-type chemistry are extremely rare or not well documented. Herein, we report on nucleophilic oxidative aldehyde deformylation reactivity by the peroxo-dicopper(II) species [Cu2II(BPMPO-)(O22-)]1+ {BPMPO-H = 2,6-bis{[(bis(2-pyridylmethyl)amino]methyl}-4-methylphenol)} and [Cu2II(XYLO-)(O22-)]1+ (XYLO- = a BPMPO- analogue possessing bis(2-{2-pyridyl}ethyl)amine chelating arms). Their dicopper(I) precursors are dioxygenase catalysts. The O2(g)-derived peroxo-dicopper(II) intermediates react rapidly with aldehydes like 2-phenylpropionaldehyde (2-PPA) and cyclohexanecarboxaldehyde (CCA) in 2-methyltetrahydrofuran at -90 °C. Warming to room temperature (RT) followed by workup results in good yields of formate (HC(O)O-) along with ketones (acetophenone or cyclohexanone). Mechanistic investigation shows that [Cu2II(BPMPO-)(O22-)]1+ species initially reacts reversibly with the aldehydes to form detectable dicopper(II) peroxyhemiacetal intermediates, for which optical titrations provide the Keq (at -90 °C) of 73.6 × 102 M-1 (2-PPA) and 10.4 × 102 M-1 (CCA). In the reaction of [Cu2II(XYLO-)(O22-)]1+ with 2-PPA, product complexes characterized by single-crystal X-ray crystallography are the anticipated dicopper(I) complex, [Cu2I(XYLO-)]1+ plus a mixed-valent Cu(I)Cu(II)-formate species. Formate was further identified and confirmed by 1H NMR spectroscopy and electrospray ionization mass spectrometry (ESI-MS) analysis. Using 18O2(g)-isotope labeling the reaction produced a high yield of 18-O incorporated acetophenone as well as formate. The overall results signify that true dioxygenase reactions have occurred, supported by a thorough mechanistic investigation.

Mechanism of the Oxidative Ring-Closure Reaction during Gliotoxin Biosynthesis by Cytochrome P450 GliF

During gliotoxin biosynthesis in fungi, the cytochrome P450 GliF enzyme catalyzes an unusual C-N ring-closure step while also an aromatic ring is hydroxylated in the same reaction cycle, which may have relevance to drug synthesis reactions in biotechnology. However, as the details of the reaction mechanism are still controversial, no applications have been developed yet. To resolve the mechanism of gliotoxin biosynthesis and gain insight into the steps leading to ring-closure, we ran a combination of molecular dynamics and density functional theory calculations on the structure and reactivity of P450 GliF and tested a range of possible reaction mechanisms, pathways and models. The calculations show that, rather than hydrogen atom transfer from the substrate to Compound I, an initial proton transfer transition state is followed by a fast electron transfer en route to the radical intermediate, and hence a non-synchronous hydrogen atom abstraction takes place. The radical intermediate then reacts by OH rebound to the aromatic ring to form a biradical in the substrate that, through ring-closure between the radical centers, gives gliotoxin products. Interestingly, the structure and energetics of the reaction mechanisms appear little affected by the addition of polar groups to the model and hence we predict that the reaction can be catalyzed by other P450 isozymes that also bind the same substrate. Alternative pathways, such as a pathway starting with an electrophilic attack on the arene to form an epoxide, are high in energy and are ruled out.

Machine learning-aided engineering of a cytochrome P450 for optimal bioconversion of lignin fragments

Using machine learning, molecular dynamics simulations, and density functional theory calculations we gain insight into the selectivity patterns of substrate activation by the cytochromes P450. In nature, the reactions catalyzed by the P450s lead to the biodegradation of xenobiotics, but recent work has shown that fungi utilize P450s for the activation of lignin fragments, such as monomer and dimer units. These fragments often are the building blocks of valuable materials, including drug molecules and fragrances, hence a highly selective biocatalyst that can produce these compounds in good yield with high selectivity would be an important step in biotechnology. In this work a detailed computational study is reported on two reaction channels of two P450 isozymes, namely the O-deethylation of guaethol by CYP255A and the O-demethylation versus aromatic hydroxylation of p-anisic acid by CYP199A4. The studies show that the second-coordination sphere plays a major role in substrate binding and positioning, heme access, and in the selectivity patterns. Moreover, the local environment affects the kinetics of the reaction through lowering or raising barrier heights. Furthermore, we predict a site-selective mutation for highly specific reaction channels for CYP199A4.

Mechanistic Insights into Oxidation of Benzaldehyde by Co-Peroxo Complexes

Transition metal-peroxide complexes play a crucial role as intermediates in oxidation reactions. To unravel the mechanism of benzaldehyde oxidation by the Co-peroxo complex, we conducted density functional theory (DFT) calculations. The identified competing mechanisms include nucleophilic attack and hydrogen atom transfer (HAT). The nucleophilic attack pathway involves Co-O cleavage and nucleophilic attack, leading to the formation of the benzoate product. And the HAT pathway comprises O-O cleavage and HAT, ultimately resulting in the benzoate product. DFT calculations revealed that the formation of the end-on Co-superoxo complex 2 through Co-O cleavage, starting from the side-on Co-peroxo complex 1, is much more favorable than the formation of the two-terminal oxyl-radical intermediate 3 through O-O cleavage. Compared with the nucleophilic attack of benzaldehyde by 2, the abstraction of a hydrogen atom from benzaldehyde by 3 requires higher energy. The nature of the nucleophilicity of 2 and 3 accounts for the reactivity of the reaction.

Harnessing the duality of bases toward controlled color and fluorescence

Bases can promote keto-enol tautomerism, a prevalent form of prototropic tautomerism, and facilitate the ring opening of anhydride ring structures. The intrinsic chemical distinctions between these processes provide an opportunity to modulate these seemingly parallel reactions. However, this potential remains largely unexplored. In this work, we report homophthalic anhydride, the first molecule exhibiting simultaneous halochromism, turn-on fluorescence (halofluorochromism), and subsequent self-destruction. Through comprehensive spectroscopic analysis and theoretical calculations, we unravel the mechanisms underlying these phenomena, revealing that the pivotal roles of the base's basicity and nucleophilicity specifically allow us to achieve controlled durations of color change and turn-on fluorescence. Capitalizing on these intriguing properties, we develop a highly dynamic CMY (cyan-magenta-yellow) palette ideal for entity encryption and anti-counterfeiting applications. Our work reshapes the understanding of the relationship between the basicity and nucleophilicity of bases, enriching the comprehension of keto-enol tautomerism and homophthalic anhydride chemistry, and unveils a spectrum of potential applications.

QM/MM Study Into the Mechanism of Oxidative C=C Double Bond Cleavage by Lignostilbene-α,β-Dioxygenase

The enzymatic biosynthesis of fragrance molecules from lignin fragments is an important reaction in biotechnology for the sustainable production of fine chemicals. In this work we investigated the biosynthesis of vanillin from lignostilbene by a nonheme iron dioxygenase using QM/MM and tested several suggested proposals via either an epoxide or dioxetane intermediate. Binding of dioxygen to the active site of the protein results in the formation of an iron(II)-superoxo species with lignostilbene cation radical. The dioxygenase mechanism starts with electrophilic attack of the terminal oxygen atom of the superoxo group on the central C=C bond of lignostilbene, and the second-coordination sphere effects in the substrate binding pocket guide the reaction towards dioxetane formation. The computed mechanism is rationalized with thermochemical cycles and valence bond schemes that explain the electron transfer processes during the reaction mechanism. Particularly, the polarity of the protein and the local electric field and dipole moments enable a facile electron transfer and an exergonic dioxetane formation pathway.

Amphoteric reactivity of a putative Cu(II)-mCPBA intermediate

In copper-based enzymes, Cu-hydroperoxo/alkylperoxo species are proposed as key intermediates for their biological activity. A vast amount of literature is available on the functional and structural mimics of enzymatic systems with heme and non-heme ligand frameworks to stabilize high valent metal intermediates, mostly at low temperatures. Herein, we report a reaction between [CuI(NCCH3)4]+ and meta-chloroperoxybenzoic acid (mCPBA) in CH3CN that produces a putative CuII(mCPBA) species (1). 1 was characterized by UV/Vis, resonance Raman, and EPR spectroscopies. 1 can catalyze both electrophilic and nucleophilic reactions, demonstrating its amphoteric behavior. Additionally, 1 can also conduct electron transfer reactions with a weak reducing agent such as diacetyl ferrocene, making it one of the reactive copper-based intermediates. One of the most important aspects of the current work is the easy synthesis of a CuII(mCPBA) adduct with no complicated ligands for stabilization. Over time, 1 decays to form a CuII paddle wheel complex (2) and is found to be unreactive towards substrate oxidation.

Catalytic divergencies in the mechanism of L-arginine hydroxylating nonheme iron enzymes

Many enzymes in nature utilize a free arginine (L-Arg) amino acid to initiate the biosynthesis of natural products. Examples include nitric oxide synthases, which generate NO from L-Arg for blood pressure control, and various arginine hydroxylases involved in antibiotic biosynthesis. Among the groups of arginine hydroxylases, several enzymes utilize a nonheme iron(II) active site and let L-Arg react with dioxygen and α-ketoglutarate to perform either C3-hydroxylation, C4-hydroxylation, C5-hydroxylation, or C4-C5-desaturation. How these seemingly similar enzymes can react with high specificity and selectivity to form different products remains unknown. Over the past few years, our groups have investigated the mechanisms of L-Arg-activating nonheme iron dioxygenases, including the viomycin biosynthesis enzyme VioC, the naphthyridinomycin biosynthesis enzyme NapI, and the streptothricin biosynthesis enzyme OrfP, using computational approaches and applied molecular dynamics, quantum mechanics on cluster models, and quantum mechanics/molecular mechanics (QM/MM) approaches. These studies not only highlight the differences in substrate and oxidant binding and positioning but also emphasize on electronic and electrostatic differences in the substrate-binding pockets of the enzymes. In particular, due to charge differences in the active site structures, there are changes in the local electric field and electric dipole moment orientations that either strengthen or weaken specific substrate C-H bonds. The local field effects, therefore, influence and guide reaction selectivity and specificity and give the enzymes their unique reactivity patterns. Computational work using either QM/MM or density functional theory (DFT) on cluster models can provide valuable insights into catalytic reaction mechanisms and produce accurate and reliable data that can be used to engineer proteins and synthetic catalysts to perform novel reaction pathways.

Mechanistic Insights into Amphoteric Reactivity of an Iron-Bispidine Complex

The reactivity of FeIII -alkylperoxido complexes has remained a riddle to inorganic chemists owing to their thermal instability and impotency towards organic substrates. These iron-oxygen adducts have been known as sluggish oxidants towards oxidative electrophilic and nucleophilic reactions. Herein, we report the synthesis and spectroscopic characterization of a relatively stable mononuclear high-spin FeIII -alkylperoxido complex supported by an engineered bispidine framework. Against the notion, this FeIII -alkylperoxido complex serves as a rare example of versatile reactivity in both electrophilic and nucleophilic reactions. Detailed mechanistic studies and computational calculations reveal a novel reaction mechanism, where a putative superoxido intermediate orchestrates the amphoteric property of the oxidant. The design of the backbone is pivotal to convey stability and reactivity to alkylperoxido and superoxido intermediates. Contrary to the well-known O-O bond cleavage that generates an FeIV -oxido species, the FeIII -alkylperoxido complex reported here undergoes O-C bond scission to generate a superoxido moiety that is responsible for the amphiphilic reactivity.

Mechanism of CO2 Reduction to Methanol with H2 on an Iron(II)-scorpionate Catalyst

CO2 utilization is an important process in the chemical industry with great environmental power. In this work we show how CO2 and H2 can be reacted to form methanol on an iron(II) center and highlight the bottlenecks for the reaction and what structural features of the catalyst are essential for efficient turnover. The calculations predict the reactions to proceed through three successive reaction cycles that start with heterolytic cleavage of H2 followed by sequential hydride and proton transfer processes. The H2 splitting process is an endergonic process and hence high pressures will be needed to overcome this step and trigger the hydrogenation reaction. Moreover, H2 cleavage into a hydride and proton requires a metal to bind hydride and a nearby source to bind the proton, such as an amide or pyrazolyl group, which the scorpionate ligand used here facilitates. As such the computations highlight the non-innocence of the ligand scaffold through proton shuttle from H2 to substrate as an important step in the reaction mechanism.

Reaction landscape of a mononuclear MnIII-hydroxo complex with hydrogen peroxide

Peroxomanganese species have been proposed as key intermediates in the catalytic cycles of both manganese enzymes and synthetic catalysts. However, many of these intermediates have yet to be observed. Here, we report the formation of a series of intermediates, each generated from the reaction of the mononuclear MnIII-hydroxo complex [MnIII(OH)(dpaq2Me)]+ with hydrogen peroxide under slightly different conditions. By changing the acidity of the reaction mixture and/or the quantity of hydrogen peroxide added, we are able to control which intermediate forms. Using a combination of electronic absorption, 1H NMR, EPR, and X-ray absorption spectroscopies, as well as density functional theory (DFT) and complete active space self-consistent field (CASSCF) calculations, we formulate these intermediates as the bis(μ-oxo)dimanganese(III,IV) complex [MnIIIMnIV(μ-O)2(dpaq2Me)2]+, the MnIII-hydroperoxo complex [MnIII(OOH)(dpaq2Me)]+, and the MnIII-peroxo complex [MnIII(O2)(dpaq2Me)]. The formation of the MnIII-hydroperoxo species from the reaction of a MnIII-hydroxo complex with hydrogen peroxide mimics an elementary reaction proposed for many synthetic manganese catalysts that activate hydrogen peroxide.

A New Thiolate-Bound Dimanganese Cluster as a Structural and Functional Model for Class Ib Ribonucleotide Reductases

In class Ib ribonucleotide reductases (RNRs) a dimanganese(II) cluster activates superoxide (O₂ ⋅^- ) rather than dioxygen (O₂ ), to access a high valent Mn^III -O₂ -Mn^IV species, responsible for the oxidation of tyrosine to tyrosyl radical. In a biomimetic approach, we report the synthesis of a thiolate-bound dimanganese complex [Mn^II ₂ (BPMT)(OAc)₂ ](ClO)₄ (BPMT=(2,6-bis{[bis(2-pyridylmethyl)amino]methyl}-4-methylthiophenolate) (1) and its reaction with O₂ ⋅^- to form a [(BPMT)MnO₂ Mn]²⁺ complex 2. Resonance Raman investigation revealed the presence of an O-O bond in 2, while EPR analysis displayed a 16-line S_t =1/2 signal at g=2 typically associated with a Mn^III Mn^IV core, as detected in class Ib RNRs. Unlike all other previously reported Mn-O₂ -Mn complexes, generated by O₂ ⋅^- activation at Mn₂ centers, 2 proved to be a capable electrophilic oxidant in aldehyde deformylation and phenol oxidation reactions, rendering it one of the best structural and functional models for class Ib RNRs.

Oxidation of Aldehydes into Carboxylic Acids by a Mononuclear Manganese(III) Iodosylbenzene Complex through Electrophilic C-H Bond Activation

The oxidation of aldehyde is one of the fundamental reactions in the biological system. Various synthetic procedures and catalysts have been developed to convert aldehydes into corresponding carboxylic acids efficiently under ambient conditions. In this work, we report the oxidation of aldehydes by a mononuclear manganese(III) iodosylbenzene complex, [Mn^III(TBDAP)(OIPh)(OH)]²⁺ (1), with kinetic and mechanistic studies in detail. The reaction of 1 with aldehydes resulted in the formation of corresponding carboxylic acids via a pre-equilibrium state. Hammett plot and reaction rates of 1 with 1°-, 2°-, and 3°-aldehydes revealed the electrophilicity of 1 in the aldehyde oxidation. A kinetic isotope effect experiment and reactivity of 1 toward cyclohexanecarboxaldehyde (CCA) analogues indicate that the reaction of 1 with aldehyde occurs through the rate-determining C-H bond activation at the formyl group. The reaction rate of 1 with CCA is correlated to the bond dissociation energy of the formyl group plotting a linear correlation with other aliphatic C-H bonds. Density functional theory calculations found that 1 electrostatically interacts with CCA at the pre-equilibrium state in which the C-H bond activation of the formyl group is performed as the most feasible pathway. Surprisingly, the rate-determining step is characterized as hydride transfer from CCA to 1, affording an (oxo)methylium intermediate. At the fundamental level, it is revealed that the hydride transfer is composed of H atom abstraction followed by a fast electron transfer. Catalytic reactions of aldehydes by 1 are also presented with a broad substrate scope. This novel mechanistic study gives better insights into the metal oxygen chemistry and would be prominently valuable for development of transition metal catalysts.

Mechanism of Melatonin Metabolism by CYP1A1: What Determines the Bifurcation Pathways of Hydroxylation versus Deformylation?

Melatonin, a widely applied cosmetic active ingredient, has a variety of uses as a skin protector through antioxidant and anti-inflammatory functions as well as giving the body UV-induced defenses and immune system support. In the body, melatonin is synthesized from a tryptophan amino acid in a cascade of reactions, but as melatonin is toxic at high concentrations, it is metabolized in the human skin by the cytochrome P450 enzymes. The P450s are diverse heme-based mono-oxygenases that catalyze oxygen atom-transfer processes that trigger metabolism and detoxification reactions in the body. In the catalytic cycle of the P450s, a short-lived high-valent iron(IV)-oxo heme cation radical is formed that has been proposed to be the active oxidant. How and why it activates melatonin in the human body and what the origin of the product distributions is, are unknown. This encouraged us to do a detailed computational study on a typical human P450 isozyme, namely CYP1A1. We initially did a series of molecular dynamics simulations with substrate docked into several orientations. These simulations reveal a number of stable substrate-bound positions in the active site, which may lead to differences in substrate activation channels. Using tunneling analysis on the full protein structures, we show that two of the four binding conformations lead to open substrate-binding pockets. As a result, in these open pockets, the substrate is not tightly bound and can escape back into the solution. In the closed conformations, in contrast, the substrate is mainly oriented with the methoxy group pointing toward the heme, although under a different angle. We then created large quantum cluster models of the enzyme and focused on the chemical reaction mechanisms for melatonin activation, leading to competitive O-demethylation and C⁶-aromatic hydroxylation pathways. The calculations show that active site positioning determines the product distributions, but the bond that is activated is not necessarily closest to the heme in the enzyme-substrate complex. As such, the docking and molecular dynamics positioning of the substrate versus oxidant can give misleading predictions on product distributions. In particular, in quantum mechanics cluster model I, we observe that through a tight hydrogen bonding network, a preferential 6-hydroxylation of melatonin is obtained. However, O-demethylation becomes possible in alternative substrate-binding orientations that have the C⁶-aromatic ring position shielded. Finally, we investigated enzymatic and non-enzymatic O-demethylation processes and show that the hydrogen bonding network in the substrate-binding pocket can assist and perform this step prior to product release from the enzyme.

Peroxide-Selective Reduction of O2 at Redox-Inactive Rare-Earth(III) Triflates Generates an Ambiphilic Peroxide

Metal peroxides are key species involved in a range of critical biological and synthetic processes. Rare-earth (group III and the lanthanides; Sc, Y, La-Lu) peroxides have been implicated as reactive intermediates in catalysis; however, reactivity studies of isolated, structurally characterized rare-earth peroxides have been limited. Herein, we report the peroxide-selective (93-99% O₂^2-) reduction of dioxygen (O₂) at redox-inactive rare-earth triflates in methanol using a mild metallocene reductant, decamethylferrocene (Fc*). The first molecular praseodymium peroxide ([Pr^III₂(O₂^2-)(18C6)₂(EG)₂][OTf]₄; 18C6 = 18-crown-6, EG = ethylene glycol, ^-OTf = ^-O₃SCF₃; 2-Pr) was isolated and characterized by single-crystal X-ray diffraction, Raman spectroscopy, and NMR spectroscopy. 2-Pr displays high thermal stability (120 °C, 50 mTorr), is protonated by mild organic acids [pK_a1(MeOH) = 5.09 ± 0.23], and engages in electrophilic (e.g., oxygen atom transfer) and nucleophilic (e.g., phosphate-ester cleavage) reactivity. Our mechanistic studies reveal that the rate of oxygen reduction is dictated by metal-ion accessibility, rather than Lewis acidity, and suggest new opportunities for differentiated reactivity of redox-inactive metal ions by leveraging weak metal-ligand binding events preceding electron transfer.

DFT Mechanistic Insights into Aldehyde Deformylations with Biomimetic Metal-Dioxygen Complexes: Distinct Mechanisms and Reaction Rules

Aldehyde deformylations occurring in organisms are catalyzed by metalloenzymes through metal-dioxygen active cores, attracting great interest to study small-molecule metal-dioxygen complexes for understanding relevant biological processes and developing biomimetic catalysts for aerobic transformations. As the known deformylation mechanisms, including nucleophilic attack, aldehyde α-H-atom abstraction, and aldehyde hydrogen atom abstraction, undergo outer-sphere pathways, we herein report a distinct inner-sphere mechanism based on density functional theory (DFT) mechanistic studies of aldehyde deformylations with a copper (II)-superoxo complex. The inner-sphere mechanism proceeds via a sequence mainly including aldehyde end-on coordination, homolytic aldehyde C-C bond cleavage, and dioxygen O-O bond cleavage, among which the C-C bond cleavage is the rate-determining step with a barrier substantially lower than those of outer-sphere pathways. The aldehyde C-C bond cleavage, enabled through the activation of the dioxygen ligand radical in a second-order nucleophilic substitution (S_N2)-like fashion, leads to an alkyl radical and facilitates the subsequent dioxygen O-O bond cleavage. Furthermore, we deduced the rules for the reactions of metal-dioxygen complexes with aldehydes and nitriles via the inner-sphere mechanism. Expectedly, our proposed inner-sphere mechanisms and the reaction rules offer another perspective to understand relevant biological processes involving metal-dioxygen cores and to discover metal-dioxygen catalysts for aerobic transformations.

Bioinspired mononuclear Mn complexes for O2 activation and biologically relevant reactions

A general interest in harnessing the oxidizing power of dioxygen (O₂) continues to motivate research efforts on bioinspired and biomimetic complexes to better understand how metalloenzymes mediate these reactions. The ubiquity of Fe- and Cu-based enzymes attracts significant attention and has resulted in many noteworthy developments for abiotic systems interested in direct O₂ reduction and small molecule activation. However, despite the existence of Mn-based metalloenzymes with important O₂-dependent activity, there has been comparatively less focus on the development of these analogues relative to Fe- and Cu-systems. In this Perspective, we summarize important contributions to the development of bioinspired mononuclear Mn complexes for O₂ activation and studies on their reactivity, emphasizing important design parameters in the primary and secondary coordination spheres and outlining mechanistic trends.

The influence of secondary interactions on the [Ni(O2)]+ mediated aldehyde oxidation reactions

A rate enhancement of one to two orders of magnitude can be obtained in the aldehyde deformylation reactions by replacing the -N(CH₃) groups of [Ni^III(O₂)(Me₄[12]aneN₄)]⁺ (Me₄[12]aneN₄ = 1,4,7,10-tetramethyl-1,4,7,10-tetraazacyclododecane) and [Ni^III(O₂)(Me₄[13]aneN₄)]⁺ (Me₄[13]aneN₄ = 1,4,7,10-tetramethyl-1,4,7,10-tetraazacyclotridecane) complexes by -NH in [Ni^III(O₂)([12]aneN₄)]⁺ (2; [12]aneN₄ = 1,4,7,10-tetraazacyclododecane) and [Ni^III(O₂)([13]aneN₄)]⁺ (4; [13]aneN₄ = 1,4,7,10-tetraazacyclotridecane). Based on detailed spectroscopic, reaction-kinetics and theoretical investigations, the higher reactivities of 2 and 4 are attributed to the changes in the secondary-sphere interactions between the [Ni^III(O₂)]⁺ and [12]aneN₄ or [13]aneN₄ moieties, which open up an alternative electrophilic pathway for the aldehyde oxidation reaction. Identification of primary kinetic isotope effects on the reactivity and stability of 2 when the -NH groups of the [12]aneN₄ ligand are deuterated may also suggest the presence of secondary interaction between the -NH groups of [12]aneN₄ and [Ni^III(O₂)]⁺ moieties, although, such interactions are not obvious in the DFT calculated optimized structure at the employed level of theory.

Spectroscopic Evidence for a Cobalt-Bound Peroxyhemiacetal Intermediate

Aldehyde deformylation reactions by metal dioxygen adducts have been proposed to involve peroxyhemiacetal species as key intermediates. However, direct evidence of such intermediates has not been obtained to date. We report the spectroscopic characterization of a mononuclear cobalt(III)-peroxyhemiacetal complex, [Co(Me₃-TPADP)(O₂CH(O)CH(CH₃)C₆H₅)]⁺ (2), in the reaction of a cobalt(III)-peroxo complex (1) with 2-phenylpropionaldehyde (2-PPA). The formation of 2 is also investigated by isotope labeling experiments and kinetic studies. The conclusion that the peroxyhemiacetalcobalt(III) intermediate is responsible for the aldehyde deformylation is supported by the product analyses. Furthermore, isotopic labeling suggests that the reactivity of the cobalt(III)-peroxo complex depends on the second reactant. The aldehyde inserts between the oxygen atoms of 1, whereas the reaction with acyl chlorides proceeds by a nucleophilic attack. The observation of the peroxyhemiacetal intermediate provides significant insight into the initial step of aldehyde deformylation by metalloenzymes.

A Mononuclear Non-heme Iron(III)-Peroxo Complex with an Unprecedented High O-O Stretch and Electrophilic Reactivity

A mononuclear non-heme iron(III)-peroxo complex, [Fe(III)(O₂)(13-TMC)]⁺ (1), was synthesized and characterized spectroscopically; the characterization with electron paramagnetic resonance, Mössbauer, X-ray absorption, and resonance Raman spectroscopies and mass spectrometry supported a high-spin S = ⁵/₂ Fe(III) species binding an O₂ unit. A notable observation was an unusually high ν_O-O at ∼1000 cm^-1 for the peroxo ligand. With regard to reactivity, 1 showed electrophilic reactivity in H atom abstraction (HAA) and O atom transfer (OAT) reactions. In the HAT reaction, a kinetic isotope effect (KIE) value of 5.8 was obtained in the oxidation of 9,10-dihydroanthracene. In the OAT reaction, a negative ρ value of -0.61 in the Hammett plot was determined in the oxidation of p-X-substituted thioanisoles. Another interesting observation was the electrophilic reactivity of 1 in the oxidation of benzaldehyde derivatives, such as a negative ρ value of -0.77 in the Hammett plot and a KIE value of 2.2. To the best of our knowledge, the present study reports the first example of a mononuclear non-heme iron(III)-peroxo complex with an unusually high ν_O-O value and unprecedented electrophilic reactivity in oxidation reactions.

Electrostatic Perturbations from the Protein Affect C-H Bond Strengths of the Substrate and Enable Negative Catalysis in the TmpA Biosynthesis Enzyme

The nonheme iron dioxygenase 2-(trimethylammonio)-ethylphosphonate dioxygenase (TmpA) is an enzyme involved in the regio- and chemoselective hydroxylation at the C¹ -position of the substrate as part of the biosynthesis of glycine betaine in bacteria and carnitine in humans. To understand how the enzyme avoids breaking the weak C² -H bond in favor of C¹ -hydroxylation, we set up a cluster model of 242 atoms representing the first and second coordination sphere of the metal center and substrate binding pocket, and investigated possible reaction mechanisms of substrate activation by an iron(IV)-oxo species by density functional theory methods. In agreement with experimental product distributions, the calculations predict a favorable C¹ -hydroxylation pathway. The calculations show that the selectivity is guided through electrostatic perturbations inside the protein from charged residues, external electric fields and electric dipole moments. In particular, charged residues influence and perturb the homolytic bond strength of the C¹ -H and C² -H bonds of the substrate, and strongly strengthens the C² -H bond in the substrate-bound orientation.

A comprehensive insight into aldehyde deformylation: mechanistic implications from biology and chemistry

Aldehyde deformylation is an important reaction in biology, organic chemistry and inorganic chemistry and the process has been widely applied and utilized. For instance, in biology, the aldehyde deformylation reaction has wide differences in biological function, whereby cyanobacteria convert aldehydes into alkanes or alkenes, which are used as natural products for, e.g., defense mechanisms. By contrast, the cytochromes P450 catalyse the biosynthesis of hormones, such as estrogen, through an aldehyde deformylation reaction step. In organic chemistry, the aldehyde deformylation reaction is a common process for replacing functional groups on a molecule, and as such, many different synthetic methods and procedures have been reported that involve an aldehyde deformylation step. In bioinorganic chemistry, a variety of metal(iii)-peroxo complexes have been synthesized as biomimetic models and shown to react efficiently with aldehydes through deformylation reactions. This review paper provides an overview of the various aldehyde deformylation reactions in organic chemistry, biology and biomimetic model systems, and shows a broad range of different chemical reaction mechanisms for this process. Although a nucleophilic attack at the carbonyl centre is the consensus reaction mechanism, several examples of an alternative electrophilic reaction mechanism starting with hydrogen atom abstraction have been reported as well. There is still much to learn and to discover on aldehyde deformylation reactions, as deciphered in this review paper.

A side-on Mn(III)-peroxo supported by a non-heme pentadentate N3Py2 ligand: synthesis, characterization and reactivity studies

A mononuclear manganese(iii)-peroxo complex [MnIII(N3Py2)(O2)]+ (1a) bearing a non-heme N,N'-dimethyl-N-(2-(methyl(pyridin-2-ylmethyl)amino)ethyl)-N'-(pyridin-2-ylmethyl)ethane-1,2-diamine (N3Py2) ligand was synthesized by the reaction of [Mn(N3Py2)(H2O)](ClO4)2 (1) with hydrogen peroxide and triethylamine in CH3CN at 25 °C. The reactivity of 1a in aldehyde deformylation using 2-phenyl propionaldehyde (2-PPA) was studied and the reaction kinetics was monitored by UV-visible spectroscopy. A kinetic isotope effect (KIE) = 1.7 was obtained in the reaction of 1a with 2-PPA and α-[D1]-PPA, suggesting nucleophilic character of 1a. The activation parameters ΔH‡ and ΔS‡ were determined using the Eyring plot while Ea was obtained from the Arrhenius equation by performing the reaction between 288 and 303 K. Hammett constants (σp) of para-substituted benzaldehydes p-X-Ph-CHO (X = Cl, F, H, and Me) were linear with a slope (ρ) = 3.0. Computational study suggested that the side-on structure of 1a is more favored over the end-on structure and facilitates the reactivity of 1a.

Mechanistic insights into dioxygen activation by a manganese corrole complex: a broken-symmetry DFT study

The Mn–oxygen species have been implicated as key intermediates in various Mn-mediated oxidation reactions. However, artificial oxidants were often used for the synthesis of the Mn–oxygen intermediates. Remarkably, the Mn(v)–oxo and Mn(iv)–peroxo species have been observed in the activation of O₂ by Mn(iii) corroles in the presence of base (OH⁻) and hydrogen donors. In this work, density functional theory methods were used to get insight into the mechanism of dioxygen activation and formation of Mn(v)–oxo. The results demonstrated that the dioxygen cannot bind to Mn without the axial OH⁻ ligand. Upon the addition of the axial OH⁻ ligand, the dioxygen can bind to Mn in an end-on fashion to give the Mn(iv)–superoxo species. The hydrogen atom transfer from the hydrogen donor (substrate) to the Mn(iv)–superoxo species is the rate-limiting step, having a high reaction barrier and a large endothermicity. Subsequently, the O–C bond formation is concerted with an electron transfer from the substrate radical to the Mn and a proton transfer from the hydroperoxo moiety to the nearby N atom of the corrole ring, generating an alkylperoxo Mn(iii) complex. The alkylperoxo O–O bond cleavage affords a Mn(v)–oxo complex and a hydroxylated substrate. This novel mechanism for the Mn(v)–oxo formation via an alkylperoxo Mn(iii) intermediate gives insight into the O–O bond activation by manganese complexes.

DFT calculations revealed a novel mechanism for the formation of Mn(v)–oxo in the dioxygen activation by a Mn(iii) corrole complex involving a Mn(iii)–alkylperoxo intermediate.

Lignin Biodegradation by a Cytochrome P450 Enzyme: A Computational Study into Syringol Activation by GcoA

A recently characterized cytochrome P450 isozyme GcoA activates lignin components through a selective O-demethylation or alternatively an acetal formation reaction. These are important reactions in biotechnology and, because lignin is readily available; it being the main component in plant cell walls. In this work we present a density functional theory study on a large active site model of GcoA to investigate syringol activation by an iron(IV)-oxo heme cation radical oxidant (Compound I) leading to hemiacetal and acetal products. Several substrate-binding positions were tested and full energy landscapes calculated. The study shows that substrate positioning determines the product distributions. Thus, with the phenol group pointing away from the heme, an O-demethylation is predicted, whereas an initial hydrogen-atom abstraction of the weak phenolic O-H group would trigger a pathway leading to ring-closure to form acetal products. Predictions on how to engineer P450 GcoA to get more selective product distributions are given.

Catalytic Mechanism of Aromatic Nitration by Cytochrome P450 TxtE: Involvement of a Ferric-Peroxynitrite Intermediate

The cytochromes P450 are heme-dependent enzymes that catalyze many vital reaction processes in the human body related to biodegradation and biosynthesis. They typically act as mono-oxygenases; however, the recently discovered P450 subfamily TxtE utilizes O₂ and NO to nitrate aromatic substrates such as L-tryptophan. A direct and selective aromatic nitration reaction may be useful in biotechnology for the synthesis of drugs or small molecules. Details of the catalytic mechanism are unknown, and it has been suggested that the reaction should proceed through either an iron(III)-superoxo or an iron(II)-nitrosyl intermediate. To resolve this controversy, we used stopped-flow kinetics to provide evidence for a catalytic cycle where dioxygen binds prior to NO to generate an active iron(III)-peroxynitrite species that is able to nitrate l-Trp efficiently. We show that the rate of binding of O₂ is faster than that of NO and also leads to l-Trp nitration, while little evidence of product formation is observed from the iron(II)-nitrosyl complex. To support the experimental studies, we performed density functional theory studies on large active site cluster models. The studies suggest a mechanism involving an iron(III)-peroxynitrite that splits homolytically to form an iron(IV)-oxo heme (Compound II) and a free NO₂ radical via a small free energy of activation. The latter activates the substrate on the aromatic ring, while compound II picks up the ipso-hydrogen to form the product. The calculations give small reaction barriers for most steps in the catalytic cycle and, therefore, predict fast product formation from the iron(III)-peroxynitrite complex. These findings provide the first detailed insight into the mechanism of nitration by a member of the TxtE subfamily and highlight how the enzyme facilitates this novel reaction chemistry.

Bioengineering of Cytochrome P450 OleTJE: How Does Substrate Positioning Affect the Product Distributions?

The cytochromes P450 are versatile enzymes found in all forms of life. Most P450s use dioxygen on a heme center to activate substrates, but one class of P450s utilizes hydrogen peroxide instead. Within the class of P450 peroxygenases, the P450 OleT_JE isozyme binds fatty acid substrates and converts them into a range of products through the α-hydroxylation, β-hydroxylation and decarboxylation of the substrate. The latter produces hydrocarbon products and hence can be used as biofuels. The origin of these product distributions is unclear, and, as such, we decided to investigate substrate positioning in the active site and find out what the effect is on the chemoselectivity of the reaction. In this work we present a detailed computational study on the wild-type and engineered structures of P450 OleT_JE using a combination of density functional theory and quantum mechanics/molecular mechanics methods. We initially explore the wild-type structure with a variety of methods and models and show that various substrate activation transition states are close in energy and hence small perturbations as through the protein may affect product distributions. We then engineered the protein by generating an in silico model of the double mutant Asn242Arg/Arg245Asn that moves the position of an active site Arg residue in the substrate-binding pocket that is known to form a salt-bridge with the substrate. The substrate activation by the iron(IV)-oxo heme cation radical species (Compound I) was again studied using quantum mechanics/molecular mechanics (QM/MM) methods. Dramatic differences in reactivity patterns, barrier heights and structure are seen, which shows the importance of correct substrate positioning in the protein and the effect of the second-coordination sphere on the selectivity and activity of enzymes.

Rhodoxanthin synthase from honeysuckle; a membrane diiron enzyme catalyzes the multistep conversation of β-carotene to rhodoxanthin

Rhodoxanthin is a vibrant red carotenoid found across the plant kingdom and in certain birds and fish. It is a member of the atypical retro class of carotenoids, which contain an additional double bond and a concerted shift of the conjugated double bonds relative to the more widely occurring carotenoid pigments, and whose biosynthetic origins have long remained elusive. Here, we identify LHRS (Lonicera hydroxylase rhodoxanthin synthase), a variant β-carotene hydroxylase (BCH)-type integral membrane diiron enzyme that mediates the conversion of β-carotene into rhodoxanthin. We identify residues that are critical to rhodoxanthin formation by LHRS. Substitution of only three residues converts a typical BCH into a multifunctional enzyme that mediates a multistep pathway from β-carotene to rhodoxanthin via a series of distinct oxidation steps in which the product of each step becomes the substrate for the next catalytic cycle. We propose a biosynthetic pathway from β-carotene to rhodoxanthin.

Peroxocobalt(iii) species activates nitriles via a superoxocobalt(ii) diradical state

The dioxygenation of nitriles by [CoIII(TBDAP)(O2)]+ (TBDAP = N,N-di-tert-butyl-2,11-diaza[3.3](2,6)-pyridinophane) is investigated using DFT-calculations. The mechanism proposed previously based on experimental observations, which invoked an outer-sphere cycloaddition, was found to be unreasonable. Instead, calculations suggest that an inner-sphere mechanism involving the cleavage of one of the Co-O bonds assisted by substrate uptake is much more likely. The reactively competent species is a triplet consisting of a Co(ii)-superoxo functionality, which can undergo O-C bond formation and O-O bond cleavage traversing low energy transition states. The role of the structurally rigid TBDAP ligand is to prevent the participation of the pyridyl ligand in the delocalization of the unpaired electron density.

Insights into plasmon induced keto-enol isomerization

Chemical reactions that are driven by plasmon-induced hot carriers are a timely topic of interest to chemists and material scientists as they provide catalytic alternatives that may reduce cost and/or waste. Herein, we monitored the localized surface plasmon resonance-induced keto-enol isomerization process of 2-mercapto-4(3H)-quinazolinone (MQ) by time-dependent surface enhanced Raman scattering (SERS), where the MQ molecules are adsorbed on gold nanoparticles (GNP) surface by Au-S bonds. The mechanism of keto-enol isomerization has been successfully investigated, and it is found that the isomerization is induced by hot hole transfer from GNPs to the adsorbed molecules. The present investigation could provide significant insights into hot hole catalyzed chemical reactions via SERS spectra and theoretical calculations.

Mechanistic dichotomies in redox reactions of mononuclear metal–oxygen intermediates

This review article focuses on various mechanistic dichotomies in redox reactions of metal–oxygen intermediates with the emphasis on understanding and controlling their redox reactivity from experimental and theoretical points of view.

Selective Hydrogen Atom Abstraction from Dihydroflavonol by a Nonheme Iron Center Is the Key Step in the Enzymatic Flavonol Synthesis and Avoids Byproducts

The plant non-heme iron dioxygenase flavonol synthase performs a regioselective desaturation reaction as part of the biosynthesis of the signaling molecule flavonol that triggers the growing of leaves and flowers. These compounds also have health benefits for humans. Desaturation of aliphatic compounds generally proceeds through two consecutive hydrogen atom abstraction steps from two adjacent carbon atoms and in nature often is performed by a high-valent iron(IV)-oxo species. We show that the order of the hydrogen atom abstraction steps, however, is opposite of those expected from the C-H bond strengths in the substrate and determines the product distributions. As such, flavonol synthase follows a negative catalysis mechanism. Using density functional theory methods on large active-site model complexes, we investigated pathways for desaturation and hydroxylation by an iron(IV)-oxo active-site model. Contrary to thermochemical predictions, we find that the oxidant abstracts the hydrogen atom from the strong C²-H bond rather than the weaker C³-H bond of the substrate first. We analyze the origin of this unexpected selective hydrogen atom abstraction pathway and find that the alternative C³-H hydrogen atom abstraction would be followed by a low-energy and competitive substrate hydroxylation mechanism hence, should give considerable amount of byproducts. Our computational modeling studies show that substrate positioning in flavonol synthase is essential, as it guides the reactivity to a chemo- and regioselective substrate desaturation from the C²-H group, leading to desaturation products efficiently.

Properties and reactivity of μ-nitrido-bridged dimetal porphyrinoid complexes: how does ruthenium compare to iron?

Methane hydroxylation by metal-oxo oxidants is one of the Holy Grails in biomimetic and biotechnological chemistry. The only enzymes known to perform this reaction in Nature are iron-containing soluble methane monooxygenase and copper-containing particulate methane monooxygenase. Furthermore, few biomimetic iron-containing oxidants have been designed that can hydroxylate methane efficiently. Recent studies reported that μ-nitrido-bridged diiron(IV)-oxo porphyrin and phthalocyanine complexes hydroxylate methane to methanol efficiently. To find out whether the reaction rates are enhanced by replacing iron by ruthenium, we performed a detailed computational study. Our work shows that the μ-nitrido-bridged diruthenium(IV)-oxo reacts with methane via hydrogen atom abstraction barriers that are considerably lower in energy (by about 5 kcal mol^‒1) as compared to the analogous diiron(IV)-oxo complex. An analysis of the electronic structure implicates similar spin and charge distributions for the diiron(IV)-oxo and diruthenium(IV)-oxo complexes, but the strength of the O‒H bond formed during the reaction is much stronger for the latter. As such a larger hydrogen atom abstraction driving force for the Ru complex than for the Fe complex is found, which should result in higher reactivity in the oxidation of methane.

Mechanism of Oxidative Activation of Fluorinated Aromatic Compounds by N-Bridged Diiron-Phthalocyanine: What Determines the Reactivity?

The biodegradation of compounds with C-F bonds is challenging due to the fact that these bonds are stronger than the C-H bond in methane. In this work, results on the unprecedented reactivity of a biomimetic model complex that contains an N-bridged diiron-phthalocyanine are presented; this model complex is shown to react with perfluorinated arenes under addition of H₂ O₂ effectively. To get mechanistic insight into this unusual reactivity, detailed density functional theory calculations on the mechanism of C₆ F₆ activation by an iron(IV)-oxo active species of the N-bridged diiron phthalocyanine system were performed. Our studies show that the reaction proceeds through a rate-determining electrophilic C-O addition reaction followed by a 1,2-fluoride shift to give the ketone product, which can further rearrange to the phenol. A thermochemical analysis shows that the weakest C-F bond is the aliphatic C-F bond in the ketone intermediate. The oxidative defluorination of perfluoroaromatics is demonstrated to proceed through a completely different mechanism compared to that of aromatic C-H hydroxylation by iron(IV)-oxo intermediates such as cytochrome P450 Compound I.

Hydrogen by Deuterium Substitution in an Aldehyde Tunes the Regioselectivity by a Nonheme Manganese(III)-Peroxo Complex

Mononuclear nonheme Mn^III -peroxo complexes are important intermediates in biology, and take part in oxygen activation by photosystem II. Herein, we present work on two isomeric biomimetic side-on Mn^III -peroxo intermediates with bispidine ligand system and reactivity patterns with aldehydes. The complexes are characterized with UV/Vis and mass spectrometric techniques and reaction rates with cyclohexane carboxaldehyde (CCA) are measured. The reaction gives an unusual regioselectivity switch from aliphatic to aldehyde hydrogen atom abstraction upon deuteration of the substrate, leading to the corresponding carboxylic acid product for the latter, while the former gives a deformylation reaction. Mechanistic details are established from kinetic isotope effect studies and density functional theory calculations. Thus, replacement of C-H by C-D raises the hydrogen atom abstraction barriers and enables a regioselectivity switch to a competitive pathway that is slightly higher in energy.

Flavonol biosynthesis by nonheme iron dioxygenases: A computational study into the structure and mechanism

Plants produce flavonol compounds for vital functions regarding plant growth, fruit and flower colouring as well as fruit ripening processes. Several of these biosynthesis steps are stereo- and regioselective and are being carried out by nonheme iron enzymes. Using density functional theory calculations on a large active site model complex of flavanone-3β-hydroxylase (FHT), we established the mechanism for conversion of naringenin to its dihydroflavonol, which is a key step in the mechanism of flavonol biosynthesis. The reaction starts with dioxygen binding to the iron(II) centre and a reaction with α-ketoglutarate co-substrate gives succinate, an iron(IV)-oxo species and CO₂ with large exothermicity and small reaction barriers. The rate-determining reaction step in the mechanism; however, is hydrogen atom abstraction of an aliphatic CH bond by the iron(IV)-oxo species. We identify a large kinetic isotope effect for the replacement of the transferring hydrogen atom by deuterium. In a final step the OH and substrate radicals combine to form the alcohol product with a barrier of several kcal mol^-1. We show that the latter is the result of geometric constraints in the active site pocket. Furthermore, the calculations show that a weak tertiary CH bond is shielded from the iron(IV)-oxo species in the substrate binding position and therefore the enzyme is able to activate a stronger CH bond. As such, the flavanone-3β-hydroxylase enzyme reacts regioselectively with one specific CH bond of naringenin by avoiding activation of weaker bonds through tight substrate and oxidant positioning.