Engineering the Reaction Pathway of a Non-heme Iron Oxygenase Using Ancestral Sequence Reconstruction

Non-heme iron (FeII), α-ketoglutarate (α-KG)-dependent oxygenases are a family of enzymes that catalyze an array of transformations that cascade forward after the formation of radical intermediates. Achieving control over the reaction pathway is highly valuable and a necessary step toward broadening the applications of these biocatalysts. Numerous approaches have been used to engineer the reaction pathway of FeII/α-KG-dependent enzymes, including site-directed mutagenesis, DNA shuffling, and site-saturation mutagenesis, among others. Herein, we showcase a novel ancestral sequence reconstruction (ASR)-guided strategy in which evolutionary information is used to pinpoint the residues critical for controlling different reaction pathways. Following this, a combinatorial site-directed mutagenesis approach was used to quickly evaluate the importance of each residue. These results were validated using a DNA shuffling strategy and through quantum mechanical/molecular mechanical (QM/MM) simulations. Using this approach, we identified a set of active site residues together with a key hydrogen bond between the substrate and an active site residue, which are crucial for dictating the dominant reaction pathway. Ultimately, we successfully converted both extant and ancestral enzymes that perform benzylic hydroxylation into variants that can catalyze an oxidative ring-expansion reaction, showcasing the potential of utilizing ASR to accelerate the reaction pathway engineering within enzyme families that share common structural and mechanistic features.

Secondary Sphere Lewis Acid Activated Heme Superoxo Adducts Mimic Crucial Non-Covalent Interactions in IDO/TDO Heme Dioxygenases

Heme enzymes play a central role in a medley of reactivities within a wide variety of crucial biological systems. Their active sites are highly decorated with pivotal evolutionarily optimized non-covalent interactions that precisely choreograph their biological functionalities with specific regio-, stereo-, and chemo-selectivities. Gaining a clear comprehension of how such weak interactions within the active sites control reactivity offers powerful information to be implemented into the design of future therapeutic agents that target these heme enzymes. To shed light on such critical details pertaining to tryptophan dioxygenating heme enzymes, this study investigates the indole dioxygenation reactivities of Lewis acid-activated heme superoxo model systems, wherein an unprecedented kinetic behavior is revealed. In that, the activated heme superoxo adduct is observed to undergo indole dioxygenation with the intermediacy of a non-covalently organized precursor complex, which forms prior to the rate-limiting step of the overall reaction landscape. Spectroscopic and theoretical characterization of this precursor complex draws close parallels to the ternary complex of heme dioxygenases, which has been postulated to be of crucial importance for successful 2,3-dioxygenative cleavage of indole moieties.

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

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

Stereoelectronic Tuning of Bioinspired Nonheme Iron(IV)-Oxo Species by Amide Groups in Primary and Secondary Coordination Spheres for Selective Oxygenation Reactions

Two mononuclear iron(II) complexes, [(6-amide2-BPMEN)FeII](OTf)2 (1) and [(6-amide-Me-BPMEN)FeII(OTf)](OTf) (2), supported by two BPMEN-derived (BPMEN = N1,N2-dimethyl-N1,N2-bis(pyridine-2-yl-methyl)ethane-1,2-diamine) ligands bearing one or two amide functionalities have been isolated to study their reactivity in the oxygenation of C-H and C═C bonds using isopropyl 2-iodoxybenzoate (iPr-IBX ester) as the oxidant. Both 1 and 2 contain six-coordinate high-spin iron(II) centers in the solid state and in solution. The 6-amide2-BPMEN ligand stabilizes an S = 1 iron(IV)-oxo intermediate, [(6-amide2-BPMEN)FeIV(O)]2+ (1A). The oxidant (1A) oxygenates the C-H and C═C bonds with a high selectivity. Oxidant 1A, upon treatment with 2,6-lutidine, is transformed into another oxidant [{(6-amide2-BPMEN)-(H)}FeIV(O)]+ (1B) through deprotonation of an amide group, resulting in a stronger equatorial ligand field and subsequent stabilization of the triplet ground state. In contrast, no iron-oxo species could be observed from complex 2 and [(6-Me2-BPMEN)FeII(OTf)2] (3) under similar experimental conditions. The iron(IV)-oxo oxidant 1A shows the highest A/K selectivity in cyclohexane oxidation and 3°/2° selectivity in adamantane oxidation reported for any synthetic nonheme iron(IV)-oxo complexes. Theoretical investigation reveals that the hydrogen bonding interaction between the -NH group of the noncoordinating amide group and Fe═O core smears out the equatorial charge density, reducing the triplet-quintet splitting, and thus helping complex 1A to achieve better reactivity.

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

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

Electrostatically Dominated Pre-Organization in Cyclodextrin Metal-Organic Frameworks

Electrostatic interactions between oppositely charged entities play a key role in pre-organizing substrates and stabilizing transition states of reactions in enzymes. The use of electrostatic interactions to pre-organize ions in nanoconfined pores, however, has not been investigated to its full potential. Herein, we describe how carboxylate anions can be pre-organized at the behest of their electrostatic interactions with K+ cations in nanoconfined tunnels present in γ-cyclodextrin metal-organic frameworks, i.e., CD-MOFs. Several carboxylate anions, which are all much smaller than the cavities of the tunnels, were visualized by X-ray crystallography when nanoconfined in CD-MOFs, despite the large voids present in the tunnels. These anions were found to be aligned within a planar array defined by four K+ cations, positioned around the periphery of the tunnels. The strong electrostatic interactions between the carboxylate anions and the K+ cations dictate the orientation of the anions and override the influence of all other possible noncovalent bonding interactions between them and the tunnels. Consequently, the aligned pairs of γ-cyclodextrin rings constituting the tunnels become distorted, resulting in them having lower symmetry and fewer disordered carboxylate anions in the solid state. Our findings offer a transformative strategy for controlling the packing and orientation of ions in nanoconfined environments.

