Molecular Characterization of Specifically Active Recombinant Fused Enzymes Consisting of CYP3A4, NADPH-Cytochrome P450 Oxidoreductase, and Cytochrome b5

Functional Expression and Construction of a Self-Sufficient Cytochrome P450 Chimera for Efficient Steroidal C14α Hydroxylation in Escherichia coli

C14-functionalized steroids enabled diverse biological activities in anti-gonadotropin and anticancer therapy. However, access to C14-functionalized steroids was impeded by the deficiency of chemical synthetic methods. Recently, several membrane-bound fungal cytochrome P450s (CYPs) have been identified with steroid C14α-hydroxylation activity. However, the lack of efficient heterologous overexpression strategy hampered their further characterization and molecular engineering. In the present study, sequences of fungi-derived CYP genes encoding putative 14α-hydroxylase were selected and bioinformatically analyzed. Substitution of the N-terminal hydrophobic helix by a soluble maltose binding protein tag significantly enhanced the soluble expression level in Escherichia coli. A novel CYP originated from Bipolaris oryzae was discovered with high steroidal C14α-hydroxylation activity when coupled with the redox partner CPRlun. A catalytically self-sufficient chimeric CYP-CPR was built by intramolecular fusion, and the electronic transfer rate was improved. A coenzyme NADPH regeneration system was finally constructed by the co-expression of glucose dehydrogenase. The developed soluble multi-enzyme cascade biotransformation system supported the selective C14α-hydroxylation toward progesterone with a final titer of 34.54 mg/L, the highest level achieved in E. coli-based heterologous expression system. This study provides insightful ideas on the functional expression of fungi-derived CYPs and promises an efficient C14α-hydroxylation system for steroidal drugs through protein engineering.

Revealing evolution of tropane alkaloid biosynthesis by analyzing two genomes in the Solanaceae family

Tropane alkaloids (TAs) are widely distributed in the Solanaceae, while some important medicinal tropane alkaloids (mTAs), such as hyoscyamine and scopolamine, are restricted to certain species/tribes in this family. Little is known about the genomic basis and evolution of TAs biosynthesis and specialization in the Solanaceae. Here, we present chromosome-level genomes of two representative mTAs-producing species: Atropa belladonna and Datura stramonium. Our results reveal that the two species employ a conserved biosynthetic pathway to produce mTAs despite being distantly related within the nightshade family. A conserved gene cluster combined with gene duplication underlies the wide distribution of TAs in this family. We also provide evidence that branching genes leading to mTAs likely have evolved in early ancestral Solanaceae species but have been lost in most of the lineages, with A. belladonna and D. stramonium being exceptions. Furthermore, we identify a cytochrome P450 that modifies hyoscyamine into norhyoscyamine. Our results provide a genomic basis for evolutionary insights into the biosynthesis of TAs in the Solanaceae and will be useful for biotechnological production of mTAs via synthetic biology approaches.

Metabolic Engineering of Saccharomyces cerevisiae for Heterologous Carnosic Acid Production

Carnosic acid (CA), a phenolic tricyclic diterpene, has many biological effects, including anti-inflammatory, anticancer, antiobesity, and antidiabetic activities. In this study, an efficient biosynthetic pathway was constructed to produce CA in Saccharomyces cerevisiae. First, the CA precursor miltiradiene was synthesized, after which the CA production strain was constructed by integrating the genes encoding cytochrome P450 enzymes (P450s) and cytochrome P450 reductase (CPR) SmCPR. The CA titer was further increased by the coexpression of CYP76AH1 and SmCPR ∼t28SpCytb5 fusion proteins and the overexpression of different catalases to detoxify the hydrogen peroxide (H₂O₂). Finally, engineering of the endoplasmic reticulum and cofactor supply increased the CA titer to 24.65 mg/L in shake flasks and 75.18 mg/L in 5 L fed-batch fermentation. This study demonstrates that the ability of engineered yeast cells to synthesize CA can be improved through metabolic engineering and synthetic biology strategies, providing a theoretical basis for microbial synthesis of other diterpenoids.

Construction and optimization of Saccharomyces cerevisiae for synthesizing forskolin

