Properties of "Stable" Mosquito Cytochrome P450 Enzymes

The use of insecticides is widespread in the control of debilitating mosquito-borne diseases. P450 enzymes (CYPs) play essential roles in mosquito physiological function but also in the enzymatic detoxification of xenobiotics. Broadly speaking, CYPs can be classified as "stable", meaning those that have no or very few paralogs, and "labile", constituting gene families with many paralogous members. The evolutionary dichotomy between "stable" and "labile" P450 genes is fuzzy and there is not a clear phylogenetic demarcation between P450s involved in detoxification and P450s involved in essential metabolic processes. In this study, bioinformatic methods were used to explore differences in the sequences of "stable" and "labile" P450s that may facilitate their functional classification. Genomic and sequence data of Anopheles gambiae (Agam), Aedes aegypti (Aaeg), and Culex quinquefasciatus (Cqui) CYPs were obtained from public databases. The results of this study show that "stable" CYPs are encoded by longer genes, have longer introns and more exons, and contain a higher proportion of hydrophobic amino acids than "labile" CYPs. Compared to "labile" CYPs, a significantly higher proportion of "stable" CYPs are associated with biosynthetic and developmental processes.

Deep mutational scanning of CYP2C19 in human cells reveals a substrate specificity-abundance tradeoff

The cytochrome P450s enzyme family metabolizes ∼80% of small molecule drugs. Variants in cytochrome P450s can substantially alter drug metabolism, leading to improper dosing and severe adverse drug reactions. Due to low sequence conservation, predicting variant effects across cytochrome P450s is challenging. Even closely related cytochrome P450s like CYP2C9 and CYP2C19, which share 92% amino acid sequence identity, display distinct phenotypic properties. Using variant abundance by massively parallel sequencing, we measured the steady-state protein abundance of 7,660 single amino acid variants in CYP2C19 expressed in cultured human cells. Our findings confirmed critical positions and structural features essential for cytochrome P450 function, and revealed how variants at conserved positions influence abundance. We jointly analyzed 4,670 variants whose abundance was measured in both CYP2C19 and CYP2C9, finding that the homologs have different variant abundances in substrate recognition sites within the hydrophobic core. We also measured the abundance of all single and some multiple wild type amino acid exchanges between CYP2C19 and CYP2C9. While most exchanges had no effect, substitutions in substrate recognition site 4 reduced abundance in CYP2C19. Double and triple mutants showed distinct interactions, highlighting a region that points to differing thermodynamic properties between the 2 homologs. These positions are known contributors to substrate specificity, suggesting an evolutionary tradeoff between stability and enzymatic function. Finally, we analyzed 368 previously unannotated human variants, finding that 43% had decreased abundance. By comparing variant effects between these homologs, we uncovered regions underlying their functional differences, advancing our understanding of this versatile family of enzymes.

Carving the Active Site of CYP153A7 Monooxygenase for Improving Terminal Hydroxylation of Medium-Chain Fatty Acids

The P450-mediated terminal hydroxylation of non-activated C-H bonds is a chemically challenging reaction. CYP153A7 monooxygenase discovered in Sphingomonas sp. HXN200 belongs to the CYP153A subfamily and shows a pronounced terminal selectivity. Herein, we report the significantly improved terminal hydroxylation activity of CYP153A7 by redesign of the substrate binding pocket based on molecular docking of CYP153A7-C 8:0 and sequence alignments. Some of the resultant single mutants were advantageous over the wild-type enzyme with higher reaction rates, achieving a complete conversion of n- octanoic acid (C 8:0. 1 mM) in a shorter period. Especially, a single-mutation variant, D258E, showed 3.8-fold higher catalytic efficiency than the wild type toward the terminal hydroxylation of medium-chain fatty acid C 8:0 into the high value-added product 8-hydroxyoctanoic acid.

Metabolism analysis of 17α-ethynylestradiol by Pseudomonas citronellolis SJTE-3 and identification of the functional genes

