Structure and Mechanism Leading to Formation of the Cysteine Sulfinate Product Complex of a Biomimetic Cysteine Dioxygenase Model

A density functional theory analysis of the C-H activation reactivity of iron(IV)-oxo complexes with an 'O' substituted tetramethylcyclam macrocycle

In this article, we present a meticulous computational study to foresee the effect of an oxygen-rich macrocycle on the reactivity for C-H activation. For this study, a widely studied nonheme Fe(IV)O molecule with a TMC (1,4,8,11-tetramethyl 1,4,8,11-tetraazacyclotetradecane) macrocycle that is equatorially attached to four nitrogen atoms (designated as N4) and acetonitrile as an axial ligand has been taken into account. For the goal of hetero-substitution, step-by-step replacement of the N4 framework with O atoms, i.e., N4, N3O1, N2O2, N1O3, and O4 systems, has been considered, and dihydroanthracene (DHA) has been used as the substrate. In order to neutralise the system and prevent the self-interaction error in DFT, triflate counterions have also been included in the calculations. The study of the energetics of these C-H bond activation reactions and the potential energy surfaces mapped therefore reveal that the initial hydrogen abstraction, which is the rate-determining step, follows the two-state reactivity (TSR) pattern, which means that the originally excited quintet state falls lower in the transition state and the product. The reaction follows the hydrogen atom transfer (HAT) mechanism, as indicated by the spin density studies. The results revealed a fascinating reactivity order, in which the reactivity increases with the enrichment of the oxygen atom in the equatorial position, namely the order follows N4 < N3O1 < N2O2 < N1O3 < O4. The impacts of oxygen substitution on quantum mechanical tunneling and the H/D kinetic isotope effect have also been investigated. When analysing the causes of this reactivity pattern, a number of variables have been identified, including the reactant-like transition structure, spin density distribution, distortion energy, and energies of the electron acceptor orbital, i.e., the energy of the LUMO (σ*z2), which validate the obtained outcome. Our results also show very good agreement with earlier combined experimental and theoretical studies considering TMC and TMCO-type complexes. The DFT predictions reported here will undoubtedly encourage experimental research in this biomimetic field, as they provide an alternative with higher reactivity in which heteroatoms can be substituted for the traditional nitrogen atom.

Enhanced Reactivity through Equatorial Sulfur Coordination in Nonheme Iron(IV)-Oxo Complexes: Insights from Experiment and Theory

Sulfur ligation in metalloenzymes often gives the active site unique properties, whether it is the axial cysteinate ligand in the cytochrome P450s or the equatorial sulfur/thiol ligation in nonheme iron enzymes. To understand sulfur ligation to iron complexes and how it affects the structural, spectroscopic, and intrinsic properties of the active species and the catalysis of substrates, we pursued a systematic study and compared sulfur with amine-ligated iron(IV)-oxo complexes. We synthesized and characterized a biomimetic N4S-ligated iron(IV)-oxo complex and compared the obtained results with an analogous N5-ligated iron(IV)-oxo complex. Our work shows that the amine for sulfur replacement in the equatorial ligand framework leads to a rate enhancement for oxygen atom and hydrogen atom transfer reactions. Moreover, the sulfur-ligated iron(IV)-oxo complex reacts through a different reaction mechanism as compared to the N5-ligated iron(IV)-oxo complex, where the former reacts through hydride transfer with the latter reacting via radical pathways. We show that the reactivity differences are caused by a dramatic change in redox potential between the two complexes. Our studies highlight the importance of implementing a sulfur ligand into the equatorial ligand framework of nonheme iron(IV)-oxo complexes and how it affects the physicochemical properties of the oxidant and its reactivity.

DFT insights into the mechanism of O2 activation catalyzed by a structural and functional model of cysteine dioxygenase with tris(2-pyridyl)methane-based ligand architecture