Forskolin, one of the primary active metabolites of labdane-type diterpenoids, exhibits significant medicinal value, such as anticancer, antiasthmatic, and antihypertensive activities. In this study, we constructed a Saccharomyces cerevisiae cell factory that efficiently produced forskolin. First, a chassis strain that can accumulate 145.8 mg/L 13R-manoyl oxide (13R-MO), the critical precursor of forskolin, was constructed. Then, forskolin was produced by integrating CfCYP76AH15, CfCYP76AH11, CfCYP76AH16, ATR1, and CfACT1-8 into the 13R-MO chassis with a titer of 76.25 μg/L. We confirmed that cytochrome P450 enzymes (P450s) are the rate-limiting step by detecting intermediate metabolite accumulation. Forskolin production reached 759.42 μg/L by optimizing the adaptations between CfCYP76AHs, t66CfCPR, and t30AaCYB5. Moreover, multiple metabolic engineering strategies, including regulation of the target genes' copy numbers, amplification of the endoplasmic reticulum (ER) area, and cofactor metabolism enhancement, were implemented to enhance the metabolic flow to forskolin from 13R-MO, resulting in a final forskolin yield of 21.47 mg/L in shake flasks and 79.33 mg/L in a 5 L bioreactor. These promising results provide guidance for the synthesis of other natural terpenoids in S. cerevisiae, especially for those containing multiple P450s in their synthetic pathways. KEY POINTS: • The forskolin biosynthesis pathway was optimized from the perspective of system metabolism for the first time in S. cerevisiae. • The adaptation and optimization of CYP76AHs, t66CfCPR, and t30AaCYB5 promote forskolin accumulation, which can provide a reference for diterpenoids containing complex pathways, especially multiple P450s pathways. • The forskolin titer of 79.33 mg/L is the highest production currently reported and was achieved by fed-batch fermentation in a 5 L bioreactor.

Overproduction of medicinal ergot alkaloids based on a fungal platform

Privileged ergot alkaloids (EAs) produced by the fungal genus Claviceps are used to treat a wide range of diseases. However, their use and research have been hampered by the challenging genetic engineering of Claviceps. Here we systematically refactored and rationally engineered the EA biosynthetic pathway in heterologous host Aspergillus nidulans by using a Fungal-Yeast-Shuttle-Vector protocol. The obtained strains allowed the production of diverse EAs and related intermediates, including prechanoclavine (PCC, 333.8 mg/L), chanoclavine (CC, 241.0 mg/L), agroclavine (AC, 78.7 mg/L), and festuclavine (FC, 99.2 mg/L), etc. This fungal platform also enabled the access to the methyl-oxidized EAs (MOEAs), including elymoclavine (EC), lysergic acid (LA), dihydroelysergol (DHLG), and dihydrolysergic acid (DHLA), by overexpressing a P450 enzyme CloA. Furthermore, by optimizing the P450 electron transfer (ET) pathway and using multi-copy of cloA, the titers of EC and DHLG have been improved by 17.3- and 9.4-fold, respectively. Beyond our demonstration of A. nidulans as a robust platform for EA overproduction, our study offers a proof of concept for engineering the eukaryotic P450s-contained biosynthetic pathways in a filamentous fungal host.

Engineering cellular metabolite transport for biosynthesis of computationally predicted tropane alkaloid derivatives in yeast

Microbial biosynthesis of plant natural products (PNPs) can facilitate access to valuable medicinal compounds and derivatives. Such efforts are challenged by metabolite transport limitations, which arise when complex plant pathways distributed across organelles and tissues are reconstructed in unicellular hosts without concomitant transport machinery. We recently reported an engineered yeast platform for production of the tropane alkaloid (TA) drugs hyoscyamine and scopolamine, in which product accumulation is limited by vacuolar transport. Here, we demonstrate that alleviation of transport limitations at multiple steps in an engineered pathway enables increased production of TAs and screening of useful derivatives. We first show that supervised classifier models trained on a tissue-delineated transcriptome from the TA-producing plant Atropa belladonna can predict TA transporters with greater efficacy than conventional regression- and clustering-based approaches. We demonstrate that two of the identified transporters, AbPUP1 and AbLP1, increase TA production in engineered yeast by facilitating vacuolar export and cellular reuptake of littorine and hyoscyamine. We incorporate four different plant transporters, cofactor regeneration mechanisms, and optimized growth conditions into our yeast platform to achieve improvements in de novo hyoscyamine and scopolamine production of over 100-fold (480 μg/L) and 7-fold (172 μg/L). Finally, we leverage computational tools for biosynthetic pathway prediction to produce two different classes of TA derivatives, nortropane alkaloids and tropane N-oxides, from simple precursors. Our work highlights the importance of cellular transport optimization in recapitulating complex PNP biosyntheses in microbial hosts and illustrates the utility of computational methods for gene discovery and expansion of heterologous biosynthetic diversity.