Probing substrate binding inside a paramagnetic cavity: a NMR spectroscopy toolbox for combined experimental and theoretical investigation

Protein cavities often rely on the paramagnetic metal present in their active site in order to catalyse various chemical transformations in biology. The selective detection and identification of the substrate is of fundamental importance in environmental monitoring and biological studies. Herein, a covalently linked Fe(iii)porphyrin dimer-based paramagnetic sensory cavity has been devised for the accurate detection and simultaneous identification of phenol (substrate) binding within the cavity that provides a unique spectroscopic signature with valuable structural and environmental information. These substrates within the paramagnetic cavity leave the fingerprints of the specific binding modes (exo vs. endo) which are well distinguished with the help of various spectroscopic studies viz. UV-vis, 1H, and 19F NMR and in their respective crystal structures also. The theoretical 19F NMR analysis plays a pivotal role in replicating the observed NMR trends with large chemical shifts of the phenolato species which in turn helps in deciphering the selective binding modes of the phenols and thereby recognizing the chemical environment within the cavity. These findings will help develop an excellent diagnostic tool for in situ monitoring of subtle conformational changes and transient interactions.

Unraveling the Mechanism of the Oxidative C-C Bond Coupling Reaction Catalyzed by Deoxypodophyllotoxin Synthase

Deoxypodophyllotoxin synthase (DPS), a nonheme Fe(II)/2-oxoglutarate (2OG)-dependent oxygenase, is a key enzyme that is involved in the construction of the fused-ring system in (-)-podophyllotoxin biosynthesis by catalyzing the C-C coupling reaction. However, the mechanistic details of DPS-catalyzed ring formation remain unclear. Herein, our quantum mechanics/molecular mechanics (QM/MM) calculations reveal a novel mechanism that involves the recycling of CO2 (a product of decarboxylation of 2OG) to prevent the formation of hydroxylated byproducts. Our results show that CO2 can react with the FeIII-OH species to generate an unusual FeIII-bicarbonate species. In this way, hydroxylation is avoided by consuming the OH group. Then, the C-C coupling followed by desaturation yields the final product, deoxypodophyllotoxin. This work highlights the crucial role of the CO2 molecule, generated in the crevice between the iron active site and the substrate, in controlling the reaction selectivity.

The Unique Role of the Second Coordination Sphere to Unlock and Control Catalysis in Nonheme Fe(II)/2-Oxoglutarate Histone Demethylase KDM2A

Nonheme Fe(II) and 2-oxoglutarate (2OG)-dependent histone lysine demethylases 2A (KDM2A) catalyze the demethylation of the mono- or dimethylated lysine 36 residue in the histone H3 peptide (H3K36me1/me2), which plays a crucial role in epigenetic regulation and can be involved in many cancers. Although the overall catalytic mechanism of KDMs has been studied, how KDM2 catalysis takes place in contrast to other KDMs remains unknown. Understanding such differences is vital for enzyme redesign and can help in enzyme-selective drug design. Herein, we employed molecular dynamics (MD) and combined quantum mechanics/molecular mechanics (QM/MM) to explore the complete catalytic mechanism of KDM2A, including dioxygen diffusion and binding, dioxygen activation, and substrate oxidation. Our study demonstrates that the catalysis of KDM2A is controlled by the conformational change of the second coordination sphere (SCS), specifically by a change in the orientation of Y222, which unlocks the 2OG rearrangement from off-line to in-line mode. The study demonstrates that the variant Y222A makes the 2OG rearrangement more favorable. Furthermore, the study reveals that it is the size of H3K36me3 that prevents the 2OG rearrangement, thus rendering the enzyme inactivity with trimethylated lysine. Calculations show that the SCS and long-range interacting residues that stabilize the HAT transition state in KDM2A differ from those in KDM4A, KDM7B, and KDM6A, thus providing the basics for the enzyme-selective redesign and modulation of KDM2A without influencing other KDMs.

Spin-Flip via Subtle Electronic Perturbation in Axially Ligated Diiron(III) Porphyrin Dimer

Spin state switching in the metal center is a crucial phenomenon in many enzymatic reactions in biology. The spin state alteration, a critical step in cytochrome P450 catalysis, is driven most likely through a weak perturbation upon substrate binding in the enzyme, which is still not well clarified. In the current work, the spin state transition of iron(III) from high to intermediate via an admixed state is observed upon a subtle electronic perturbation to the sulphonate moieties coordinated axially to a diiron(III)porphyrin dimer. While electron-donating substituents stabilize the high-spin state of iron(III), strongly electron-withdrawing groups stabilize an intermediate-spin state, whereas the moderate electron-withdrawing nature of axial ligands resulted in an admixed state. Confirmation of the molecular structures and their spin states have been made utilizing single-crystal X-ray structure analysis, Mössbauer, magnetic, EPR, and 1H NMR spectroscopic investigations. The position of the signals of the porphyrin macrocycle in the paramagnetic 1H NMR is found to be very characteristic of the spin state of the iron center in solution. The Curie plot for the pure high-spin complexes shows the signals' temperature dependency in line with the Curie law. Conversely, the pure intermediate-spin state of iron exhibits an anti-Curie temperature dependence, whereas the admixed-spin state of iron displays significant curvature of the lines in the Curie plot. An extensive DFT analysis displays a linear dependence between the energy difference between dx 2 - y 2 ${{_{x{^{2}}- y{^{2}}}}}$ and dz 2 ${{_{z{^{2}}}}}$ orbital versus Fe-Npor distance for the complexes reported here. Furthermore, a strong linear correlation between the Fe-O distance and the spin density over the oxygen atom, as well as the Fe-Npor distance for the complexes, has been observed. Thus, a slight electronic perturbation at the axial ligand of the diheme resulted in a large change in the electronic structures with a spin-flip. This is at par with the metalloenzymes, which employ minute perturbations around the periphery of the active sites, leading to spin state transitions.