Synthetic estrogens are the most hazardous and persistent environmental estrogenic contaminants, with few reports on their biodegradation. Pseudomonas citronellolis SJTE-3 degraded natural steroids efficiently and metabolized 17α-ethynylestradiol (EE2) with the addition of different easily used energy sources (glucose, peptone, ethanol, yeast extract, fulvic acid and ammonia). Over 92% of EE2 (1 mg/L) and 55% of EE2 (10 mg/L) in culture were removed in seven days with the addition of 0.1% ethanol, and the EE2-biotransforming efficiency increased with the increasing ethanol concentrations. Two novel intermediate metabolites of EE2 (C₂₂H₂₂O and C₁₈H₃₄O₂) were identified with high-performance liquid chromatography (HPLC) and GC-Orbitrap/MS. Comparative analysis and genome mining revealed strain SJTE-3 contained a unique genetic basis for EE2 metabolism, and the putative EE2-degrading genes exhibited dispersed distribution. The EE2 metabolism of strain SJTE-3 was inducible and the transcription of eight genes were significantly induced by EE2. Three genes (sdr3, yjcH and cyp2) encoding a short-chain dehydrogenase, a membrane transporter and a cytochrome P450 hydroxylase, respectively, were vital for EE2 metabolism in strain SJTE-3; their over-expression accelerated EE2 metabolic processes and advanced the generation of intermediate metabolites. This work could promote the study of bacterial EE2 metabolism mechanisms and facilitate efficient bioremediation for EE2 pollution.

Cytochromes P450 (P450s): A review of the class system with a focus on prokaryotic P450s

Cytochromes P450 (P450s) are a large superfamily of heme-containing monooxygenases. P450s are found in all Kingdoms of life and exhibit incredible diversity, both at sequence level and also on a biochemical basis. In the majority of cases, P450s can be assigned into one of ten classes based on their associated redox partners, domain architecture and cellular localization. Prokaryotic P450s now represent a large diverse collection of annotated/known enzymes, of which many have great potential biocatalytic potential. The self-sufficient P450 classes (Class VII/VIII) have been explored significantly over the past decade, with many annotated and biochemically characterized members. It is clear that the prokaryotic P450 world is expanding rapidly, as the number of published genomes and metagenome studies increases, and more P450 families are identified and annotated (CYP families).

Blooming of Unusual Cytochrome P450s by Tandem Duplication in the Pathogenic Fungus Conidiobolus coronatus

While the Zygomycete fungus Conidiobolus coronatus primarily infects insects, it can be pathogenic to mammals as well, including humans. High variability in the treatment of this fungal infection with currently available drugs, including azole drugs is a very common phenomenon. Azoles bind to the cytochrome P450 monooxygenases (P450s/CYP) including CYP51, a sterol 14-α-demethylase, inhibiting the synthesis of cell membrane ergosterol and thus leading to the elimination of infecting fungi. Despite P450’s role as a drug target, to date, no information on C. coronatus P450s has been reported. Genome-wide data mining has revealed the presence of 142 P450s grouped into 12 families and 21 subfamilies in C. coronatus. Except for CYP51, the remaining 11 P450 families are new (CYP5854-CYP5864). Despite having a large number of P450s among entomopathogenic fungi, C. coronatus has the lowest number of P450 families, which suggests blooming P450s. Further analysis has revealed that 79% of the same family P450s is tandemly positioned, suggesting that P450 tandem duplication led to the blooming of P450s. The results of this study; i.e., unravelling the C. coronatus P450 content, will certainly help in designing experiments to understand P450s’ role in C. coronatus physiology, including a highly variable response to azole drugs with respect to P450s.

Comparative analyses and structural insights of the novel cytochrome P450 fusion protein family CYP5619 in Oomycetes

Phylogenetic and structural analysis of P450 proteins fused to peroxidase/dioxygenase has not been reported yet. We present phylogenetic and in silico structural analysis of the novel P450 fusion family CYP5619 from the deadliest fish pathogenic oomycete, Saprolegnia diclina. Data-mining and annotation of CYP5619 members revealed their unique presence in oomycetes. CYP5619 members have the highest number of conserved amino acids among eukaryotic P450s. The highest number of conserved amino acids (78%) occurred in the peroxidase/dioxygenase domain compared to the P450 domain (22%). In silico structural analysis using a high-quality CYP5619A1 model revealed that CYP5619A1 has characteristic P450 structural motifs including EXXR and CXG. However, the heme-binding domain (CXG) in CYP5619 members was found to be highly degenerated. The in silico substrate binding pattern revealed that CYP5619A1 have a high affinity to medium chain fatty acids. Interestingly, the controlling agent of S. diclina malachite green was predicted to have the highest binding affinity, along with linoleic acid. However, unlike fatty acids, none of the active site amino acids formed hydrogen bonds with malachite green. The study’s results will pave the way for assessing CYP5619A1’s role in S. diclina physiology, including the nature of malachite green binding.

