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.

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

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

An Active Site Tyr Residue Guides the Regioselectivity of Lysine Hydroxylation by Nonheme Iron Lysine-4-hydroxylase Enzymes through Proton-Coupled Electron Transfer

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 C4-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 C4-H bond of the substrate and, therefore, the electrostatic and electric field perturbations of the protein assist the enzyme in creating C4-H hydroxylation selectivity. Furthermore, an active site Tyr233 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 Tyr233 phenol loses a proton to the nearby Asp179 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 C4-H bond of the substrate and guides the selectivity to the C4-hydroxylation of the substrate.

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.

%VBur index and steric maps: from predictive catalysis to machine learning

Steric indices are parameters used in chemistry to describe the spatial arrangement of atoms or groups of atoms in molecules. They are important in determining the reactivity, stability, and physical properties of chemical compounds. One commonly used steric index is the steric hindrance, which refers to the obstruction or hindrance of movement in a molecule caused by bulky substituents or functional groups. Steric hindrance can affect the reactivity of a molecule by altering the accessibility of its reactive sites and influencing the geometry of its transition states. Notably, the Tolman cone angle and %VBur are prominent among these indices. Actually, steric effects can also be described using the concept of steric bulk, which refers to the space occupied by a molecule or functional group. Steric bulk can affect the solubility, melting point, boiling point, and viscosity of a substance. Even though electronic indices are more widely used, they have certain drawbacks that might shift preferences towards others. They present a higher computational cost, and often, the weight of electronics in correlation with chemical properties, e.g. binding energies, falls short in comparison to %VBur. However, it is worth noting that this may be because the steric index inherently captures part of the electronic content. Overall, steric indices play an important role in understanding the behaviour of chemical compounds and can be used to predict their reactivity, stability, and physical properties. Predictive chemistry is an approach to chemical research that uses computational methods to anticipate the properties and behaviour of these compounds and reactions, facilitating the design of new compounds and reactivities. Within this domain, predictive catalysis specifically targets the prediction of the performance and behaviour of catalysts. Ultimately, the goal is to identify new catalysts with optimal properties, leading to chemical processes that are both more efficient and sustainable. In this framework, %VBur can be a key metric for deepening our understanding of catalysis, emphasizing predictive catalysis and sustainability. Those latter concepts are needed to direct our efforts toward identifying the optimal catalyst for any reaction, minimizing waste, and reducing experimental efforts while maximizing the efficacy of the computational methods.

Evidence of a Distinctive Enantioselective Binding Mode for the Photoinduced Radical Cyclization of α-Chloroamides in Ene-Reductases

We demonstrate here through molecular simulations and mutational studies the origin of the enantioselectivity in the photoinduced radical cyclization of α-chloroacetamides catalyzed by ene-reductases, in particular the Gluconobacter oxidans ene-reductase and the Old Yellow Enzyme 1, which show opposite enantioselectivity. Our results reveal that neither the π-facial selectivity model nor a protein-induced selective stabilization of the transition states is able to explain the enantioselectivity of the radical cyclization in the studied flavoenzymes. We propose a new enantioinduction scenario according to which enantioselectivity is indeed controlled by transition-state stability; however, the relative stability of the prochiral transition states is not determined by direct interaction with the protein but is rather dependent on an inherent degree of freedom within the substrate itself. This intrinsic degree of freedom, distinct from the traditional π-facial exposure mode, can be controlled by the substrate conformational selection upon binding to the enzyme.

Computation-Aided Engineering of Cytochrome P450 for the Production of Pravastatin

CYP105AS1 is a cytochrome P450 from Amycolatopsis orientalis that catalyzes monooxygenation of compactin to 6-epi-pravastatin. For fermentative production of the cholesterol-lowering drug pravastatin, the stereoselectivity of the enzyme needs to be inverted, which has been partially achieved by error-prone PCR mutagenesis and screening. In the current study, we report further optimization of the stereoselectivity by a computationally aided approach. Using the CoupledMoves protocol of Rosetta, a virtual library of mutants was designed to bind compactin in a pro-pravastatin orientation. By examining the frequency of occurrence of beneficial substitutions and rational inspection of their interactions, a small set of eight mutants was predicted to show the desired selectivity and these variants were tested experimentally. The best CYP105AS1 variant gave >99% stereoselective hydroxylation of compactin to pravastatin, with complete elimination of the unwanted 6-epi-pravastatin diastereomer. The enzyme-substrate complexes were also examined by ultrashort molecular dynamics simulations of 50 × 100 ps and 5 × 22 ns, which revealed that the frequency of occurrence of near-attack conformations agreed with the experimentally observed stereoselectivity. These results show that a combination of computational methods and rational inspection could improve CYP105AS1 stereoselectivity beyond what was obtained by directed evolution. Moreover, the work lays out a general in silico framework for specificity engineering of enzymes of known structure.

