Genome-wide identification, annotation and characterization of novel thermostable cytochrome P450 monooxygenases from the thermophilic biomass-degrading fungi Thielavia terrestris and Myceliophthora thermophila

Influence of Heat Treatment on Tea Polyphenols and Their Impact on Improving Heat Tolerance in Drosophila melanogaster

This study investigated the potential mechanism of action of tea polyphenols (TPs), one of the major active ingredients in tea, to enhance heat resistance in Drosophila and the attenuating effect of heat treatment of TPs on their efficacy. The results showed that TPs were able to prolong the average survival time of Drosophila under high-temperature stress (p < 0.05), but the effect of TPs in prolonging the survival time of Drosophila melanogaster was significantly reduced (p < 0.05) with increasing TP heat-treatment time until it disappeared. The composition of TPs changed after heat treatment. It was also shown that the weakening of the effect of TPs in improving the heat tolerance of Drosophila was related to the decrease in the content of catechins and phenolic acids in their fractions as well as with the increase in the content of laccase. Transcriptomic analysis showed that the effect of TPs on heat tolerance in Drosophila melanogaster was closely related to the longevity regulation pathway, the neuroactive ligand-receptor interaction signaling pathway, and the drug metabolism-cytochrome P450 pathway. Metabolomics analysis showed that the effect of TP intervention in improving the body's heat tolerance was mainly related to amino acid metabolism and energy metabolism. However, thermal processing weakened the relevance of these transcriptomes and metabolomes. The present study reveals the mechanism of action by which heat-treated TPs affect the body's heat tolerance, which is important for the development and utilization of the heat-protection function of tea.

A Novel Thermostable Cytochrome P450 from Sequence-Based Metagenomics of Binh Chau Hot Spring as a Promising Catalyst for Testosterone Conversion

Biotechnological applications of cytochromes P450 show difficulties, such as low activity, thermal and/or solvent instability, narrow substrate specificity and redox partner dependence. In an attempt to overcome these limitations, an exploitation of novel thermophilic P450 enzymes from nature via uncultured approaches is desirable due to their great advantages that can resolve nearly all mentioned impediments. From the metagenomics library of the Binh Chau hot spring, an open reading frame (ORF) encoding a thermostable cytochrome P450—designated as P450-T3—which shared 66.6% amino acid sequence identity with CYP109C2 of Sorangium cellulosum So ce56 was selected for further identification and characterization. The ORF was synthesized artificially and heterologously expressed in Escherichia coli C43(DE3) using the pET17b system. The purified enzyme had a molecular weight of approximately 43 kDa. The melting temperature of the purified enzyme was 76.2 °C and its apparent half-life at 60 °C was 38.7 min. Redox partner screening revealed that P450-T3 was reduced well by the mammalian AdR-Adx4-108 and the yeast Arh1-Etp1 redox partners. Lauric acid, palmitic acid, embelin, retinoic acid (all-trans) and retinoic acid (13-cis) demonstrated binding to P450-T3. Interestingly, P450-T3 also bound and converted testosterone. Overall, P450-T3 might become a good candidate for biocatalytic applications on a larger scale.

Computer-aided search for a cold-active cellobiose 2-epimerase

Cellobiose 2-epimerase (CE) is a promising industrial enzyme that can catalyze bioconversion of lactose to its high-value derivatives, namely epilactose and lactulose. A need exists in the dairy industry to catalyze lactose bioconversions at low temperatures to avoid microbial growth. We focused on the discovery of cold-active CE in this study. A genome mining method based on computational prediction was used to screen the potential genes encoding cold-active enzymes. The CE-encoding gene from Roseburia intestinalis, with a predicted high structural flexibility, was expressed heterologously in Escherichia coli. The catalytic property of the recombinant enzyme was extensively studied. The optimum temperature and pH of the enzyme were 45°C and 7.0, respectively. The specific activity of this enzyme to catalyze conversion of lactose to epilactose was measured to be 77.3 ± 1.6 U/mg. The kinetic parameters, including turnover number (k_cat), Michaelis constant (K_m), and catalytic efficiency (k_cat/K_m) using lactose as a substrate were 117.0 ± 7.7 s^-1, 429.9 ± 57.3 mM, and 0.27 mM^-1s^-1, respectively. In situ production of epilactose was carried out at 8°C: 20.9% of 68.4 g/L lactose was converted into epilactose in 4 h using 0.02 mg/mL (1.5 U/mL, measured at 45°C) of recombinant enzyme. The enzyme discovered by this in silico method is suitable for low-temperature applications.

