Variations on a (t)heme—novel mechanisms, redox partners and catalytic functions in the cytochrome P450 superfamily

Effect of CO binding to P450 BM3 F393 mutants on electron density distribution in the heme cofactor

Resonance Raman spectroscopy has been performed on a set of cytochrome P450 BM3 heme domains in which mutation of the highly conserved Phe393 induces significant variation in heme iron reduction potential. In previous work [Chen, Z., Ost, T.W.B., and Schelvis, J.P.M. (2004) Biochemistry 43, 1798-1808], a correlation between heme vinyl conformation and the heme iron reduction potential indicated a steric control by the protein over the distribution of electron density in the reduced heme cofactor. The current study aims to monitor changes in electron density on the ferrous heme cofactor following CO binding. In addition, ferric-NO complexes have been studied to investigate potential changes to the proximal Cys400 thiolate. We find that binding of CO to the ferrous heme domains results in a reorientation of the vinyl groups to a largely out-of-plane conformation, the extent of which correlates with the size of the residue at position 393. We conclude that FeII dπ back bonding to the CO ligand largely takes away the need for conjugation of the vinyl groups with the porphyrin ring to accommodate FeII dπ back bonding to the porphyrin ligand. The ferrous-CO and ferric-NO data are consistent with a small decrease in σ-electron donation from the proximal Cys400 thiolate in the F393A mutant and, to a lesser extent, the F393H mutant, potentially due to a small increase in hydrogen bonding to the proximal ligand. Phe393 seems strategically placed to preserve robust σ-electron donation to the heme iron and to fine-tune its electron density by limiting vinyl group rotation.

Choose Your Own Adventure: A Comprehensive Database of Reactions Catalyzed by Cytochrome P450 BM3 Variants

Cytochrome P450 BM3 monooxygenase is the topic of extensive research as many researchers have evolved this enzyme to generate a variety of products. However, the abundance of information on increasingly diversified variants of P450 BM3 that catalyze a broad array of chemistry is not in a format that enables easy extraction and interpretation. We present a database that categorizes variants by their catalyzed reactions and includes details about substrates to provide reaction context. This database of >1500 P450 BM3 variants is downloadable and machine-readable and includes instructions to maximize ease of gathering information. The database allows rapid identification of commonly reported substitutions, aiding researchers who are unfamiliar with the enzyme in identifying starting points for enzyme engineering. For those actively engaged in engineering P450 BM3, the database, along with this review, provides a powerful and user-friendly platform to understand, predict, and identify the attributes of P450 BM3 variants, encouraging the further engineering of this enzyme.

Multiscale Study on the Intramolecular C-S Bond Formation Catalyzed by P450 Monooxygenase CxnD Involved in the Biosynthesis of Chuangxinmycin: The Critical Roles of Noncrystal Water Molecule and Conformational Change

Cytochrome P450 monooxygenase CxnD catalyzes intramolecular C-S bond formation in the biosynthesis of chuangxinmycin, which is representative of the synthesis of sulfur-containing natural heterocyclic compounds. The intramolecular cyclization usually requires the activation of two reaction sites and a large conformational change; thus, illuminating its detailed reaction mechanism remains challengeable. Here, the reaction pathway of CxnD-catalyzed C-S bond formation was clarified by a series of calculations, including Gaussian accelerated molecular dynamics simulations and quantum mechanical-molecular mechanical calculations. Our results revealed that the C-S formation follows a diradical coupling mechanism. CxnD first employs Cpd I to abstract the hydrogen atom from the imino group of the indole ring, and then, the resulted Cpd II further extracts another hydrogen atom from the thiol group of the side chain to afford a diradical intermediate, in which a noncrystal water molecule entering into the active site after the formation of Cpd I was proved to play an indispensable role. Moreover, the diradical intermediate cannot directly perform the coupling reaction. It should first undergo a series of conformational changes leading to the proximity of two reaction sites. It is the flexibility of the active site of the enzyme and the side chain of the substrate that makes the diradical coupling to be successful.

Cholesterol and oxysterols in retinal neuron-glia interactions: relevance for glaucoma

