Cloning, functional expression, and subcellular localization of multiple NADPH-cytochrome P450 reductases from hybrid poplar

Mining and functional characterization of NADPH-cytochrome P450 reductases of the DNJ biosynthetic pathway in mulberry leaves

BACKGROUND

1-Deoxynojirimycin (DNJ), the main active ingredient in mulberry leaves, with wide applications in the medicine and food industries due to its significant functions in lowering blood sugar, and lipids, and combating viral infections. Cytochrome P450 is a key enzyme for DNJ biosynthesis, its activity depends on the electron supply of NADPH-cytochrome P450 reductases (CPRs). However, the gene for MaCPRs in mulberry leaves remains unknown.

RESULTS

In this study, we successfully cloned and functionally characterized two key genes, MaCPR1 and MaCPR2, based on the transcriptional profile of mulberry leaves. The MaCPR1 gene comprised 2064 bp, with its open reading frame (ORF) encoding 687 amino acids. The MaCPR2 gene comprised 2148 bp, and its ORF encoding 715 amino acids. The phylogenetic tree indicates that MaCPR1 and MaCPR2 belong to Class I and Class II, respectively. In vitro, we found that the recombinant enzymes MaCPR2 protein could reduce cytochrome c and ferricyanide using NADPH as an electron donor, while MaCPR1 did not. In yeast, heterologous co-expression indicates that MaCPR2 delivers electrons to MaC3'H hydroxylase, a key enzyme catalyzing the production of chlorogenic acid from 3-O-p-coumaroylquinic acid.

CONCLUSIONS

These findings highlight the orchestration of hydroxylation process mediated by MaCPR2 during the biosynthesis of secondary metabolite biosynthesis in mulberry leaves. These results provided a foundational understanding for fully elucidating the DNJ biosynthetic pathway within mulberry leaves.

Cloning and Functional Characterization of NADPH-Cytochrome P450 Reductases in Aconitum vilmorinianum

Diterpenoid alkaloids (DAs) are major pharmacologically active ingredients of Aconitum vilmorinianum, an important medicinal plant. Cytochrome P450 monooxygenases (P450s) are involved in the DA biosynthetic pathway, and the electron transfer reaction of NADPH-cytochrome P450 reductase (CPR) with P450 is the rate-limiting step of the P450 redox reaction. Here, we identified and characterized two homologs of CPR from Aconitum vilmorinianum. The open reading frames of AvCPR1 and AvCPR2 were found to be 2103 and 2100 bp, encoding 700 and 699 amino acid residues, respectively. Phylogenetic analysis characterized both AvCPR1 and AvCPR2 as class II CPRs. Cytochrome c and ferricyanide could be reduced with the recombinant proteins of AvCPR1 and AvCPR2. Both AvCPR1 and AvCPR2 were expressed in the roots, stems, leaves, and flowers of A. vilmorinianum. The expression levels of AvCPR1 and AvCPR2 were significantly increased in response to methyl jasmonate (MeJA) treatment. The yeasts co-expressing AvCPR1/AvCPR2/SmCPR1 and CYP76AH1 all produced ferruginol, indicating that AvCPR1 and AvCPR2 can transfer electrons to CYP76AH1 in the same manner as SmCPR1. Docking analysis confirmed the experimentally deduced functional activities of AvCPR1 and AvCPR2 for FMN, FAD, and NADPH. The functional characterization of AvCPRs will be helpful in disclosing molecular mechanisms relating to the biosynthesis of diterpene alkaloids in A. vilmorinianum.

Biochemical Alterations in Triticale Seedlings Pretreated with Selective Herbicide and Subjected to Drought or Waterlogging Stress

Waterlogging and drought disrupt crop development and productivity. Triticale is known to be relatively tolerant to different stress factors. In natural conditions, plants are rather subjected to multiple environmental factors. Serrate® (Syngenta) is a systemic selective herbicide suitable for cereal crops such as triticale and wheat to restrain annual grass and broadleaf weeds. Triticale (×Triticosecale Wittm., cv. Rozhen) was grown as soil culture under controlled conditions. Seventeen-day-old plantlets were leaf sprayed with Serrate®. The water stress (drought or waterlogging) was applied after 72 h for 7 days, and then the seedlings were left for recovery. The herbicide does not provoke sharp alterations in the antioxidant state (stress markers level, and antioxidant and xenobiotic-detoxifying enzymes activity). The water stresses and combined treatments enhanced significantly the content of stress markers (malondialdehyde, proline, hydrogen peroxide), non-enzymatic (total phenolics and thiol groups-containing compounds), and enzymatic (activities of superoxide dismutase, catalase, guaiacol peroxidase, glutathione reductase) antioxidants, and xenobiotic-detoxifying enzymes (activities of glutathione S-transferase, NADPH:cytochrome P450 reductase, NADH:cytochrome b5 reductase). These effects were more severely expressed after the drought stress, suggesting that this cultivar is more tolerant to waterlogging than to drought stress.

An alternative angiosperm DGAT1 topology and potential motifs in the N-terminus

The highly variable cytoplasmic N-terminus of the plant diacylglycerol acyltransferase 1 (DGAT1) has been shown to have roles in oligomerization as well as allostery; however, the biological significance of the variation within this region is not understood. Comparing the coding sequences over the variable N-termini revealed the Poaceae DGAT1s contain relatively high GC compositional gradients as well as numerous direct and inverted repeats in this region. Using a variety of reciprocal chimeric DGAT1s from angiosperms we show that related N-termini had similar effects (positive or negative) on the accumulation of the recombinant protein in Saccharomyces cerevisiae. When expressed in Camelina sativa seeds the recombinant proteins of specific chimeras elevated total lipid content of the seeds as well as increased seed size. In addition, we combine N- and C-terminal as well as internal tags with high pH membrane reformation, protease protection and differential permeabilization. This led us to conclude the C-terminus is in the ER lumen; this contradicts earlier reports of the cytoplasmic location of plant DGAT1 C-termini.

Subcellular Localization Prediction of Human Proteins Using Multifeature Selection Methods

Subcellular localization attempts to assign proteins to one of the cell compartments that performs specific biological functions. Finding the link between proteins, biological functions, and subcellular localization is an effective way to investigate the general organization of living cells in a systematic manner. However, determining the subcellular localization of proteins by traditional experimental approaches is difficult. Here, protein–protein interaction networks, functional enrichment on gene ontology and pathway, and a set of proteins having confirmed subcellular localization were applied to build prediction models for human protein subcellular localizations. To build an effective predictive model, we employed a variety of robust machine learning algorithms, including Boruta feature selection, minimum redundancy maximum relevance, Monte Carlo feature selection, and LightGBM. Then, the incremental feature selection method with random forest and support vector machine was used to discover the essential features. Furthermore, 38 key features were determined by integrating results of different feature selection methods, which may provide critical insights into the subcellular location of proteins. Their biological functions of subcellular localizations were discussed according to recent publications. In summary, our computational framework can help advance the understanding of subcellular localization prediction techniques and provide a new perspective to investigate the patterns of protein subcellular localization and their biological importance.

Expanding the terpene biosynthetic code with non-canonical 16 carbon atom building blocks

Humankind relies on specialized metabolites for medicines, flavors, fragrances, and numerous other valuable biomaterials. However, the chemical space occupied by specialized metabolites, and, thus, their application potential, is limited because their biosynthesis is based on only a handful of building blocks. Engineering organisms to synthesize alternative building blocks will bypass this limitation and enable the sustainable production of molecules with non-canonical chemical structures, expanding the possible applications. Herein, we focus on isoprenoids and combine synthetic biology with protein engineering to construct yeast cells that synthesize 10 non-canonical isoprenoid building blocks with 16 carbon atoms. We identify suitable terpene synthases to convert these building blocks into C₁₆ scaffolds and a cytochrome P450 to decorate the terpene scaffolds and produce different oxygenated compounds. Thus, we reconstruct the modular structure of terpene biosynthesis on 16-carbon backbones, synthesizing 28 different non-canonical terpenes, some of which have interesting odorant properties.