Cysteine dioxygenation is an important step in the metabolism of toxic L-cysteine (Cys) in the human body, carried out by cysteine dioxygenase enzyme (CDO). The disruption of this process is found to elicit neurological health issues. This work reports a computational investigation of mechanistic aspects of this reaction, using a recently reported tris(2-pyridyl)methane-based biomimetic model complex of CDO. The computed results indicate that, the initial SO2 bond formation process is the slowest step in the S-dioxygenation process, possessing an activation barrier of 12.7 kcal/mol. The remaining steps were found to be downhill requiring very small activation energies. The transition states were found to undergo spin crossover between triplet and quintet states, while the singlet surface remained unstable throughout the entire reaction. In essence, the mechanistic scheme and multistate reactivity pattern together with the relatively small computed rate-limiting activation barrier as well as the exothermic formation energy demonstrate that the model complex is an efficient biomimetic CDO model. In addition, the study also substantiates the involvement of Fe(IV)oxido intermediates in the mechanism of S-dioxygenation by the chosen model complex. The insights derived from the O₂ activation process might pave way for development of more accurate CDO model catalysts that might be capable of even more efficiently mimicking the geometric, spectroscopic and functional features of the CDO enzyme.

Role of "S" Substitution on C-H Activation Reactivity of Iron(IV)-Oxo Cyclam Complexes: a Computational Investigation

A comprehensive density functional theory (DFT) investigation has been presented in this article to address the role of equatorial sulfur ligation in C-H activation. A non-heme iron-oxo compound with four nitrogen atoms constituting the equatorially connected macrocyclic framework (represented as N₄) [Fe(IV)═O(THC)(CH₃CN)]²⁺(THC = 1,4,8,11-tetrahydro1,4,8,11-tetraazacyclotetradecane) has been considered as the base compound. Other complexes have been anticipated by the sequential replacement of this nitrogen by sulfur, that is, N₄, N₃S₁, N₂S₂, N₁S₃, and S₄. Counterions, as always, have been considered to avoid the self-interaction error in DFT. Generally, the anti-conformers (with respect to equatorial N-H and Fe═O) turned out to be the most stable. It was found that with the enrichment of the equatorial sulfur atom, reactivity increases successively, that is, we get the trend N₄ < N₃S₁ < N₂S₂ < N₁S₃ < S_4. Our investigations have also verified the available experimental results where it has been reported that N₂S₂ is more reactive than N₄ in their mixed conformation. In search of insights into this typical pattern of reactivity, the interplay of several factors has been recognized, such as the distortion energy which decreases for the transition states with the addition of sulfur; the spin density on the oxygen atom which increases implying that the radical character of abstractor increases on sulfur ligation; the energy of the electron acceptor orbital (the lowest unoccupied molecular orbital (σ_z²*)) which decreases continuously with the sulfur substitution; and the triplet-quintet oxidant energy gap which decreases consistently with S enrichment in the equatorial position. The computational predictions reported here, if further validated by experiments, will definitely encourage the synthesis of sulfur-ligated bio-inspired complexes instead of the ones constituting nitrogen exclusively.

Proton-Coupled Electron-Transfer Reactivity Controls Iron versus Sulfur Oxidation in Nonheme Iron-Thiolate Complexes

Reaction of the five-coordinate Fe^II(N₄S) complexes, [Fe^II(iPr₃TACN)(abt^X)](OTf) (abt = aminobenzenethiolate, X = H, CF₃), with a one-electron oxidant and an appropriate base leads to net H atom loss, generating new Fe^III(iminobenzenethiolate) complexes that were characterized by single-crystal X-ray diffraction (XRD), as well as UV-vis, EPR, and Mössbauer spectroscopies. The spectroscopic data indicate that the iminobenzenethiolate complexes have S = 3/2 ground states. In the absence of a base, oxidation of the Fe^II(abt) complexes leads to disulfide formation instead of oxidation at the metal center. Bracketing studies with separated proton-coupled electron-transfer (PCET) reagents show that the Fe^II(aminobenzenethiolate) and Fe^III(iminobenzenethiolate) forms are readily interconvertible by H⁺/e^- transfer and provide a measure of the bond dissociation free energy (BDFE) for the coordinated N-H bond between 64 and 69 kcal mol^-1. This work shows that coordination to the iron center causes a dramatic weakening of the N-H bond and that Fe- versus S-oxidation in a nonheme iron complex can be controlled by the protonation state of an ancillary amino donor.