Cholesterol is an essential component of cellular membranes, crucial for maintaining their structural and functional integrity. It is especially important for nervous tissues, including the retina, which rely on high amounts of plasma membranes for the transmission of the nervous signal. While cholesterol is by far the most abundant sterol, the retina also contains cholesterol precursors and metabolites, especially oxysterols, which are bioactive molecules. Cholesterol lack or excess is deleterious and some oxysterols are known for their effect on neuron survival. Cholesterol homeostasis must therefore be maintained. Retinal glial cells, especially Müller cells, the principal glial cells of the vertebrate retina, provide mechanical, nutritional, and metabolic support for the neighboring neurons. Several pieces of evidence indicate that Müller cells are major actors of cholesterol homeostasis in the retina, as it is known for other glial cells in the brain. This process is based on a close cooperation with neurons, and sterols can be signaling molecules participating in glia-neuron interactions. While some implication of cholesterol in age-related macular degeneration is now recognized, based on epidemiological and laboratory data, evidence for its role in glaucoma is still scarce. The association between cholesterolemia and glaucoma is controversial, but experimental data suggest that sterols could take part in the pathological processes. It has been demonstrated that Müller glial cells are implicated in the development of glaucoma through an ambivalent reactive retinal gliosis process. The early steps contribute to maintaining retinal homeostasis and favor the survival of ganglion cells, which are targeted during glaucoma. If gliosis persists, dysregulation of the neuroprotective functions, cytotoxic effects of gliotic Müller cells and disruption of glia-neuron interactions lead to an acceleration of ganglion cell death. Sterols could play a role in the glial cell response to glaucomatous injury. This represents an understudied but attractive topic to better understand glaucoma and conceive novel preventive or curative strategies. The present review describes the current knowledge on i) sterol metabolism in retinal glial cells, ii) the potential role of cholesterol in glaucoma, and iii) the possible relationships between cholesterol and oxysterols, glial cells and glaucoma. Focus is put on glia-neuron interactions.

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.

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.

Design of a H2 O2 -generating P450SPα fusion protein for high yield fatty acid conversion

Sphingomonas paucimobilis' P450SPα (CYP152B1) is a good candidate as industrial biocatalyst. This enzyme is able to use hydrogen peroxide as unique cofactor to catalyze the fatty acids conversion to α-hydroxy fatty acids, thus avoiding the use of expensive electron-donor(s) and redox partner(s). Nevertheless, the toxicity of exogenous H2 O2 toward proteins and cells often results in the failure of the reaction scale-up when it is directly added as co-substrate. In order to bypass this problem, we designed a H2 O2 self-producing enzyme by fusing the P450SPα to the monomeric sarcosine oxidase (MSOX), as H2 O2 donor system, in a unique polypeptide chain, obtaining the P450SPα -polyG-MSOX fusion protein. The purified P450SPα -polyG-MSOX protein displayed high purity (A417 /A280  = 0.6) and H2 O2 -tolerance (kdecay  = 0.0021 ± 0.000055 min-1 ; ΔA417  = 0.018 ± 0.001) as well as good thermal stability (Tm : 59.3 ± 0.3°C and 63.2 ± 0.02°C for P450SPα and MSOX domains, respectively). The data show how the catalytic interplay between the two domains can be finely regulated by using 500 mM sarcosine as sacrificial substrate to generate H2 O2 . Indeed, the fusion protein resulted in a high conversion yield toward fat waste biomass-representative fatty acids, that is, lauric acid (TON = 6,800 compared to the isolated P450SPα TON = 2,307); myristic acid (TON = 6,750); and palmitic acid (TON = 1,962).

Toxicity of chronic waterborne zinc exposure in the hepatopancreas of white shrimp Litopenaeus vannamei

Zinc (Zn) is necessary for the survival of aquatic organisms; nevertheless, the accumulation of Zn in excessive amounts may have toxic consequences. Few studies focusing on the biochemical, morphological, and transcriptional effects of aqueous Zn in Litopenaeus vannamei have been reported, and the underlying toxic mechanism remains largely unknown. The present study was performed to investigate the growth performance, morphological alterations, physiological changes, and transcriptional responses after Zn exposure at 0 (control), 0.01, 0.1, and 1 mg/L concentrations for 30 days in white shrimp L. vannamei hepatopancreas. The results found that survival rate (SR) and growth performance were significantly reduced in 1 mg/L Zn group. Significant structural damage and significant Zn accumulation in hepatopancreas were observed. The activities of trypsin and amylase (AMS), and the total antioxidant capacity (T-AOC) were attenuated, while the production of reactive oxygen species (ROS) and malondialdehyde (MDA) content were significantly increased after Zn exposure. Many differentially expressed genes (DEGs) were obtained after Zn exposure, and the majority of these DEGs were downregulated. Ten DEGs involved in oxidative stress, immunological response, apoptosis, and other processes were selected for qRT-PCR validation and the expression profiles of these DEGs kept well consistent with the transcriptome data, which confirmed the accuracy and reliability of the transcriptome results. Subsequently, we screened 12 genes to examine the changes of expression in different concentrations in more detail. All the results implying that Zn exposure caused severe histopathological changes and increased Zn accumulation in hepatopancreas, altered immune, antioxidant and detoxifying response by regulating the gene expressions of related genes, and eventually might trigger apoptosis. These findings provide valuable information and a new perspective on the molecular toxicity of crustaceans in response to environmental heavy metal exposure.

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.

Recent Advances in Molecular Fluorescent Probes for CYP450 Sensing and Imaging

