How external perturbations affect the chemoselectivity of substrate activation by cytochrome P450 OleTJE

Simultaneous production of linear α-olefins and 2,5-furandicarboxylic acid by combining two recombinant enzymes OleT-ELP and HMFO-ELP

The enzyme OleT can utilize H2O2 as the co-substrate, and this biocatalysis is an H2O2-driven enzymatic catalysis. In this work, OleT was recombinated by being fused to an elastin-like polypeptide (ELP). The recombinant enzyme OleT-ELP exhibits higher stability and resistance to H2O2 interference compared to native OleT. OleT-ELP showed improved catalytic efficiency in producing α-olefins via fatty acid decarboxylation. The recombinant 5-hydroxymethylfurfural oxidase (HMFO-ELP) catalyzes the oxidation of 5-hydroxymethylfurfural (HMF) to 2,5-furandicarboxylic acid (FDCA), generating H2O2 as a byproduct. Combining OleT-ELP with HMFO-ELP enabled simultaneous conversion of fatty acids and HMF. The in situ H2O2 generated by HMFO-ELP was transferred to OleT-ELP, enhancing catalytic efficiencies for both α-olefins and FDCA production.

Regulation of ferryl reactivity by the cytochrome P450 decarboxylase OleT

The cytochrome P450 OleT catalyzes the decarboxylation of long-chain fatty acid substrates to produce terminal alkenes using hydrogen peroxide as a co-substrate. The facile activation of peroxide to form Compound I in the first step of the reaction, and subsequent CC bond cleavage mediated by Compound II, provides a unique opportunity to visualize both ferryl intermediates using transient kinetic approaches. Analysis of the Arrhenius behavior yields activation barriers of ∼6 kcal/mol and ∼ 18 kcal/mol for the decay of Compound I and Compound II respectively. The influence of the secondary coordination sphere, probed through site-directed mutagenesis approaches, suggests that restriction of the donor-acceptor distance contributes to the reactivity of Compound I. The reactivity of Compound II was further probed using kinetic solvent isotope effect approaches, confirming that the large barrier owes to a proton-gated mechanism in the decarboxylation reaction coordinate. Hydrogen-bonding to an active-site histidine (H85) in the distal pocket plays a key role in this process.

Microbial P450 repertoire (P450ome) and its application feasibility in pharmaceutical industry, chemical industry, and environmental protection

Cytochrome P450s (CYPs) are heme-thiolated enzymes that catalyze the oxidation of C-H bonds in a regio- and stereo-selective manner. CYPs are widely present in the biological world. With the completion of more biological genome sequencing, the number and types of P450 enzymes have increased rapidly. P450 in microorganisms is easy to clone and express, rich in catalytic types, and strong in substrate adaptability, which has good application potential. Although the number of P450 enzymes found in microorganisms is huge, the function of most of the microorganism P450s has not been studied, and it contains a large number of excellent biocatalysts to be developed. This review is based on the P450 groups in microorganisms. First, it reviews the distribution of P450 groups in different microbial species, and then studies the application of microbial P450 enzymes in the pharmaceutical industry, chemical industry and environmental pollutant treatment in recent years. And focused on the application fields of P450 enzymes of different families to guide the selection of suitable P450s from the huge P450 library. In view of the current shortcomings of microbial P450 in the application process, the final solution is the most likely to assist the application of P450 enzymes in large-scale, that is, whole cell transformation combined with engineering, fusion P450 combined with immobilization technology.

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

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.

Phenylalanine Residues in the Active Site of CYP2E1 Participate in Determining the Binding Orientation and Metabolism-Dependent Genotoxicity of Aromatic Compounds