Alternative splicing of Arabidopsis G6PD5 recruits NADPH-producing OPPP reactions to the endoplasmic reticulum

Glucose-6-phosphate dehydrogenase is the rate-limiting enzyme of the oxidative pentose-phosphate pathway (OPPP). The OPPP mainly provides NADPH and sugar-phosphate building blocks for anabolic pathways and is present in all eukaryotes. In plant cells, the irreversible part of the OPPP is found in several compartments. Among the isoforms catalyzing the first OPPP step in Arabidopsis, G6PD1 to G6PD4 target plastids (with G6PD1 being also directed to peroxisomes), whereas G6PD5 and G6PD6 operate in the cytosol. We noticed that alternative splice forms G6PD5.4 and G6PD5.5 encode N-terminally extended proteoforms. Compared to G6PD5.1, RT-PCR signals differed and fluorescent reporter fusions expressed in Arabidopsis protoplasts accumulated in distinct intracellular sites. Co-expression with organelle-specific markers revealed that the G6PD5.4 and G6PD5.5 proteoforms label different subdomains of the endoplasmic reticulum (ER), and analysis of C-terminal roGFP fusions showed that their catalytic domains face the cytosol. In g6pd5-1 g6pd6-2 mutant protoplasts lacking cytosolic G6PDH activity, the ER-bound proteoforms were both active and thus able to form homomers. Among the Arabidopsis 6-phosphogluconolactonases (catalyzing the second OPPP step), we noticed that isoform PGL2 carries a C-terminal CaaX motif that may be prenylated for membrane attachment. Reporter-PGL2 fusions co-localized with G6PD5.4 in ER subdomains, which was abolished by Cys-to-Ser exchange in the ²⁵⁶CSIL motif. Among the Arabidopsis 6-phosphogluconate dehydrogenases (catalyzing the third OPPP step), S-acylated peptides were detected for all three isoforms in a recent palmitoylome, with dual cytosolic/peroxisomal PGD2 displaying three sites. Co-expression of GFP-PGD2 diminished crowding of OFP-G6PD5.4 at the ER, independent of PGL2's presence. Upon pull-down of GFP-G6PD5.4, not only unlabeled PGD2 and PGL2 were enriched, but also enzymes that depend on NADPH provision at the ER, indicative of physical interaction with the OPPP enzymes. When membrane-bound G6PD5.5 and 5.4 variants were co-expressed with KCR1 (ketoacyl-CoA reductase, involved in fatty acid elongation), ATR1 (NADPH:cytochrome-P450 oxidoreductase), or pulled C4H/CYP73A5 (cinnamate 4-hydroxylase) as indirectly (via ATR) NADPH-dependent cytochrome P450 enzyme, co-localization in ER subdomains was observed. Thus, alternative splicing of G6PD5 can direct the NADPH-producing OPPP reactions to the cytosolic face of the ER, where they may operate as membrane-bound metabolon to support several important biosynthetic pathways of plant cells.

Phylogeny-aware chemoinformatic analysis of chemical diversity in the Lamiaceae enables iridoid pathway assembly and discovery of aucubin synthase

Countless reports describe the isolation and structural characterization of natural products, yet this information remains disconnected and under-utilized. Using a cheminformatics approach, we leverage the reported observations of iridoid glucosides with the known phylogeny of a large iridoid producing plant family (Lamiaceae), to generate a set of biosynthetic pathways that best explain the extant iridoid chemical diversity. We developed a pathway reconstruction algorithm that connects iridoid reports via reactions, and prunes this solution space by considering phylogenetic relationships between genera. We formulate a model that emulates the evolution of iridoid glucosides to create a synthetic dataset, used to select the parameters that would best reconstruct the pathways, and apply them to the iridoid dataset to generate Pathway Hypotheses. These computationally generated pathways were then used as the basis by which to select and screen biosynthetic enzyme candidates. Our model was successfully applied to discover a cytochrome P450 enzyme from Callicarpa americana that catalyzes the oxidation of bartsioside to aucubin, predicted by our model despite neither molecule having been observed in the genus. We also demonstrate aucubin synthase activity in orthologues of Vitex agnus-castus, and the outgroup Paulownia tomentosa, further strengthening the hypothesis, enabled by our model, that the reaction was present in the ancestral biosynthetic pathway. This is the first systematic hypothesis on the epi-iridoid glucosides biosynthesis in 25 years, and sets the stage for streamlined work on the iridoid pathway. This work highlights how curation and computational analysis of widely available structural data can facilitate hypothesis-based gene discovery.

Functional and structural insight into the flexibility of cytochrome P450 reductases from Sorghum bicolor and its implications for lignin composition

Plant NADPH-dependent cytochrome P450 reductase (CPR) is a multidomain enzyme that donates electrons for hydroxylation reactions catalyzed by class II cytochrome P450 monooxygenases involved in the synthesis of many primary and secondary metabolites. These P450 enzymes include trans-cinnamate-4-hydroxylase, p-coumarate-3′-hydroxylase, and ferulate-5-hydroxylase involved in monolignol biosynthesis. Because of its role in monolignol biosynthesis, alterations in CPR activity could change the composition and overall output of lignin. Therefore, to understand the structure and function of three CPR subunits from sorghum, recombinant subunits SbCPR2a, SbCPR2b, and SbCPR2c were subjected to X-ray crystallography and kinetic assays. Steady-state kinetic analyses demonstrated that all three CPR subunits supported the oxidation reactions catalyzed by SbC4H1 (CYP73A33) and SbC3′H (CYP98A1). Furthermore, comparing the SbCPR2b structure with the well-investigated CPRs from mammals enabled us to identify critical residues of functional importance and suggested that the plant flavin mononucleotide–binding domain might be more flexible than mammalian homologs. In addition, the elucidated structure of SbCPR2b included the first observation of NADP⁺ in a native CPR. Overall, we conclude that the connecting domain of SbCPR2, especially its hinge region, could serve as a target to alter biomass composition in bioenergy and forage sorghums through protein engineering.

Isolation and Functional Characterization of Multiple NADPH-Cytochrome P450 Reductase Genes from Camellia sinensis in View of Catechin Biosynthesis

Catechins are critical constituents for the sensory quality and health-promoting benefits of tea. Cytochrome P450 monooxygenases are required for catechin biosynthesis and are dependent on NADPH-cytochrome P450 reductases (CPRs) to provide reducing equivalents for their activities. However, CPRs have not been identified in tea, and their relationship to catechin accumulation also remains unknown. Thus, three CsCPR genes were identified in this study, all of which had five CPR-related conserved domains and were targeted to the endoplasmic reticulum. These three recombinant CsCPR proteins could reduce cytochrome c using NADPH as an electron donor. Heterologous co-expression in yeast demonstrated that all the three CsCPRs could support the enzyme activities of CsC4H and CsF3'H. Correlation analysis indicated that the expression level of CsCPR1 (or CsCPR2 or CsCPR3) was positively correlated with 3',4',5'-catechin (or total catechins) content. Our results indicate that the CsCPRs are involved in the biosynthesis of catechins in tea leaves.

Molecular Cloning and Identification of NADPH Cytochrome P450 Reductase from Panax ginseng