Cytochrome P450 (CYP450) is a major drug-metabolizing enzyme system mainly distributed in liver microsomes and involved in the metabolism of many endogenous substances (such as fatty acids and arachidonic acids), and exogenous compounds (such as drugs, toxicants, carcinogens, and procarcinogens). Due to the similarity in structures and catalytic functions between CYP450 isoforms, the lack of effective selective detection tools greatly limits the understanding and the research of their respective physiological roles in living organisms. Until now, several small-molecular fluorescent probes have been employed for selective detection and monitoring of CYP450s (Cytochrome P450 enzymes) in vitro or in vivo owing to the tailored properties, biodegradability, and high temporal and spatial resolution imaging in situ. In this review, we summarize the recent advances in fluorescent probes for CYP450s (including CYP1, CYP2, and CYP3 families), and we discuss and focus on their identification mechanisms, general probe design strategies, and bioimaging applications. We also highlight the potential challenges and prospects of designing new generations of fluorescent probes in CYP450 studies, which will further enhance the diversity, practicality, and clinical feasibility of research into CYP450.

Rational Design of P450 aMOx for Improving Anti-Markovnikov Selectivity Based on the “Butterfly” Model

Aromatic aldehydes are important industrial raw materials mainly synthesized by anti-Markovnikov (AM) oxidation of corresponding aromatic olefins. The AM product selectivity remains a big challenge. P450 aMOx is the first reported enzyme that could catalyze AM oxidation of aromatic olefins. Here, we reported a rational design strategy based on the “butterfly” model of the active site of P450 aMOx. Constrained molecular dynamic simulations and a binding energy analysis of key residuals combined with an experimental alanine scan were applied. As a result, the mutant A275G showed high AM selectivity of >99%. The results also proved that the “butterfly” model is an effective design strategy for enzymes.

First purified recombinant CYP75B including transmembrane helix with unexpected high substrate specificity to (2R)-naringenin

Anthochlor pigments (chalcones and aurones) play an important role in yellow flower colourization, the formation of UV-honey guides and show numerous health benefits. The B-ring hydroxylation of chalcones is performed by membrane bound cytochrome P450 enzymes. It was assumed that usual flavonoid 3′-hydroxlases (F3′Hs) are responsible for the 3,4- dihydroxy pattern of chalcones, however, we previously showed that a specialized F3′H, namely chalcone 3-hydroxylase (CH3H), is necessary for the hydroxylation of chalcones. In this study, a sequence encoding membrane bound CH3H from Dahlia variabilis was recombinantly expressed in yeast and a purification procedure was developed. The optimized purification procedure led to an overall recovery of 30% recombinant DvCH3H with a purity of more than 84%. The enzyme was biochemically characterized with regard to its kinetic parameters on various substrates, including racemic naringenin, as well as its enantiomers (2S)-, and (2R)-naringenin, apigenin and kaempferol. We report for the first time the characterization of a purified Cytochrome P450 enzyme from the flavonoid biosynthesis pathway, including the transmembrane helix. Further, we show for the first time that recombinant DvCH3H displays a higher affinity for (2R)-naringenin than for (2S)-naringenin, although (2R)-flavanones are not naturally formed by chalcone isomerase.

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.

Peroxide-Mediated Oxygenation of Organic Compounds by Fungal Peroxygenases

Unspecific peroxygenases (UPOs), whose sequences can be found in the genomes of thousands of filamentous fungi, many yeasts and certain fungus-like protists, are fascinating biocatalysts that transfer peroxide-borne oxygen (from H2O2 or R-OOH) with high efficiency to a wide range of organic substrates, including less or unactivated carbons and heteroatoms. A twice-proline-flanked cysteine (PCP motif) typically ligates the heme that forms the heart of the active site of UPOs and enables various types of relevant oxygenation reactions (hydroxylation, epoxidation, subsequent dealkylations, deacylation, or aromatization) together with less specific one-electron oxidations (e.g., phenoxy radical formation). In consequence, the substrate portfolio of a UPO enzyme always combines prototypical monooxygenase and peroxidase activities. Here, we briefly review nearly 20 years of peroxygenase research, considering basic mechanistic, molecular, phylogenetic, and biotechnological aspects.

pH Changes That Induce an Axial Ligand Effect on Nonheme Iron(IV) Oxo Complexes with an Appended Aminopropyl Functionality

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

Mechanism-Based Inactivation of Cytochrome P450 Enzymes: Computational Insights

