Role of Perferryl–Oxo Oxidant in Alkane Hydroxylation Catalyzed by Cytochrome P450: A Hybrid Density Functional Study

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.

The Nature of the Chemical Bonds of High-Valent Transition-Metal Oxo (M=O) and Peroxo (MOO) Compounds: A Historical Perspective of the Metal Oxyl-Radical Character by the Classical to Quantum Computations

This review article describes a historical perspective of elucidation of the nature of the chemical bonds of the high-valent transition metal oxo (M=O) and peroxo (M-O-O) compounds in chemistry and biology. The basic concepts and theoretical backgrounds of the broken-symmetry (BS) method are revisited to explain orbital symmetry conservation and orbital symmetry breaking for the theoretical characterization of four different mechanisms of chemical reactions. Beyond BS methods using the natural orbitals (UNO) of the BS solutions, such as UNO CI (CC), are also revisited for the elucidation of the scope and applicability of the BS methods. Several chemical indices have been derived as the conceptual bridges between the BS and beyond BS methods. The BS molecular orbital models have been employed to explain the metal oxyl-radical character of the M=O and M-O-O bonds, which respond to their radical reactivity. The isolobal and isospin analogy between carbonyl oxide R2C-O-O and metal peroxide LFe-O-O has been applied to understand and explain the chameleonic chemical reactivity of these compounds. The isolobal and isospin analogy among Fe=O, O=O, and O have also provided the triplet atomic oxygen (3O) model for non-heme Fe(IV)=O species with strong radical reactivity. The chameleonic reactivity of the compounds I (Cpd I) and II (Cpd II) is also explained by this analogy. The early proposals obtained by these theoretical models have been examined based on recent computational results by hybrid DFT (UHDFT), DLPNO CCSD(T0), CASPT2, and UNO CI (CC) methods and quantum computing (QC).

Genome-Wide Identification and Expression Profiling of Cytochrome P450 Monooxygenase Superfamily in Foxtail Millet

The cytochrome P450 monooxygenases (CYP450) are the largest enzyme family in plant metabolism and widely involved in the biosynthesis of primary and secondary metabolites. Foxtail millet (Setaria italica (L.) P. Beauv) can respond to abiotic stress through a highly complex polygene regulatory network, in which the SiCYP450 family is also involved. Although the CYP450 superfamily has been systematically studied in a few species, the research on the CYP450 superfamily in foxtail millet has not been completed. In this study, three hundred and thirty-one SiCYP450 genes were identified in the foxtail millet genome by bioinformatics methods, which were divided into four groups, including forty-six subgroups. One hundred and sixteen genes were distributed in thirty-three tandem duplicated gene clusters. Chromosome mapping showed that SiCYP450 was distributed on seven chromosomes. In the SiCYP450 family of foxtail millet, 20 conserved motifs were identified. Cis-acting elements in the promoter region of SiCYP450 genes showed that hormone response elements were found in all SiCYP450 genes. Of the three hundred and thirty-one SiCYP450 genes, nine genes were colinear with the Arabidopsis thaliana genes. Two hundred SiCYP450 genes were colinear with the Setaria viridis genes, including two hundred and forty-five gene duplication events. The expression profiles of SiCYP450 genes in different organs and developmental stages showed that SiCYP450 was preferentially expressed in specific tissues, and many tissue-specific genes were identified, such as SiCYP75B6, SiCYP96A7, SiCYP71A55, SiCYP71A61, and SiCYP71A62 in the root, SiCYP78A1 and SiCYP94D9 in leaves, and SiCYP78A6 in the ear. The RT-PCR data showed that SiCYP450 could respond to abiotic stresses, ABA, and herbicides in foxtail millet. Among them, the expression levels of SiCYP709B4, SiCYP71A11, SiCYP71A14, SiCYP78A1, SiCYP94C3, and SiCYP94C4 were significantly increased under the treatment of mesotrione, florasulam, nicosulfuron, fluroxypyr, and sethoxydim, indicating that the same gene might respond to multiple herbicides. The results of this study will help reveal the biological functions of the SiCYP450 family in development regulation and stress response and provide a basis for molecular breeding of foxtail millet.

Mechanism of Melatonin Metabolism by CYP1A1: What Determines the Bifurcation Pathways of Hydroxylation versus Deformylation?