Enhancing Tris(pyrazolyl)borate-Based Models of Cysteine/Cysteamine Dioxygenases through Steric Effects: Increased Reactivities, Full Product Characterization and Hints to Initial Superoxide Formation

The design of biomimetic model complexes for the cysteine dioxygenase (CDO) and cysteamine dioxygenase (ADO) is reported, where the 3-His coordination of the iron ion is simulated by three pyrazole donors of a trispyrazolyl borate ligand (Tp) and protected cysteine and cysteamine represent substrate ligands. It is found that the replacement of phenyl groups-attached at the 3-positions of the pyrazole units in a previous model-by mesityl residues has massive consequences, as the latter arrange to a more spacious reaction pocket. Thus, the reaction with O2 proceeds much faster and afterwards the first structural characterization of an iron(II) η2 -O,O-sulfinate product became possible. If one of the three Tp-mesityl groups is placed in the 5-position, an even larger reaction pocket results, which leads to yet faster rates and accumulation of a reaction intermediate at low temperatures, as shown by UV/Vis and Mössbauer spectroscopy. After comparison with the results of investigations on the cobalt analogues this intermediate is tentatively assigned to an iron(III) superoxide species.

A synthetic model of the nonheme iron-superoxo intermediate of cysteine dioxygenase

A nonheme Fe(ii) complex (1) that models substrate-bound cysteine dioxygenase (CDO) reacts with O2 at -80 °C to yield a purple intermediate (2). Analysis with spectroscopic and computational methods determined that 2 features a thiolate-ligated Fe(iii) center bound to a superoxide radical, mimicking the putative structure of a key CDO intermediate.

A Comparative Review on the Catalytic Mechanism of Nonheme Iron Hydroxylases and Halogenases

Enzymatic halogenation and haloperoxidation are unusual processes in biology; however, a range of halogenases and haloperoxidases exist that are able to transfer an aliphatic or aromatic C–H bond into C–Cl/C–Br. Haloperoxidases utilize hydrogen peroxide, and in a reaction with halides (Cl−/Br−), they react to form hypohalides (OCl−/OBr−) that subsequently react with substrate by halide transfer. There are three types of haloperoxidases, namely the iron-heme, nonheme vanadium, and flavin-dependent haloperoxidases that are reviewed here. In addition, there are the nonheme iron halogenases that show structural and functional similarity to the nonheme iron hydroxylases and form an iron(IV)-oxo active species from a reaction of molecular oxygen with α-ketoglutarate on an iron(II) center. They subsequently transfer a halide (Cl−/Br−) to an aliphatic C–H bond. We review the mechanism and function of nonheme iron halogenases and hydroxylases and show recent computational modelling studies of our group on the hectochlorin biosynthesis enzyme and prolyl-4-hydroxylase as examples of nonheme iron halogenases and hydroxylases. These studies have established the catalytic mechanism of these enzymes and show the importance of substrate and oxidant positioning on the stereo-, chemo- and regioselectivity of the reaction that takes place.

Understanding How Prolyl-4-hydroxylase Structure Steers a Ferryl Oxidant toward Scission of a Strong C-H Bond