Mechanism-based inactivation (MBI) refers to the metabolic bioactivation of a xenobiotic by cytochrome P450s to a highly reactive intermediate which subsequently binds to the enzyme and leads to the quasi-irreversible or irreversible inhibition. Xenobiotics, mainly drugs with specific functional units, are the major sources of MBI. Two possible consequences of MBI by medicinal compounds are drug-drug interaction and severe toxicity that are observed and highlighted by clinical experiments. Today almost all of these latent functional groups (e.g., thiophene, furan, alkylamines, etc.) are known, and their features and mechanisms of action, owing to the vast experimental and theoretical studies, are determined. In the past decade, molecular modeling techniques, mostly density functional theory, have revealed the most feasible mechanism that a drug undergoes by P450 enzymes to generate a highly reactive intermediate. In this review, we provide a comprehensive and detailed picture of computational advances toward the elucidation of the activation mechanisms of various known groups with MBI activity. To this aim, we briefly describe the computational concepts to carry out and analyze the mechanistic investigations, and then, we summarize the studies on compounds with known inhibition activity including thiophene, furan, alkylamines, terminal acetylene, etc. This study can be reference literature for both theoretical and experimental (bio)chemists in several different fields including rational drug design, the process of toxicity prevention, and the discovery of novel inhibitors and catalysts.

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.

Exploring the Potential of Cytochrome P450 CYP109B1 Catalyzed Regio-and Stereoselective Steroid Hydroxylation

Cytochrome P450 enzyme CYP109B1 is a versatile biocatalyst exhibiting hydroxylation activities toward various substrates. However, the regio- and stereoselective steroid hydroxylation by CYP109B1 is far less explored. In this study, the oxidizing activity of CYP109B1 is reconstituted by coupling redox pairs from different sources, or by fusing it to the reductase domain of two self-sufficient P450 enzymes P450RhF and P450BM3 to generate the fused enzyme. The recombinant Escherichia coli expressing necessary proteins are individually constructed and compared in steroid hydroxylation. The ferredoxin reductase (Fdr_0978) and ferredoxin (Fdx_1499) from Synechococcus elongates is found to be the best redox pair for CYP109B1, which gives above 99% conversion with 73% 15β selectivity for testosterone. By contrast, the rest ones and the fused enzymes show much less or negligible activity. With the aid of redox pair of Fdr_0978/Fdx_1499, CYP109B1 is used for hydroxylating different steroids. The results show that CYP109B1 displayed good to excellent activity and selectivity toward four testosterone derivatives, giving all 15β-hydroxylated steroids as main products except for 9 (10)-dehydronandrolone, for which the selectivity is shifted to 16β. While for substrates bearing bulky substitutions at C17 position, the activity is essentially lost. Finally, the origin of activity and selectivity for CYP109B1 catalyzed steroid hydroxylation is revealed by computational analysis, thus providing theoretical basis for directed evolution to further improve its catalytic properties.

Cytochrome P450 enzymes in fungal natural product biosynthesis

Covering: 2015 to the end of 2020 Fungal-derived polyketides, non-ribosomal peptides, terpenoids and their hybrids contribute significantly to the chemical space of total natural products. Cytochrome P450 enzymes play essential roles in fungal natural product biosynthesis with their broad substrate scope, great catalytic versatility and high frequency of involvement. Due to the membrane-bound nature, the functional and mechanistic understandings for fungal P450s have been limited for quite a long time. However, recent technical advances, such as the efficient and precise genome editing techniques and the development of several filamentous fungal strains as heterologous P450 expression hosts, have led to remarkable achievements in fungal P450 studies. Here, we provide a comprehensive review to cover the most recent progresses from 2015 to 2020 on catalytic functions and mechanisms, research methodologies and remaining challenges in the fast-growing field of fungal natural product biosynthetic P450s.

Catalytic approach to in vivo metabolism of atractylenolide III using biomimetic iron–porphyrin complexes

In vivo oxygenation of atractylenolide III during the metabolism and in situ oxidation mechanism by an iron–porphyrin iron(iv)-oxo porphyrin π-cation–radical complex.

Advanced Understanding of the Electron Transfer Pathway of Cytochrome P450s

Cytochrome P450s are heme-thiolate enzymes that participate in carbon source assimilation, natural compound biosynthesis and xenobiotic metabolism in all kingdoms of life. P450s can catalyze various reactions by using a wide range of organic compounds, thus exhibiting great potential in biotechnological applications. The catalytic reactions of P450s are driven by electron equivalents that are sourced from pyridine nucleotides and delivered by cognate or matching redox partners (RPs). The electron transfer (ET) route from RPs to P450s involves one or more redox center-containing domains. As the rate of ET is one of the main determinants of P450 efficacy, an in-depth understanding of the P450 ET pathway should increase our knowledge of these important enzymes and benefit their further applications. Here, the various P450 RP systems along with current understanding of their ET routes will be reviewed. Notably, state-of-the-art structural studies of the two main types of self-sufficient P450 will also be summarized.

Assessing the Covalent Attachment and Energy Transfer Capabilities of Upconverting Phosphors With Cofactor Containing Bioactive Enzymes

