Search

The cytochrome P450 homolog, TxtE, efficiently catalyzes the direct and regioselective aromatic nitration of the indolyl moiety of L-tryptophan to 4-nitro-L-tryptophan, using nitric oxide (NO) and dioxygen (O 2 ) as co-substrates. Pathways for such direct and selective nitration of heteroaromatic motifs present platforms for engineering new nitration biocatalysts for pharmacologically beneficial targets, among a medley of other pivotal industrial applications. Precise mechanistic details concerning this pathway are only weakly understood, albeit a heme iron(III)-peroxynitrite active species has been postulated. To shed light on this unique reaction landscape, we investigated the indole nitration pathway of a series of biomimetic ferric heme superoxide mimics, [(Por)Fe III (O 2 - ⋅)], in the presence of NO. Therein, our model systems gave rise to three distinct nitroindole products, including 4-nitroindole, the product analogous to that obtained with TxtE. Moreover, 15 N and 18 O isotope labeling studies, along with meticulously designed control experiments lend credence to a heme peroxynitrite active nitrating agent, drawing close similarities to the tryptophan nitration mechanism of TxtE. All organic and inorganic reaction components have been fully characterized using spectroscopic methods. Theoretical investigation into several mechanistic possibilities deem a unique indolyl radical based reaction pathway as the most energetically favorable, products of which, are in excellent agreement with experimental findings., (© 2024 Wiley-VCH GmbH.)

The nonheme iron dioxygenase capreomycin C (CmnC) hydroxylates a free L-arginine amino acid regio- and stereospecifically at the C 3 -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 C 3 -H, C 4 -H and C 5 -H hydroxylation and C 3 -C 4 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 C 3 -H activation, whereby the polar groups in the substrate binding pocket induce an electric field effect that guides the selectivity., (© 2024 The Author(s). Chemistry - A European Journal published by Wiley-VCH GmbH.)

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., (© 2024 The Author(s). Chemistry - A European Journal published by Wiley-VCH GmbH.)

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.

The ability to introduce noncanonical amino acids as axial ligands in heme enzymes has provided a powerful experimental tool for studying the structure and reactivity of their Fe IV =O ("ferryl") intermediates. Here, we show that a similar approach can be used to perturb the conserved Fe coordination environment of 2-oxoglutarate (2OG) dependent oxygenases, a versatile class of enzymes that employ highly-reactive ferryl intermediates to mediate challenging C-H functionalizations. Replacement of one of the cis-disposed histidine ligands in the oxygenase VioC with a less electron donating N δ -methyl-histidine (MeHis) preserves both catalytic function and reaction selectivity. Significantly, the key ferryl intermediate responsible for C-H activation can be accumulated in both the wildtype and the modified protein. In contrast to heme enzymes, where metal-oxo reactivity is extremely sensitive to the nature of the proximal ligand, the rates of C-H activation and the observed large kinetic isotope effects are only minimally affected by axial ligand replacement in VioC. This study showcases a powerful tool for modulating the coordination sphere of nonheme iron enzymes that will enhance our understanding of the factors governing their divergent activities., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

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.

The nature of the axial ligand in high-valent iron-oxo heme enzyme intermediates and related synthetic catalysts is a critical structural element for controlling proton-coupled electron-transfer (PCET) reactivity of these species. Herein, we describe the generation and characterization of three new 6-coordinate, iron(IV)-oxo porphyrinoid-π-cation-radical complexes and report their PCET reactivity together with a previously published 5-coordinate analogue, Fe IV (O)(TBP 8 Cz +• ) (TBP 8 Cz = octakis( p - tert -butylphenyl)corrolazinato 3- ) ( 2 ) (Cho, K. A high-valent iron-oxo corrolazine activates C-H bonds via hydrogen-atom transfer. J. Am. Chem. Soc. 2012, 134, 7392-7399). The new complexes Fe IV (O)(TBP 8 Cz +• )(L) (L = 1-methyl imidazole (1-MeIm) ( 4a ), 4-dimethylaminopyridine (DMAP) ( 4b ), cyanide (CN - )( 4c )) can be generated from either oxidation of the ferric precursors or by addition of L to the Compound-I (Cpd-I) analogue at low temperatures. These complexes were characterized by UV-vis, electron paramagnetic resonance (EPR), and Mössbauer spectroscopies, and cryospray ionization mass spectrometry (CSI-MS). Kinetic studies using 4-OMe-TEMPOH as a test substrate indicate that coordination of a sixth axial ligand dramatically lowers the PCET reactivity of the Cpd-I analogue (rates up to 7000 times slower ). Extensive density functional theory (DFT) calculations together with the experimental data show that the trend in reactivity with the axial ligands does not correlate with the thermodynamic driving force for these reactions or the calculated strengths of the O-H bonds being formed in the Fe IV (O-H) products, pointing to non-Bell-Evans-Polanyi behavior. However, the PCET reactivity does follow a trend with the bracketed reduction potential of Cpd-I analogues and calculated electron affinities. The combined data suggest a concerted mechanism (a concerted proton electron transfer (CPET)) and an asynchronous movement of the electron/proton pair in the transition state.