Ginseng (Panax ginseng C.A. Mey.) is a precious Chinese traditional medicine, for which ginsenosides are the most important medicinal ingredients. Cytochrome P450 enzymes (CYP450) and their primary redox molecular companion NADPH cytochrome P450 reductase (CPR) play a key role in ginsenoside biosynthesis pathway. However, systematic studies of CPR genes in ginseng have not been reported. Numerous studies on ginsenoside synthesis biology still use Arabidopsis CPR (AtCPR1) as a reductase. In this study, we isolated two CPR genes (PgCPR1, PgCPR2) from ginseng adventitious roots. Phylogenetic tree analysis showed that both PgCPR1 and PgCPR2 are grouped in classⅡ of dicotyledonous CPR. Enzyme experiments showed that recombinant proteins PgCPR1, PgCPR2 and AtCPR1 can reduce cytochrome c and ferricyanide with NADPH as the electron donor, and PgCPR1 had the highest enzymatic activities. Quantitative real-time PCR analysis showed that PgCPR1 and PgCPR2 transcripts were detected in all examined tissues of Panax ginseng and both showed higher expression in stem and main root. Expression levels of the PgCPR1 and PgCPR2s were both induced after a methyl jasmonate (MeJA) treatment and its pattern matched with ginsenoside accumulation. The present investigation suggested PgCPR1 and PgCPR2 are associated with the biosynthesis of ginsenoside. This report will assist in future CPR family studies and ultimately improving ginsenoside production through transgenic engineering and synthetic biology.

Matrix metalloproteinases are involved in eclosion and wing expansion in the American cockroach, Periplaneta americana

Matrix metalloproteinases (MMPs) are the major proteinases that process or degrade numerous extracellular matrix (ECM) components and are evolutionarily conserved from nematodes to humans. During molting in insects, the old cuticle is removed and replaced by a new counterpart. Although the regulatory mechanisms of hormones and nutrients in molting have been well studied, very little is known about the roles of ECM-modifying enzymes in this process. Here, we found that MMPs are necessary for imaginal molting of the American cockroach, Periplaneta americana. Inhibition of Mmp activity via inhibitor treatment led to the failure of eclosion and wing expansion. Five Mmps genes were identified from the P. americana genome, and PaMmp2 played the dominant roles during molting. Further microscopic investigations showed that newly formed adult cuticles were attenuated and that then chitin content was reduced upon Mmp inhibition. Transcriptomic analysis of the integument demonstrated that multiple signaling and metabolic pathways were changed. Microscopic investigation of the wings showed that epithelial cells were restrained together because they were incapable of degrading the ECM upon Mmp inhibition. Transcriptomic analysis of the wing identified dozens of possible genes functioned in wing expansion. This is the first study to show the essential roles of Mmps in the nymph-adult transition of hemimetabolous insects.

Limonoid biosynthesis 3: Functional characterization of crucial genes involved in neem limonoid biosynthesis

Neem (Azadirachta indica L.) is well known for its medicinal, agricultural, and pesticidal applications since ages. The secondary metabolites, limonoids, confer these biological properties, wherein over 150 different limonoids have been reported from neem. To understand limonoid biosynthesis, we analyzed tissue-specific (kernel, pericarp, leaves, and flower) transcriptome that resulted in the identification of one farnesyl diphosphate synthase (AiFDS), one squalene synthase (AiSQS), three squalene epoxidases (AiSQE1, AiSQE2, and AiSQE3), two triterpene synthases (AiTTS1 and AiTTS2), cycloartenol synthase (AiCAS), two cytochrome P450 reductases, and ten cytochrome P450 systems. Comparative tissue-expression analysis indicated that AiFDS, AiSQS, AiSQE3, and AiTTS1 are expressed higher in the kernel than in the other tissues. Heterologously expressed recombinant AiTTS1 produced tirucalla-7,24-dien-3β-ol as the sole product. Expression profile data, phylogeny with triterpene synthases from Meliaceae and Rutaceae families, real-time PCR of different tissues, and transient transformation revealed the involvement of tirucalla-7,24-dien-3β-ol synthase (AiTTS1) in limonoid biosynthesis. Further, mutagenesis studies of AiTTS1 indicated that Y125 and F260 are probably involved in stabilization of dammarenyl cation. A 2.6-fold increase in production of tirucalla-7,24-dien-3β-ol was observed when AiSQE1 was co-expressed with mutant AiTTS1 in a yeast system. Furthermore, we functionally characterized the highly expressed cytochrome P450 reductases and cycloartenol synthase. This study helps in further analysis and identification of genes involved in limonoid biosynthesis in Meliaceae/Rutaceae and their production in a metabolically tractable heterologous system.

Two bi-functional cytochrome P450 CYP72 enzymes from olive (Olea europaea) catalyze the oxidative C-C bond cleavage in the biosynthesis of secoxy-iridoids - flavor and quality determinants in olive oil

Olive (Olea europaea) is an important crop in Europe, with high cultural, economic and nutritional significance. Olive oil flavor and quality depend on phenolic secoiridoids, but the biosynthetic pathway of these iridoids remains largely uncharacterized. We discovered two bifunctional cytochrome P450 enzymes, catalyzing the rare oxidative C-C bond cleavage of 7-epi-loganin to produce oleoside methyl ester (OeOMES) and secoxyloganin (OeSXS), both through a ketologanin intermediary. Although these enzymes are homologous to the previously reported Catharanthus roseus secologanin synthase (CrSLS), the substrate and product profiles differ. Biochemical assays provided mechanistic insights into the two-step OeOMES and CrSLS reactions. Model-guided mutations of OeOMES changed the product profile in a predictable manner, revealing insights into the molecular basis for this change in product specificity. Our results suggest that, in contrast to published hypotheses, in planta production of secoxy-iridoids is secologanin-independent. Notably, sequence data of cultivated and wild olives point to a relation between domestication and OeOMES expression. Thus, the discovery of this key biosynthetic gene suggests a link between domestication and secondary metabolism, and could potentially be used as a genetic marker to guide next-generation breeding programs.

Comparative Analysis of NADPH-Cytochrome P450 Reductases From Legumes for Heterologous Production of Triterpenoids in Transgenic Saccharomyces cerevisiae

Triterpenoids are plant specialized metabolites with various pharmacological activities. They are widely distributed in higher plants, such as legumes. Because of their low accumulation in plants, there is a need for improving triterpenoid production. Cytochrome P450 monooxygenases (CYPs) play critical roles in the structural diversification of triterpenoids. To perform site-specific oxidations, CYPs require the electrons that are transferred by NADPH-cytochrome P450 reductase (CPR). Plants possess two main CPR classes, class I and class II. CPR classes I and II have been reported to be responsible for primary and specialized (secondary) metabolism, respectively. In this study, we first analyzed the CPR expression level of three legumes species, Medicago truncatula, Lotus japonicus, and Glycyrrhiza uralensis, showing that the expression level of CPR class I was lower and more stable, while that of CPR class II was higher in almost all the samples. We then co-expressed different combinations of CYP716As and CYP72As with different CPR classes from these three legumes in transgenic yeast. We found that CYP716As worked better with CPR-I from the same species, while CYP72As worked better with any CPR-IIs. Using engineered yeast strains, CYP88D6 paired with class II GuCPR produced the highest level of 11-oxo-β-amyrin, the important precursor of high-value metabolites glycyrrhizin. This study provides insight into co-expressing genes from legumes for heterologous production of triterpenoids in yeast.

Functional characterization of cinnamate 4-hydroxylase from Helianthus annuus Linn using a fusion protein method

Sunflower (Helianthus annuus L.) is an important oil crop, the secondary metabolites of it include many compounds such as flavonoids and lignin. However, the research on the biosynthesis of phenolic compounds in sunflowers is still scarce. Cinnamate 4-hydroxylase (C4H) belongs to the cytochrome P450-dependent monooxygenase family and is involved in the synthesis of many phenolic compounds, but C4H in sunflowers has not yet been cloned and functionally characterized. In this study, we screened three C4H genes from the sunflower transcriptome and genomic databases, named HaC4H1, HaC4H2, and, HaC4H3, respectively. In heterologous expression experiments, we had improved a method from previous studies by the addition of restriction sites to make it easier to express multiple C4H functions and suitable for in vitro activity verification. HaC4Hs without the N-terminal membrane anchor region was fused with a redox partner of Arabidopsis thaliana cytochrome P450 enzyme (CYP450) by the method and functionally expressed in E. coli and the results showed that these three enzymes catalyzed the formation of p-coumaric acid. To further investigate whether our fusion protein approach is applicable to other C4Hs, we used this method to explore the functions of C4H from Peucedanum praeruptorum and Angelica decursiva, and they can also convert trans-cinnamic acid to p-coumaric acid. The gene expression profile showed that all three HaC4H genes showed the highest transcription levels in the roots and might be up-regulated by MeJA. In summary, these results reveal the function of HaC4Hs in sunflower and provide a simpler way to explore C4H and even other cytochrome P450 enzymes in prokaryotic expression systems.