Upconverting phosphors (UCPs) convert multiple low energy photons into higher energy emission via the process of photon upconversion and offer an attractive alternative to organic fluorophores for use as luminescent probes. Examples of biosensors utilizing the apparent energy transfer of UCPs and nanophosphors (UCNPs) with biomolecules have started to appear in the literature but very few exploit the covalent anchoring of the biomolecule to the surface of the UCP to improve the sensitivity of the systems. Here, we demonstrate a robust and versatile method for the covalent attachment of biomolecules to the surface of a variety of UCPs and UCNPs in which the UCPs were capped with functionalized silica in order to provide a surface to covalently conjugate biomolecules with surface-accessible cysteines. Variants of BM3Heme, cytochrome C, glucose oxidase, and glutathione reductase were then attached via maleimide-thiol coupling. BM3Heme, glucose oxidase, and glutathione reductase were shown to retain their activity when coupled to the UCPs potentially opening up opportunities for biosensing applications.

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.

Catalytic Mechanism of Aromatic Nitration by Cytochrome P450 TxtE: Involvement of a Ferric-Peroxynitrite Intermediate

The cytochromes P450 are heme-dependent enzymes that catalyze many vital reaction processes in the human body related to biodegradation and biosynthesis. They typically act as mono-oxygenases; however, the recently discovered P450 subfamily TxtE utilizes O₂ and NO to nitrate aromatic substrates such as L-tryptophan. A direct and selective aromatic nitration reaction may be useful in biotechnology for the synthesis of drugs or small molecules. Details of the catalytic mechanism are unknown, and it has been suggested that the reaction should proceed through either an iron(III)-superoxo or an iron(II)-nitrosyl intermediate. To resolve this controversy, we used stopped-flow kinetics to provide evidence for a catalytic cycle where dioxygen binds prior to NO to generate an active iron(III)-peroxynitrite species that is able to nitrate l-Trp efficiently. We show that the rate of binding of O₂ is faster than that of NO and also leads to l-Trp nitration, while little evidence of product formation is observed from the iron(II)-nitrosyl complex. To support the experimental studies, we performed density functional theory studies on large active site cluster models. The studies suggest a mechanism involving an iron(III)-peroxynitrite that splits homolytically to form an iron(IV)-oxo heme (Compound II) and a free NO₂ radical via a small free energy of activation. The latter activates the substrate on the aromatic ring, while compound II picks up the ipso-hydrogen to form the product. The calculations give small reaction barriers for most steps in the catalytic cycle and, therefore, predict fast product formation from the iron(III)-peroxynitrite complex. These findings provide the first detailed insight into the mechanism of nitration by a member of the TxtE subfamily and highlight how the enzyme facilitates this novel reaction chemistry.

Late-Stage Diversification of Natural Products

Late-stage diversification of natural products is an efficient way to generate natural product derivatives for drug discovery and chemical biology. Benefiting from the development of site-selective synthetic methodologies, late-stage diversification of natural products has achieved notable success. This outlook will outline selected examples of novel methodologies for site-selective transformations of reactive functional groups and inert C-H bonds that enable late-stage diversification of complex natural products. Accordingly, late-stage diversification provides an opportunity to rapidly access various derivatives for modifying lead compounds, identifying cellular targets, probing protein-protein interactions, and elucidating natural product biosynthetic relationships.

Artificial control of the multistep oxidation reactions catalyzed by the cytochrome P450 enzyme RosC

The cytochrome P450 monooxygenase RosC catalyzes the three-step oxidation reactions, which leads to the formation of a hydroxy, formyl, and carboxy group at C-20 during rosamicin biosynthesis in Micromonospora rosaria IFO13697. To determine if amino acid substitutions in RosC could allow for the control of the multistep oxidation reactions, we screened RosC random mutants. The RosC mutant RM30, with five amino acid substitutions (P107S, L176Q, S254N, V277A, and I319N), catalyzed only the first step of the oxidation reaction. Whole-cell assays using Escherichia coli cells expressing RosC mutants with single and double amino acid substitutions derived from RM30 indicated that P107S/L176Q, P107S/V277A, P107S/I319N, L176Q/V277A, L176Q/I319N, and S254N/V277A significantly reduced the catalytic activity of the second reaction, which is alcohol oxidation. Of the previously mentioned mutants, double mutants containing L176Q, which was presumed to occur in the FG loop region, lost the total catalytic activity of the third reaction (aldehyde oxidation). Additionally, an engineered M. rosaria strain with rosC disruption, which introduced the gene encoding the RosC mutants P107S/L176Q and P107S/V277A preferentially produced 20-dihydrorosamicin, which is formed after the first oxidation reaction of RosC.

P450-BM3-Catalyzed Sulfoxidation versus Hydroxylation: A Common or Two Different Catalytically Active Species?