Cytochrome b5 plays a dual role in the reaction cycle of cytochrome P450 3A4 during oxidation of the anticancer drug ellipticine

Abstract

Ellipticine is an anticancer agent that forms covalent DNA adducts after enzymatic activation by cytochrome P450 (CYP) enzymes, mainly by CYP3A4. This process is one of the most important ellipticine DNA-damaging mechanisms for its antitumor action. Here, we investigated the efficiencies of human hepatic microsomes and human recombinant CYP3A4 expressed with its reductase, NADPH:CYP oxidoreductase (POR), NADH:cytochrome b₅ reductase and/or cytochrome b₅ in Supersomes™ to oxidize this drug. We also evaluated the effectiveness of coenzymes of two of the microsomal reductases, NADPH as a coenzyme of POR, and NADH as a coenzyme of NADH:cytochrome b₅ reductase, to mediate ellipticine oxidation in these enzyme systems. Using HPLC analysis we detected up to five ellipticine metabolites, which were formed by human hepatic microsomes and human CYP3A4 in the presence of NADPH or NADH. Among ellipticine metabolites, 9-hydroxy-, 12-hydroxy-, and 13-hydroxyellipticine were formed by hepatic microsomes as the major metabolites, while 7-hydroxyellipticine and the ellipticine N²-oxide were the minor ones. Human CYP3A4 in Supersomes™ generated only three metabolic products, 9-hydroxy-, 12-hydroxy-, and 13-hydroxyellipticine. Using the ³²P-postlabeling method two ellipticine-derived DNA adducts were generated by microsomes and the CYP3A4-Supersome system, both in the presence of NADPH and NADH. These adducts were derived from the reaction of 13-hydroxy- and 12-hydroxyellipticine with deoxyguanosine in DNA. In the presence of NADPH or NADH, cytochrome b₅ stimulated the CYP3A4-mediated oxidation of ellipticine, but the stimulation effect differed for individual ellipticine metabolites. This heme protein also stimulated the formation of both ellipticine-DNA adducts. The results demonstrate that cytochrome b₅ plays a dual role in the CYP3A4-catalyzed oxidation of ellipticine: (1) cytochrome b₅ mediates CYP3A4 catalytic activities by donating the first and second electron to this enzyme in its catalytic cycle, indicating that NADH:cytochrome b₅ reductase can substitute NADPH-dependent POR in this enzymatic reaction and (2) cytochrome b₅ can act as an allosteric modifier of the CYP3A4 oxygenase.

Graphical abstract

Mechanistic basis of electron transfer to cytochromes p450 by natural redox partners and artificial donor constructs

Cytochromes P450 (P450s) are hemoproteins catalyzing oxidative biotransformation of a vast array of natural and xenobiotic compounds. Reducing equivalents required for dioxygen cleavage and substrate hydroxylation originate from different redox partners including diflavin reductases, flavodoxins, ferredoxins and phthalate dioxygenase reductase (PDR)-type proteins. Accordingly, circumstantial analysis of structural and physicochemical features governing donor-acceptor recognition and electron transfer poses an intriguing challenge. Thus, conformational flexibility reflected by togging between closed and open states of solvent exposed patches on the redox components was shown to be instrumental to steered electron transmission. Here, the membrane-interactive tails of the P450 enzymes and donor proteins were recognized to be crucial to proper orientation toward each other of surface sites on the redox modules steering functional coupling. Also, mobile electron shuttling may come into play. While charge-pairing mechanisms are of primary importance in attraction and complexation of the redox partners, hydrophobic and van der Waals cohesion forces play a minor role in docking events. Due to catalytic plasticity of P450 enzymes, there is considerable promise in biotechnological applications. Here, deeper insight into the mechanistic basis of the redox machinery will permit optimization of redox processes via directed evolution and DNA shuffling. Thus, creation of hybrid systems by fusion of the modified heme domain of P450s with proteinaceous electron carriers helps obviate the tedious reconstitution procedure and induces novel activities. Also, P450-based amperometric biosensors may open new vistas in pharmaceutical and clinical implementation and environmental monitoring.