Functional characterization of NADPH-cytochrome P450 reductase and cinnamic acid 4-hydroxylase encoding genes from Scoparia dulcis L

BACKGROUND

Most plant cytochrome P450 (P450) proteins need to be supplied with electrons from a redox partner, e.g. an NADPH-cytochrome P450 reductase (CPR), for the activation of oxygen molecules via heme. CPR is a flavoprotein with an N-terminal transmembrane domain, which transfers electrons from NADPH to the P450 via coenzymes flavin adenine dinucleotide and flavin mononucleotide.

RESULTS

In this study, a novel CPR (SdCPR) was isolated from a tropical medicinal plant Scoparia dulcis L. The deduced amino acid of SdCPR showed high homology of > 76% with CPR from higher plants and belonged to the class II CPRs of dicots. Recombinant SdCPR protein reduced cytochrome c, ferricyanide (K₃Fe(CN)₆), and dichlorophenolindophenol in an NADPH-dependent manner. To elucidate the P450 monooxygenase activity of SdCPR, we isolated a cinnamic acid 4-hydroxylase (SdC4H, CYP73A111) gene from S. dulcis. Biochemical characterization of SdCPR/SdC4H demonstrated that SdCPR supports the oxidation step of SdC4H. Real-time qPCR results showed that expression levels of SdCPR and SdC4H were inducible by mechanical wounding treatment and phytohormone elicitation (methyl jasmonate, salicylic acid), which were consistent with the results of promotor analyses.

CONCLUSIONS

Our results showed that the SdCPR and SdC4H are related to defense reactions, including the biosynthesis of secondary metabolites.

Orthogonal monoterpenoid biosynthesis in yeast constructed on an isomeric substrate

Synthetic biology efforts for the production of valuable chemicals are frequently hindered by the structure and regulation of the native metabolic pathways of the chassis. This is particularly evident in the case of monoterpenoid production in Saccharomyces cerevisiae, where the canonical terpene precursor geranyl diphosphate is tightly coupled to the biosynthesis of isoprenoid compounds essential for yeast viability. Here, we establish a synthetic orthogonal monoterpenoid pathway based on an alternative precursor, neryl diphosphate. We identify structural determinants of isomeric substrate selectivity in monoterpene synthases and engineer five different enzymes to accept the alternative substrate with improved efficiency and specificity. We combine the engineered enzymes with dynamic regulation of metabolic flux to harness the potential of the orthogonal substrate and improve the production of industrially-relevant monoterpenes by several-fold compared to the canonical pathway. This approach highlights the introduction of synthetic metabolism as an effective strategy for high-value compound production.

Functional characterization of squalene epoxidase and NADPH-cytochrome P450 reductase in Dioscorea zingiberensis

Dioscorea zingiberensis is a perennial medicinal herb rich in a variety of pharmaceutical steroidal saponins. Squalene epoxidase (SE) is the key enzyme in the biosynthesis pathways of triterpenoids and sterols, and catalyzes the epoxidation of squalene in coordination with NADPH-cytochrome P450 reductase (CPR). In this study, we cloned DzSE and DzCPR gene sequences from D. zingiberensis leaves, encoding proteins with 514 and 692 amino acids, respectively. Recombinant proteins were successfully expressed in vitro, and enzymatic analysis indicated that, when SE and CPR were incubated with the substrates squalene and NADPH, 2,3-oxidosqualene was formed as the product. Subcellular localization revealed that both the DzSE and DzCPR proteins are localized to the endoplasmic reticulum. The changes in transcription of DzSE and DzCPR were similar in several tissues. DzSE expression was enhanced in a time-dependent manner after methyl jasmonate (MeJA) treatments, while DzCPR expression was not inducible.

De Novo Biosynthesis of p-Coumaric Acid in E. coli with a trans-Cinnamic Acid 4-Hydroxylase from the Amaryllidaceae Plant Lycoris aurea

p-Coumaric acid is a commercially available phenolcarboxylic acid with a great number of important applications in the nutraceutical, pharmaceutical, material and chemical industries. p-Coumaric acid has been biosynthesized in some engineered microbes, but the potential of the plant CYP450-involved biosynthetic route has not investigated in Escherichia coli. In the present study, a novel trans-cinnamic acid 4-hydroxylase (C4H) encoding the LauC4H gene was isolated from Lycoris aurea (L’ Hér.) Herb via rapid amplification of cDNA ends. Then, N-terminal 28 amino acids of LauC4H were characterized, for the subcellular localization, at the endoplasmic reticulum membrane in protoplasts of Arabidopsis thaliana. In E. coli, LauC4H without the N-terminal membrane anchor region was functionally expressed when fused with the redox partner of A. thaliana cytochrome P450 enzyme (CYP450), and was verified to catalyze the trans-cinnamic acid to p-coumaric acid transformation by whole-cell bioconversion, HPLC detection and LC-MS analysis as well. Further, with phenylalanine ammonia-lyase 1 of A. thaliana, p-coumaric acid was de novo biosynthesized from glucose as the sole carbon source via the phenylalanine route in the recombinant E. coli cells. By regulating the level of intracellular NADPH, the production of p-coumaric acid was dramatically improved by 9.18-fold, and achieved with a titer of 156.09 μM in shake flasks. The recombinant cells harboring functional LauC4H afforded a promising chassis for biological production of p-coumaric acid, even other derivatives, via a plant CYP450-involved pathway.

The scaffold proteins of lignin biosynthetic cytochrome P450 enzymes

Lignin is a complex and irregular biopolymer of crosslinked phenylpropanoid units in plant secondary cell walls. Its biosynthesis requires three endoplasmic reticulum (ER)-resident cytochrome P450 monooxygenases, C4H, C3'H and F5H, to establish the structural characteristics of its monomeric precursors. These P450 enzymes were reported to associate with each other or potentially with other soluble monolignol biosynthetic enzymes to form an enzyme complex or a metabolon. However, the molecular basis governing such enzyme or pathway organization remains elusive. Here, we show that Arabidopsis membrane steroid-binding proteins (MSBPs) serve as a scaffold to physically organize monolignol P450 monooxygenases, thereby regulating the lignin biosynthetic process. We find that although C4H, C3'H and F5H are in spatial proximity to each other on the ER membrane in vivo, they do not appear to directly interact with each other. Instead, two MSBP proteins physically interact with all three P450 enzymes and, moreover, MSBPs themselves associate as homomers and heteromers on the ER membrane, thereby organizing P450 clusters. Downregulation of MSBP genes does not affect the transcription levels of monolignol biosynthetic P450 genes but substantially impairs the stability and activity of the MSBP-interacting P450 enzymes and, consequently, lignin deposition, and the accumulation of soluble phenolics in the monolignol branch but not in the flavonoid pathway. Our study suggests that MSBP proteins are essential structural components in the ER membrane that physically organize and stabilize the monolignol biosynthetic P450 enzyme complex, thereby specifically controlling phenylpropanoid-monolignol branch biosynthesis.

Biosynthesis of bioactive diterpenoids in the medicinal plant Vitex agnus-castus