Unraveling the mechanism of enantio-controlling switches of an alcohol dehydrogenase toward sterically small ketone

Efficient synthesis of chiral compounds under mild conditions is highly desirable in the chemical and pharmaceutical communities, but it often faces difficulties. Although various enzymes have been harnessed as biocatalysts...

Computational studies of DNA base repair mechanisms by nonheme iron dioxygenases: selective epoxidation and hydroxylation pathways

DNA base repair mechanisms of alkylated DNA bases is an important reaction in chemical biology and particularly in the human body. It is typically catalyzed by an α-ketoglutarate-dependent nonheme iron dioxygenase named the AlkB repair enzyme. In this work we report a detailed computational study into the structure and reactivity of AlkB repair enzymes with alkylated DNA bases. In particular, we investigate the aliphatic hydroxylation and C[double bond, length as m-dash]C epoxidation mechanisms of alkylated DNA bases by a high-valent iron(iv)-oxo intermediate. Our computational studies use quantum mechanics/molecular mechanics methods on full enzymatic structures as well as cluster models on active site systems. The work shows that the iron(iv)-oxo species is rapidly formed after dioxygen binding to an iron(ii) center and passes a bicyclic ring structure as intermediate. Subsequent cluster models explore the mechanism of substrate hydroxylation and epoxidation of alkylated DNA bases. The work shows low energy barriers for substrate activation and consequently energetically feasible pathways are predicted. Overall, the work shows that a high-valent iron(iv)-oxo species can efficiently dealkylate alkylated DNA bases and return them into their original form.

How Does Replacement of the Axial Histidine Ligand in Cytochrome c Peroxidase by Nδ-Methyl Histidine Affect Its Properties and Functions? A Computational Study

Heme peroxidases have important functions in nature related to the detoxification of H₂O₂. They generally undergo a catalytic cycle where, in the first stage, the iron(III)-heme-H₂O₂ complex is converted into an iron(IV)-oxo-heme cation radical species called Compound I. Cytochrome c peroxidase Compound I has a unique electronic configuration among heme enzymes where a metal-based biradical is coupled to a protein radical on a nearby Trp residue. Recent work using the engineered N_δ-methyl histidine-ligated cytochrome c peroxidase highlighted changes in spectroscopic and catalytic properties upon axial ligand substitution. To understand the axial ligand effect on structure and reactivity of peroxidases and their axially N_δ-methyl histidine engineered forms, we did a computational study. We created active site cluster models of various sizes as mimics of horseradish peroxidase and cytochrome c peroxidase Compound I. Subsequently, we performed density functional theory studies on the structure and reactivity of these complexes with a model substrate (styrene). Thus, the work shows that the N_δ-methyl histidine group has little effect on the electronic configuration and structure of Compound I and little changes in bond lengths and the same orbital occupation is obtained. However, the N_δ-methyl histidine modification impacts electron transfer processes due to a change in the reduction potential and thereby influences reactivity patterns for oxygen atom transfer. As such, the substitution of the axial histidine by N_δ-methyl histidine in peroxidases slows down oxygen atom transfer to substrates and makes Compound I a weaker oxidant. These studies are in line with experimental work on N_δ-methyl histidine-ligated cytochrome c peroxidases and highlight how the hydrogen bonding network in the second coordination sphere has a major impact on the function and properties of the enzyme.

Rapid prototyping of microbial production strains for the biomanufacture of potential materials monomers