Enhanced Reactivities of Iron(IV)-Oxo Porphyrin Species in Oxidation Reactions Promoted by Intramolecular Hydrogen-Bonding

High-valent iron-oxo species are one of the common intermediates in both biological and biomimetic catalytic oxidation reactions. Recently, hydrogen-bonding (H-bonding) has been proved to be critical in determining the selectivity and reactivity. However, few examples have been established for mechanistic insights into the H-bonding effect. Moreover, intramolecular H-bonding effect on both C-H activation and oxygen atom transfer (OAT) reactions in synthetic porphyrin model system has not been investigated yet. In this study, a series of heme-containing iron(IV)-oxo porphyrin species with or without intramolecular H-bonding are synthesized and characterized. Kinetic studies revealed that intramolecular H-bonding can significantly enhance the reactivity of iron(IV)-oxo species in OAT, C-H activation, and electron-transfer reactions. This unprecedented unified H-bonding effect is elucidated by theoretical calculations, which showed that intramolecular H-bonding interactions lower the energy of the anti-bonding orbital of iron(IV)-oxo porphyrin species, resulting in the enhanced reactivities in oxidation reactions irrespective of the reaction type. To the best of the knowledge, this is the first extensive investigation on the intramolecular H-bonding effect in heme system. The results show that H-bonding interactions have a unified effect with iron(IV)-oxo porphyrin species in all three investigated reactions.

The thin line between monooxygenases and peroxygenases. P450s, UPOs, MMOs, and LPMOs: A brick to bridge fields of expertise

Many scientific fields, although driven by similar purposes and dealing with similar technologies, often appear so isolated and far from each other that even the vocabularies to describe the very same phenomenon might differ. Concerning the vast field of biocatalysis, a special role is played by those redox enzymes that employ oxygen-based chemistry to unlock transformations otherwise possible only with metal-based catalysts. As such, greener chemical synthesis methods and environmentally-driven biotechnological approaches were enabled over the last decades by the use of several enzymes and ultimately resulted in the first industrial applications. Among what can be called today the environmental biorefinery sector, biomass transformation, greenhouse gas reduction, bio-gas/fuels production, bioremediation, as well as bulk or fine chemicals and even pharmaceuticals manufacturing are all examples of fields in which successful prototypes have been demonstrated employing redox enzymes. In this review we decided to focus on the most prominent enzymes (MMOs, LPMO, P450 and UPO) capable of overcoming the ∼100 kcal mol-1 barrier of inactivated CH bonds for the oxyfunctionalization of organic compounds. Harnessing the enormous potential that lies within these enzymes is of extreme value to develop sustainable industrial schemes and it is still deeply coveted by many within the aforementioned fields of application. Hence, the ambitious scope of this account is to bridge the current cutting-edge knowledge gathered upon each enzyme. By creating a broad comparison, scientists belonging to the different fields may find inspiration and might overcome obstacles already solved by the others. This work is organised in three major parts: a first section will be serving as an introduction to each one of the enzymes regarding their structural and activity diversity, whereas a second one will be encompassing the mechanistic aspects of their catalysis. In this regard, the machineries that lead to analogous catalytic outcomes are depicted, highlighting the major differences and similarities. Finally, a third section will be focusing on the elements that allow the oxyfunctionalization chemistry to occur by delivering redox equivalents to the enzyme by the action of diverse redox partners. Redox partners are often overlooked in comparison to the catalytic counterparts, yet they represent fundamental elements to better understand and further develop practical applications based on mono- and peroxygenases.

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

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

Structural, Spectroscopic, and Computational Insights from Canavanine-Bound and Two Catalytically Compromised Variants of the Ethylene-Forming Enzyme

The ethylene-forming enzyme (EFE) is an Fe(II), 2-oxoglutarate (2OG), and l-arginine (l-Arg)-dependent oxygenase that either forms ethylene and three CO2/bicarbonate from 2OG or couples the decarboxylation of 2OG to C5 hydroxylation of l-Arg. l-Arg binds with C5 toward the metal center, causing 2OG to change from monodentate to chelate metal interaction and OD1 to OD2 switch of D191 metal coordination. We applied anaerobic UV-visible spectroscopy, X-ray crystallography, and computational approaches to three EFE systems with high-resolution structures. The ineffective l-Arg analogue l-canavanine binds to the EFE with O5 pointing away from the metal center while promoting chelate formation by 2OG but fails to switch the D191 metal coordination from OD1 to OD2. Substituting alanine for R171 that interacts with 2OG and l-Arg inactivates the protein, prevents metal chelation by 2OG, and weakens l-Arg binding. The R171A EFE had electron density at the 2OG binding site that was identified by mass spectrometry as benzoic acid. The substitution by alanine of Y306 in the EFE, a residue 12 Å away from the catalytic metal center, generates an interior cavity that leads to multiple local and distal structural changes that reduce l-Arg binding and significantly reduce the enzyme activity. Flexibility analyses revealed correlated and anticorrelated motions in each system, with important distinctions from the wild-type enzyme. In combination, the results are congruent with the currently proposed enzyme mechanism, reinforce the importance of metal coordination by OD2 of D191, and highlight the importance of the second coordination sphere and longer range interactions in promoting EFE activity.

Computational Study Into the Oxidative Ring-Closure Mechanism During the Biosynthesis of Deoxypodophyllotoxin

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 C7'-H bond closest to the iron(IV)-oxo species. The reaction proceeds through a radical mechanism starting with hydrogen atom abstraction from the C7'-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.