Prolyl-4-hydroxylase (P4H) is a non-heme iron hydroxylase that regio- and stereospecifically hydroxylates proline residues in a peptide chain into R-4-hydroxyproline, which is essential for collagen cross-linking purposes in the human body. Surprisingly, in P4H, a strong aliphatic C-H bond is activated, while thermodynamically much weaker aliphatic C-H groups, that is, at the C³ and C⁵ positions, are untouched. Little is known on the origins of the high regio- and stereoselectivity of P4H and many non-heme and heme enzymes in general, and insight into this matter may be relevant to Biotechnology as well as Drug Development. The active site of the protein contains two aromatic residues (Tyr₁₄₀ and Trp₂₄₃) that we expected to be crucial for guiding the regioselectivity of the reaction. We performed a detailed quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics (MD) study on wild-type and mutant structures. The work shows that Trp₂₄₃ is involved in key protein loop-loop interactions that affect the shape and size of the substrate binding pocket and its mutation has major long-range effects. By contrast, the Tyr₁₄₀ residue is shown to guide the regio- and stereoselectivity by holding the substrate and ferryl oxidant in a specific orientation through hydrogen bonding and π-stacking interactions. Compelling evidence is found that the Tyr₁₄₀ residue is involved in expelling the product from the binding pocket after the reaction is complete. It is shown that mutations where the hydrogen bonding network that involves the Tyr₁₄₀ and Trp₂₄₃ residues is disrupted lead to major changes in folding of the protein and the size and shape of the substrate binding pocket. Specifically, the Trp₂₄₃ residue positions the amino acid side chains of Arg₁₆₁ and Glu₁₂₇ in specific orientations with substrate. As such, the P4H enzyme is a carefully designed protein with a subtle and rigid secondary structure that enables the binding of substrate, guides the regioselectivity, and expels product efficiently.

Sulfoxide Synthase versus Cysteine Dioxygenase Reactivity in a Nonheme Iron Enzyme

The sulfoxide synthase EgtB represents a unique family of nonheme iron enzymes that catalyze the formation of a C-S bond between N-α-trimethyl histidine and γ-glutamyl cysteine, which is the key step in the biosynthesis of ergothioneine, an important amino acid related to aging. A controversy has arisen regarding its catalytic mechanism related to the function of the active-site Tyr₃₇₇ residue. The biosynthesis of ergothioneine in EgtB shows structural similarities to cysteine dioxygenase which transfers two oxygen atoms to the thiolate group of cysteine. The question, therefore, is how do EgtB enzymes catalyze the C-S bond-formation reaction, while also preventing a dioxygenation of its cysteinate substrate? In this work we present a quantum mechanics/molecular mechanics study into the mechanism of sulfoxide synthase enzymes as compared to cysteine dioxygenase enzymes and present pathways for both reaction channels in EgtB. We show that EgtB contains a conserved tyrosine residue that reacts via proton-coupled electron transfer with the iron(III)-superoxo species and creates an iron(III)-hydroperoxo intermediate, thereby preventing the possible thiolate dioxygenation side reaction. The nucleophilic C-S bond-formation step happens subsequently concomitant to relay of the proton of the iron(II)-hydroperoxo back to Tyr₃₇₇. This is the rate-determining step in the reaction cycle and is followed by hydrogen-atom transfer from the CE1-H group of trimethyl histidine substrate to iron(II)-superoxo. In the final step, a quick and almost barrierless sulfoxidation leads to the sulfoxide product complexes. The work highlights a unique machinery and active-site setup of the enzyme that drives the sulfoxide synthase reaction.

Spectroscopic and computational studies of reversible O2 binding by a cobalt complex of relevance to cysteine dioxygenase

Spectroscopic and computational studies of reversible O₂ binding by a cobalt active-site mimic shed light on the catalytic mechanism of cysteine dioxygenases.

Synthesis, X-ray Structures, Electronic Properties, and O2/NO Reactivities of Thiol Dioxygenase Active-Site Models