The composition of amino acids forming the active site of a CYP enzyme is impactful in its substrate selectivity. For CYP2E1, the role of PHE residues in the formation of effective binding orientations for its aromatic substrates remains unclear. In this study, molecular docking and molecular dynamics analysis were performed to reflect the interactions between PHEs in the active site of human CYP2E1 and various aromatic compounds known as its substrates. The results indicated that the orientation of 1-methylpyrene (1-MP) in the active site was highly determined by the presence of PHEs, PHE478 contributing to the binding free energy most significantly. Moreover, by building a random forest model the relationship between each of 19 molecular descriptors of polychlorinated biphenyl (PCB) compounds (from molecular docking, quantum mechanics, and physicochemical properties) and their human CYP2E1-dependent mutagenicityas established mostly in our lab, was investigated. The presence of PHEs did not appear to significantly modify the electronic or structural feature of each bound ligand (PCB), instead, the flexibility of the conformation of PHEs contributed substantially to the effective binding energy and orientation. It is supposed that PHE residues adjust their own conformation to permit a suitablly shaped cavity for holding the ligand and forming its orientation as favorable for a biochemical reaction. This study has provided some insights into the role of PHEs in guiding the interactive adaptation of the active site of human CYP2E1 for the binding and metabolism of aromatic substrates.

Comparative oxidative ability of mononuclear and dinuclear high-valent iron-oxo species towards the activation of methane: does the axial/bridge atom modulate the reactivity?

Over the years, mononuclear Fe^IVO species have been extensively studied, but the presence of dinuclear Fe^IVO species in soluble methane monooxygenase (sMMO) has inspired the development of biomimic models that could activate inert substrates such as methane. There are some successful attempts; particularly the [(Por)(m-CBA) Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- species has been reported to activate methane and yield decent catalytic turnover numbers and therefore regarded as the closest to the sMMO enzyme functional model, as no mononuclear Fe^IVO analogues could achieve this feat. In this work, we have studied a series of mono and dinuclear models using DFT and ab initio DLPNO-CCSD(T) calculations to probe the importance of nuclearity in enhancing the reactivity. We have probed the catalytic activities of four complexes: [(HO)Fe^IV(O)(Por)]^- (1), [(HO)Fe^IV(O)(Por˙⁺)] (2), μ-oxo dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-O)Fe^IV(O) (Por˙⁺)]^- (3) and N-bridged dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- (4) towards the activation of methane. Additionally, calculations were performed on the mononuclear models [(X)Fe^IV(O)(Por˙⁺)]ⁿ {X = N 4a (n = -2), NH 4b (n = -1) and NH₂4c (n = 0)} to understand the role of nuclearity in the reactivity. DFT calculations performed on species 1-4 suggest an interesting variation among them, with species 1-3 possessing an intermediate spin (S = 1) as a ground state and species 4 possessing a high-spin (S = 2) as a ground state. Furthermore, the two Fe^IV centres in species 3 and 4 are antiferromagnetically coupled, yielding a singlet state with a distinct difference in their electronic structure. On the other hand, species 2 exhibits a ferromagnetic coupling between the Fe^IV and the Por˙⁺ moiety. Our calculations suggest that the higher barriers for the C-H bond activation of methane and the rebound step for species 1 and 3 are very high in energy, rendering them unreactive towards methane, while species 2 and 4 have lower barriers, suggesting their reactivity towards methane. Studies on the system reveal that model 4a has multiple FeN bonds facilitating greater reactivity, whereas the other two models have longer Fe-N bonds and less radical character with steeper barriers. Strong electronic cooperativity is found to be facilitated by the bridging nitride atom, and this cooperativity is suppressed by substituents such as oxygen, rendering them inactive. Thus, our study unravels that apart from enhancing the nuclearity, bridging atoms that facilitate strong cooperation between the metals are required to activate very inert substrates such as methane, and our results are broadly in agreement with earlier experimental findings.

Local Electric Fields Dictate Function: The Different Product Selectivities Observed for Fatty Acid Oxidation by Two Deceptively Very Similar P450-Peroxygenases OleT and BSβ