While the mechanism of the P450-catalyzed oxidative hydroxylation of organic compounds has been studied in detail for many years, less is known about sulfoxidation. Depending upon the structure of the respective substrate, heme-Fe═O (Cpd I), heme-Fe(III)-OOH (Cpd 0), and heme-Fe(III)-H2O2 (protonated Cpd 0) have been proposed as reactive intermediates. In the present study, we consider the transformation of isosteric substrates via sulfoxidation and oxidative hydroxylation, respectively, catalyzed by regio- and enantioselective mutants of P450-BM3 which were constructed by directed evolution. 1-Thiochromanone and 1-tetralone were used as the isosteric substrates because, unlike previous studies involving fully flexible compounds such as thia-fatty acids and fatty acids, respectively, these compounds are rigid and cannot occur in a multitude of different conformations and binding modes in the large P450-BM3 binding pocket. The experimental results comprising activity and regio- and enantioselectivity, flanked by molecular dynamics computations within a time scale of 300 ns and QM/MM calculations of transition-state energies, unequivocally show that heme-Fe═O (Cpd I) is the common catalytically active intermediate in both sulfoxidation and oxidative hydroxylation.

Engineering Plant Cytochrome P450s for Enhanced Synthesis of Natural Products: Past Achievements and Future Perspectives

Cytochrome P450s (P450s) are the most versatile catalysts and are widely used by plants to synthesize a vast array of structurally diverse specialized metabolites that not only play essential ecological roles but also constitute a valuable resource for the development of new drugs. To accelerate the metabolic engineering of these high-value metabolites, it is imperative to identify and characterize pathway P450s, and to further improve their activities through protein engineering. In this review, we focus on P450 engineering and summarize the major strategies for enhancing the stability and activity of P450s and successful cases of P450 engineering. Studies in which the functions of P450s were altered to create de novo metabolic pathways or novel compounds are discussed as well. We also overview emerging tools, specifically DNA synthesis, machine learning, and de novo protein design, as well as the evolutionary patterns of P450s unveiled from a massive number of DNA sequences that could be integrated to accelerate the engineering of these enzymes. These approaches would greatly aid in the exploitation of plant-specialized metabolites or derivatives for various uses including medical applications.

Bacterial cytochrome P450-catalyzed regio- and stereoselective steroid hydroxylation enabled by directed evolution and rational design

Steroids are the most widely marketed products by the pharmaceutical industry after antibiotics. Steroid hydroxylation is one of the most important functionalizations because their derivatives enable a higher biological activity compared to their less polar non-hydroxylated analogs. Bacterial cytochrome P450s constitute promising biocatalysts for steroid hydroxylation due to their high expression level in common workhorses like Escherichia coli. However, they often suffer from wrong or insufficient regio- and/or stereoselectivity, low activity, narrow substrate range as well as insufficient thermostability, which hampers their industrial application. Fortunately, these problems can be generally solved by protein engineering based on directed evolution and rational design. In this work, an overview of recent developments on the engineering of bacterial cytochrome P450s for steroid hydroxylation is presented.

Cross-linking of aromatic phenolate groups by cytochrome P450 enzymes: a model for the biosynthesis of vancomycin by OxyB

Aromatic cross-linking by cytochrome P450 enzymes was studied computationally. P450 Compound I rapidly abstracts two weak phenolic H-atoms that link up via a rate-determining C–O bond formation.

Iron- and cobalt-catalyzed C(sp3)–H bond functionalization reactions and their application in organic synthesis

This review highlights the developments in iron and cobalt catalyzed C(sp³)–H bond functionalization reactions with emphasis on their applications in organic synthesis, i.e. natural products and pharmaceuticals synthesis and/or modification.

Structural Perspectives of the CYP3A Family and Their Small Molecule Modulators in Drug Metabolism

Cytochrome P450 enzymes function to catalyze a wide range of reactions, many of which are critically important for drug response. Members of the human cytochrome P450 3A (CYP3A) family are particularly important in drug clearance, and they collectively metabolize more than half of all currently prescribed medications. The ability of these enzymes to bind a large and structurally diverse set of compounds increases the chances of their modulating or facilitating drug metabolism in unfavorable ways. Emerging evidence suggests that individual enzymes in the CYP3A family play discrete and important roles in catalysis and disease progression. Here we review the similarities and differences among CYP3A enzymes with regard to substrate recognition, metabolism, modulation by small molecules, and biological consequence, highlighting some of those with clinical significance. We also present structural perspectives to further characterize the basis of these comparisons.

Comparative Sequence Analysis of TRI1 of Fusarium

Trichothecene mycotoxins are a class of secondary metabolites produced by multiple genera of fungi, including certain plant pathogenic Fusarium species. Functional variation in the TRI1 gene produces a novel Type A trichothecene called NX-2 in strains of F. graminearum. Using a bioinformatics approach, a systematic analysis of 52 translated TRI1 sequences of Fusarium species, including five F. graminearum NX-2 producers and four F. graminearum non-NX-2 producers, was conducted to explain the functional difference of TRI1p of FGNX-2. An assessment of several signature motifs of fungal P450s revealed amino acid substitutions in addition to the post-translational N-X-S/T sequons motif, which is indicative of N-linked glycosylation of this TRI1-encoded protein characteristic of NX-2 producers. There was evidence of selection bias, where TRI1 gene sequences were found to be under positive selection and, therefore, under functional constraints. The cumulative amino acid changes in the TRI1p sequences were reflected in the phylogenetic analyses which revealed species-specific clustering with a distinct separation of FGNX-2 from FG-non-NX-2 producers with high bootstrap support. Together, our findings provide insight into the amino acid sequence features responsible for the functional diversification of this TRI1p.