Bio-based production of industrial chemicals using synthetic biology can provide alternative green routes from renewable resources, allowing for cleaner production processes. To efficiently produce chemicals on-demand through microbial strain engineering, biomanufacturing foundries have developed automated pipelines that are largely compound agnostic in their time to delivery. Here we benchmark the capabilities of a biomanufacturing pipeline to enable rapid prototyping of microbial cell factories for the production of chemically diverse industrially relevant material building blocks. Over 85 days the pipeline was able to produce 17 potential material monomers and key intermediates by combining 160 genetic parts into 115 unique biosynthetic pathways. To explore the scale-up potential of our prototype production strains, we optimized the enantioselective production of mandelic acid and hydroxymandelic acid, achieving gram-scale production in fed-batch fermenters. The high success rate in the rapid design and prototyping of microbially-produced material building blocks reveals the potential role of biofoundries in leading the transition to sustainable materials production.

Generation of amine dehydrogenases with increased catalytic performance and substrate scope from ε-deaminating L-Lysine dehydrogenase

Amine dehydrogenases (AmDHs) catalyse the conversion of ketones into enantiomerically pure amines at the sole expense of ammonia and hydride source. Guided by structural information from computational models, we create AmDHs that can convert pharmaceutically relevant aromatic ketones with conversions up to quantitative and perfect chemical and optical purities. These AmDHs are created from an unconventional enzyme scaffold that apparently does not operate any asymmetric transformation in its natural reaction. Additionally, the best variant (LE-AmDH-v1) displays a unique substrate-dependent switch of enantioselectivity, affording S- or R-configured amine products with up to >99.9% enantiomeric excess. These findings are explained by in silico studies. LE-AmDH-v1 is highly thermostable (Tm of 69 °C), retains almost entirely its catalytic activity upon incubation up to 50 °C for several days, and operates preferentially at 50 °C and pH 9.0. This study also demonstrates that product inhibition can be a critical factor in AmDH-catalysed reductive amination.

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

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

Does Substrate Positioning Affect the Selectivity and Reactivity in the Hectochlorin Biosynthesis Halogenase?

In this work we present the first computational study on the hectochlorin biosynthesis enzyme HctB, which is a unique three-domain halogenase that activates non-amino acid moieties tethered to an acyl-carrier, and as such may have biotechnological relevance beyond other halogenases. We use a combination of small cluster models and full enzyme structures calculated with quantum mechanics/molecular mechanics methods. Our work reveals that the reaction is initiated with a rate-determining hydrogen atom abstraction from substrate by an iron (IV)-oxo species, which creates an iron (III)-hydroxo intermediate. In a subsequent step the reaction can bifurcate to either halogenation or hydroxylation of substrate, but substrate binding and positioning drives the reaction to optimal substrate halogenation. Furthermore, several key residues in the protein have been identified for their involvement in charge-dipole interactions and induced electric field effects. In particular, two charged second coordination sphere amino acid residues (Glu223 and Arg245) appear to influence the charge density on the Cl ligand and push the mechanism toward halogenation. Our studies, therefore, conclude that nonheme iron halogenases have a chemical structure that induces an electric field on the active site that affects the halide and iron charge distributions and enable efficient halogenation. As such, HctB is intricately designed for a substrate halogenation and operates distinctly different from other nonheme iron halogenases.

Bacterial bifunctional chorismate mutase-prephenate dehydratase PheA increases flux into the yeast phenylalanine pathway and improves mandelic acid production

Mandelic acid is an important aromatic fine chemical and is currently mainly produced via chemical synthesis. Recently, mandelic acid production was achieved by microbial fermentations using engineered Escherichia coli and Saccharomyces cerevisiae expressing heterologous hydroxymandelate synthases (hmaS). The best-performing strains carried a deletion of the gene encoding the first enzyme of the tyrosine biosynthetic pathway and therefore were auxotrophic for tyrosine. This was necessary to avoid formation of the competing intermediate hydroxyphenylpyruvate, the preferred substrate for HmaS, which would have resulted in the predominant production of hydroxymandelic acid. However, feeding tyrosine to the medium would increase fermentation costs. In order to engineer a tyrosine prototrophic mandelic acid-producing S. cerevisiae strain, we tested three strategies: (1) rational engineering of the HmaS active site for reduced binding of hydroxyphenylpyruvate, (2) compartmentalization of the mandelic acid biosynthesis pathway by relocating HmaS together with the two upstream enzymes chorismate mutase Aro7 and prephenate dehydratase Pha2 into mitochondria or peroxisomes, and (3) utilizing a feedback-resistant version of the bifunctional E. coli enzyme PheA (PheA^fbr) in an aro7 deletion strain. PheA has both chorismate mutase and prephenate dehydratase activity. Whereas the enzyme engineering approaches were only successful in respect to reducing the preference of HmaS for hydroxyphenylpyruvate but not in increasing mandelic acid titers, we could show that strategies (2) and (3) significantly reduced hydroxymandelic acid production in favor of increased mandelic acid production, without causing tyrosine auxotrophy. Using the bifunctional enzyme PheA^fbr turned out to be the most promising strategy, and mandelic acid production could be increased 12-fold, yielding titers up to 120 mg/L. Moreover, our results indicate that utilizing PheA^fbr also shows promise for other industrial applications with S. cerevisiae that depend on a strong flux into the phenylalanine biosynthetic pathway.