Phylogenetic re-evaluation of Thielavia with the introduction of a new family Podosporaceae

The genus Thielavia is morphologically defined by having non-ostiolate ascomata with a thin peridium composed of textura epidermoidea, and smooth, single-celled, pigmented ascospores with one germ pore. Thielavia is typified with Th. basicola that grows in close association with a hyphomycete which was traditionally identified as Thielaviopsis basicola. Besides Th. basicola exhibiting the mycoparasitic nature, the majority of the described Thielavia species are from soil, and some have economic and ecological importance. Unfortunately, no living type material of Th. basicola exists, hindering a proper understanding of the classification of Thielavia. Therefore, Thielavia basicola was neotypified by material of a mycoparasite presenting the same ecology and morphology as described in the original description. We subsequently performed a multi-gene phylogenetic analyses (rpb2, tub2, ITS and LSU) to resolve the phylogenetic relationships of the species currently recognised in Thielavia. Our results demonstrate that Thielavia is highly polyphyletic, being related to three family-level lineages in two orders. The redefined genus Thielavia is restricted to its type species, Th. basicola, which belongs to the Ceratostomataceae (Melanosporales) and its host is demonstrated to be Berkeleyomyces rouxiae, one of the two species in the "Thielaviopsis basicola" species complex. The new family Podosporaceae is sister to the Chaetomiaceae in the Sordariales and accommodates the re-defined genera Podospora, Trangularia and Cladorrhinum, with the last genus including two former Thielavia species (Th. hyalocarpa and Th. intermedia). This family also includes the genetic model species Podospora anserina, which was combined in Triangularia (as Triangularia anserina). The remaining Thielavia species fall in ten unrelated clades in the Chaetomiaceae, leading to the proposal of nine new genera (Carteria, Chrysanthotrichum, Condenascus, Hyalosphaerella, Microthielavia, Parathielavia, Pseudothielavia, Stolonocarpus and Thermothielavioides). The genus Canariomyces is transferred from Microascaceae (Microascales) to Chaetomiaceae based on its type species Can. notabilis. Canariomyces is closely related to the human-pathogenic genus Madurella, and includes three thielavia-like species and one novel species. Three monotypic genera with a chaetomium-like morph (Brachychaeta, Chrysocorona and Floropilus) are introduced to better resolve the Chaetomiaceae and the thielavia-like species in the family. Chrysocorona lucknowensis and Brachychaeta variospora are closely related to Acrophialophora and three newly introduced genera containing thielavia-like species; Floropilus chiversii is closely related to the industrially important and thermophilic species Thermothielavioides terrestris (syn. Th. terrestris). This study shows that the thielavia-like morph is a homoplastic form that originates from several separate evolutionary events. Furthermore, our results provide new insights into the taxonomy of Sordariales and the polyphyletic Lasiosphaeriaceae.

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.

In silico analysis of cytochrome P450 monooxygenases in chronic granulomatous infectious fungus Sporothrix schenckii: Special focus on CYP51