De novo biosynthesis of 8-hydroxyoctanoic acid via a medium-chain length specific fatty acid synthase and cytochrome P450 in Saccharomyces cerevisiae

Terminally hydroxylated fatty acids or dicarboxylic acids are industrially relevant compounds with broad applications. Here, we present the proof of principle for the de novo biosynthesis of 8-hydroxyoctanoic acid from glucose and ethanol in the yeast Saccharomyces cerevisiae. Toxicity tests with medium-chain length ω-hydroxy fatty acids and dicarboxylic acids revealed little or no growth impairments on yeast cultures even at higher concentrations. The ability of various heterologous cytochrome P450 enzymes in combination with their cognate reductases for ω-hydroxylation of externally fed octanoic acid were compared. Finally, the most efficient P450 enzyme system was expressed in a yeast strain, whose fatty acid synthase was engineered for octanoic acid production, resulting in de novo biosynthesis of 8-hydroxyoctanoic acid up to 3 ​mg/l. Accumulation of octanoic acid revealed that cytochromes P450 activities were limiting 8-hydroxyoctanoic acid synthesis. The hydroxylation of both externally added and intracellularly produced octanoic acid was strongly dependent on the carbon source used, with ethanol being preferred. We further identified the availability of heme, a cofactor needed for P450 activity, as a limiting factor of 8-hydroxyoctanoic acid biosynthesis.

•

Low toxic effects of medium-chain ω-hydroxy fatty acids on yeast cells .

•

Systematic comparison of cytochrome P450 enzyme activities on octanoic acid .

•

De novo biosynthesis of 8-hydroxyoctanoic acid .

•

Improvement of cytochrome P450 activity with ethanol or by addition of hemin .

Hydroxylation of Steroids by a Microbial Substrate-Promiscuous P450 Cytochrome (CYP105D7): Key Arginine Residues for Rational Design

Our previous study showed that CYP105D7, a substrate-promiscuous P450, catalyzes the hydroxylation of 1-deoxypentalenic acid, diclofenac, naringenin, and compactin. In this study, 14 steroid compounds were screened using recombinant Escherichia coli cells harboring genes encoding CYP105D7 and redox partners (Pdx/Pdr, RhFRED, and FdxH/FprD), and the screening identified steroid A-ring 2β- and D-ring 16β-hydroxylation activity. Wild-type CYP105D7 was able to catalyze the hydroxylation of five steroids (testosterone, progesterone, 4-androstene-3,17-dione, adrenosterone, and cortisone) with low (<10%) conversion rates. Structure-guided site-directed mutagenesis of arginine residues around the substrate entrance and active site showed that the R70A and R190A single mutants and an R70A/R190A double mutant exhibited greatly enhanced conversion rates for steroid hydroxylation. For the conversion of testosterone in particular, the R70A/R190A mutant's k _cat/K_m values increased 1.35-fold and the in vivo conversion rates increased significantly by almost 9-fold with high regio- and stereoselectivity. Molecular docking analysis revealed that when Arg70 and Arg190 were replaced with alanine, the volume of the substrate access and binding pocket increased 1.08-fold, which might facilitate improvement of the hydroxylation efficiency of steroids.IMPORTANCE Cytochrome P450 monooxygenases (P450s) are able to introduce oxygen atoms into nonreactive hydrocarbon compounds under mild conditions, thereby offering significant advantages compared to chemical catalysts. Promiscuous P450s with broad substrate specificity and reaction diversity have significant potential for applications in various fields, including synthetic biology. The study of the function, molecular mechanisms, and rational engineering of substrate-promiscuous P450s from microbial sources is important to fulfill this potential. Here, we present a microbial substrate-promiscuous P450, CYP105D7, which can catalyze hydroxylation of steroids. The loss of the bulky side chains of Arg70 and Arg190 in the active site and substrate entrance resulted in an up to 9-fold increase in the substrate conversion rate. These findings will support future rational and semirational engineering of P450s for applications as biocatalysts.

A Continuing Career in Biocatalysis: Frances H. Arnold

On the occasion of Professor Frances H. Arnold's recent acceptance of the 2018 Nobel Prize in Chemistry, we honor her numerous contributions to the fields of directed evolution and biocatalysis. Arnold pioneered the development of directed evolution methods for engineering enzymes as biocatalysts. Her highly interdisciplinary research has provided a ground not only for understanding the mechanisms of enzyme evolution but also for developing commercially viable enzyme biocatalysts and biocatalytic processes. In this Account, we highlight some of her notable contributions in the past three decades in the development of foundational directed evolution methods and their applications in the design and engineering of enzymes with desired functions for biocatalysis. Her work has created a paradigm shift in the broad catalysis field.