Unusual catalytic strategy by non-heme Fe(ii)/2-oxoglutarate-dependent aspartyl hydroxylase AspH

Biocatalytic C-H oxidation reactions are of important synthetic utility, provide a sustainable route for selective synthesis of important organic molecules, and are an integral part of fundamental cell processes. The multidomain non-heme Fe(ii)/2-oxoglutarate (2OG) dependent oxygenase AspH catalyzes stereoselective (3R)-hydroxylation of aspartyl- and asparaginyl-residues. Unusually, compared to other 2OG hydroxylases, crystallography has shown that AspH lacks the carboxylate residue of the characteristic two-His-one-Asp/Glu Fe-binding triad. Instead, AspH has a water molecule that coordinates Fe(ii) in the coordination position usually occupied by the Asp/Glu carboxylate. Molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) studies reveal that the iron coordinating water is stabilized by hydrogen bonding with a second coordination sphere (SCS) carboxylate residue Asp721, an arrangement that helps maintain the six coordinated Fe(ii) distorted octahedral coordination geometry and enable catalysis. AspH catalysis follows a dioxygen activation-hydrogen atom transfer (HAT)-rebound hydroxylation mechanism, unusually exhibiting higher activation energy for rebound hydroxylation than for HAT, indicating that the rebound step may be rate-limiting. The HAT step, along with substrate positioning modulated by the non-covalent interactions with SCS residues (Arg688, Arg686, Lys666, Asp721, and Gln664), are essential in determining stereoselectivity, which likely proceeds with retention of configuration. The tetratricopeptide repeat (TPR) domain of AspH influences substrate binding and manifests dynamic motions during catalysis, an observation of interest with respect to other 2OG oxygenases with TPR domains. The results provide unique insights into how non-heme Fe(ii) oxygenases can effectively catalyze stereoselective hydroxylation using only two enzyme-derived Fe-ligating residues, potentially guiding enzyme engineering for stereoselective biocatalysis, thus advancing the development of non-heme Fe(ii) based biomimetic C-H oxidation catalysts, and supporting the proposal that the 2OG oxygenase superfamily may be larger than once perceived.

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

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

From Biomimicking to Bioinspired Design of Electrocatalysts for CO2 Reduction to C1 Products

The electrochemical reduction of CO2 (CO2 RR) is a promising approach to maintain a carbon cycle balance and produce value-added chemicals. However, CO2 RR technology is far from mature, since the conventional CO2 RR electrocatalysts suffer from low activity (leading to currents <10 mA cm-2 in an H-cell), stability (<120 h), and selectivity. Hence, they cannot meet the requirements for commercial applications (>200 mA cm-2 , >8000 h, >90 % selectivity). Significant improvements are possible by taking inspiration from nature, considering biological organisms that efficiently catalyze the CO2 to various products. In this minireview, we present recent examples of enzyme-inspired and enzyme-mimicking CO2 RR electrocatalysts enabling the production of C1 products with high faradaic efficiency (FE). At present, these designs do not typically follow a methodical approach, but rather focus on isolated features of biological systems. To achieve disruptive change, we advocate a systematic design methodology that leverages fundamental mechanisms associated with desired properties in nature and adapts them to the context of engineering applications.

Monovalent metal ion binding promotes the first transesterification reaction in the spliceosome

Cleavage and formation of phosphodiester bonds in nucleic acids is accomplished by large cellular machineries composed of both protein and RNA. Long thought to rely on a two-metal-ion mechanism for catalysis, structure comparisons revealed many contain highly spatially conserved second-shell monovalent cations, whose precise function remains elusive. A recent high-resolution structure of the spliceosome, essential for pre-mRNA splicing in eukaryotes, revealed a potassium ion in the active site. Here, we employ biased quantum mechanics/ molecular mechanics molecular dynamics to elucidate the function of this monovalent ion in splicing. We discover that the K+ ion regulates the kinetics and thermodynamics of the first splicing step by rigidifying the active site and stabilizing the substrate in the pre- and post-catalytic state via formation of key hydrogen bonds. Our work supports a direct role for the K+ ion during catalysis and provides a mechanistic hypothesis likely shared by other nucleic acid processing enzymes.

Computationally guided bioengineering of the active site, substrate access pathway, and water channels of thermostable cytochrome P450, CYP175A1, for catalyzing the alkane hydroxylation reaction

Understanding structure-function relationships in proteins is pivotal in their development as industrial biocatalysts. In this regard, rational engineering of protein active site access pathways and various tunnels and channels plays a central role in designing competent enzymes with high stability and enhanced efficiency. Here, we report the rational evolution of a thermostable cytochrome P450, CYP175A1, to catalyze the C-H activation reaction of longer-chain alkanes. A strategy combining computational tools with experiments has shown that the substrate scope and enzymatic activity can be enhanced by rational engineering of certain important channels such as the substrate entry and water channels along with the active site of the enzyme. The evolved enzymes showed an improved catalytic rate for hexadecane hydroxylation with high regioselectivity. The Q67L/Y68F mutation showed binding of the substrate in the active site, water channel mutation L80F/V220T showed improved catalytic activity through the peroxide shunt pathway and substrate entry channel mutation W269F/I270A showed better substrate accessibility to the active pocket. All-atom MD simulations provided the rationale for the inactivity of the wild-type CYP175A1 for hexadecane hydroxylation and predicted the above hot-spot residues to enhance the activity. The reaction mechanism was studied by QM/MM calculations for enzyme-substrate complexes and reaction intermediates. Detailed thermal and thermodynamic stability of all the mutants were analyzed and the results showed that the evolved enzymes were thermally stable. The present strategy showed promising results, and insights gained from this work can be applied to the general enzymatic system to expand substrate scope and improve catalytic activity.

Mechanisms and Opportunities for Rational In Silico Design of Enzymes to Degrade Per- and Polyfluoroalkyl Substances (PFAS)