Rational design of a tripartite fusion protein of heparinase I enables one-step affinity purification and real-time activity detection

Enzymatic degradation of heparin has great potential as an ecological and specific way to produce low molecular weight heparin. However, the commercial use of heparinase I (HepA), one of the most important heparin lyases, has been hampered by low productivity and poor thermostability. Fusion with green fluorescent protein (GFP) or maltose-binding protein (MBP) has shown potential in facilitating the industrial use of HepA. Thus, tripartite fusion of GFP, MBP and HepA would be a promising approach. Therefore, in the present study, the tripartite fusion strategy was systematically studied, mainly focusing on the fusion order and the linker sequence, to obtain a fusion protein offering one-step purification and real-time detection of HepA activity by fluorescence as well as high HepA activity and thermostability. Our results show that fusion order is important for MBP binding affinity and HepA activity, while the linker sequences at domain junctions have significant effects on protein expression level, HepA activity and thermostability as well as GFP fluorescence. The best tripartite fusion was identified as MBP-(EAAAK)(3)-GFP-(GGGGS)(3)-HepA, which shows potential to facilitate the production of HepA and its application in industrial preparation of low molecular weight heparin.

Two novel alkane hydroxylase-rubredoxin fusion genes isolated from a Dietzia bacterium and the functions of fused rubredoxin domains in long-chain n-alkane degradation

Two alkane hydroxylase-rubredoxin fusion gene homologs (alkW1 and alkW2) were cloned from a Dietzia strain, designated DQ12-45-1b, which can grow on crude oil and n-alkanes ranging in length from 6 to 40 carbon atoms as sole carbon sources. Both AlkW1 and AlkW2 have an integral-membrane alkane monooxygenase (AlkB) conserved domain and a rubredoxin (Rd) conserved domain which are fused together. Phylogenetic analysis showed that these two AlkB-fused Rd domains formed a novel third cluster with all the Rds from the alkane hydroxylase-rubredoxin fusion gene clusters in Gram-positive bacteria and that this third cluster was distant from the known AlkG1- and AlkG2-type Rds. Expression of the alkW1 gene in DQ12-45-1b was induced when cells were grown on C(8) to C(32) n-alkanes as sole carbon sources, but expression of the alkW2 gene was not detected. Functional heterologous expression in an alkB deletion mutant of Pseudomonas fluorescens KOB2Δ1 suggested the alkW1 could restore the growth of KOB2Δ1 on C(14) and C(16) n-alkanes and induce faster growth on C(18) to C(32) n-alkanes than alkW1ΔRd, the Rd domain deletion mutant gene of alkW1, which also caused faster growth than KOB2Δ1 itself. In addition, the artificial fusion of AlkB from the Gram-negative P. fluorescens CHA0 and the Rds from both Gram-negative P. fluorescens CHA0 and Gram-positive Dietzia sp. DQ12-45-1b significantly increased the degradation of C(32) alkane compared to that seen with AlkB itself. In conclusion, the alkW1 gene cloned from Dietzia species encoded an alkane hydroxylase which increased growth on and degradation of n-alkanes up to C(32) in length, with its fused rubredoxin domain being necessary to maintain the functions. In addition, the fusion of alkane hydroxylase and rubredoxin genes from both Gram-positive and -negative bacteria can increase the degradation of long-chain n-alkanes (such as C(32)) in the Gram-negative bacterium.

Structure of the open conformation of a functional chimeric NADPH cytochrome P450 reductase

Two catalytic domains, bearing FMN and FAD cofactors, joined by a connecting domain, compose the core of the NADPH cytochrome P450 reductase (CPR). The FMN domain of CPR mediates electron shuttling from the FAD domain to cytochromes P450. Together, both enzymes form the main mixed-function oxidase system that participates in the metabolism of endo- and xenobiotic compounds in mammals. Available CPR structures show a closed conformation, with the two cofactors in tight proximity, which is consistent with FAD-to-FMN, but not FMN-to-P450, electron transfer. Here, we report the 2.5 A resolution crystal structure of a functionally competent yeast-human chimeric CPR in an open conformation, compatible with FMN-to-P450 electron transfer. Comparison with closed structures shows a major conformational change separating the FMN and FAD cofactors from 86 A.