A Structural and Functional Model for the Tris-Histidine Motif in Cysteine Dioxygenase

The iron(II) complexes [Fe(L)(MeCN)₃ ](SO₃ CF₃ )₂ (L are two derivatives of tris(2-pyridyl)-based ligands) have been synthesized as models for cysteine dioxygenase (CDO). The molecular structure of one of the complexes exhibits octahedral coordination geometry and the Fe-N_py bond lengths [1.953(4)-1.972(4) Å] are similar to those in the Cys-bound Fe^II -CDO; Fe-N_His : 1.893-2.199 Å. The iron(II) centers of the model complexes exhibit relatively high Fe^III/II redox potentials (E_1/2 =0.988-1.380 V vs. ferrocene/ferrocenium electrode, Fc/Fc⁺ ), within the range for O₂ activation and typical for the corresponding nonheme iron enzymes. The reaction of in situ generated [Fe(L)(MeCN)(SPh)]⁺ with excess O₂ in acetonitrile (MeCN) yields selectively the doubly oxygenated phenylsulfinic acid product. Isotopic labeling studies using ¹⁸ O₂ confirm the incorporation of both oxygen atoms of O₂ into the product. Kinetic and preliminary DFT studies reveal the involvement of an Fe^III peroxido intermediate with a rhombic S= 1 / 2 Fe^III center (687-696 nm; g≈2.46-2.48, 2.13-2.15, 1.92-1.94), similar to the spectroscopic signature of the low-spin Cys-bound Fe^III CDO (650 nm, g≈2.47, 2.29, 1.90). The proposed Fe^III peroxido intermediates have been trapped, and the O-O stretching frequencies are in the expected range (approximately 920 and 820 cm^-1 for the alkyl- and hydroperoxido species, respectively). The model complexes have a structure similar to that of the enzyme and structural aspects as well as the reactivity are discussed.

One small step for cytochrome P450 in its catalytic cycle, one giant leap for enzymology

The intermediates operating in the cytochrome P450 catalytic cycle have been investigated for more than half a century, fascinating many enzymologists. Each intermediate has its unique role to carry out diverse oxidations. Natural time course of the catalytic cycle is quite fast, hence, not all of the reactive intermediates could be isolated during physiological catalysis. Different high-valent iron intermediates have been proposed as primary oxidants: the candidates are compound 0 (Cpd 0, [FeOOH][Formula: see text]P450) and compound I (Cpd I, Fe(IV)[Formula: see text]O por[Formula: see text]P450). Among them, the role of Cpd I in hydroxylation is fairly well understood due the discovery of the peroxide shunt. This review endeavors to put the outstanding research efforts conducted to isolate and characterize the intermediates together. In addition to spectral features of each intermediate in the catalytic cycle, the oxidizing powers of Cpd 0 and Cpd I will be discussed along with most recent scientific findings.

Advances in Sustainable Catalysis: A Computational Perspective

The enormous challenge of moving our societies to a more sustainable future offers several exciting opportunities for computational chemists. The first principles approach to "catalysis by design" will enable new and much greener chemical routes to produce vital fuels and fine chemicals. This prospective outlines a wide variety of case studies to underscore how the use of theoretical techniques, from QM/MM to unrestricted DFT and periodic boundary conditions, can be applied to biocatalysis and to both homogeneous and heterogenous catalysts of all sizes and morphologies to provide invaluable insights into the reaction mechanisms they catalyze.

Novel insights into P450 BM3 interactions with FDA-approved antifungal azole drugs

Flavocytochrome P450 BM3 is a natural fusion protein constructed of cytochrome P450 and NADPH-cytochrome P450 reductase domains. P450 BM3 binds and oxidizes several mid- to long-chain fatty acids, typically hydroxylating these lipids at the ω-1, ω-2 and ω-3 positions. However, protein engineering has led to variants of this enzyme that are able to bind and oxidize diverse compounds, including steroids, terpenes and various human drugs. The wild-type P450 BM3 enzyme binds inefficiently to many azole antifungal drugs. However, we show that the BM3 A82F/F87V double mutant (DM) variant binds substantially tighter to numerous azole drugs than does the wild-type BM3, and that their binding occurs with more extensive heme spectral shifts indicative of complete binding of several azoles to the BM3 DM heme iron. We report here the first crystal structures of P450 BM3 bound to azole antifungal drugs - with the BM3 DM heme domain bound to the imidazole drugs clotrimazole and tioconazole, and to the triazole drugs fluconazole and voriconazole. This is the first report of any protein structure bound to the azole drug tioconazole, as well as the first example of voriconazole heme iron ligation through a pyrimidine nitrogen from its 5-fluoropyrimidine ring.