Per and polyfluoroalkyl substances (PFAS) present a unique challenge to remediation techniques because their strong carbon-fluorine bonds make them difficult to degrade. This review explores the use of in silico enzymatic design as a potential PFAS degradation technique. The scope of the enzymes included is based on currently known PFAS degradation techniques, including chemical redox systems that have been studied for perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) defluorination, such as those that incorporate hydrated electrons, sulfate, peroxide, and metal catalysts. Bioremediation techniques are also discussed, namely the laccase and horseradish peroxidase systems. The redox potential of known reactants and enzymatic radicals/metal-complexes are then considered and compared to potential enzymes for degrading PFAS. The molecular structure and reaction cycle of prospective enzymes are explored. Current knowledge and techniques of enzyme design, particularly radical-generating enzymes, and application are also discussed. Finally, potential routes for bioengineering enzymes to enable or enhance PFAS remediation are considered as well as the future outlook for computational exploration of enzymatic in situ bioremediation routes for these highly persistent and globally distributed contaminants.

From random to rational: improving enzyme design through electric fields, second coordination sphere interactions, and conformational dynamics

Enzymes are versatile and efficient biological catalysts that drive numerous cellular processes, motivating the development of enzyme design approaches to tailor catalysts for diverse applications. In this perspective, we investigate the unique properties of natural, evolved, and designed enzymes, recognizing their strengths and shortcomings. We highlight the challenges and limitations of current enzyme design protocols, with a particular focus on their limited consideration of long-range electrostatic and dynamic effects. We then delve deeper into the impact of the protein environment on enzyme catalysis and explore the roles of preorganized electric fields, second coordination sphere interactions, and protein dynamics for enzyme function. Furthermore, we present several case studies illustrating successful enzyme-design efforts incorporating enzyme strategies mentioned above to achieve improved catalytic properties. Finally, we envision the future of enzyme design research, spotlighting the challenges yet to be overcome and the synergy of intrinsic electric fields, second coordination sphere interactions, and conformational dynamics to push the state-of-the-art boundaries.

Insights into Cytochrome P450 Enzyme Catalyzed Defluorination of Aromatic Fluorides

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.

Melatonin Activation by Human Cytochrome P450 Enzymes: A Comparison between Different Isozymes

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

Catalysis by KDM6 Histone Demethylases - A Synergy between the Non-Heme Iron(II) Center, Second Coordination Sphere, and Long-Range Interactions

KDM6A (UTX) and KDM6B (JMJD3) are human non-heme Fe(II) and 2-oxoglutarate (2OG) dependent JmjC oxygenases that catalyze the demethylation of trimethylated lysine 27 in the N-terminal tail of histone H3, a post-translational modification that regulates transcription. A Combined Quantum Mechanics/ Molecular Mechanics (QM/MM) and Molecular Dynamics (MD) study on the catalytic mechanism of KDM6A/B reveals that the transition state for the rate-limiting hydrogen atom transfer (HAT) reaction in KDM6A catalysis is stabilized by polar (Asn217) and aromatic (Trp369)/non-polar (Pro274) residues in contrast to KDM4, KDM6B and KDM7 demethylases where charged residues (Glu, Arg, Asp) are involved. KDM6A employs both σ- and π-electron transfer pathways for HAT, whereas KDM6B employs the σ-electron pathway. Differences in hydrogen bonding of the Fe-chelating Glu252(KDM6B) contribute to the lower energy barriers in KDM6B vs. KDM6A. The study reveals a dependence of the activation barrier of the rebound hydroxylation on the Fe-O-C angle in the transition state of KDM6A. Anti-correlation of the Zn-binding domain with the active site residues is a key factor distinguishing KDM6A/B from KDM7/4s. The results reveal the importance of communication between the Fe center, second coordination sphere, and long-range interactions in catalysis by KDMs and, by implication, other 2OG oxygenases.

Regioselective oxidation of heterocyclic aromatic hydrocarbons catalyzed by cytochrome P450: A case study of carbazole

Recently there are increasing interests in accurately evaluating the health effects of heterocyclic PAHs. However, the activation mechanism and possible metabolites of heterocyclic PAHs catalyzed by human CYP1A1 is still elusive to a great extent. Here, leveraged to high level QM/MM calculations, the corresponding activation pathways of a representative heterocyclic PAHs, carbazole, were systematically explored. The first stage is electrophilic addition or hydrogen abstraction from N-H group. Electrophilic addition was evidenced to be more feasible and regioselectivity at C3 and C4 sites were identified. Correlations between energy barriers and key structural/electrostatic parameters reveal that O-Cα distance and Fe-O-Cα angle are the main origin for the catalytic regioselectivity. Electrophilic addition was determined as the rate-determining step and the subsequent possible reactions include epoxidation, NIH shift (the hydrogen migration from the site of hydroxylation to the adjacent carbon) and proton shuttle. The corresponding products are epoxides, ketones and hydroxylated carbazoles, respectively. The main metabolites (hydroxylated carbazoles) are estimated to be more toxic than carbazole. The regioselectivity of carbazole activated by CYP1A1 is different from the environmental processes (gas and aqueous phase). Collectively, these results will inform the in-depth understanding the metabolic processes of heterocyclic PAHs and aid the accurate evaluation of their health effects.

Dioxygen Binding Is Controlled by the Protein Environment in Non-heme FeII and 2-Oxoglutarate Oxygenases: A Study on Histone Demethylase PHF8 and an Ethylene-Forming Enzyme