Vitex agnus-castus L. (Lamiaceae) is a medicinal plant historically used throughout the Mediterranean region to treat menstrual cycle disorders, and is still used today as a clinically effective treatment for premenstrual syndrome. The pharmaceutical activity of the plant extract is linked to its ability to lower prolactin levels. This feature has been attributed to the presence of dopaminergic diterpenoids that can bind to dopamine receptors in the pituitary gland. Phytochemical analyses of V. agnus-castus show that it contains an enormous array of structurally related diterpenoids and, as such, holds potential as a rich source of new dopaminergic drugs. The present work investigated the localisation and biosynthesis of diterpenoids in V. agnus-castus. With the assistance of matrix-assisted laser desorption ionisation-mass spectrometry imaging (MALDI-MSI), diterpenoids were localised to trichomes on the surface of fruit and leaves. Analysis of a trichome-specific transcriptome database, coupled with expression studies, identified seven candidate genes involved in diterpenoid biosynthesis: three class II diterpene synthases (diTPSs); three class I diTPSs; and a cytochrome P450 (CYP). Combinatorial assays of the diTPSs resulted in the formation of a range of different diterpenes that can account for several of the backbones of bioactive diterpenoids observed in V. agnus-castus. The identified CYP, VacCYP76BK1, was found to catalyse 16-hydroxylation of the diol-diterpene, peregrinol, to labd-13Z-ene-9,15,16-triol when expressed in Saccharomyces cerevisiae. Notably, this product is a potential intermediate in the biosynthetic pathway towards bioactive furan- and lactone-containing diterpenoids that are present in this species.

Gene Duplication Leads to Altered Membrane Topology of a Cytochrome P450 Enzyme in Seed Plants

Evolution of the phenolic metabolism was critical for the transition of plants from water to land. A cytochrome P450, CYP73, with cinnamate 4-hydroxylase (C4H) activity, catalyzes the first plant-specific and rate-limiting step in this pathway. The CYP73 gene is absent from green algae, and first detected in bryophytes. A CYP73 duplication occurred in the ancestor of seed plants and was retained in Taxaceae and most angiosperms. In spite of a clear divergence in primary sequence, both paralogs can fulfill comparable cinnamate hydroxylase roles both in vitro and in vivo. One of them seems dedicated to the biosynthesis of lignin precursors. Its N-terminus forms a single membrane spanning helix and its properties and length are highly constrained. The second is characterized by an elongated and variable N-terminus, reminiscent of ancestral CYP73s. Using as proxies the Brachypodium distachyon proteins, we show that the elongation of the N-terminus does not result in an altered subcellular localization, but in a distinct membrane topology. Insertion in the membrane of endoplasmic reticulum via a double-spanning open hairpin structure allows reorientation to the lumen of the catalytic domain of the protein. In agreement with participation to a different functional unit and supramolecular organization, the protein displays modified heme proximal surface. These data suggest the evolution of divergent C4H enzymes feeding different branches of the phenolic network in seed plants. It shows that specialization required for retention of gene duplicates may result from altered protein topology rather than change in enzyme activity.

GFDP: the gene family database in poplar

A gene family is formed by duplication of a single original gene. Poplar trees (genus Populus) are important, principally because of their ecological and economic benefits, and are one of the most widely distributed and adaptable trees in the world. Systematic identification and annotation of gene family members are primary steps in studying the function and evolution of poplar genomes. Here, we describe the construction of the Gene Family Database in Poplar (GFDP), which contains information that systematically describes 6551 genes distributed in 145 gene families. GFDP is designed to present important biological information, such as gene structure, protein length, isoelectric point and functional and evolutionary information, using highly visual displays. Data and graphs are visualized by a web-based interface. Users can browse and download data through all the major browsers. GFDP provides a comprehensive platform with a solid foundation for further study of poplar gene families. GFDP is free available.

Heterologous production of a ginsenoside saponin (compound K) and its precursors in transgenic tobacco impairs the vegetative and reproductive growth

MAIN CONCLUSION

Production of compound K (a ginsenoside saponin) and its precursors in transgenic tobacco resulted in stunted growth and seed set failure, which may be caused by strong autotoxicity of heterologously produced phytochemicals against the tobacco itself. Panax ginseng roots contain various saponins (ginsenosides), which are major bioactive compounds. A monoglucosylated saponin, compound K (20-O-(β-D-glucopyranosyl)-20(S)-protopanaxadiol), has high medicinal and cosmetic values but is present in undetectable amounts in naturally grown ginseng roots. The production of compound K (CK) requires complicated deglycosylation of ginsenosides using physicochemical and/or enzymatic degradation. In this work, we report the production of CK in transgenic tobacco by co-overexpressing three genes (PgDDS, CYP716A47 and UGT71A28) isolated from P. ginseng. Introduction and expression of the transgenes in tobacco lines were confirmed by genomic PCR and RT-PCR. All the lines of transgenic tobacco produced CK including its precursors, protopanaxadiol and dammarenediol-II (DD). The concentrations of CK in the leaves ranged from 1.55 to 2.64 µg/g dry weight, depending on the transgenic line. Interestingly, production of CK in tobacco brought stunted plant growth and gave rise to seed set failure. This seed set failure was caused by both long-styled flowers and abnormal pollen development in transgenic tobacco. Both CK and DD treatments highly suppressed in vitro germination and tube growth in wild-type pollens. Based on these results, metabolic engineering for CK production in transgenic tobacco was successfully achieved, but the production of CK and its precursors in tobacco severely affects vegetative and reproductive growth due to the cytotoxicity of phytochemicals that are heterologously produced in transgenic tobacco.

Assembly of Dynamic P450-Mediated Metabolons-Order Versus Chaos

Purpose of Review

We provide an overview of the current knowledge on cytochrome P450-mediated metabolism organized as metabolons and factors that facilitate their stabilization. Essential parameters will be discussed including those that are commonly disregarded using the dhurrin metabolon from Sorghum bicolor as a case study.

Recent Findings

Sessile plants control their metabolism to prioritize their resources between growth and development, or defense. This requires fine-tuned complex dynamic regulation of the metabolic networks involved. Within the recent years, numerous studies point to the formation of dynamic metabolons playing a major role in controlling the metabolic fluxes within such networks.

Summary

We propose that P450s and their partners interact and associate dynamically with POR, which acts as a charging station possibly in concert with Cytb5. Solvent environment, lipid composition, and non-catalytic proteins guide metabolon formation and thereby activity, which have important implications for synthetic biology approaches aiming to produce high-value specialized metabolites in heterologous hosts.

Class II Cytochrome P450 Reductase Governs the Biosynthesis of Alkaloids

Expansion of the biosynthesis of plant specialized metabolites notably results from the massive recruitment of cytochrome P450s that catalyze multiple types of conversion of biosynthetic intermediates. For catalysis, P450s require a two-electron transfer catalyzed by shared cytochrome P450 oxidoreductases (CPRs), making these auxiliary proteins an essential component of specialized metabolism. CPR isoforms usually group into two distinct classes with different proposed roles, namely involvement in primary and basal specialized metabolisms for class I and inducible specialized metabolism for class II. By studying the role of CPRs in the biosynthesis of monoterpene indole alkaloids, we provide compelling evidence of an operational specialization of CPR isoforms in Catharanthus roseus (Madagascar periwinkle). Global analyses of gene expression correlation combined with transcript localization in specific leaf tissues and gene-silencing experiments of both classes of CPR all point to the strict requirement of class II CPRs for monoterpene indole alkaloid biosynthesis with a minimal or null role of class I. Direct assays of interaction and reduction of P450s in vitro, however, showed that both classes of CPR performed equally well. Such high specialization of class II CPRs in planta highlights the evolutionary strategy that ensures an efficient reduction of P450s in specialized metabolism.

Yeast factories for the production of aromatic compounds: from building blocks to plant secondary metabolites