•

Increased mandelic acid production in tyrosine prototrophic S. cerevisiae.

•

Bifunctional E. coli enzyme PheA increases flux into yeast phenylalanine branch.

•

PheA allows mandelic acid production in prototrophic S. cerevisiae.

•

Compartmentalized mandelic acid production in yeast mitochondria/peroxisomes.

Inhibition of the Exoglucanase Cel7A by a Douglas-Fir-Condensed Tannin

Douglas-fir forestry residues are a potential feedstock for saccharification-based biofuels, and condensed tannins are expected to make up ∼3% of the dry mass of this feedstock. Condensed tannins are well-known for their ability to interact with proteins and can bind and inhibit cellulase enzymes used in saccharification. In this study, we use molecular docking and classical molecular dynamics simulations to investigate how a characterized condensed tannin from Douglas-fir bark binds to the exoglucanase Cel7A from Trichoderma reesei. Through looking at the "occupancy" and "residency" of specific amino acid residue-tannin interactions, we find that the binding sites are characterized by many simultaneous tannin-enzyme interactions with the strongest occurring on the catalytic module as opposed to the carbohydrate-binding module. The simulations indicate that tannin inhibition can result from binding at or near the catalytic tunnel's entrance and exit. The analyzed tannin further prefers to bind to loops around the catalytic region and has affinity for aromatic and charged amino acid residues. These insights provide direction for the rational design of tannin-resistant cellulases.

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

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

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

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

Mechanistic Insight on the Activity and Substrate Selectivity of Nonheme Iron Dioxygenases

Nonheme iron dioxygenases catalyze vital reactions for human health particularly related to aging processes. They are involved in the biosynthesis of amino acids, but also the biodegradation of toxic compounds. Typically they react with their substrate(s) through oxygen atom transfer, although often with the assistance of a co-substrate like α-ketoglutarate that is converted to succinate and CO₂ . Many reaction processes catalyzed by the nonheme iron dioxygenases are stereoselective or regiospecific and hence understanding the mechanism and protein involvement in the selectivity is important for the design of biotechnological applications of these enzymes. To this end, I will review recent work of our group on nonheme iron dioxygenases and include background information on their general structure and catalytic cycle. Examples of stereoselective and regiospecific reaction mechanisms we elucidated are for the AlkB repair enzyme, prolyl-4-hydroxylase and the ergothioneine biosynthesis enzyme. Finally, I cover an example where we bioengineered S-p-hydroxymandelate synthase into the R-p-hydroxymandelate synthase.

How Are Substrate Binding and Catalysis Affected by Mutating Glu127 and Arg161 in Prolyl-4-hydroxylase? A QM/MM and MD Study

Prolyl-4-hydroxylase is a vital enzyme for human physiology involved in the biosynthesis of 4-hydroxyproline, an essential component for collagen formation. The enzyme performs a unique stereo- and regioselective hydroxylation at the C4 position of proline despite the fact that the C5 hydrogen atoms should be thermodynamically easier to abstract. To gain insight into the mechanism and find the origin of this regioselectivity, we have done a quantum mechanics/molecular mechanics (QM/MM) study on wildtype and mutant structures. In a previous study (Timmins et al., 2017) we identified several active site residues critical for substrate binding and positioning. In particular, the Glu127 and Arg161 were shown to form multiple hydrogen bonding and ion-dipole interactions with substrate and could thereby affect the regio- and stereoselectivity of the reaction. In this work, we decided to test that hypothesis and report a QM/MM and molecular dynamics (MD) study on prolyl-4-hydroxylase and several active site mutants where Glu127 or Arg161 are mutated for Asp, Gln, or Lys. Thus, the R161D and R161Q mutants give very high barriers for hydrogen atom abstraction from any proline C-H bond and therefore will be inactive. The R161K mutant, by contrast, sees the regio- and stereoselectivity of the reaction change but still is expected to hydroxylate proline at room temperature. By contrast, the Glu127 mutants E127D and E127Q show possible changes in regioselectivity with the former being more probable to react compared to the latter.