This study investigates dioxygen binding and 2-oxoglutarate (2OG) coordination by two model non-heme FeII /2OG enzymes: a class 7 histone demethylase (PHF8) that catalyzes the hydroxylation of its H3K9me2 histone substrate leading to demethylation reactivity and the ethylene-forming enzyme (EFE), which catalyzes two competing reactions of ethylene generation and substrate l-Arg hydroxylation. Although both enzymes initially bind 2OG by using an off-line 2OG coordination mode, in PHF8, the substrate oxidation requires a transition to an in-line mode, whereas EFE is catalytically productive for ethylene production from 2OG in the off-line mode. We used classical molecular dynamics (MD), quantum mechanics/molecular mechanics (QM/MM) MD and QM/MM metadynamics (QM/MM-MetD) simulations to reveal that it is the dioxygen binding process and, ultimately, the protein environment that control the formation of the in-line FeIII -OO⋅- intermediate in PHF8 and the off-line FeIII -OO⋅- intermediate in EFE.

Emerging Metabolic Profiles of Sulfonamide Antibiotics by Cytochromes P450: A Computational-Experimental Synergy Study on Emerging Pollutants

Metabolism, especially by CYP450 enzymes, is the main reason for mediating the toxification and detoxification of xenobiotics in humans, while some uncommon metabolic pathways, especially for emerging pollutants, probably causing idiosyncratic toxicity are easily overlooked. The pollution of sulfonamide antibiotics in aqueous system has attracted increasing public attention. Hydroxylation of the central amine group can trigger a series of metabolic processes of sulfonamide antibiotics in humans; however, this work parallelly reported the coupling and fragmenting initiated by amino H-abstraction of sulfamethoxazole (SMX) catalyzed by human CYP450 enzymes. Elucidation of the emerging metabolic profiles was mapped via a multistep synergy between computations and experiments, involving preliminary DFT computations and in vitro and in vivo assays, profiling adverse effects, and rationalizing the fundamental factors via targeted computations. Especially, the confirmed SMX dimer was shown to potentially act as a metabolism disruptor in humans, while spin aromatic delocalization resulting in the low electron donor ability of amino radicals was revealed as the fundamental factor to enable coupling of sulfonamide antibiotics by CYP450 through the nonconventional nonrebound pathway. This work may further strengthen the synergistic use of computations prior to experiments to avoid wasteful experimental screening efforts in environmental chemistry and toxicology.

Hydrogen Bonding-Assisted and Nonheme Manganese-Catalyzed Remote Hydroxylation of C-H Bonds in Nitrogen-Containing Molecules

The development of catalytic systems capable of oxygenating unactivated C-H bonds with excellent site-selectivity and functional group tolerance under mild conditions remains a challenge. Inspired by the secondary coordination sphere (SCS) hydrogen bonding in metallooxygenases, reported herein is an SCS solvent hydrogen bonding strategy that employs 1,1,1,3,3,3-hexafluoroisopropanol (HFIP) as a strong hydrogen bond donor solvent to enable remote C-H hydroxylation in the presence of basic aza-heteroaromatic rings with a low loading of a readily available and inexpensive manganese complex as a catalyst and hydrogen peroxide as a terminal oxidant. We demonstrate that this strategy represents a promising compliment to the current state-of-the-art protection approaches that rely on precomplexation with strong Lewis and/or Brønsted acids. Mechanistic studies with experimental and theoretical approaches reveal the existence of a strong hydrogen bonding between the nitrogen-containing substrate and HFIP, which prevents the catalyst deactivation by nitrogen binding and deactivates the basic nitrogen atom toward oxygen atom transfer and the α-C-H bonds adjacent to the nitrogen center toward H-atom abstraction. Moreover, the hydrogen bonding exerted by HFIP has also been demonstrated not only to facilitate the O-O bond heterolytic cleavage of a putative Mn^III-OOH precursor to generate Mn^V(O)(OC(O)CH₂Br) as an active oxidant but also to affect the stability and the activity of Mn^V(O)(OC(O)CH₂Br).

Melatonin Activation by Cytochrome P450 Isozymes: How Does CYP1A2 Compare to CYP1A1?

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.

Biotransformation of Bisphenol by Human Cytochrome P450 2C9 Enzymes: A Density Functional Theory Study

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

Can Second Coordination Sphere and Long-Range Interactions Modulate Hydrogen Atom Transfer in a Non-Heme Fe(II)-Dependent Histone Demethylase?

Fe(II)-dependent oxygenases employ hydrogen atom transfer (HAT) to produce a myriad of products. Understanding how such enzymes use dynamic processes beyond the immediate vicinity of the active site to control the selectivity and efficiency of HAT is important for metalloenzyme engineering; however, obtaining such knowledge by experiments is challenging. This study develops a computational framework for identifying second coordination sphere (SCS) and especially long-range (LR) residues relevant for catalysis through dynamic cross-correlation analysis (DCCA) using the human histone demethylase PHF8 (KDM7B) as a model oxygenase. Furthermore, the study explores the mechanistic pathways of influence of the SCS and LR residues on the HAT reaction. To demonstrate the plausibility of the approach, we investigated the effect of a PHF8 F279S clinical mutation associated with X-linked intellectual disability, which has been experimentally shown to ablate PHF8-catalyzed demethylation. In agreement, the molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) studies showed a change in the H3_1-14K9me2 substrate orientation and an increased HAT barrier. We systematically analyzed the pathways by which the identified SCS and LR residues may influence HAT by exploring changes in H3K9me2 substrate orientation, interdomain correlated motions, HAT transition state stabilization, reaction energetics, electron transfer mechanism, and alterations in the intrinsic electric field of PHF8. Importantly, SCS and LR variations decrease key motions of α9-α12 of the JmjC domain toward the Fe(IV)-center that are associated with tighter binding of the H3_1-14K9me2 substrate. SCS and LR residues alter the intrinsic electric field of the enzyme along the reaction coordinate and change the individual energetic contributions of residues toward TS stabilization. The overall results suggest that DCCA can indeed identify non-active-site residues relevant for catalysis. The substitutions of such dynamically correlated residues might be used as a tool to tune HAT in non-heme Fe(II)- and 2OG-dependent enzymes.