Lysine dioxygenase (KDO) is an important enzyme in human physiology involved in bioprocesses that trigger collagen cross-linking and blood pressure control. There are several KDOs in nature; however, little is known about the factors that govern the regio- and stereoselectivity of these enzymes. To understand how KDOs can selectively hydroxylate their substrate, we did a comprehensive computational study into the mechanisms and features of 4-lysine dioxygenase. In particular, we selected a snapshot from the MD simulation on KDO5 and created large QM cluster models ( A , B , and C ) containing 297, 312, and 407 atoms, respectively. The largest model predicts regioselectivity that matches experimental observation with rate-determining hydrogen atom abstraction from the C 4 -H position, followed by fast OH rebound to form 4-hydroxylysine products. The calculations show that in model C , the dipole moment is positioned along the C 4 -H bond of the substrate and, therefore, the electrostatic and electric field perturbations of the protein assist the enzyme in creating C 4 -H hydroxylation selectivity. Furthermore, an active site Tyr 233 residue is identified that reacts through proton-coupled electron transfer akin to the axial Trp residue in cytochrome c peroxidase. Thus, upon formation of the iron(IV)-oxo species in the catalytic cycle, the Tyr 233 phenol loses a proton to the nearby Asp 179 residue, while at the same time, an electron is transferred to the iron to create an iron(III)-oxo active species. This charged tyrosyl residue directs the dipole moment along the C 4 -H bond of the substrate and guides the selectivity to the C 4 -hydroxylation of the substrate.

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., (© 2024 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

The nonheme iron dioxygenase deoxypodophyllotoxin synthase performs an oxidative ring-closure reaction as part of natural product synthesis in plants. How the enzyme enables the oxidative ring-closure reaction of (-)-yatein and avoids substrate hydroxylation remains unknown. To gain insight into the reaction mechanism and understand the details of the pathways leading to products and by-products we performed a comprehensive computational study. The work shows that substrate is bound tightly into the substrate binding pocket with the C 7' -H bond closest to the iron(IV)-oxo species. The reaction proceeds through a radical mechanism starting with hydrogen atom abstraction from the C 7' -H position followed by ring-closure and a final hydrogen transfer to form iron(II)-water and deoxypodophyllotoxin. Alternative mechanisms including substrate hydroxylation and an electron transfer pathway were explored but found to be higher in energy. The mechanism is guided by electrostatic perturbations of charged residues in the second-coordination sphere that prevent alternative pathways., (© 2024 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

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 N 4 S-ligated iron(IV)-oxo complex and compared the obtained results with an analogous N 5 -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 N 5 -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.

Transforming CO 2 into valuable materials is an important reaction in catalysis, especially because CO 2 concentrations in the atmosphere have been growing steadily due to extensive fossil fuel usage. From an environmental perspective, reduction of CO 2 to valuable materials should be catalyzed by an environmentally benign catalyst and avoid the use of heavy transition-metal ions. In this work, we present a computational study into a novel iron(I) porphyrin catalyst for CO 2 reduction, namely, with a tetraphenylporphyrin ligand and analogues. In particular, we investigated iron(I) tetraphenylporphyrin with one of the meso -phenyl groups substituted with o -urea, p -urea, or o -2-amide groups. These substituents can provide hydrogen-bonding interactions in the second coordination sphere with bound ligands and assist with proton relay. Furthermore, our studies investigated bicarbonate and phenol as stabilizers and proton donors in the reaction mechanism. Potential energy landscapes for double protonation of iron(I) porphyrinate with bound CO 2 are reported. The work shows that the bicarbonate bridges the urea/amide groups to the CO 2 and iron center and provides a tight bonding pattern with strong hydrogen-bonding interactions that facilitates easy proton delivery and reduction of CO 2 . Specifically, bicarbonate provides a low-energy proton shuttle mechanism to form CO and water efficiently. Furthermore, the o -urea group locks bicarbonate and CO 2 in a tight orientation and helps with ideal proton transfer, while there is more mobility and lesser stability with an o -amide group in that position instead. Our calculations show that the o -urea group leads to reduction in proton-transfer barriers, in line with experimental observation. We then applied electric-field-effect calculations to estimate the environmental effects on the two proton-transfer steps in the reaction. These calculations describe the perturbations that enhance the driving forces for the proton-transfer steps and have been used to make predictions about how the catalysts can be further engineered for more enhanced CO 2 reduction processes.

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 C 3 -hydroxylation, C 4 -hydroxylation, C 5 -hydroxylation, or C 4 -C 5 -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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declare that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Ali and de Visser.)

CO 2 utilization is an important process in the chemical industry with great environmental power. In this work we show how CO 2 and H 2 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 H 2 followed by sequential hydride and proton transfer processes. The H 2 splitting process is an endergonic process and hence high pressures will be needed to overcome this step and trigger the hydrogenation reaction. Moreover, H 2 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 H 2 to substrate as an important step in the reaction mechanism., (© 2023 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

Density functional calculations establish a novel mechanism of aromatic defluorination by P450 Compound I. This is achieved via either an initial epoxide intermediate or through a 1,2-fluorine shift in an electrophilic intermediate, which highlights that the P450s can defluorinate fluoroarenes. However, in the absence of a proton donor a strong Fe-F bond can be obtained as shown from the calculations., (© 2023 Wiley-VCH GmbH.)

[MFe 3 S 4 ] cubanes have for some time been of interest for their ability to mimic the electronic and geometric structure of the active site of nitrogenase, the enzyme responsible for fixing N 2 to NH 3 . Nitrogenase naturally occurs in three forms, with the major difference being that the metal ion present in the cofactor active site is either molybdenum (FeMoco), vanadium (FeVco), or iron. The molybdenum and vanadium versions of these cofactors are more closely studied, owing to their larger abundance and rate of catalysis. In this study, we compare free energy profiles and electronic properties of the Mo/V cubanes at various stages during the reduction of N 2 H 4 to NH 3 . Our findings highlight the differences in how the complexes facilitate the reaction, in particular, vanadium's comparatively weaker ability to interact with the Fe/S network and stabilize reducing electrons prior to N-N bond cleavage, which may have implications when considering the lower efficiency of the vanadium-dependent nitrogenase.

Cytochrome P450 enzymes in the human body play a pivotal role in both the biosynthesis and the degradation of the hormone melatonin. Melatonin plays a key role in circadian rhythms in the body, but its concentration is also linked to mood fluctuations as well as emotional well-being. In the present study, we present a computational analysis of the binding and activation of melatonin by various P450 isozymes that are known to yield different products and product distributions. In particular, the P450 isozymes 1A1, 1A2, and 1B1 generally react with melatonin to provide dominant aromatic hydroxylation at the C 6 -position, whereas the P450 2C19 isozyme mostly provides O -demethylation products. To gain insight into the origin of these product distributions of the P450 isozymes, we performed a comprehensive computational study of P450 2C19 isozymes and compared our work with previous studies on alternative isozymes. The work covers molecular mechanics, molecular dynamics and quantum mechanics approaches. Our work highlights major differences in the size and shape of the substrate binding pocket amongst the different P450 isozymes. Consequently, substrate binding and positioning in the active site varies substantially within the P450 isozymes. Thus, in P450 2C19, the substrate is oriented with its methoxy group pointing towards the heme, and therefore reacts favorably through hydrogen atom abstraction, leading to the production of O -demethylation products. On the other hand, the substrate-binding pockets in P450 1A1, 1A2, and 1B1 are tighter, direct the methoxy group away from the heme, and consequently activate an alternative site and lead to aromatic hydroxylation instead.

Nitrogenase is a fascinating enzyme in biology that reduces dinitrogen from air to ammonia through stepwise reduction and protonation. Despite it being studied in detail by experimental and computational groups, there are still many unknown factors in the catalytic cycle of nitrogenase, especially related to the addition of protons and electrons and their order. A recent biomimetic study characterized a potential dinitrogen-bridged diiron cluster as a synthetic model of nitrogenase. Using strong acid and reductants, the dinitrogen was converted into ammonia molecules, but details of the mechanism remains unknown. In particular, it was unclear from the experimental studies whether the proton and electron transfer steps are sequential or alternating. Moreover, the work failed to establish what the function of the diiron core is and whether it split into mononuclear iron fragments during the reaction. To understand the structure and reactivity of the biomimetic dinitrogen-bridged diiron complex [(P 2 P'Ph FeH) 2 (μ-N 2 )] with triphenylphosphine ligands, we performed a density functional theory study. Our computational methods were validated against experimental crystal structure coordinates, Mössbauer parameters, and vibrational frequencies and show excellent agreement. Subsequently, we investigated the alternating and consecutive addition of electrons and protons to the system. The calculations identify a number of possible reaction channels, namely, same-site protonation, alternating protonation, and complex dissociation into mononuclear iron centers. The calculations show that the overall mechanism is not a pure sequential set of electron and proton transfers but a mixture of alternating and consecutive steps. In particular, the first reaction steps will start with double proton transfer followed by an electron transfer, while thereafter, there is another proton transfer and a second electron transfer to give a complex whereby ammonia can split off with a low energetic barrier. The second channel starts with alternating protonation of the two nitrogen atoms, whereafter the initial double proton transfer, electrons and protons are added sequentially to form a hydrazine-bound complex. The latter split off ammonia spontaneously after further protonation. The various reaction channels are analyzed with valence bond and orbital diagrams. We anticipate the nitrogenase enzyme to operate with mixed alternating and consecutive protonation and electron transfer steps.

Heterogenisation of homogeneous catalysts onto solid supports represents a potential strategy to make the homogeneous catalytic function recyclable and reuseable. Yet, it is usually the case that immobilised catalysts have much lower catalytic activity than their homogeneous counterpart. In addition, the presence of a solid interface introduces a higher degree of complexity by modulating solid/fluid interactions, which can often influence adsorption properties of solvents and reactive species and, ultimately, catalytic activity. In this work, the influence of support and solvent in the H-transfer reduction of propionaldehyde over Al(O i Pr) 3 -SiO 2 , Al(O i Pr) 3 -TiO 2 and Al(O i Pr) 3 -Al 2 O 3 heterogenised catalysts has been studied. Reaction studies are coupled with both NMR relaxation measurements as well as molecular dynamics (MD) simulations in order to unravel surface and solvation effects during the reaction. The results show that, whilst the choice of the support does not influence significantly catalytic activity, reactions carried out in solvents with high affinity for the catalyst surface, or able to hinder access to active sites due to solvation effects, have a lower activity. MD calculations provide key insights into bulk solvation effects involved in such reactions, which are thought to play an important role in determining the catalytic behaviour. The activity of the heterogenised catalysts was found to be comparable with that of the homogeneous Al(O i Pr) 3 catalysts for all supports used, showing that for the type of reaction studied immobilisation of the homogeneous catalyst onto solid supports is a viable, robust and effective strategy.

High-valent metal-oxo species play critical roles in enzymatic catalysis yet their properties are still poorly understood. In this work we report a combined experimental and computational study into biomimetic iron(IV)-oxo and iron(III)-oxo complexes with tight second-coordination sphere environments that restrict substrate access. The work shows that the second-coordination sphere slows the hydrogen atom abstraction step from toluene dramatically and the kinetics is zeroth order in substrate. However, the iron(II)-hydroxo that is formed has a low reduction potential and hence cannot do OH rebound favorably. The tolyl radical in solution then reacts further with alternative reaction partners. By contrast, the iron(IV)-oxo species reacts predominantly through OH rebound to form alcohol products. Our studies show that the oxidation state of the metal influences reactivities and selectivities with substrate dramatically and that enzymes will likely need an iron(IV) center to catalyze C-H hydroxylation reactions., (© 2023 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

High-valent iron(IV)-oxo intermediates are versatile oxidants in the biotransformation of various substrates by metalloenzymes and catalyze essential reactions for human health as well as in the biodegradation of toxic organic pollutants in the environment. Herein, we report a biomimetic system that efficiently reacts with halophenols through defluorination reactions and characterize various short-lived intermediates along the reaction mechanism. We study the reactivity pattern of a nonheme iron(IV)-oxo species with a series of trihalophenols (X=F, Cl, Br). A combined experimental and computational study reveals that the oxidative dehalogenation of 2,4,6-trifluorophenol is initiated with an H-atom abstraction from the phenolic group by the iron(IV)-oxo species resulting in the formation of a phenolate radical and an iron(III)-hydroxo species. This iron(III)-hydroxo species forms an adduct with the oxidized substrate with λ max at 558 nm which subsequently decays to give quinones as products., (© 2023 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

Caffeine is a natural compound found in plant seeds that after consumption by humans effects the central nervous system as well as the cardiovascular system. In general, the cytochrome P450 enzymes in the liver are involved in the biodegradation of caffeine, which gives paraxanthine, theobromine and theophylline products. There has been debate for many years why multiple products are obtained and how their distributions are determined. To this end we performed a high-level computational study using a combination of molecular dynamics and quantum mechanical approaches. A series of quantum chemical cluster models on the mechanism of caffeine activation by P450 model complexes give hydrogen atom abstraction barriers that predicts the correct ordering and statistical distribution of products. Our studies highlight that second-coordination sphere effects and thermochemical properties of the substrate determine the product distributions., (© 2023 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

Vanadium haloperoxidases (VHPOs) are unique enzymes in biology that catalyze a challenging halogen transfer reaction and convert a strong aromatic C-H bond into C-X (X = Cl, Br, I) with the use of a vanadium cofactor and H 2 O 2 . The VHPO catalytic cycle starts with the conversion of hydrogen peroxide and halide (X = Cl, Br, I) into hypohalide on the vanadate cofactor, and the hypohalide subsequently reacts with a substrate. However, it is unclear whether the hypohalide is released from the enzyme or otherwise trapped within the enzyme structure for the halogenation of organic substrates. A substrate-binding pocket has never been identified for the VHPO enzyme, which questions the role of the protein in the overall reaction mechanism. Probing its role in the halogenation of small molecules will enable further engineering of the enzyme and expand its substrate scope and selectivity further for use in biotechnological applications as an environmentally benign alternative to current organic chemistry synthesis. Using a combined experimental and computational approach, we elucidate the role of the vanadium haloperoxidase protein in substrate halogenation. Activity studies show that binding of the substrate to the enzyme is essential for the reaction of the hypohalide with substrate. Stopped-flow measurements demonstrate that the rate-determining step is not dependent on substrate binding but partially on hypohalide formation. Using a combination of molecular mechanics (MM) and molecular dynamics (MD) simulations, the substrate binding area in the protein is identified and even though the selected substrates (methylphenylindole and 2-phenylindole) have limited hydrogen-bonding abilities, they are found to bind relatively strongly and remain stable in a binding tunnel. A subsequent analysis of the MD snapshots characterizes two small tunnels leading from the vanadate active site to the surface that could fit small molecules such as hypohalide, halide, and hydrogen peroxide. Density functional theory studies using electric field effects show that a polarized environment in a specific direction can substantially lower barriers for halogen transfer. A further analysis of the protein structure indeed shows a large dipole orientation in the substrate-binding pocket that could enable halogen transfer through an applied local electric field. These findings highlight the importance of the enzyme in catalyzing substrate halogenation by providing an optimal environment to lower the energy barrier for this challenging aromatic halide insertion reaction., Competing Interests: The authors declare no competing financial interest., (© 2023 The Authors. Published by American Chemical Society.)

The catalytic functions of metalloenzymes are often strongly correlated with metal elements in the active sites. However, dioxygen-activating nonheme quercetin dioxygenases (QueD) are found with various first-row transition-metal ions when metal swapping inactivates their innate catalytic activity. To unveil the molecular basis of this seemingly promiscuous yet metal-specific enzyme, we transformed manganese-dependent QueD into a nickel-dependent enzyme by sequence- and structure-based directed evolution. Although the net effect of acquired mutations was primarily to rearrange hydrophobic residues in the active site pocket, biochemical, kinetic, X-ray crystallographic, spectroscopic, and computational studies suggest that these modifications in the secondary coordination spheres can adjust the electronic structure of the enzyme-substrate complex to counteract the effects induced by the metal substitution. These results explicitly demonstrate that such noncovalent interactions encrypt metal specificity in a finely modulated manner, revealing the underestimated chemical power of the hydrophobic sequence network in enzyme catalysis.

Cytochrome P450 enzymes are versatile enzymes found in most biosystems that catalyze mono-oxygenation reactions as a means of biosynthesis and biodegradation steps. In the liver, they metabolize xenobiotics, but there are a range of isozymes with differences in three-dimensional structure and protein chain. Consequently, the various P450 isozymes react with substrates differently and give varying product distributions. To understand how melatonin is activated by the P450s in the liver, we did a thorough molecular dynamics and quantum mechanics study on cytochrome P450 1A2 activation of melatonin forming 6-hydroxymelatonin and N -acetylserotonin products through aromatic hydroxylation and O -demethylation pathways, respectively. We started from crystal structure coordinates and docked substrate into the model, and obtained ten strong binding conformations with the substrate in the active site. Subsequently, for each of the ten substrate orientations, long (up to 1 μs) molecular dynamics simulations were run. We then analyzed the orientations of the substrate with respect to the heme for all snapshots. Interestingly, the shortest distance does not correspond to the group that is expected to be activated. However, the substrate positioning gives insight into the protein residues it interacts with. Thereafter, quantum chemical cluster models were created and the substrate hydroxylation pathways calculated with density functional theory. These relative barrier heights confirm the experimental product distributions and highlight why certain products are obtained. We make a detailed comparison with previous results on CYP1A1 and identify their reactivity differences with melatonin.

Bisphenol A (BPA, 2,2-bis-(4-hydroxyphenyl)propane) is used as a precursor in the synthesis of polycarbonate and epoxy plastics; however, its availability in the environment is causing toxicity as an endocrine-disrupting chemical. Metabolism of BPA and their analogues (substitutes) is generally performed by liver cytochrome P450 enzymes and often leads to a mixture of products, and some of those are toxic. To understand the product distributions of P450 activation of BPA, we have performed a computational study into the mechanisms and reactivities using large model structures of a human P450 isozyme (P450 2C9) with BPA bound. Density functional theory (DFT) calculations on mechanisms of BPA activation by a P450 compound I model were investigated, leading to a number of possible products. The substrate-binding pocket is tight, and as a consequence, aliphatic hydroxylation is not feasible as the methyl substituents of BPA cannot reach compound I well due to constraints of the substrate-binding pocket. Instead, we find low-energy pathways that are initiated with phenol hydrogen atom abstraction followed by OH rebound to the phenolic ortho - or para -position. The barriers of para -rebound are well lower in energy than those for ortho -rebound, and consequently, our P450 2C9 model predicts dominant hydroxycumyl alcohol products. The reactions proceed through two-state reactivity on competing doublet and quartet spin state surfaces. The calculations show fast and efficient substrate activation on a doublet spin state surface with a rate-determining electrophilic addition step, while the quartet spin state surface has multiple high-energy barriers that can also lead to various side products including C 4 -aromatic hydroxylation. This work shows that product formation is more feasible on the low spin state, while the physicochemical properties of the substrate govern barrier heights of the rate-determining step of the reaction. Finally, the importance of the second-coordination sphere is highlighted that determines the product distributions and guides the bifurcation pathways.

Bacterial cytochromes P450 BM3 (CYP450 BM3) catalyze reactions of industrial importance. Despite many successful biotransformations, robust (re)design for novel applications remains challenging. Rational design and evolutionary approaches are not always successful, highlighting a lack of complete understanding of the mechanisms of electron transfer (ET) modulations. Thus, the full potential of CYP450 reactions remains under-exploited. In this work, we report the first molecular dynamics (MD)-based explicit prediction of BM3 ET parameters (reorganization energies; λ and ET free energies; Δ G °), and log ET rates (log k ET ) using the Marcus theory. Overall, the calculated ET rates for the BM3 wild-type (WT), mutants (F393 and L86), ligand-bound state, and ion concentrations agree well with experimental data. In ligand-free (LF) BM3, mutations modulate k ET via ET Δ G °. Simulations show that the experimental ET rate enhancement is due to increased driving force (more negative Δ G °) upon ligation. This increase is related to the protein reorganization required to accommodate the ligand in the binding pocket rather than binding interactions with the ligand. Our methodology (CYPWare 1.0) automates all the stages of the MD simulation step-up, energy calculations, and estimation of ET parameters. CYPWare 1.0 and this work thus represent an important advancement in the CYP450 ET rate predictions, which has the potential to guide the redesign of ET enzymes. This program and a Web tool are available on GitHub for academic research.

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 6 -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 6 -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.

Many drug molecules contain isonitrile substituents; however, synthesizing these compounds remains challenging in organic chemistry. The isonitrile synthesizing enzyme ScoE utilizes a substrate with the γ-Gly substituent, and using two molecules of dioxygen and α-ketoglutarate converts it to an isonitrile group through an oxidative decarboxylation reaction. To explore its substrate scope and whether this process could be used for the biosynthesis of isonitrile-containing drug molecules, we performed a predictive computational study. We started with the recent crystal structure coordinates of ScoE, removed the substrate and inserted two potential precursor molecules of the drug molecules axisonitrile-1 and xanthocillin into the structure, whereby both molecules have their isonitrile groups replaced by γ-Gly. Both substrates fit into the substrate binding pocket of the enzyme well and position them in the correct orientation for catalysis on the iron center. Based on a molecular dynamics simulation, we created a quantum chemical cluster model of the enzyme active site with γ-Gly-substituted axisonitrile-1 and studied the oxidative decarboxylation reaction to form axisonitrile-1 products. The calculations give similar barriers to wildtype substrate for either the initial C-H or N-H hydrogen atom abstraction, which leads to a radical intermediate and form desaturated reactants. We then took the desaturated substrate and created another iron(IV)-oxo model complex to study the subsequent hydrogen atom abstraction and decarboxylation and found this to be feasible as well although we predict to see by-products for hydroxylation in the second cycle. Nevertheless, we believe that the ScoE enzyme can be utilized for the biosynthesis of isonitrile substituents in substrates with γ-Gly components as an environmentally benign alternative to organic chemistry approaches for the synthesis of isonitrile groups. We hope that experimental studies will be able to confirm our hypothesis.

Ferritins are ubiquitous diiron enzymes involved in iron(II) detoxification and oxidative stress responses and can act as metabolic iron stores. The overall reaction mechanisms of ferritin enzymes are still unclear, particularly concerning the role of the conserved, near catalytic center Tyr residue. Thus, we carried out a computational study of a ferritin using a large cluster model of well over 300 atoms including its first- and second-coordination sphere. The calculations reveal important insight into the structure and reactivity of ferritins. Specifically, the active site Tyr residue delivers a proton and electron in the catalytic cycle prior to iron(II) oxidation. In addition, the calculations highlight a likely cation binding site at Asp 65 , which through long-range electrostatic interactions, influences the electronic configuration and charge distributions of the metal center. The results are consistent with experimental observations but reveal novel detail of early mechanistic steps that lead to an unusual mixed-valent iron(III)-iron(II) center., (© 2022 The Authors. ChemBioChem published by Wiley-VCH GmbH.)

Nonheme iron halogenases are unique enzymes in nature that selectively activate an aliphatic C-H bond of a substrate to convert it into C-X (X = Cl/Br, but not F/I). It is proposed that they generate an Fe III (OH)(X) intermediate in their catalytic cycle. The analogous Fe III (OH) intermediate in nonheme iron hydroxylases transfers OH • to give alcohol product, whereas the halogenases transfer X • to the carbon radical substrate. There remains significant debate regarding what factors control their remarkable selectivity of the halogenases. The reactivity of the complexes Fe III (BNPA Ph2 O)(OH)(X) (X = Cl, Br) with a secondary carbon radical (R • ) is described. It is found that X • transfer occurs with a secondary carbon radical, as opposed to OH • transfer with tertiary radicals. Comprehensive computational studies involving density functional theory were carried out to examine the possible origins of this selectivity. The calculations reproduce the experimental findings, which indicate that halogen transfer is not observed for the tertiary radicals because of a nonproductive equilibrium that results from the endergonic nature of these reactions, despite a potentially lower reaction barrier for the halogenation pathway. In contrast, halogen transfer is favored for secondary carbon radicals, for which the halogenated product complex is thermodynamically more stable than the reactant complex. These results are rationalized by considering the relative strengths of the C-X bonds that are formed for tertiary versus secondary carbon centers. The computational analysis also shows that the reaction barrier for halogen transfer is significantly dependent on secondary coordination sphere effects, including steric and H-bonding interactions.

Steady-state kinetics of cytochrome-c dependent denitrifying NO reductases (cNORs) show evidence of substrate inhibition at NO concentrations higher than 10 μM, but the mechanism of inhibition remains unclear. Here, we present low-temperature FTIR photolysis experiments carried out on the NO complex formed by addition of NO to the oxidized cNORs. A differential signal at 1261 cm -1 that downshifts with 15 NO and 15 N 18 O is assigned to a ν(NO 2 ) from a bridging diiron-nitrito complex at the heme-nonheme diron site. Theoretical calculations reproduces observed frequencies and isotope shifts. Our experimental results confirm a prior theoretical study by Blomberg and Siegbahn [Blomberg, M. R., and Siegbahn, P. E. M. Biochemistry 2012, 51, 5173-5186] that proposed substrate inhibition through a radical combination reaction between the diferric μ-oxo group and an NO molecule to form a heme Fe(III)-nitrito-Fe B (II) inhibitory complex. Stopped-flow experiments suggest that substrate inhibition also occurs after a half-reduction cycle, i.e. when fully-reduced cNOR reduces two NO molecules at the heme-nonheme diferrous active site cluster to produce one N 2 O molecule and the diferric cluster. These results support catalytic mechanisms that proceed through isomerization of a diferric-hyponitrite transient complex to produce a bridging diferric μ-oxo group and N 2 O without protonation of the putative hyponitrite intermediate., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Mononuclear high-valent iron(IV)-oxo intermediates are excellent oxidants towards oxygenation reactions by heme and nonheme metalloenzymes and their model systems. One of the most important functions of these intermediates in nature is to detoxify various environmental pollutants. Organic substrates, such as halogenated phenols, are known to be water pollutants which can be degraded to their less hazardous forms through an oxidation reaction by iron(IV)-oxo complexes. Metalloproteins in nature utilize various types of second-coordination sphere interactions to anchor the substrate in the vicinity of the active site. This concept of substrate-binding is well-known for natural enzymes, but is elusive for the relevant biomimetic model systems. Herein, we report the oxidative reactivity patterns of an iron(IV)-oxo intermediate, [Fe IV (O)(2PyN2Q)] 2+ , (2PyN2Q = 1,1-di(pyridin-2yl)- N , N -bis(quinolin-2-ylmethyl)methanamine) with a series of mono-, di- and tri-halophenols. A detailed experimental study shows that the dehalogenation reactions of the halophenols by such iron(IV)-oxo intermediates proceed via an initial hydrogen atom abstraction from the phenolic O-H group. Furthermore, based on the size and nucleophilicity of the halophenol, an intermediate substrate-bound species forms that is a phenolate adduct to the ferric species, which thereafter leads to the formation of the corresponding products.

Taurine/α-ketoglutarate dioxygenase is an important enzyme that takes part in the cysteine catabolism process in the human body and selectively hydroxylates taurine at the C 1 -position. Recent computational studies showed that in the gas-phase the C 2 -H bond of taurine is substantially weaker than the C 1 -H bond, yet no evidence exists of 2-hydroxytaurine products. To this end, a detailed computational study on the selectivity patterns in TauD was performed. The calculations show that the second-coordination sphere and the protonation states of residues play a major role in guiding the enzyme to the right selectivity. Specifically, a single proton on an active site histidine residue can change the regioselectivity of the reaction through its electrostatic perturbations in the active site and effectively changes the C 1 -H and C 2 -H bond strengths of taurine. This is further emphasized by many polar and hydrogen bonding interactions of the protein cage in TauD with the substrate and the oxidant that weaken the pro-R C 1 -H bond and triggers a chemoselective reaction process. The large cluster models reproduce the experimental free energy of activation excellently., (© 2021 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

Aryloxyalkanoate dioxygenases are unique herbicide biodegrading nonheme iron enzymes found in plants and hence, from environmental and agricultural point of view they are important and valuable. However, they often are substrate specific and little is known on the details of the mechanism and the substrate scope. To this end, we created enzyme models and calculate the mechanism for 2,4-dichlorophenoxyacetic acid biodegradation and 2-methyl substituted analogues by density functional theory. The work shows that the substrate binding is tight and positions the aliphatic group close to the metal center to enable a chemoselective reaction mechanism to form the C 2 -hydroxy products, whereas the aromatic hydroxylation barriers are well higher in energy. Subsequently, we investigated the metabolism of R- and S-methyl substituted inhibitors and show that these do not react as efficiently as 2,4-dichlorophenoxyacetic acid substrate due to stereochemical clashes in the active site and particularly for the R-isomer give high rebound barriers., (© 2021 Wiley-VCH GmbH.)

Non-heme iron dioxygenases catalyze vital processes for human health related to the biosynthesis of essential products and the biodegradation of toxic metabolites. Often the natural product biosyntheses by these non-heme iron dioxygenases is highly regio- and chemoselective, which are commonly assigned to tight substrate-binding and positioning. However, recent high-level computational modeling has shown that substrate-binding and positioning is only part of the story and long-range electrostatic interactions can play a major additional role.In this Account, we review and summarize computational viewpoints on the high regio- and chemoselectivity of α-ketoglutarate-dependent non-heme iron dioxygenases and how external perturbations affect the catalysis. In particular, studies from our groups have shown that often a regioselectivity in enzymes can be accomplished by stabilization of the rate-determining transition state for the reaction through external charges, electric dipole moments, or local electric field effects. Furthermore, bond dissociation energies in molecules are shown to be influenced by an electric field effect, and through targeting a specific bond in an electric field, this can lead to an unusually specific reaction. For instance, in the carbon-induced starvation protein, we studied two substrate-bound conformations and showed that regardless of what C-H bond of the substrate is closest to the iron(IV)-oxo oxidant, the lowest hydrogen atom abstraction barrier is always for the pro- S C 2 -H abstraction due to an induced dipole moment of the protein that weakens this bond. In another example of the hygromycin biosynthesis enzyme, an oxidative ring-closure reaction in the substrate forms an ortho-δ-ester ring. Calculations on this enzyme show that the selectivity is guided by a protonated lysine residue in the active site that, through its positive charge, triggers a low energy hydrogen atom abstraction barrier. A final set of examples in this Account discuss the viomycin biosynthesis enzyme and the 2-(trimethylammonio)ethylphosphonate dioxygenase (TmpA) enzyme. Both of these enzymes are shown to possess a significant local dipole moment and local electric field effect due to charged residues surrounding the substrate and oxidant binding pockets. The protein dipole moment and local electric field strength changes the C-H bond strengths of the substrate as compared to the gas-phase triggers the regioselectivity of substrate activation. In particular, we show that in the gas phase and in a protein environment C-H bond strengths are different due to local electric dipole moments and electric field strengths. These examples show that enzymes have an intricately designed structure that enables a chemical reaction under ambient conditions through the positioning of positively and negatively charged residues that influence and enhance reaction mechanisms. These computational insights create huge possibilities in bioengineering to apply local electric field and dipole moments in proteins to achieve an unusual selectivity and specificity and trigger a fit-for-purpose biocatalyst for unique biotransformations.

Thiol dioxygenases are important enzymes for human health; they are involved in the detoxification and catabolism of toxic thiol-containing natural products such as cysteine. As such, these enzymes have relevance to the development of Alzheimer's and Parkinson's diseases in the brain. Recent crystal structure coordinates of cysteine and 3-mercaptopropionate dioxygenase (CDO and MDO) showed major differences in the second-coordination spheres of the two enzymes. To understand the difference in activity between these two analogous enzymes, we created large, active-site cluster models. We show that CDO and MDO have different iron(III)-superoxo-bound structures due to differences in ligand coordination. Furthermore, our studies show that the differences in the second-coordination sphere and particularly the position of a positively charged Arg residue results in changes in substrate positioning, mobility and enzymatic turnover. Furthermore, the substrate scope of MDO is explored with cysteinate and 2-mercaptosuccinic acid and their reactivity is predicted., (© 2021 Wiley-VCH GmbH.)

Nonheme iron enzymes often utilize a high-valent iron(IV) oxo species for the biosynthesis of natural products, but their high reactivity often precludes structural and functional studies of these complexes. In this work, a combined experimental and computational study is presented on a biomimetic nonheme iron(IV) oxo complex bearing an aminopyridine macrocyclic ligand and its reactivity toward olefin epoxidation upon changes in the identity and coordination ability of the axial ligand. Herein, we show a dramatic effect of the pH on the oxygen-atom-transfer (OAT) reaction with substrates. In particular, these changes have occurred because of protonation of the axial-bound pendant amine group, where its coordination to iron is replaced by a solvent molecule or anionic ligand. This axial ligand effect influences the catalysis, and we observe enhanced cyclooctene epoxidation yields and turnover numbers in the presence of the unbound protonated pendant amine group. Density functional theory studies were performed to support the experiments and highlight that replacement of the pendant amine with a neutral or anionic ligand dramatically lowers the rate-determining barriers of cyclooctene epoxidation. The computational work further establishes that the change in OAT is due to electrostatic interactions of the pendant amine cation that favorably affect the barrier heights.

Free energy drives a wide range of molecular processes such as solvation, binding, chemical reactions and conformational change. Given the central importance of binding, a wide range of methods exist to calculate it, whether based on scoring functions, machine-learning, classical or electronic structure methods, alchemy, or explicit evaluation of energy and entropy. Here we present a new energy-entropy (EE) method to calculate the host-guest binding free energy directly from molecular dynamics (MD) simulation. Entropy is evaluated using Multiscale Cell Correlation (MCC) which uses force and torque covariance and contacts at two different length scales. The method is tested on a series of seven host-guest complexes in the SAMPL8 (Statistical Assessment of the Modeling of Proteins and Ligands) "Drugs of Abuse" Blind Challenge. The EE-MCC binding free energies are found to agree with experiment with an average error of 0.9 kcal mol -1 . MCC makes clear the origin of the entropy changes, showing that the large loss of positional, orientational, and to a lesser extent conformational entropy of each binding guest is compensated for by a gain in orientational entropy of water released to bulk, combined with smaller decreases in vibrational entropy of the host, guest and contacting water., (© 2021. The Author(s).)

Nature employs high-energy metal-oxo intermediates embedded within enzyme active sites to perform challenging oxidative transformations with remarkable selectivity. Understanding how different local metal-oxo coordination environments control intermediate reactivity and catalytic function is a long-standing objective. However, conducting structure-activity relationships directly in active sites has proven challenging due to the limited range of amino acid substitutions achievable within the constraints of the genetic code. Here, we use an expanded genetic code to examine the impact of hydrogen bonding interactions on ferryl heme structure and reactivity, by replacing the N-H group of the active site Trp51 of cytochrome c peroxidase by an S atom. Removal of a single hydrogen bond stabilizes the porphyrin π-cation radical state of C c P W191F compound I. In contrast, this modification leads to more basic and reactive neutral ferryl heme states, as found in C c P W191F compound II and the wild-type ferryl heme-Trp191 radical pair of compound I. This increased reactivity manifests in a >60-fold activity increase toward phenolic substrates but remarkably has negligible effects on oxidation of the biological redox partner cyt c . Our data highlight how Trp51 tunes the lifetimes of key ferryl intermediates and works in synergy with the redox active Trp191 and a well-defined substrate binding site to regulate catalytic function. More broadly, this work shows how noncanonical substitutions can advance our understanding of active site features governing metal-oxo structure and reactivity., Competing Interests: The authors declare no competing financial interest., (© 2021 The Authors. Published by American Chemical Society.)

There are two types of cytochrome P450 enzymes in nature, namely, the monooxygenases and the peroxygenases. Both enzyme classes participate in substrate biodegradation or biosynthesis reactions in nature, but the P450 monooxygenases use dioxygen, while the peroxygenases take H 2 O 2 in their catalytic cycle instead. By contrast to the P450 monooxygenases, the P450 peroxygenases do not require an external redox partner to deliver electrons during the catalytic cycle, and also no external proton source is needed. Therefore, they are fully self-sufficient, which affords them opportunities in biotechnological applications. One specific P450 peroxygenase, namely, P450 OleT JE , reacts with long-chain linear fatty acids through oxidative decarboxylation to form hydrocarbons and, as such, has been implicated as a suitable source for the biosynthesis of biofuels. Unfortunately, the reactions were shown to produce a considerable amount of side products originating from C α and C β hydroxylation and desaturation. These product distributions were found to be strongly dependent on whether the substrate had substituents on the C α and/or C β atoms. To understand the bifurcation pathways of substrate activation by P450 OleT JE leading to decarboxylation, C α hydroxylation, C β hydroxylation and C α -C β desaturation, we performed a computational study using 3-phenylpropionate and 2-phenylbutyrate as substrates. We set up large cluster models containing the heme, the substrate and the key features of the substrate binding pocket and calculated (using density functional theory) the pathways leading to the four possible products. This work predicts that the two substrates will react with different reaction rates due to accessibility differences of the substrates to the active oxidant, and, as a consequence, these two substrates will also generate different products. This work explains how the substrate binding pocket of P450 OleT JE guides a reaction to a chemoselectivity.

The nonheme iron dioxygenase 2-(trimethylammonio)-ethylphosphonate dioxygenase (TmpA) is an enzyme involved in the regio- and chemoselective hydroxylation at the C 1 -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 2 -H bond in favor of C 1 -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 1 -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 1 -H and C 2 -H bonds of the substrate, and strongly strengthens the C 2 -H bond in the substrate-bound orientation., (© 2021 Wiley-VCH GmbH.)

Heme superoxides are one of the most versatile metallo-intermediates in biology, and they mediate a vast variety of oxidation and oxygenation reactions involving O 2(g) . Overall proton-coupled electron transfer (PCET) processes they facilitate may proceed via several different mechanistic pathways, attributes of which are not yet fully understood. Herein we present a detailed investigation into concerted PCET events of a series of geometrically similar, but electronically disparate synthetic heme superoxide mimics, where unprecedented, PCET feasibility-determining electronic effects of the heme center have been identified. These electronic factors firmly modulate both thermodynamic and kinetic parameters that are central to PCET, as supported by our experimental and theoretical observations. Consistently, the most electron-deficient superoxide adduct shows the strongest driving force for PCET, whereas the most electron-rich system remains unreactive. The pivotal role of these findings in understanding significant heme systems in biology, as well as in alternative energy applications is also discussed., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

The carbon starvation-induced protein D (CsiD) is a recently characterized iron(II)/α-ketoglutarate-dependent oxygenase that activates a glutarate molecule as substrate at the C 2 position to exclusively produce ( S )-2-hydroxyglutarate products. This selective hydroxylation reaction by CsiD is an important component of the lysine biodegradation pathway in Escherichia coli ; however, little is known on the details and the origin of the selectivity of the reaction. So far, experimental studies failed to trap and characterize any short-lived catalytic cycle intermediates. As no computational studies have been reported on this enzyme either, we decided to investigate the chemical reaction mechanism of glutarate activation by an iron(IV)-oxo model of the CsiD enzyme. In this work, we present a density functional theory study on a large active site cluster model of CsiD and investigate the glutarate hydroxylation pathways by a high-valent iron(IV)-oxo species leading to ( S )-2-hydroxyglutarate, ( R )-2-hydroxyglutarate, and 3-hydroxyglutarate. In agreement with experimental observation, the favorable product channel leads to pro- S C 2 -H hydrogen atom abstraction to form ( S )-2-hydroxyglutarate. The reaction is stepwise with a hydrogen atom abstraction by an iron(IV)-oxo species followed by OH rebound from a radical intermediate. The work presented in this paper shows that despite the fact that the C-H bond strengths at the C 2 and C 3 positions of glutarate are similar in the gas phase, substrate binding and positioning guide the reaction to an enantioselective reaction process by destabilizing the hydrogen atom abstraction pathways for the pro- R C 2 -H and C 3 -H positions. Our studies predict the chemical properties of the iron(IV)-oxo species and its rate constants with glutarate and deuterated-glutarate. Moreover, the work shows little protein motions during the catalytic process, while the substrate entrance into the substrate binding pocket appears to be guided by three active site arginine residues that position the substrate for pro- S C 2 -H hydrogen atom abstraction. Finally, the calculations show that irrespective of the position of the substrate and what C-H bond is closest to the metal center, the lowest energy pathway is for a selective pro- S C 2 -H hydrogen atom abstraction.

DFT calculations reveal that the iodine of cobalt(ii)-iodylarene complexes acts as a directing group via halogen bonding interaction to substrates. A transient 3c-4e bond is formed during oxidation reactions to decrease the activation energy by electron delocalization. Dehydrogenation of dihydroantharacene proceeds via a novel concerted hydride transfer/proton transfer mechanism.

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.