Structural and Biophysical Characterization of the Mycobacterium tuberculosis Protein Rv0577, a Protein Associated with Neutral Red Staining of Virulent Tuberculosis Strains and Homologue of the Streptomyces coelicolor Protein KbpA

Mycobacterium tuberculosis protein Rv0577 is a prominent antigen in tuberculosis patients, the component responsible for neutral red staining of virulent strains of M. tuberculosis, a putative component in a methylglyoxal detoxification pathway, and an agonist of toll-like receptor 2. It also has an amino acid sequence that is 36% identical to that of Streptomyces coelicolor AfsK-binding protein A (KbpA), a component in the complex secondary metabolite pathways in the Streptomyces genus. To gain insight into the biological function of Rv0577 and the family of KpbA kinase regulators, the crystal structure for Rv0577 was determined to a resolution of 1.75 Å, binding properties with neutral red and deoxyadenosine were surveyed, backbone dynamics were measured, and thermal stability was assayed by circular dichroism spectroscopy. The protein is composed of four approximate repeats with a βαβββ topology arranged radially in consecutive pairs to form two continuous eight-strand β-sheets capped on both ends with an α-helix. The two β-sheets intersect in the center at roughly a right angle and form two asymmetric deep "saddles" that may serve to bind ligands. Nuclear magnetic resonance chemical shift perturbation experiments show that neutral red and deoxyadenosine bind to Rv0577. Binding to deoxyadenosine is weaker with an estimated dissociation constants of 4.1 ± 0.3 mM for saddle 1. Heteronuclear steady-state {¹H}-¹⁵N nuclear Overhauser effect, T₁, and T₂ values were generally uniform throughout the sequence with only a few modest pockets of differences. Circular dichroism spectroscopy characterization of the thermal stability of Rv0577 indicated irreversible unfolding upon heating with an estimated melting temperature of 56 °C.

Sulfoxide Synthase versus Cysteine Dioxygenase Reactivity in a Nonheme Iron Enzyme

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

CADEE: Computer-Aided Directed Evolution of Enzymes

The tremendous interest in enzymes as biocatalysts has led to extensive work in enzyme engineering, as well as associated methodology development. Here, a new framework for computer-aided directed evolution of enzymes (CADEE) is presented which allows a drastic reduction in the time necessary to prepare and analyze in silico semi-automated directed evolution of enzymes. A pedagogical example of the application of CADEE to a real biological system is also presented in order to illustrate the CADEE workflow.

A Systematic Account on Aromatic Hydroxylation by a Cytochrome P450 Model Compound I: A Low-Pressure Mass Spectrometry and Computational Study

Cytochrome P450 enzymes are heme-containing mono-oxygenases that mainly react through oxygen-atom transfer. Specific features of substrate and oxidant that determine the reaction rate constant for oxygen atom transfer are still poorly understood and therefore, we did a systematic gas-phase study on reactions by iron(IV)-oxo porphyrin cation radical structures with arenes. We present herein the first results obtained by using Fourier transform-ion cyclotron resonance mass spectrometry and provide rate constants and product distributions for the assayed reactions. Product distributions and kinetic isotope effect studies implicate a rate-determining aromatic hydroxylation reaction that correlates with the ionization energy of the substrate and no evidence of aliphatic hydroxylation products is observed. To further understand the details of the reaction mechanism, a computational study on a model complex was performed. These studies confirm the experimental hypothesis of dominant aromatic over aliphatic hydroxylation and show that the lack of an axial ligand affects the aliphatic pathways. Moreover, a two-parabola valence bond model is used to rationalize the rate constant and identify key properties of the oxidant and substrate that drive the reaction. In particular, the work shows that aromatic hydroxylation rates correlate with the ionization energy of the substrate as well as with the electron affinity of the oxidant.