Sporotrichosis is an emerging chronic, granulomatous, subcutaneous, mycotic infection caused by Sporothrix species. Sporotrichosis is treated with the azole drug itraconazole as ketoconazole is ineffective. It is a well-known fact that azole drugs act by inhibiting cytochrome P450 monooxygenases (P450s), heme-thiolate proteins. To date, nothing is known about P450s in Sporothrix schenckii and the molecular basis of its resistance to ketoconazole. Here we present genome-wide identification, annotation, phylogenetic analysis and comprehensive P450 family-level comparative analysis of S. schenckii P450s with pathogenic fungi P450s, along with a rationale for ketoconazole resistance by S. schenckii based on in silico structural analysis of CYP51. Genome data-mining of S. schenckii revealed 40 P450s in its genome that can be grouped into 32 P450 families and 39 P450 subfamilies. Comprehensive comparative analysis of P450s revealed that S. schenckii shares 11 P450 families with plant pathogenic fungi and has three unique P450 families: CYP5077, CYP5386 and CYP5696 (novel family). Among P450s, CYP51, the main target of azole drugs was also found in S. schenckii. 3D modeling of S. schenckii CYP51 revealed the presence of characteristic P450 motifs with exceptionally large reductase interaction site 2. In silico analysis revealed number of mutations that can be associated with ketoconazole resistance, especially at the channel entrance to the active site. One of possible reason for better stabilization of itraconazole, compared to ketoconazole, is that the more extended molecule of itraconazole may form a hydrogen bond with ASN-230. This in turn may explain its effectiveness against S. schenckii vis-a-vis resistant to ketoconazole. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Expression, Purification, and Biochemical Characterization of the Flavocytochrome P450 CYP505A30 from Myceliophthora thermophila