Effect of Redox-Inactive Metal Ion-Nickel(III) Interactions on the Redox Properties and Proton-Coupled Electron Transfer Reactivity

Mononuclear nickel(II) and nickel(III) complexes of a bisamidate-bisalkoxide ligand, (NMe₄)₂[Ni^II(HMPAB)] (1) and (NMe₄)[Ni^III(HMPAB)] (2), respectively, have been synthesized and characterized by various spectroscopic techniques including X-ray crystallography. The reaction of redox-inactive metal ions (Mⁿ⁺ = Ca²⁺, Mg²⁺, Zn²⁺, Y³⁺, and Sc³⁺) with 2 resulted in 2-Mⁿ⁺ adducts, which was assessed by an array of spectroscopic techniques including X-ray absorption spectroscopy (XAS), electron paramagnetic resonance (EPR), and reactivity studies. The X-ray structure of Ca²⁺ coordinated to Ni(III) complexes, 2-Ca²⁺T, was determined and exhibited an average Ni-Ca distance of 3.1253 Å, close to the metal ions' covalent radius. XAS analysis of 2-Ca²⁺ and 2-Y³⁺ in solution further revealed an additional coordination to Ca and Y in the 2-Mⁿ⁺ adducts with shortened Ni-M distances of 2.15 and 2.11 Å, respectively, implying direct bonding interactions between Ni and Lewis acids (LAs). Such a short interatomic distance between Ni(III) and M is unprecedented and was not observed before. EPR analysis of 2 and 2-Mⁿ⁺ species, moreover, displayed rhombic signals with g_av > 2.12 for all complexes, supporting the +III oxidation state of Ni. The Ni^III/Ni^II redox potential of 2 and 2-Mⁿ⁺ species was determined, and a plot of E_1/2 of 2-Mⁿ⁺ versus pK_a of [M(H₂O)_n]^m+ exhibited a linear relationship, implying that the Ni^III/Ni^II potential of 2 can be tuned with different redox-inactive metal ions. Reactivity studies of 2 and 2-Mⁿ⁺ with different 4-X-2,6-ditert-butylphenol (4-X-DTBP) and other phenol derivatives were performed, and based on kinetic studies, we propose the involvement of a proton-coupled electron transfer (PCET) pathway. Analysis of the reaction products after the reaction of 2 with 4-OMe-DTBP showed the formation of a Ni(II) complex (1a) where one of the alkoxide arms of the ligand is protonated. A pK_a value of 24.2 was estimated for 1a. The reaction of 2-Mⁿ⁺ species was examined with 4-OMe-DTBP, and it was observed that the k₂ values of 2-Mⁿ⁺ species increase by increasing the Lewis acidity of redox-inactive metal ions. However, the obtained k₂ values for 2-Mⁿ⁺ species are much lower compared to the k₂ value for 2. Such a variation of PCET reactivity between 2 and 2-Mⁿ⁺ species may be attributed to the interactions between Ni(III) and LAs. Our findings show the significance of the secondary coordination sphere effect on the PCET reactivity of Ni(III) complexes and furnish important insights into the reaction mechanism involving high-valent nickel species, which are frequently invoked as key intermediates in Ni-mediated enzymatic reactions, solar-fuel catalysis, and biomimetic/synthetic transformation reactions.

Charge Maintenance during Catalysis in Nonheme Iron Oxygenases

Here, the choice of the first coordination shell of the metal center is analyzed from the perspective of charge maintenance in a binary enzyme-substrate complex and an O₂-bound ternary complex in the nonheme iron oxygenases. Comparing homogentisate 1,2-dioxygenase and gentisate dioxygenase highlights the significance of charge maintenance after substrate binding as an important factor that drives the reaction coordinate. We then extend the charge analysis to several common types of nonheme iron oxygenases containing either a 2-His-1-carboxylate facial triad or a 3-His or 4-His ligand motif, including extradiol and intradiol ring-cleavage dioxygenases, thiol dioxygenases, α-ketoglutarate-dependent oxygenases, and carotenoid cleavage oxygenases. After forming the productive enzyme-substrate complex, the overall charge of the iron complex at the 0, +1, or +2 state is maintained in the remaining catalytic steps. Hence, maintaining a constant charge is crucial to promote the reaction of the iron center beginning from the formation of the Michaelis or ternary complex. The charge compensation to the iron ion is tuned not only by protein-derived carboxylate ligands but also by substrates. Overall, these analyses indicate that charge maintenance at the iron center is significant when all the necessary components form a productive complex. This charge maintenance concept may apply to most oxygen-activating metalloenzymes systems that do not draw electrons and protons step-by-step from a separate reactant, such as NADH, via a reductase. The charge maintenance perception may also be useful in proposing catalytic pathways or designing prototypical reactions using artificial or engineered enzymes for biotechnological applications.

Second coordination sphere effects on the mechanistic pathways for dioxygen activation by a ferritin: involvement of a Tyr radical and the identification of a cation binding site

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 Asp65, 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.

Peptidic Scaffolds Enable Rapid and Multivariate Secondary Sphere Evolution for an Abiotic Metallocatalyst