Arene activation by a nonheme iron(III)-hydroperoxo complex: pathways leading to phenol and ketone products

Iron(III)-hydroperoxo complexes are found in various nonheme iron enzymes as catalytic cycle intermediates; however, little is known on their catalytic properties. The recent work of Banse and co-workers on a biomimetic nonheme iron(III)-hydroperoxo complex provided evidence of its involvement in reactivity with arenes. This contrasts the behavior of heme iron(III)-hydroperoxo complexes that are known to be sluggish oxidants. To gain insight into the reaction mechanism of the biomimetic iron(III)-hydroperoxo complex with arenes, we performed a computational (density functional theory) study. The calculations show that iron(III)-hydroperoxo reacts with substrates via low free energies of activation that should be accessible at room temperature. Moreover, a dominant ketone reaction product is observed as primary products rather than the thermodynamically more stable phenols. These product distributions are analyzed and the calculations show that charge interaction between the iron(III)-hydroxo group and the substrate in the intermediate state pushes the transferring proton to the meta-carbon atom of the substrate and guides the selectivity of ketone formation. These studies show that the relative ratio of ketone versus phenol as primary products can be affected by external interactions of the oxidant with the substrate. Moreover, iron(III)-hydroperoxo complexes are shown to selectively give ketone products, whereas iron(IV)-oxo complexes will react with arenes to form phenols instead.

Challenging Density Functional Theory Calculations with Hemes and Porphyrins

In this paper we review recent advances in computational chemistry and specifically focus on the chemical description of heme proteins and synthetic porphyrins that act as both mimics of natural processes and technological uses. These are challenging biochemical systems involved in electron transfer as well as biocatalysis processes. In recent years computational tools have improved considerably and now can reproduce experimental spectroscopic and reactivity studies within a reasonable error margin (several kcal·mol(-1)). This paper gives recent examples from our groups, where we investigated heme and synthetic metal-porphyrin systems. The four case studies highlight how computational modelling can correctly reproduce experimental product distributions, predicted reactivity trends and guide interpretation of electronic structures of complex systems. The case studies focus on the calculations of a variety of spectroscopic features of porphyrins and show how computational modelling gives important insight that explains the experimental spectra and can lead to the design of porphyrins with tuned properties.

Conversion of a non-heme iron-dependent sulfoxide synthase into a thiol dioxygenase by a single point mutation

EgtB from Mycobacterium thermoresistibile catalyzes O₂-dependent sulfur–carbon bond formation between the side chains of Nα-trimethyl histidine and γ-glutamyl cysteine as a central step in ergothioneine biosynthesis.

Quantum Mechanics/Molecular Mechanics Modeling of Enzymatic Processes: Caveats and Breakthroughs

Nature has developed large groups of enzymatic catalysts with the aim to transfer substrates into useful products, which enables biosystems to perform all their natural functions. As such, all biochemical processes in our body (we drink, we eat, we breath, we sleep, etc.) are governed by enzymes. One of the problems associated with research on biocatalysts is that they react so fast that details of their reaction mechanisms cannot be obtained with experimental work. In recent years, major advances in computational hardware and software have been made and now large (bio)chemical systems can be studied using accurate computational techniques. One such technique is the quantum mechanics/molecular mechanics (QM/MM) technique, which has gained major momentum in recent years. Unfortunately, it is not a black-box method that is easily applied, but requires careful set-up procedures. In this work we give an overview on the technical difficulties and caveats of QM/MM and discuss work-protocols developed in our groups for running successful QM/MM calculations.

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

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

Properties and reactivities of nonheme iron(IV)-oxo versus iron(V)-oxo: long-range electron transfer versus hydrogen atom abstraction