Mononuclear non-heme iron complexes that serve as structural and functional mimics of the thiol dioxygenases (TDOs), cysteine dioxygenase (CDO) and cysteamine dioxygenase (ADO), have been prepared and characterized with crystallographic, spectroscopic, kinetic, and computational methods. The high-spin Fe(II) complexes feature the facially coordinating tris(4,5-diphenyl-1-methylimidazol-2-yl)phosphine (^Ph2TIP) ligand that replicates the three histidine (3His) triad of the TDO active sites. Further coordination with bidentate l-cysteine ethyl ester (CysOEt) or cysteamine (CysAm) anions yielded five-coordinate (5C) complexes that resemble the substrate-bound forms of CDO and ADO, respectively. Detailed electronic-structure descriptions of the [Fe(^Ph2TIP)(L_S,N)]BPh₄ complexes, where L_S,N = CysOEt (1) or CysAm (2), were generated through a combination of spectroscopic techniques [electronic absorption, magnetic circular dichroism (MCD)] and density functional theory (DFT). Complexes 1 and 2 decompose in the presence of O₂ to yield the corresponding sulfinic acid (RSO₂H) products, thereby emulating the reactivity of the TDO enzymes and related complexes. Rate constants and activation parameters for the dioxygenation reactions were measured and interpreted with the aid of DFT calculations for O₂-bound intermediates. Treatment of the TDO models with nitric oxide (NO)-a well-established surrogate of O₂-led to a mixture of high-spin and low-spin {FeNO}⁷ species at low temperature (-70 °C), as indicated by electron paramagnetic resonance (EPR) spectroscopy. At room temperature, these Fe/NO adducts convert to a common species with EPR and infrared (IR) features typical of cationic dinitrosyl iron complexes (DNICs). To complement these results, parallel spectroscopic, computational, and O₂/NO reactivity studies were carried out using previously reported TDO models that feature an anionic hydrotris(3-phenyl-5-methyl-pyrazolyl)borate (^Ph,MeTp^-) ligand. Though the O₂ reactivities of the ^Ph2TIP- and ^Ph,MeTp-based complexes are quite similar, the supporting ligand perturbs the energies of Fe 3d-based molecular orbitals and modulates Fe-S bond covalency, suggesting possible rationales for the presence of neutral 3His coordination in CDO and ADO.

Singlet versus Triplet Reactivity in an Mn(V)-Oxo Species: Testing Theoretical Predictions Against Experimental Evidence

Discerning the factors that control the reactivity of high-valent metal-oxo species is critical to both an understanding of metalloenzyme reactivity and related transition metal catalysts. Computational studies have suggested that an excited higher spin state in a number of metal-oxo species can provide a lower energy barrier for oxidation reactions, leading to the conclusion that this unobserved higher spin state complex should be considered as the active oxidant. However, testing these computational predictions by experiment is difficult and has rarely been accomplished. Herein, we describe a detailed computational study on the role of spin state in the reactivity of a high-valent manganese(V)-oxo complex with para-Z-substituted thioanisoles and utilize experimental evidence to distinguish between the theoretical results. The calculations show an unusual change in mechanism occurs for the dominant singlet spin state that correlates with the electron-donating property of the para-Z substituent, while this change is not observed on the triplet spin state. Minimum energy crossing point calculations predict small spin-orbit coupling constants making the spin state change from low spin to high spin unlikely. The trends in reactivity for the para-Z-substituted thioanisole derivatives provide an experimental measure for the spin state reactivity in manganese-oxo corrolazine complexes. Hence, the calculations show that the V-shaped Hammett plot is reproduced by the singlet surface but not by the triplet state trend. The substituent effect is explained with valence bond models, which confirm a change from an electrophilic to a nucleophilic mechanism through a change of substituent.

Utilizing the Trispyrazolyl Borate Ligand for the Mimicking of O2-Activating Mononuclear Nonheme Iron Enzymes