Melatonin, a widely applied cosmetic active ingredient, has a variety of uses as a skin protector through antioxidant and anti-inflammatory functions as well as giving the body UV-induced defenses and immune system support. In the body, melatonin is synthesized from a tryptophan amino acid in a cascade of reactions, but as melatonin is toxic at high concentrations, it is metabolized in the human skin by the cytochrome P450 enzymes. The P450s are diverse heme-based mono-oxygenases that catalyze oxygen atom-transfer processes that trigger metabolism and detoxification reactions in the body. In the catalytic cycle of the P450s, a short-lived high-valent iron(IV)-oxo heme cation radical is formed that has been proposed to be the active oxidant. How and why it activates melatonin in the human body and what the origin of the product distributions is, are unknown. This encouraged us to do a detailed computational study on a typical human P450 isozyme, namely CYP1A1. We initially did a series of molecular dynamics simulations with substrate docked into several orientations. These simulations reveal a number of stable substrate-bound positions in the active site, which may lead to differences in substrate activation channels. Using tunneling analysis on the full protein structures, we show that two of the four binding conformations lead to open substrate-binding pockets. As a result, in these open pockets, the substrate is not tightly bound and can escape back into the solution. In the closed conformations, in contrast, the substrate is mainly oriented with the methoxy group pointing toward the heme, although under a different angle. We then created large quantum cluster models of the enzyme and focused on the chemical reaction mechanisms for melatonin activation, leading to competitive O-demethylation and C⁶-aromatic hydroxylation pathways. The calculations show that active site positioning determines the product distributions, but the bond that is activated is not necessarily closest to the heme in the enzyme-substrate complex. As such, the docking and molecular dynamics positioning of the substrate versus oxidant can give misleading predictions on product distributions. In particular, in quantum mechanics cluster model I, we observe that through a tight hydrogen bonding network, a preferential 6-hydroxylation of melatonin is obtained. However, O-demethylation becomes possible in alternative substrate-binding orientations that have the C⁶-aromatic ring position shielded. Finally, we investigated enzymatic and non-enzymatic O-demethylation processes and show that the hydrogen bonding network in the substrate-binding pocket can assist and perform this step prior to product release from the enzyme.

Low-Lying Electromeric States in Chloro-Ligated Iron(IV)-Oxo Porphyrin as a Model for Compound I, Studied with Second-Order Perturbation Theory Based on Density Matrix Renormalization Group

Employing second-order perturbation theory based on the density matrix renormalization group (DMRG-CASPT2), this work aims at providing a quantitative description of the spin state energetics of a chloro-ligated iron(IV)-oxo porphyrin as a model for the cytochromes P450 active species, also known as compound I (Cpd I). We explored DMRG-CASPT2 to its full extent with an extensive active space (up to 31 active orbitals) as well as a large number of renormalized states m (up to 10000). Different flavors of DMRG-CASPT2, using either the costly exact 4-particle reduced density matrix (4-RDM) or the cheaper cumulant approximated 4-RDM (cu(4)), were analyzed. All flavors essentially converge to similar relative energies between different spin states. Including a correction for the protein environment, we found a quartet Fe^IVO ground state and, more importantly, a thermally accessible doublet Fe^VO excited state that might directly contribute to the reactivity of this iron-oxo species. Our results also showed that cheaper approaches, such as CASPT2 based on a smaller active space or the cumulant approximation DMRG-cu(4)-CASPT2, are capable of accurately describing the spin state energetics of this species.

Role of the Proximal Cysteine Hydrogen Bonding Interaction in Cytochrome P450 2B4 Studied by Cryoreduction, Electron Paramagnetic Resonance, and Electron-Nuclear Double Resonance Spectroscopy