Recent work of Nam and co-workers [J. Yoon, S. A. Wilson, Y. K. Jang, M. S. Seo, K. Nehru, B. Hedman, K. O. Hodgson, E. Bill, E. I. Solomon and W. Nam, Angew. Chem., Int. Ed., 2009, 48, 1257] on a biomimetic iron complex implicated a mixture of iron(IV)-oxo and iron(V)-oxo intermediates but the latter could not be spectroscopically characterized, hence its involvement was postulated. To gain insight into the relative activity of these iron(IV)-oxo versus iron(V)-oxo intermediates, we have performed an extensive density functional theory (DFT) study on the chemical properties of the chemical system of Nam et al., namely [Fe(O)(BQEN)(NCCH3)](2+/3+) with BQEN = N,N'-dimethyl-N,N'-bis(8-quinolyl)ethane-1,2-diamine and their reactivity in hydrogen atom abstraction from ethylbenzene. We show that the perceived iron(V)-oxo species actually is an iron(IV)-oxo ligand cation radical, similar to cytochrome P450 compound I. Moreover, this intermediate has an extremely large electron affinity and therefore can abstract electrons from substrates readily. In our particular system, this means that prior to the hydrogen atom abstraction, an electron is abstracted to form an iron(IV)-oxo species, which subsequently abstracts a hydrogen atom from the substrate. Thus, our calculations show for the first time how some nonheme iron complexes react by long-range electron transfer and others directly via hydrogen atom abstraction. We have rationalized our results with detailed thermochemical cycles that explain the observed reactivity patterns.

Neocarzinostatin-based hybrid biocatalysts for oxidation reactions

An anionic iron(III)-porphyrin-testosterone conjugate 1-Fe has been synthesized and fully characterized. It has been further associated with a neocarzinostatin variant, NCS-3.24, to generate a new artificial metalloenzyme following the so-called 'Trojan Horse' strategy. This new 1-Fe-NCS-3.24 biocatalyst showed an interesting catalytic activity as it was found able to catalyze the chemoselective and slightly enantioselective (ee = 13%) sulfoxidation of thioanisole by H2O2. Molecular modelling studies show that a synergy between the binding of the steroid moiety and that of the porphyrin macrocycle into the protein binding site can explain the experimental results, indicating a better affinity of 1-Fe for the NCS-3.24 variant than testosterone and testosterone-hemisuccinate themselves. They also show that the Fe-porphyrin complex is sandwiched between the two subdomains of the protein providing with good complementarities. However, the artificial cofactor entirely fills the cavity and its metal ion remains widely exposed to the solvent which explains the moderate enantioselectivity observed. Some possible improvements in the "Trojan Horse" strategy for obtaining better catalysts of selective oxidations are presented.

Dioxygen activation and two consecutive oxidative decarboxylations of phenylpyruvate by nonheme iron(II) complexes: functional models of hydroxymandelate synthase (HMS) and CloR

Two mononuclear iron(ii)-phenylpyruvate complexes of monoanionic facial N3 ligands are reported to react with dioxygen to undergo two consecutive oxidative decarboxylation steps via an iron-mandelate complex mimicking the function of HMS and CloR.

Various strategies for obtaining oxidative artificial hemoproteins with a catalytic oxidative activity: from "Hemoabzymes" to "Hemozymes"?

The design of artificial hemoproteins that could lead to new biocatalysts for selective oxidation reactions using clean oxidants such as O 2 or H 2 O 2 under ecocompatible conditions constitutes a really promising challenge for a wide range of industrial applications. In vivo, such reactions are performed by heme-thiolate proteins, cytochromes P450, that catalyze the oxidation of drugs by dioxygen in the presence of electrons delivered from NADPH by cytochrome P450 reductase. Several strategies were used to design new artificial hemoproteins to mimic these enzymes, that associate synthetic metalloporphyrin derivatives to a protein that is supposed to induce a selectivity in the catalyzed reaction. A first generation of artificial hemoproteins or "hemoabzymes" was obtained by the non-covalent association of synthetic hemes such as N-methyl-mesoporphyrin IX, Fe(III) -α3β-tetra-o-carboxyphenylporphyrin or microperoxidase 8 with monoclonal antibodies raised against these cofactors. The obtained antibody-metalloporphyrin complexes displayed a peroxidase activity and some of them catalyzed the regio-selective nitration of phenols by H 2 O 2/ NO 2 and the stereo-selective oxidation of sulphides by H 2 O 2. A second generation of artificial hemoproteins or "hemozymes", was obtained by the non-covalent association of non-relevant proteins with metalloporphyrin derivatives. Several strategies were used, the most successful of which, named "host-guest" strategy involved the non-covalent incorporation of metalloporphyrin derivatives into easily affordable proteins. The artificial hemoproteins obtained were found to be able to perform efficiently the stereoselective oxidation of organic compounds such as sulphides and alkenes by H 2 O 2 and KHSO 5.