The cytochrome P450/P450 reductase fusion enzyme CYP505A30 from the thermophilic fungus Myceliophthora thermophila and its heme (P450) domain were expressed in Escherichia coli and purified using affinity, ion exchange, and size exclusion chromatography. CYP505A30 binds straight chain fatty acids (from ∼C10 to C20), with highest affinity for tridecanoic acid (K_D = 2.7 μM). Reduced nicotinamide adenine dinucleotide phosphate is the preferred reductant for CYP505A30 (K_M = 3.1 μM compared to 330 μM for reduced nicotinamide adenine dinucleotide in cytochrome c reduction). Electron paramagnetic resonance confirmed cysteine thiolate coordination of heme iron in CYP505A30 and its heme domain. Redox potentiometry revealed an unusually positive midpoint potential for reduction of the flavin adenine dinucleotide and flavin mononucleotide cofactors (E⁰' ∼ -118 mV), and a large increase in the CYP505A30 heme domain Fe^III/Fe^II redox couple (ca. 230 mV) on binding arachidonic acid substrate. This switch brings the ferric heme iron potential into the same range as that of the reductase flavins. Multiangle laser light scattering analysis revealed CYP505A30's ability to dimerize, whereas the heme domain is monomeric. These data suggest CYP505A30 may function catalytically as a dimer (as described for Bacillus megaterium P450 BM3), and that binding interactions between CYP505A30 heme domains are not required for dimer formation. CYP505A30 catalyzed hydroxylation of straight chain fatty acids at the ω-1 to ω-3 positions, with a strong preference for ω-1 over ω-3 hydroxylation in the oxidation of dodecanoic and tetradecanoic acids (88 vs 2% products and 63 vs 9% products, respectively). CYP505A30 has important structural and catalytic similarities to P450 BM3 but distinct regioselectivity of lipid substrate oxidation with potential biotechnological applications.

Determinants of thermostability in the cytochrome P450 fold

Cytochromes P450 are found throughout the biosphere in a wide range of environments, serving a multitude of physiological functions. The ubiquity of the P450 fold suggests that it has been co-opted by evolution many times, and likely presents a useful compromise between structural stability and conformational flexibility. The diversity of substrates metabolized and reactions catalyzed by P450s makes them attractive starting materials for use as biocatalysts of commercially useful reactions. However, process conditions impose different requirements on enzymes to those in which they have evolved naturally. Most natural environments are relatively mild, and therefore most P450s have not been selected in Nature for the ability to withstand temperatures above ~40°C, yet industrial processes frequently require extended incubations at much higher temperatures. Thus, there has been considerable interest and effort invested in finding or engineering thermostable P450 systems. Numerous P450s have now been identified in thermophilic organisms and analysis of their structures provides information as to mechanisms by which the P450 fold can be stabilized. In addition, protein engineering, particularly by directed or artificial evolution, has revealed mutations that serve to stabilize particular mesophilic enzymes of interest. Here we review the current understanding of thermostability as it applies to the P450 fold, gleaned from the analysis of P450s characterized from thermophilic organisms and the parallel engineering of mesophilic forms for greater thermostability. We then present a perspective on how this information might be used to design stable P450 enzymes for industrial application. This article is part of a Special Issue entitled: Cytochrome P450 biodiversity and biotechnology, edited by Erika Plettner, Gianfranco Gilardi, Luet Wong, Vlada Urlacher, Jared Goldstone.

Molecular evolutionary dynamics of cytochrome P450 monooxygenases across kingdoms: Special focus on mycobacterial P450s

Since the initial identification of cytochrome P450 monooxygenases (CYPs/P450s), great progress has been made in understanding their structure-function relationship, diversity and application in producing compounds beneficial to humans. However, the molecular evolution of P450s in terms of their dynamics both at protein and DNA levels and functional conservation across kingdoms still needs investigation. In this study, we analyzed 17 598 P450s belonging to 113 P450 families (bacteria -42; fungi -19; plant -28; animal -22; plant and animal -1 and common P450 family -1) and found highly conserved and rapidly evolving P450 families. Results suggested that bacterial P450s, particularly P450s belonging to mycobacteria, are highly conserved both at protein and DNA levels. Mycobacteria possess the highest P450 diversity percentage compared to other microbes and have a high coverage of P450s (≥1%) in their genomes, as found in fungi and plants. Phylogenetic and functional analyses revealed the functional conservation of P450s despite belonging to different biological kingdoms, suggesting the adherence of P450s to their innate function such as their involvement in either generation or oxidation of steroids and structurally related molecules, fatty acids and terpenoids. This study's results offer new understanding of the dynamic structural nature of P450s.

Versatile biocatalysis of fungal cytochrome P450 monooxygenases

Cytochrome P450 (CYP) monooxygenases, the nature’s most versatile biological catalysts have unique ability to catalyse regio-, chemo-, and stereospecific oxidation of a wide range of substrates under mild reaction conditions, thereby addressing a significant challenge in chemocatalysis. Though CYP enzymes are ubiquitous in all biological kingdoms, the divergence of CYPs in fungal kingdom is manifold. The CYP enzymes play pivotal roles in various fungal metabolisms starting from housekeeping biochemical reactions, detoxification of chemicals, and adaptation to hostile surroundings. Considering the versatile catalytic potentials, fungal CYPs has gained wide range of attraction among researchers and various remarkable strategies have been accomplished to enhance their biocatalytic properties. Numerous fungal CYPs with multispecialty features have been identified and the number of characterized fungal CYPs is constantly increasing. Literature reveals ample reviews on mammalian, plant and bacterial CYPs, however, modest reports on fungal CYPs urges a comprehensive review highlighting their novel catalytic potentials and functional significances. In this review, we focus on the diversification and functional diversity of fungal CYPs and recapitulate their unique and versatile biocatalytic properties. As such, this review emphasizes the crucial issues of fungal CYP systems, and the factors influencing efficient biocatalysis.

Genome-Wide Annotation and Comparative Analysis of Cytochrome P450 Monooxygenases in Basidiomycete Biotrophic Plant Pathogens

Fungi are an exceptional source of diverse and novel cytochrome P450 monooxygenases (P450s), heme-thiolate proteins, with catalytic versatility. Agaricomycotina saprophytes have yielded most of the available information on basidiomycete P450s. This resulted in observing similar P450 family types in basidiomycetes with few differences in P450 families among Agaricomycotina saprophytes. The present study demonstrated the presence of unique P450 family patterns in basidiomycete biotrophic plant pathogens that could possibly have originated from the adaptation of these species to different ecological niches (host influence). Systematic analysis of P450s in basidiomycete biotrophic plant pathogens belonging to three different orders, Agaricomycotina (Armillaria mellea), Pucciniomycotina (Melampsora laricis-populina, M. lini, Mixia osmundae and Puccinia graminis) and Ustilaginomycotina (Ustilago maydis, Sporisorium reilianum and Tilletiaria anomala), revealed the presence of numerous putative P450s ranging from 267 (A. mellea) to 14 (M. osmundae). Analysis of P450 families revealed the presence of 41 new P450 families and 27 new P450 subfamilies in these biotrophic plant pathogens. Order-level comparison of P450 families between biotrophic plant pathogens revealed the presence of unique P450 family patterns in these organisms, possibly reflecting the characteristics of their order. Further comparison of P450 families with basidiomycete non-pathogens confirmed that biotrophic plant pathogens harbour the unique P450 families in their genomes. The CYP63, CYP5037, CYP5136, CYP5137 and CYP5341 P450 families were expanded in A. mellea when compared to other Agaricomycotina saprophytes and the CYP5221 and CYP5233 P450 families in P. graminis and M. laricis-populina. The present study revealed that expansion of these P450 families is due to paralogous evolution of member P450s. The presence of unique P450 families in these organisms serves as evidence of how a host/ecological niche can influence shaping the P450 content of an organism. The present study initiates our understanding of P450 family patterns in basidiomycete biotrophic plant pathogens.

Diversity and evolution of cytochrome P450 monooxygenases in Oomycetes

Cytochrome P450 monooxygenases (P450s) are heme-thiolate proteins whose role as drug targets against pathogens, as well as in valuable chemical production and bioremediation, has been explored. In this study we performed comprehensive comparative analysis of P450s in 13 newly explored oomycete pathogens. Three hundred and fifty-six P450s were found in oomycetes. These P450s were grouped into 15 P450 families and 84 P450 subfamilies. Among these, nine P450 families and 31 P450 subfamilies were newly found in oomycetes. Research revealed that oomycetes belonging to different orders contain distinct P450 families and subfamilies in their genomes. Evolutionary analysis and sequence homology data revealed P450 family blooms in oomycetes. Tandem arrangement of a large number of P450s belonging to the same family indicated that P450 family blooming is possibly due to its members’ duplications. A unique combination of amino acid patterns was observed at EXXR and CXG motifs for the P450 families CYP5014, CYP5015 and CYP5017. A novel P450 fusion protein (CYP5619 family) with an N-terminal P450 domain fused to a heme peroxidase/dioxygenase domain was discovered in Saprolegnia declina. Oomycete P450 patterns suggested host influence in shaping their P450 content. This manuscript serves as reference for future P450 annotations in newly explored oomycetes.

Cytochrome P450 monooxygenase CYP53 family in fungi: comparative structural and evolutionary analysis and its role as a common alternative anti-fungal drug target

Cytochrome P450 monooxygenases (CYPs/P450s) are heme-thiolate proteins whose role as a drug target against pathogenic microbes has been explored because of their stereo- and regio-specific oxidation activity. We aimed to assess the CYP53 family's role as a common alternative drug target against animal (including human) and plant pathogenic fungi and its role in fungal-mediated wood degradation. Genome-wide analysis of fungal species revealed the presence of CYP53 members in ascomycetes and basidiomycetes. Basidiomycetes had a higher number of CYP53 members in their genomes than ascomycetes. Only two CYP53 subfamilies were found in ascomycetes and six subfamilies in basidiomycetes, suggesting that during the divergence of phyla ascomycetes lost CYP53 P450s. According to phylogenetic and gene-structure analysis, enrichment of CYP53 P450s in basidiomycetes occurred due to the extensive duplication of CYP53 P450s in their genomes. Numerous amino acids (103) were found to be conserved in the ascomycetes CYP53 P450s, against only seven in basidiomycetes CYP53 P450s. 3D-modelling and active-site cavity mapping data revealed that the ascomycetes CYP53 P450s have a highly conserved protein structure whereby 78% amino acids in the active-site cavity were found to be conserved. Because of this rigid nature of ascomycetes CYP53 P450s' active site cavity, any inhibitor directed against this P450 family can serve as a common anti-fungal drug target, particularly toward pathogenic ascomycetes. The dynamic nature of basidiomycetes CYP53 P450s at a gene and protein level indicates that these P450s are destined to acquire novel functions. Functional analysis of CYP53 P450s strongly supported our hypothesis that the ascomycetes CYP53 P450s ability is limited for detoxification of toxic molecules, whereas basidiomycetes CYP53 P450s play an additional role, i.e. involvement in degradation of wood and its derived components. This study is the first report on genome-wide comparative structural (gene and protein structure-level) and evolutionary analysis of a fungal P450 family.

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.