Crystallographic studies have shown that the F429H mutation of cytochrome P450 2B4 introduces an H-bond between His429 and the proximal thiolate ligand, Cys436, without altering the protein fold but sharply decreases the enzymatic activity and stabilizes the oxyferrous P450 2B4 complex. To characterize the influence of this hydrogen bond on the states of the catalytic cycle, we have used radiolytic cryoreduction combined with electron paramagnetic resonance (EPR) and (electron-nuclear double resonance (ENDOR) spectroscopy to study and compare their characteristics for wild-type (WT) P450 2B4 and the F429H mutant. (i) The addition of an H-bond to the axial Cys436 thiolate significantly changes the EPR signals of both low-spin and high-spin heme-iron(III) and the hyperfine couplings of the heme-pyrrole (14)N but has relatively little effect on the (1)H ENDOR spectra of the water ligand in the six-coordinate low-spin ferriheme state. These changes indicate that the H-bond introduced between His and the proximal cysteine decreases the extent of S → Fe electron donation and weakens the Fe(III)-S bond. (ii) The added H-bond changes the primary product of cryoreduction of the Fe(II) enzyme, which is trapped in the conformation of the parent Fe(II) state. In the wild-type enzyme, the added electron localizes on the porphyrin, generating an S = (3)/2 state with the anion radical exchange-coupled to the Fe(II). In the mutant, it localizes on the iron, generating an S = (1)/2 Fe(I) state. (iii) The additional H-bond has little effect on g values and (1)H-(14)N hyperfine couplings of the cryogenerated, ferric hydroperoxo intermediate but noticeably slows its decay during cryoannealing. (iv) In both the WT and the mutant enzyme, this decay shows a significant solvent kinetic isotope effect, indicating that the decay reflects a proton-assisted conversion to Compound I (Cpd I). (v) We confirm that Cpd I formed during the annealing of the cryogenerated hydroperoxy intermediate and that it is the active hydroxylating species in both WT P450 2B4 and the F429H mutant. (vi) Our data also indicate that the added H-bond of the mutation diminishes the reactivity of Cpd I.

Monooxygenase, peroxidase and peroxygenase properties and reaction mechanisms of cytochrome P450 enzymes

This review examines the monooxygenase, peroxidase and peroxygenase properties and reaction mechanisms of cytochrome P450 (CYP) enzymes in bacterial, archaeal and mammalian systems. CYP enzymes catalyze monooxygenation reactions by inserting one oxygen atom from O2 into an enormous number and variety of substrates. The catalytic versatility of CYP stems from its ability to functionalize unactivated carbon-hydrogen (C-H) bonds of substrates through monooxygenation. The oxidative prowess of CYP in catalyzing monooxygenation reactions is attributed primarily to a porphyrin π radical ferryl intermediate known as Compound I (CpdI) (Por•+FeIV=O), or its ferryl radical resonance form (FeIV-O•). CYP-mediated hydroxylations occur via a consensus H atom abstraction/oxygen rebound mechanism involving an initial abstraction by CpdI of a H atom from the substrate, generating a highly-reactive protonated Compound II (CpdII) intermediate (FeIV-OH) and a carbon-centered alkyl radical that rebounds onto the ferryl hydroxyl moiety to yield the hydroxylated substrate. CYP enzymes utilize hydroperoxides, peracids, perborate, percarbonate, periodate, chlorite, iodosobenzene and N-oxides as surrogate oxygen atom donors to oxygenate substrates via the shunt pathway in the absence of NAD(P)H/O2 and reduction-oxidation (redox) auxiliary proteins. It has been difficult to isolate the historically elusive CpdI intermediate in the native NAD(P)H/O2-supported monooxygenase pathway and to determine its precise electronic structure and kinetic and physicochemical properties because of its high reactivity, unstable nature (t½~2 ms) and short life cycle, prompting suggestions for participation in monooxygenation reactions of alternative CYP iron-oxygen intermediates such as the ferric-peroxo anion species (FeIII-OO-), ferric-hydroperoxo species (FeIII-OOH) and FeIII-(H2O2) complex.

Applications of density functional theory to iron-containing molecules of bioinorganic interest

The past decades have seen an explosive growth in the application of density functional theory (DFT) methods to molecular systems that are of interest in a variety of scientific fields. Owing to its balanced accuracy and efficiency, DFT plays particularly useful roles in the theoretical investigation of large molecules. Even for biological molecules such as proteins, DFT finds application in the form of, e.g., hybrid quantum mechanics and molecular mechanics (QM/MM), in which DFT may be used as a QM method to describe a higher prioritized region in the system, while a MM force field may be used to describe remaining atoms. Iron-containing molecules are particularly important targets of DFT calculations. From the viewpoint of chemistry, this is mainly because iron is abundant on earth, iron plays powerful (and often enigmatic) roles in enzyme catalysis, and iron thus has the great potential for biomimetic catalysis of chemically difficult transformations. In this paper, we present a brief overview of several recent applications of DFT to iron-containing non-heme synthetic complexes, heme-type cytochrome P450 enzymes, and non-heme iron enzymes, all of which are of particular interest in the field of bioinorganic chemistry. Emphasis will be placed on our own work.