De Novo Sequencing of a Sparassis latifolia Genome and Its Associated Comparative Analyses

Known to be rich in β-glucan, Sparassis latifolia (S. latifolia) is a valuable edible fungus cultivated in East Asia. A few studies have suggested that S. latifolia is effective on antidiabetic, antihypertension, antitumor, and antiallergen medications. However, it is still unclear genetically why the fungus has these medical effects, which has become a key bottleneck for its further applications. To provide a better understanding of this fungus, we sequenced its whole genome, which has a total size of 48.13 megabases (Mb) and contains 12,471 predicted gene models. We then performed comparative and phylogenetic analyses, which indicate that S. latifolia is closely related to a few species in the antrodia clade including Fomitopsis pinicola, Wolfiporia cocos, Postia placenta, and Antrodia sinuosa. Finally, we annotated the predicted genes. Interestingly, the S. latifolia genome encodes most enzymes involved in carbohydrate and glycoconjugate metabolism and is also enriched in genes encoding enzymes critical to secondary metabolite biosynthesis and involved in indole, terpene, and type I polyketide pathways. As a conclusion, the genome content of S. latifolia sheds light on its genetic basis of the reported medicinal properties and could also be used as a reference genome for comparative studies on fungi.

Molecular cloning and 3D model of first cytochrome P450 from CYP3A subfamily in saltwater crocodile (Crocodylus porosus)

Cytochrome P450s (CYPs) play critical role in oxidative metabolism of numerous xenobiotics and endogenous compounds. The first CYP3A subfamily member in saltwater crocodile has been cloned and modelled for three-dimensional (3D) structure. The full-length cDNA was obtained employing reverse transcription polymerase chain reaction (RT-PCR) strategy and rapid amplification of cDNA ends (RACE). The cDNA sequence of 1659 nucleotides includes 132 nucleotides from 5' untranslated region (UTR), an open reading frame of 1527 nucleotides encoding 509 amino acids designated as CYP3A163. The alignment of CYP3A163 sequence with CYP3A subfamily across the lineages exhibit the loss of 1 residue in birds and 7 residues in mammals in comparison to reptiles suggesting the adaptation processes during evolution. The amino acid identity of CYP3A163 with Alligator mississippiensis CYP3A77 and Homo sapiens CYP3A4 is 91% and 62% respectively. The 3D structure of CYP3A163 modelled using human CYP3A4 structure as a template with Phyre² software, represents high similarity with its functionally important motifs and catalytic domain. Both sequence and structure of CYP3A163 display the common and conserved features of CYP3A subfamily. Overall, this study provides primary molecular and structural data of CYP3A163 required to investigate the xenobiotic metabolism in saltwater crocodiles.

Modulating proposed electron transfer pathways in P450BM3 led to improved activity and coupling efficiency

Electrochemical in vitro reduction of P450 enzymes is a promising alternative to in vivo applications. Previously we presented three engineered P450_BM3 variants for aniline hydroxylation, equipped with a carbon nanotube binding-peptide (CNT-tag) for self-assembly on CNT electrodes. Compared to wildtype P450_BM3 the NADPH-dependent activity was enhanced, but the coupling efficiency remained low. For P450_BM3 Verma, Schwaneberg and Roccatano (2014, Biopolymers 101, 197-209) calculated putative electron transfer pathways (eTPs) by MD simulations. We hypothesised that knockouts of these transfer pathways would alter the coupling efficiency of the system. The results revealed no improved system for the electrically-driven P450s. For the NADPH-driven P450s, however, the most active eTP-mutant showed a 13-fold increased activity and a 32-fold elevated coupling efficiency using NADPH as reducing equivalent. This suggests an alternative principle of electron transport for the reduction by NADPH and an electrode, respectively. The work presents moreover a tool to improve the coupling and activity of P450s with non-natural substrates.

Elucidation of terpenoid metabolism in Scoparia dulcis by RNA-seq analysis