The aromatic amino acid biosynthesis pathway is a source to a plethora of commercially relevant chemicals with very diverse industrial applications. Tremendous efforts in microbial engineering have led to the production of compounds ranging from small aromatic molecular building blocks all the way to intricate plant secondary metabolites. Particularly, the yeast Saccharomyces cerevisiae has been a great model organism given its superior capability to heterologously express long metabolic pathways, especially the ones containing cytochrome P450 enzymes. This review contains a collection of state-of-the-art metabolic engineering work devoted towards unraveling the mechanisms for enhancing the flux of carbon into the aromatic pathway. Some of the molecules discussed include the polymer precursor muconic acid, as well as important nutraceuticals (flavonoids and stilbenoids), and opium-derived drugs (benzylisoquinoline alkaloids).

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.

Carnosic acid biosynthesis elucidated by a synthetic biology platform

Synthetic biology approaches achieving the reconstruction of specific plant natural product biosynthetic pathways in dedicated microbial "chassis" have provided access to important industrial compounds (e.g., artemisinin, resveratrol, vanillin). However, the potential of such production systems to facilitate elucidation of plant biosynthetic pathways has been underexplored. Here we report on the application of a modular terpene production platform in the characterization of the biosynthetic pathway leading to the potent antioxidant carnosic acid and related diterpenes in Salvia pomifera and Rosmarinus officinalis.Four cytochrome P450 enzymes are identified (CYP76AH24, CYP71BE52, CYP76AK6, and CYP76AK8), the combined activities of which account for all of the oxidation events leading to the biosynthesis of the major diterpenes produced in these plants. This approach develops yeast as an efficient tool to harness the biotechnological potential of the numerous sequencing datasets that are increasingly becoming available through transcriptomic or genomic studies.

Production of the forskolin precursor 11β-hydroxy-manoyl oxide in yeast using surrogate enzymatic activities

Background

Several plant diterpenes have important biological properties. Among them, forskolin is a complex labdane-type diterpene whose biological activity stems from its ability to activate adenylyl cyclase and to elevate intracellular cAMP levels. As such, it is used in the control of blood pressure, in the protection from congestive heart failure, and in weight-loss supplements. Chemical synthesis of forskolin is challenging, and production of forskolin in engineered microbes could provide a sustainable source. To this end, we set out to establish a platform for the production of forskolin and related epoxy-labdanes in yeast.

Results

Since the forskolin biosynthetic pathway has only been partially elucidated, and enzymes involved in terpene biosynthesis frequently exhibit relaxed substrate specificity, we explored the possibility of reconstructing missing steps of this pathway employing surrogate enzymes. Using CYP76AH24, a Salvia pomifera cytochrome P450 responsible for the oxidation of C-12 and C-11 of the abietane skeleton en route to carnosic acid, we were able to produce the forskolin precursor 11β-hydroxy-manoyl oxide in yeast. To improve 11β-hydroxy-manoyl oxide production, we undertook a chassis engineering effort involving the combination of three heterozygous yeast gene deletions (mct1/MCT1, whi2/WHI2, gdh1/GDH1) and obtained a 9.5-fold increase in 11β-hydroxy-manoyl oxide titers, reaching 21.2 mg L⁻¹.

Conclusions

In this study, we identify a surrogate enzyme for the specific and efficient hydroxylation of manoyl oxide at position C-11β and establish a platform that will facilitate the synthesis of a broad range of tricyclic (8,13)-epoxy-labdanes in yeast. This platform forms a basis for the heterologous production of forskolin and will facilitate the elucidation of subsequent steps of forskolin biosynthesis. In addition, this study highlights the usefulness of using surrogate enzymes for the production of intermediates of complex biosynthetic pathways. The combination of heterozygous deletions and the improved yeast strain reported here will provide a useful tool for the production of numerous other isoprenoids.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-016-0440-8) contains supplementary material, which is available to authorized users.

Evolution of NADPH-cytochrome P450 oxidoreductases (POR) in Apiales - POR 1 is missing

The NADPH-dependent cytochrome P450 oxidoreductase (POR) is the obligate electron donor to eukaryotic microsomal cytochromes P450 enzymes. The number of PORs within plant species is limited to one to four isoforms, with the most common being two PORs per plant. These enzymes provide electrons to a huge number of different cytochromes P450s (from 50 to several hundred within one plant). Within the eudicotyledons, PORs can be divided into two major clades, POR 1 and POR 2. Based on our own sequencing analysis and publicly available data, we have identified 45 PORs from the angiosperm order Apiales. These were subjected to a phylogenetic analysis along with 237 other publicly available (NCBI and oneKP) POR sequences found within the clade Asterids. Here, we show that the order Apiales only harbor members of the POR 2 clade, which are further divided into two distinct subclades. This is in contrast to most other eudicotyledon orders that have both POR 1 and POR 2. This suggests that through gene duplications and one gene deletion, Apiales only contain members of the POR 2 clade. Three POR 2 isoforms from Thapsia garganica L., Apiaceae, were all full-length in an Illumina root transcriptome dataset (available from the SRA at NCBI). All three genes were shown to be functional upon reconstitution into nanodiscs, confirming that none of the isoforms are pseudogenes.

A Decade of Molecular Understanding of Withanolide Biosynthesis and In vitro Studies in Withania somnifera (L.) Dunal: Prospects and Perspectives for Pathway Engineering

Withania somnifera, a multipurpose medicinal plant is a rich reservoir of pharmaceutically active triterpenoids that are steroidal lactones known as withanolides. Though the plant has been well-characterized in terms of phytochemical profiles as well as pharmaceutical activities, limited attempts have been made to decipher the biosynthetic route and identification of key regulatory genes involved in withanolide biosynthesis. This scenario limits biotechnological interventions for enhanced production of bioactive compounds. Nevertheless, recent emergent trends vis-à-vis, the exploration of genomic, transcriptomic, proteomic, metabolomics, and in vitro studies have opened new vistas regarding pathway engineering of withanolide production. During recent years, various strategic pathway genes have been characterized with significant amount of regulatory studies which allude toward development of molecular circuitries for production of key intermediates or end products in heterologous hosts. Another pivotal aspect covering redirection of metabolic flux for channelizing the precursor pool toward enhanced withanolide production has also been attained by deciphering decisive branch point(s) as robust targets for pathway modulation. With these perspectives, the current review provides a detailed overview of various studies undertaken by the authors and collated literature related to molecular and in vitro approaches employed in W. somnifera for understanding various molecular network interactions in entirety.

Accumulation of cynaropicrin in globe artichoke and localization of enzymes involved in its biosynthesis

Globe artichoke (Cynara cardunculus var. scolymus) belongs to the Asteraceae family, in which one of the most biologically significant class of secondary metabolites are sesquiterpene lactones (STLs). In globe artichoke the principal STL is the cynaropicrin, which contributes to approximately 80% of its characteristic bitter taste. Cynaropicrin content was assessed in globe artichoke tissues and was observed to accumulate in leaves of different developmental stages. In the receptacle, a progressive decrease was observed during inflorescence development, while the STL could not be detected in the inflorescence bracts. Almost undetectable amounts were found in the roots and inflorescence stems at the commercial stage. Cynaropicrin content was found to correlate with expression of genes encoding CcGAS, CcGAO and CcCOS, which are involved in the STL biosynthesis. A more detailed study of leaf material revealed that cynaropicrin predominantly accumulates in the trichomes, and not in the apoplastic cavity fluids. Analysis of the promoter regions of CcGAO and CcCOS revealed the presence of L1-box motifs, which confers trichome-specific expression in Arabidopsis, suggesting that cynaropicrin is not only stored but also synthesized in trichomes. A transient expression of GFP fusion proteins was performed in Nicotiana benthamiana plants: the CcGAS fluorescence signal was located in the cytoplasm while the CcGAO and CcCOS localized to the endoplasmatic reticulum.

Production of the dammarene sapogenin (protopanaxadiol) in transgenic tobacco plants and cultured cells by heterologous expression of PgDDS and CYP716A47

KEY MESSAGE