Mononuclear, O2-activating nonheme iron enzymes are a fascinating class of metalloproteines, capable of realizing the most different reactions, ranging from C-H activation, via O atom transfer to C-C bond cleavage, in the course of O2 activation. They can lead us the way to achieve similar reactions with comparable efficiency and selectivity in chemical laboratories, which would be highly desirable aiming at accessing value-added products or to achieve degradation of unwanted compounds. Hence, these enyzmes motivate attempts to construct artificial low-molecular weight analogues, mimicking structural or functional characteristics. Such models can, for instance, provide insights about which of the features inherent to an active site are essential and guarantee the enzyme function, and from this kind of information the minimal requirements for a biomimetic or bioinspired complex that may be applied in catalysis can be derived. On the other hand, they can contribute to an understanding of the enzyme functioning. In order to create such replicates, it is important to faithfully mimic the surroundings of the iron centers in their active sites. Most of them feature two histidine residues and one carboxylate donor, while a few exhibit a deceptively simple (His)3Fe active site. For the simulation of these, the trispyrazolyl borate ligand (Tp) particularly offers itself, as the facial arrangement of three pyrazole donors is reminiscent of the three histidine-derived imidazole donors. The focus of this Account will be on bioinorganic/biomimetic research from our laboratory utilizing Tp ligands to develop molecular models for (i) two representatives of the (His)3Fe-enzyme family, namely, the cysteine dioxygenase (CDO) and acetyl acetone dioxygenase (Dke1), (ii) a related but less well-explored variant of the CDO-the 2-aminoethanethiol dioxygenase-as well as (iii) the 2-His-1-carboxylate representative 1-aminocyclopropane-1-carboxylic acid oxidase (ACCO). The CDO catalyzes the dioxygenation of cysteine with O2 to give cysteine sulfinic acid, which could be mimicked at TpFe units in a realistic manner. Furthermore, the successful dioxygenation of 2-aminoethanethiol at the same complex metal fragments lends further support to the hypothesis that the active sites of CDO and the one of 2-aminoethanethiol dioxygenase, whose structure is unknown, are quite similar. Dke1 is capable of cleaving diketones and ketoesters to give the corresponding carboxylic acids and α-keto aldehydes, and Tp-based models have achieved comparable C-C bond cleavage reactions. The ACCO develops ethylene from ACC in the course of oxidation, and recently this has been achieved the first time for a TpFe model, too.

Alkyl Chain Growth on a Transition Metal Center: How Does Iron Compare to Ruthenium and Osmium?

Industrial Fischer-Tropsch processes involve the synthesis of hydrocarbons usually on metal surface catalysts. On the other hand, very few homogeneous catalysts are known to perform a Fischer-Tropsch style of reaction. In recent work, we established the catalytic properties of a diruthenium-platinum carbene complex, [(CpRu)₂(μ²-H)(μ²-NHCH₃)(μ³-C)PtCH₃(P(CH₃)₃)₂](CO)_n⁺ with n = 0, 2 and Cp = η⁵-C₅(CH₃)₅, and showed it to react efficiently by initial hydrogen atom transfer followed by methyl transfer to form an alkyl chain on the Ru-center. In particular, the catalytic efficiency was shown to increase after the addition of two CO molecules. As such, this system could be viewed as a potential homogeneous Fischer-Tropsch catalyst. Herein, we have engineered the catalytic center of the catalyst and investigated the reactivity of trimetal carbene complexes of the same type using iron, ruthenium and osmium at the central metal scaffold. The work shows that the reactivity should increase from diosmium to diruthenium to diiron; however, a non-linear trend is observed due to multiple factors contributing to the individual barrier heights. We identified all individual components of these reaction steps in detail and established the difference in reactivity of the various complexes.

A structural and functional model for the 1-aminocyclopropane-1-carboxylic acid oxidase

The hitherto most realistic low-molecular-weight analogue for the 1-aminocyclopropane-1-carboxylic acid oxidase (ACCO) is reported. The ACCOs 2-His-1-carboxylate iron(II) active site was mimicked by a TpFe moiety, to which the natural substrate ACC could be bound. The resulting complex [Tp(Me,Ph) FeACC] (1), according to X-ray diffraction analysis performed for the nickel analogue, represents an excellent structural model, featuring ACC coordinated in a bidentate fashion-as proposed for the enzymatic substrate complex-as well as a vacant coordination site that forms the basis for the first successful replication also of the ACCO function: 1 is the first known ACC complex that reacts with O2 to produce ethylene. As a FeOOH species had been suggested as intermediate in the catalytic cycle, H2 O2 was tested as the oxidant, too, and indeed evolution of ethylene proceeded even more rapidly to give 65 % yield.