Scoparia dulcis biosynthesize bioactive diterpenes, such as scopadulcic acid B (SDB), which are known for their unique molecular skeleton. Although the biosynthesis of bioactive diterpenes is catalyzed by a sequence of class II and class I diterpene synthases (diTPSs), the mechanisms underlying this process are yet to be fully identified. To elucidate these biosynthetic machinery, we performed a high-throughput RNA-seq analysis, and de novo assembly of clean reads revealed 46,332 unique transcripts and 40,503 two unigenes. We found diTPSs genes including a putative syn-copalyl diphosphate synthase (SdCPS2) and two kaurene synthase-like (SdKSLs) genes. Besides them, total 79 full-length of cytochrome P450 (CYP450) genes were also discovered. The expression analyses showed selected CYP450s associated with their expression pattern of SdCPS2 and SdKSL1, suggesting that CYP450 candidates involved diterpene modification. SdCPS2 represents the first predicted gene to produce syn-copalyl diphosphate in dicots. In addition, SdKSL1 potentially contributes to the SDB biosynthetic pathway. Therefore, these identified genes associated with diterpene biosynthesis lead to the development of genetic engineering focus on diterpene metabolism in S. dulcis.

Semirational Protein Engineering of CYP153AM.aq. -CPRBM3 for Efficient Terminal Hydroxylation of Short- to Long-Chain Fatty Acids

The regioselective terminal hydroxylation of alkanes and fatty acids is of great interest in a variety of industrial applications, such as in cosmetics, in fine chemicals, and in the fragrance industry. The chemically challenging activation and oxidation of non-activated C-H bonds can be achieved with cytochrome P450 enzymes. CYP153AM.aq. -CPRBM3 is an artificial fusion construct consisting of the heme domain from Marinobacter aquaeolei and the reductase domain of CYP102A1 from Bacillus megaterium. It has the ability to hydroxylate medium- and long-chain fatty acids selectively at their terminal positions. However, the activity of this interesting P450 construct needs to be improved for applications in industrial processes. For this purpose, the design of mutant libraries including two consecutive steps of mutagenesis is demonstrated. Targeted positions and residues chosen for substitution were based on semi-rational protein design after creation of a homology model of the heme domain of CYP153AM.aq. , sequence alignments, and docking studies. Site-directed mutagenesis was the preferred method employed to address positions within the binding pocket, whereas diversity was created with the aid of a degenerate codon for amino acids located at the substrate entrance channel. Combining the successful variants led to the identification of a double variant-G307A/S233G-that showed alterations of one position within the binding pocket and one position located in the substrate access channel. This double variant showed twofold increased activity relative to the wild type for the terminal hydroxylation of medium-chain-length fatty acids. This variant furthermore showed improved activity towards short- and long-chain fatty acids and enhanced stability in the presence of higher concentrations of fatty acids.

Bifunctional CYP81AA proteins catalyse identical hydroxylations but alternative regioselective phenol couplings in plant xanthone biosynthesis

Xanthones are natural products present in plants and microorganisms. In plants, their biosynthesis starts with regioselective cyclization of 2,3′,4,6-tetrahydroxybenzophenone to either 1,3,5- or 1,3,7-trihydroxyxanthones, catalysed by cytochrome P450 (CYP) enzymes. Here we isolate and express CYP81AA-coding sequences from Hypericum calycinum and H. perforatum in yeast. Microsomes catalyse two consecutive reactions, that is, 3′-hydroxylation of 2,4,6-trihydroxybenzophenone and C–O phenol coupling of the resulting 2,3′,4,6-tetrahydroxybenzophenone. Relative to the inserted 3′-hydroxyl, the orthologues Hc/HpCYP81AA1 cyclize via the para position to form 1,3,7-trihydroxyxanthone, whereas the paralogue HpCYP81AA2 directs cyclization to the ortho position, yielding the isomeric 1,3,5-trihydroxyxanthone. Homology modelling and reciprocal mutagenesis reveal the impact of S375, L378 and A483 on controlling the regioselectivity of HpCYP81AA2, which is converted into HpCYP81AA1 by sextuple mutation. However, the reciprocal mutations in HpCYP81AA1 barely affect its regiospecificity. Product docking rationalizes the alternative C–O phenol coupling reactions. Our results help understand the machinery of bifunctional CYPs.

Xanthones are pharmacologically and biosynthetically intriguing compounds. Here, the authors identify two cytochrome P450 enzymes, which hydroxylate and cyclize the benzophenone precursor to either 1,3,7- or 1,3,5-trihydroxyxanthones, and pinpoint residues that determine the alternative regioselectivities.