Cytochrome P450 peroxygenases use hydrogen peroxide to hydroxylate long-chain fatty acids by bypassing the use of O₂ and a redox partner. Among the peroxygenases, P450_OleT uniquely performs decarboxylation of fatty acids and production of terminal olefins. This route taken by P450_OleT is intriguing, and its importance is augmented by the practical importance of olefin production. As such, this mechanistic choice merits elucidation. To address this puzzle, we use hybrid QM/MM calculations and MD simulations for the OleT enzyme as well as for the structurally analogous enzyme, P450_BSβ. The study of P450_OleT reveals that the protonated His85 in the wild-type P450_OleT plays a crucial role in steering decarboxylation activity by stabilizing the corresponding hydroxoiron(IV) intermediate (Cpd II). In contrast, for P450_BSβ in which Q85 replaces H85, the respective Cpd II species is unstable and it reacts readily with the substrate radical by rebound, producing hydroxylation products. As shown, this single-site difference creates in P450_OleT a local electric field (LEF), which is significantly higher than that in P450_BSβ. In turn, these LEF differences are responsible for the different stabilities of the respective Cpd II/radical intermediates and hence for different functions of the two enzymes. P450_BSβ uses the common rebound mechanism and leads to hydroxylation, whereas P450_OleT proceeds via decarboxylation and generates terminal olefins. Olefin production projects the power of a single residue to alter the LEF and the enzyme's function.

Unexpected Reactions of α,β-Unsaturated Fatty Acids Provide Insight into the Mechanisms of CYP152 Peroxygenases

CYP152 peroxygenases catalyze decarboxylation and hydroxylation of fatty acids using H₂ O₂ as cofactor. To understand the molecular basis for the chemo- and regioselectivity of these unique P450 enzymes, we analyze the activities of three CYP152 peroxygenases (OleT_JE , P450_SPα , P450_BSβ ) towards cis- and trans-dodecenoic acids as substrate probes. The unexpected 6S-hydroxylation of the trans-isomer and 4R-hydroxylation of the cis-isomer by OleT_JE , and molecular docking results suggest that the unprecedented selectivity is due to OleT_JE 's preference of C2-C3 cis-configuration. In addition to the common epoxide products, undecanal is the unexpected major product of P450_SPα and P450_BSβ regardless of the cis/trans-configuration of substrates. The combined H₂ ¹⁸ O₂ tracing experiments, MD simulations, and QM/MM calculations unravel an unusual mechanism for Compound I-mediated aldehyde formation in which the active site water derived from H₂ O₂ activation is involved in the generation of a four-membered ring lactone intermediate. These findings provide new insights into the unusual mechanisms of CYP152 peroxygenases.

Product Distributions of Cytochrome P450 OleTJE with Phenyl-Substituted Fatty Acids: A Computational Study

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

Density Functional Theory Study into the Reaction Mechanism of Isonitrile Biosynthesis by the Nonheme Iron Enzyme ScoE

The nonheme iron enzyme ScoE catalyzes the biosynthesis of an isonitrile substituent in a peptide chain. To understand details of the reaction mechanism we created a large active site cluster model of 212 atoms that contains substrate, the active oxidant and the first- and second-coordination sphere of the protein and solvent. Several possible reaction mechanisms were tested and it is shown that isonitrile can only be formed through two consecutive catalytic cycles that both use one molecule of dioxygen and α-ketoglutarate. In both cycles the active species is an iron(IV)-oxo species that in the first reaction cycle reacts through two consecutive hydrogen atom abstraction steps: first from the N–H group and thereafter from the C–H group to desaturate the NH-CH2 bond. The alternative ordering of hydrogen atom abstraction steps was also tested but found to be higher in energy. Moreover, the electronic configurations along that pathway implicate an initial hydride transfer followed by proton transfer. We highlight an active site Lys residue that is shown to donate charge in the transition states and influences the relative barrier heights and bifurcation pathways. A second catalytic cycle of the reaction of iron(IV)-oxo with desaturated substrate starts with hydrogen atom abstraction followed by decarboxylation to give isonitrile directly. The catalytic cycle is completed with a proton transfer to iron(II)-hydroxo to generate the iron(II)-water resting state. The work is compared with experimental observation and previous computational studies on this system and put in a larger perspective of nonheme iron chemistry.

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

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