Protopanaxadiol (PPD) is an aglycone of dammarene-type ginsenoside and has high medicinal values. In this work, we reported the PPD production in transgenic tobacco co-overexpressing PgDDS and CYP716A47. PPD is an aglycone of ginsenosides produced by Panax species and has a wide range of pharmacological activities. PPD is synthesized via the hydroxylation of dammarenediol-II (DD) by CYP716A47 enzyme. Here, we established a PPD production system via cell suspension culture of transgenic tobacco co-overexpressing the genes for PgDDS and CYP716A47. The concentration of PPD in transgenic tobacco leaves was 2.3-5.7 µg/g dry weight (DW), depending on the transgenic line. Leaf segments were cultured on medium with various types of hormones to induce callus. Auxin treatment, particularly 2,4-D, strongly enhanced the production of DD (783.8 µg g(-1) DW) and PPD (125.9 µg g(-1) DW). Treatment with 2,4-D enhanced the transcription of the HMG-Co reductase (HMGR) and squalene epoxidase genes. PPD production reached 166.9 and 980.9 µg g(-1) DW in a 250-ml shake flask culture and in 5-l airlift bioreactor culture, respectively.

A collection of cytochrome P450 monooxygenase genes involved in modification and detoxification of herbicide atrazine in rice (Oryza sativa) plants

Plant cytochrome P450 monooxygenases constitute one of the largest families of protein genes involved in plant growth, development and acclimation to biotic and abiotic stresses. However, whether these genes respond to organic toxic compounds and their biological functions for detoxifying toxic compounds such as herbicides in rice are poorly understood. The present study identified 201 genes encoding cytochrome P450s from an atrazine-exposed rice transcriptome through high-throughput sequencing. Of these, 69 cytochrome P450 genes were validated by microarray and some of them were confirmed by real time PCR. Activities of NADPH-cytochrome P450 reductase (CPR) and p-nitroanisole O-demethylase (PNOD) related to toxicity were determined and significantly induced by atrazine exposure. To dissect the mechanism underlying atrazine modification and detoxification by P450, metabolites (or derivatives) of atrazine in plants were analyzed by ultra performance liquid chromatography mass spectrometry (UPLC/MS). Major metabolites comprised desmethylatrazine (DMA), desethylatrazine (DEA), desisopropylatrazine (DIA), hydroxyatrazine (HA), hydroxyethylatrazine (HEA) and hydroxyisopropylatrazine (HIA). All of them were chemically modified by P450s. Furthermore, two specific inhibitors of piperonyl butoxide (PBO) and malathion (MAL) were used to assess the correlation between the P450s activity and rice responses including accumulation of atrazine in tissues, shoot and root growth and detoxification.

Molecular Cloning, Heterologous Expression, and Functional Characterization of an NADPH-Cytochrome P450 Reductase Gene from Camptotheca acuminata, a Camptothecin-Producing Plant

Camptothecin (CAM), a complex pentacyclic pyrroloqinoline alkaloid, is the starting material for CAM-type drugs that are well-known antitumor plant drugs. Although many chemical and biological research efforts have been performed to produce CAM, a few attempts have been made to uncover the enzymatic mechanism involved in the biosynthesis of CAM. Enzyme-catalyzed oxidoreduction reactions are ubiquitously presented in living organisms, especially in the biosynthetic pathway of most secondary metabolites such as CAM. Due to a lack of its reduction partner, most catalytic oxidation steps involved in the biosynthesis of CAM have not been established. In the present study, an NADPH-cytochrome P450 reductase (CPR) encoding gene CamCPR was cloned from Camptotheca acuminata, a CAM-producing plant. The full length of CamCPR cDNA contained an open reading frame of 2127-bp nucleotides, corresponding to 708-amino acid residues. CamCPR showed 70 ~ 85% identities to other characterized plant CPRs and it was categorized to the group II of CPRs on the basis of the results of multiple sequence alignment of the N-terminal hydrophobic regions. The intact and truncate CamCPRs with N- or C-terminal His₆-tag were heterologously overexpressed in Escherichia coli. The recombinant enzymes showed NADPH-dependent reductase activity toward a chemical substrate ferricyanide and a protein substrate cytochrome c. The N-terminal His₆-tagged CamCPR showed 18- ~ 30-fold reduction activity higher than the C-terminal His₆-tagged CamCPR, which supported a reported conclusion, i.e., the last C-terminal tryptophan of CPRs plays an important role in the discrimination between NADPH and NADH. Co-expression of CamCPR and a P450 monooxygenase, CYP73A25, a cinnamate 4-hydroxylase from cotton, and the following catalytic formation of p-coumaric acid suggested that CamCPR transforms electrons from NADPH to the heme center of P450 to support its oxidation reaction. Quantitative real-time PCR analysis showed that CamCPR was expressed in the roots, stems, and leaves of C. acuminata seedlings. The relative transcript level of CamCPR in leaves was 2.2-fold higher than that of roots and the stems showed 1.5-fold transcript level higher than the roots. The functional characterization of CamCPR will be helpful to disclose the mysterious mechanisms of the biosynthesis of CAM. The present study established a platform to characterize the P450 enzymes involved in the growth, development, and metabolism of eukaryotic organisms.

Biphenyl 4-Hydroxylases Involved in Aucuparin Biosynthesis in Rowan and Apple Are Cytochrome P450 736A Proteins

Upon pathogen attack, fruit trees such as apple (Malus spp.) and pear (Pyrus spp.) accumulate biphenyl and dibenzofuran phytoalexins, with aucuparin as a major biphenyl compound. 4-Hydroxylation of the biphenyl scaffold, formed by biphenyl synthase (BIS), is catalyzed by a cytochrome P450 (CYP). The biphenyl 4-hydroxylase (B4H) coding sequence of rowan (Sorbus aucuparia) was isolated and functionally expressed in yeast (Saccharomyces cerevisiae). SaB4H was named CYP736A107. No catalytic function of CYP736 was known previously. SaB4H exhibited absolute specificity for 3-hydroxy-5-methoxybiphenyl. In rowan cell cultures treated with elicitor from the scab fungus, transient increases in the SaB4H, SaBIS, and phenylalanine ammonia lyase transcript levels preceded phytoalexin accumulation. Transient expression of a carboxyl-terminal reporter gene construct directed SaB4H to the endoplasmic reticulum. A construct lacking the amino-terminal leader and transmembrane domain caused cytoplasmic localization. Functional B4H coding sequences were also isolated from two apple (Malus × domestica) cultivars. The MdB4Hs were named CYP736A163. When stems of cv Golden Delicious were infected with the fire blight bacterium, highest MdB4H transcript levels were observed in the transition zone. In a phylogenetic tree, the three B4Hs were closest to coniferaldehyde 5-hydroxylases involved in lignin biosynthesis, suggesting a common ancestor. Coniferaldehyde and related compounds were not converted by SaB4H.

Towards Elucidating Carnosic Acid Biosynthesis in Lamiaceae: Functional Characterization of the Three First Steps of the Pathway in Salvia fruticosa and Rosmarinus officinalis

Carnosic acid (CA) is a phenolic diterpene with anti-tumour, anti-diabetic, antibacterial and neuroprotective properties that is produced by a number of species from several genera of the Lamiaceae family, including Salvia fruticosa (Cretan sage) and Rosmarinus officinalis (Rosemary). To elucidate CA biosynthesis, glandular trichome transcriptome data of S. fruticosa were mined for terpene synthase genes. Two putative diterpene synthase genes, namely SfCPS and SfKSL, showing similarities to copalyl diphosphate synthase and kaurene synthase-like genes, respectively, were isolated and functionally characterized. Recombinant expression in Escherichia coli followed by in vitro enzyme activity assays confirmed that SfCPS is a copalyl diphosphate synthase. Coupling of SfCPS with SfKSL, both in vitro and in yeast, resulted in the synthesis miltiradiene, as confirmed by 1D and 2D NMR analyses (1H, 13C, DEPT, COSY H-H, HMQC and HMBC). Coupled transient in vivo assays of SfCPS and SfKSL in Nicotiana benthamiana further confirmed production of miltiradiene in planta. To elucidate the subsequent biosynthetic step, RNA-Seq data of S. fruticosa and R. officinalis were searched for cytochrome P450 (CYP) encoding genes potentially involved in the synthesis of the first phenolic compound in the CA pathway, ferruginol. Three candidate genes were selected, SfFS, RoFS1 and RoFS2. Using yeast and N. benthamiana expression systems, all three where confirmed to be coding for ferruginol synthases, thus revealing the enzymatic activities responsible for the first three steps leading to CA in two Lamiaceae genera.