Whole-cell microtiter plate screening assay for terminal hydroxylation of fatty acids by P450s

A readily available galactose oxidase (GOase) variant was used to develop a whole cell screening assay. This endpoint detection system was applied in a proof-of-concept approach by screening a focussed mutant library. This led to the discovery of the thus far most active P450 Marinobacter aquaeolei mutant catalysing the terminal hydroxylation of fatty acids.

Generate a bioactive natural product library by mining bacterial cytochrome P450 patterns

The increased number of annotated bacterial genomes provides a vast resource for genome mining. Several bacterial natural products with epoxide groups have been identified as pre-mRNA spliceosome inhibitors and antitumor compounds through genome mining. These epoxide-containing natural products feature a common biosynthetic characteristic that cytochrome P450s (CYPs) and its patterns such as epoxidases are employed in the tailoring reactions. The tailoring enzyme patterns are essential to both biological activities and structural diversity of natural products, and can be used for enzyme pattern-based genome mining. Recent development of direct cloning, heterologous expression, manipulation of the biosynthetic pathways and the CRISPR-CAS9 system have provided molecular biology tools to turn on or pull out nascent biosynthetic gene clusters to generate a microbial natural product library. This review focuses on a library of epoxide-containing natural products and their associated CYPs, with the intention to provide strategies on diversifying the structures of CYP-catalyzed bioactive natural products. It is conceivable that a library of diversified bioactive natural products will be created by pattern-based genome mining, direct cloning and heterologous expression as well as the genomic manipulation.

A Comprehensive Review of Aliphatic Hydrocarbon Biodegradation by Bacteria

Hydrocarbons are relatively recalcitrant compounds and are classified as high-priority pollutants. However, these compounds are slowly degraded by a large variety of microorganisms. Bacteria are able to degrade aliphatic saturated and unsaturated hydrocarbons via both aerobic and anaerobic pathways. Branched hydrocarbons and cyclic hydrocarbons are also degraded by bacteria. The aerobic bacteria use different types of oxygenases, including monooxygenase, cytochrome-dependent oxygenase and dioxygenase, to insert one or two atoms of oxygen into their targets. Anaerobic bacteria, on the other hand, employ a variety of simple organic and inorganic molecules, including sulphate, nitrate, carbonate and metals, for hydrocarbon oxidation.

Expanding the toolbox of organic chemists: directed evolution of P450 monooxygenases as catalysts in regio- and stereoselective oxidative hydroxylation

Cytochrome P450 enzymes (CYPs) have been used for more than six decades as catalysts for the CH-activating oxidative hydroxylation of organic compounds with formation of added-value products.

Comparative analysis of P450 signature motifs EXXR and CXG in the large and diverse kingdom of fungi: identification of evolutionarily conserved amino acid patterns characteristic of P450 family

Cytochrome P450 monooxygenases (P450s) are heme-thiolate proteins distributed across the biological kingdoms. P450s are catalytically versatile and play key roles in organisms primary and secondary metabolism. Identification of P450s across the biological kingdoms depends largely on the identification of two P450 signature motifs, EXXR and CXG, in the protein sequence. Once a putative protein has been identified as P450, it will be assigned to a family and subfamily based on the criteria that P450s within a family share more than 40% homology and members of subfamilies share more than 55% homology. However, to date, no evidence has been presented that can distinguish members of a P450 family. Here, for the first time we report the identification of EXXR- and CXG-motifs-based amino acid patterns that are characteristic of the P450 family. Analysis of P450 signature motifs in the under-explored fungal P450s from four different phyla, ascomycota, basidiomycota, zygomycota and chytridiomycota, indicated that the EXXR motif is highly variable and the CXG motif is somewhat variable. The amino acids threonine and leucine are preferred as second and third amino acids in the EXXR motif and proline and glycine are preferred as second and third amino acids in the CXG motif in fungal P450s. Analysis of 67 P450 families from biological kingdoms such as plants, animals, bacteria and fungi showed conservation of a set of amino acid patterns characteristic of a particular P450 family in EXXR and CXG motifs. This suggests that during the divergence of P450 families from a common ancestor these amino acids patterns evolve and are retained in each P450 family as a signature of that family. The role of amino acid patterns characteristic of a P450 family in the structural and/or functional aspects of members of the P450 family is a topic for future research.