Metalloenzymes have benefited from the iterative process of evolution to achieve the precise arrangements of secondary sphere non-covalent interactions that enhance metal-centered catalysis. Iterative synthesis of scaffolds that display complex secondary sphere elements in abiotic systems can be highly challenging and time-intensive. To overcome this synthetic bottleneck, we developed a highly modular and rapid synthetic strategy, leveraging the efficiency of solid-phase peptide synthesis and conformational control afforded by non-canonical residues to construct a ligand platform displaying up to four unique residues of varying electronics and sterics in the secondary coordination sphere. As a proof-of-concept that peptidic secondary sphere can cooperate with the metal complex, we applied this scaffold to a well-known, modestly active C-H oxidizing Fe catalyst to evolve specific non-covalent interactions that is more than double its catalytic activity. Solution-state NMR structures of several catalyst variants suggest that higher activity is correlated with a hydrophobic pocket above the Fe center that may enhance the formation of a catalyst-substrate complex. Above all, we show that peptides are a convenient, highly modular, and structurally defined ligand platform for creating secondary coordination spheres that comprise multiple, diverse functional groups.

On the roles of methionine and the importance of its microenvironments in redox metalloproteins

The amino acid residue methionine (Met) is commonly thought of as a ligand in redox metalloproteins, for example in cytochromes c and in blue copper proteins. However, the roles of Met can go beyond a simple ligand. The thioether functional group of Met allows it to be considered as a hydrophobic residue as well as one that is capable of weak dipolar interactions. In addition, the lone pairs on sulphur allow Met to interact with other groups, inluding the aforementioned metal ions. Because of its properties, Met can play diverse roles in metal coordination, fine tuning of redox reactions, or supporting protein structures. These roles are strongly influenced by the nature of the surrounding medium. Herein, we describe several common interactions between Met and surrounding aromatic amino acids and how they affect the physical properties of both copper and iron metalloproteins. While the importance of interactions between Met and other groups is established in biological systems, less is known about their roles in redox metalloproteins and our view is that this is an area that is ready for greater attention.

Properties of the Reactants and Their Interactions within and with the Enzyme Binding Cavity Determine Reaction Selectivities. The Case of Fe(II)/2-Oxoglutarate Dependent Enzymes

Fe(II)/2-oxoglutarate dependent dioxygenases (ODDs) share a double stranded beta helix (DSBH) fold and utilise a common reactive intermediate, ferryl species, to catalyse oxidative transformations of substrates. Despite the structural similarities, ODDs accept a variety of substrates and facilitate a wide range of reactions, that is hydroxylations, desaturations, (oxa)cyclisations and ring rearrangements. In this review we present and discuss the factors contributing to the observed (regio)selectivities of ODDs. They span from inherent properties of the reactants, that is, substrate molecule and iron cofactor, to the interactions between the substrate and the enzyme's binding cavity; the latter can counterbalance the effect of the former. Based on results of both experimental and computational studies dedicated to ODDs, we also line out the properties of the reactants which promote reaction outcomes other than the "default" hydroxylation. It turns out that the reaction selectivity depends on a delicate balance of interactions between the components of the investigated system.

Structural Insight into the Catalytic Mechanism of the Endoperoxide Synthase FtmOx1

The 2-oxoglutarate (2OG)-dependent non-heme enzyme FtmOx1 catalyzes the endoperoxide biosynthesis of verruculogen. Although several mechanistic studies have been carried out, the catalytic mechanism of FtmOx1 is not well determined owing to the lack of a reliable complex structure of FtmOx1 with fumitremorgin B. Herein we provide the X-ray crystal structure of the ternary complex FtmOx1⋅2OG⋅fumitremorgin B at a resolution of 1.22 Å. Our structures show that the binding of fumitremorgin B induces significant compression of the active pocket and that Y68 is in close proximity to C26 of the substrate. Further MD simulation and QM/MM calculations support a CarC-like mechanism, in which Y68 acts as the H atom donor for quenching the C26-centered substrate radical. Our results are consistent with all available experimental data and highlight the importance of accurate complex structures in the mechanistic study of enzymatic catalysis.

Following Nature's Footprint: Mimicking the High-Valent Heme-Oxo Mediated Indole Monooxygenation Reaction Landscape of Heme Enzymes

Pathways for direct conversion of indoles to oxindoles have accumulated considerable interest in recent years due to their significance in the clear comprehension of various pathogenic processes in humans and the multipotent therapeutic value of oxindole pharmacophores. Heme enzymes are predominantly responsible for this conversion in biology and are thought to proceed with a compound-I active oxidant. These heme-enzyme-mediated indole monooxygenation pathways are rapidly emerging therapeutic targets; however, a clear mechanistic understanding is still lacking. Additionally, such knowledge holds promise in the rational design of highly specific indole monooxygenation synthetic protocols that are also cost-effective and environmentally benign. We herein report the first examples of synthetic compound-I and activated compound-II species that can effectively monooxygenate a diverse array of indoles with varied electronic and steric properties to exclusively produce the corresponding 2-oxindole products in good to excellent yields. Rigorous kinetic, thermodynamic, and mechanistic interrogations clearly illustrate an initial rate-limiting epoxidation step that takes place between the heme oxidant and indole substrate, and the resulting indole epoxide intermediate undergoes rearrangement driven by a 2,3-hydride shift on indole ring to ultimately produce 2-oxindole. The complete elucidation of the indole monooxygenation mechanism of these synthetic heme models will help reveal crucial insights into analogous biological systems, directly reinforcing drug design attempts targeting those heme enzymes. Moreover, these bioinspired model compounds are promising candidates for the future development of better synthetic protocols for the selective, efficient, and sustainable generation of 2-oxindole motifs, which are already known for a plethora of pharmacological benefits.

Local Charge Distributions, Electric Dipole Moments, and Local Electric Fields Influence Reactivity Patterns and Guide Regioselectivities in α-Ketoglutarate-Dependent Non-heme Iron Dioxygenases

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

Can the isonitrile biosynthesis enzyme ScoE assist with the biosynthesis of isonitrile groups in drug molecules? A computational study

Computational studies show that the isonitrile synthesizing enzyme ScoE can catalyse the conversion of γ-Gly substituents in substrates to isonitrile. This enables efficient isonitrile substitution into target molecules such as axisonitrile-1.