Small glycosylated lignin oligomers are stored in Arabidopsis leaf vacuoles

Lignin is an aromatic polymer derived from the combinatorial coupling of monolignol radicals in the cell wall. Recently, various glycosylated lignin oligomers have been revealed in Arabidopsis thaliana. Given that monolignol oxidation and monolignol radical coupling are known to occur in the apoplast, and glycosylation in the cytoplasm, it raises questions about the subcellular localization of glycosylated lignin oligomer biosynthesis and their storage. By metabolite profiling of Arabidopsis leaf vacuoles, we show that the leaf vacuole stores a large number of these small glycosylated lignin oligomers. Their structural variety and the incorporation of alternative monomers, as observed in Arabidopsis mutants with altered monolignol biosynthesis, indicate that they are all formed by combinatorial radical coupling. In contrast to the common believe that combinatorial coupling is restricted to the apoplast, we hypothesized that the aglycones of these compounds are made within the cell. To investigate this, leaf protoplast cultures were cofed with 13C6-labeled coniferyl alcohol and a 13C4-labeled dimer of coniferyl alcohol. Metabolite profiling of the cofed protoplasts provided strong support for the occurrence of intracellular monolignol coupling. We therefore propose a metabolic pathway involving intracellular combinatorial coupling of monolignol radicals, followed by oligomer glycosylation and vacuolar import, which shares characteristics with both lignin and lignan biosynthesis.

RNA interference of NADPH-cytochrome P450 reductase of the rice brown planthopper, Nilaparvata lugens, increases susceptibility to insecticides

BACKGROUND

NADPH-cytochrome P450 reductase (CPR) is essential for numerous biological reactions catalysed by microsomal cytochrome P450 monooxygenases (P450s). Knockdown of CPR in several insects leads to developmental defects and increased susceptibility to insecticides. However, information about the role of CPR in the brown planthopper, Nilaparvata lugens, is still unavailable.

RESULTS

A full-length cDNA encoding CPR was cloned from N. lugens (NlCPR). The deduced amino acid sequence showed marked features of classical CPRs, such as an N-terminus membrane anchor, conserved domains for flavin mononucleotide, flavin adenine dinucleotide (FAD) and nicotinamide adenine dinucleotide phosphate binding, as well as an FAD-binding motif and catalytic residues. Phylogenetic analysis revealed that NlCPR was located in a branch along with bed bug and pea aphid hemipteran insects. NlCPR mRNA was detectable in all tissues and developmental stages of N. lugens, as determined by real-time quantitative PCR. NlCPR transcripts were most abundant in the abdomen in adults, and in first-instar nymphs. Injection of N. lugens with double-strand RNA (dsRNA) against NlCPR significantly reduced the transcription level of the mRNA, and silencing of NlCPR resulted in increased susceptibility in N. lugens to beta-cypermethrin and imidacloprid.

CONCLUSION

The results provide first evidence that NlCPR contributes to the susceptibility to beta-cypermethrin and imidacloprid in N. lugens.

Molecular cloning, bacterial expression and functional characterisation of cytochrome P450 monooxygenase, CYP97C27, and NADPH-cytochrome P450 reductase, CPR I, from Croton stellatopilosus Ohba

The cDNAs for cytochrome P450 monooxygenase (designated as CYP97C27 by D. Nelson's group) and NADPH-cytochrome P450 reductase (designated as CPR I based on its classification) were isolated from Croton stellatopilosus leaves, which actively biosynthesise plaunotol (18-OH geranylgeraniol). CYP97C27 and CPR I contain open reading frames encoding proteins of 471 and 711 amino acids with predicted molecular masses of 53 and 79kDa, respectively. By aligning the deduced sequences of CYP97C27 and CPR I with other plant species, all functional domains of CYP97C27 (heme and oxygen binding) and CPR I (CYP- and FMN, FAD, and NADPH cofactor binding) were identified. Amino acid sequence comparison indicated that both CYP97C27 (85-93%) and CPR I (79-83%) share high sequence identities with homologous proteins in other plant species, suggesting that CYP97C27 belongs to the CYP97C subfamily and that CPR I belongs to class I of the dicotyledonous CPR. Functional characterisation of both enzymes, produced in Escherichia coli (pET32a/BL21(DE3)) as recombinant proteins, showed that simultaneous incubation of CYP97C27 and CPR I with the substrate geranylgeraniol (GGOH) and coenzyme NADPH led to formation of the product plaunotol. In C. stellatopilosus, the levels of the CYP97C27 and CPR I transcripts were highly correlated with those of several mRNAs involved in the plaunotol biosynthetic pathway, suggesting that CYP97C27 and CPR I are the enzymes that catalyse the last hydroxylation step of the pathway.

Mutation of the inducible ARABIDOPSIS THALIANA CYTOCHROME P450 REDUCTASE2 alters lignin composition and improves saccharification

ARABIDOPSIS THALIANA CYTOCHROME P450 REDUCTASE1 (ATR1) and ATR2 provide electrons from NADPH to a large number of CYTOCHROME P450 (CYP450) enzymes in Arabidopsis (Arabidopsis thaliana). Whereas ATR1 is constitutively expressed, the expression of ATR2 appears to be induced during lignin biosynthesis and upon stresses. Therefore, ATR2 was hypothesized to be preferentially involved in providing electrons to the three CYP450s involved in lignin biosynthesis: CINNAMATE 4-HYDROXYLASE (C4H), p-COUMARATE 3-HYDROXYLASE1 (C3H1), and FERULATE 5-HYDROXYLASE1 (F5H1). Here, we show that the atr2 mutation resulted in a 6% reduction in total lignin amount in the main inflorescence stem and a compositional shift of the remaining lignin to a 10-fold higher fraction of p-hydroxyphenyl units at the expense of syringyl units. Phenolic profiling revealed shifts in lignin-related phenolic metabolites, in particular with the substrates of C4H, C3H1 and F5H1 accumulating in atr2 mutants. Glucosinolate and flavonol glycoside biosynthesis, both of which also rely on CYP450 activities, appeared less affected. The cellulose in the atr2 inflorescence stems was more susceptible to enzymatic hydrolysis after alkaline pretreatment, making ATR2 a potential target for engineering plant cell walls for biofuel production.

Heterologous expression and functional characterization of the NADPH-cytochrome P450 reductase from Capsicum annuum

Two NADPH-cytochrome P450 reductase (CPR) genes (CaCPR1 and CaCPR2) were isolated from hot pepper (Capsicum annuum L. cv. Bukang). At the red ripe stage, the expression level of CaCPR1 was more than 6-fold greater than that in leaves or flowers. It gradually increased during fruit ripening. The CaCPR2 gene seemed to be expressed constitutively in all of the tested tissues. To investigate the enzymatic properties of CaCPR1, the cDNA of CaCPR1 was heterologously expressed in Escherichia coli without any modification of amino acid sequences, and CaCPR1 was purified. The enzymatic properties of CaCPR1 were confirmed using cytochrome c and cytochrome b5 as protein substrates. The CaCPR1 could support human CYP1A2-catalyzed reaction. It also reduced tetrazolium salts and ferricyanide. These results show that CaCPR1 is the major CPR in most hot pepper tissues. It is suggested that the CaCPR1 can be used a prototype for studying biological functions and biotechnological applications of plant CPRs.