Next Article in Journal
Evaluation of the Potential Release Tendency of Metals and Metalloids from the Estuarine Sediments: Case Study of Raša Bay
Previous Article in Journal
Protein Hydrolysis by Subcritical Water: A New Perspective on Obtaining Bioactive Peptides
Previous Article in Special Issue
Monoterpene Indole Alkaloids with Cav3.1 T-Type Calcium Channel Inhibitory Activity from Catharanthus roseus
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 26
Issue 21
10.3390/molecules26216654
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Xian Zou
,
Yue Zhang
,
Xu Zeng
,
Tuo Liu
,
Gui Li
,
Yuxin Dai
,
Yuanzhu Xie
and
Zhiyong Luo
*
Department of Biochemistry and Molecular Biology, School of Life Sciences, Central South University, Changsha 410008, China
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(21), 6654; https://doi.org/10.3390/molecules26216654
Submission received: 31 August 2021 / Revised: 29 September 2021 / Accepted: 8 October 2021 / Published: 3 November 2021
(This article belongs to the Special Issue A Themed Issue in Memory of Academician Xing Qiyi (1911–2002))
Browse Figures
Review Reports Versions Notes

Abstract

:
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.

1. Introduction

Panax ginseng C. A. Meyer is a traditional Chinese herbal medicine [1,2]. Ginsenoside is the most important medicinal active component in ginseng [3]. More and more research has found that ginsenoside can reduce oxidative stress response, protect the heart [4], slow down depression [5], prevent and treat Alzheimer’s disease [6,7,8], inhibit the growth of tumor cells [9,10,11,12] and other pharmacological effects. However, ginseng, a crop with a long growth period and a continuous cropping obstacle are seriously limited by the growing environment in yield and quality. Meanwhile, the contents of ginsenosides with unique pharmacological activities are very few [13]. So, the output of ginsenoside is far from meeting the market demand. Therefore, many researchers attempt to use functional gene identification and synthetic biology to improve the yield of ginsenoside.
In recent years, research into the ginsenoside biosynthesis pathway and its key enzymes have made continuous progress. Its biosynthesis pathway can be generally divided into three stages: Step one is to synthesize isoprene pyrophosphate (IPP) through the MVA pathway (Mevalonate pathway) and the MEP pathway (Methylerythritol 4-phosphate pathway); step two is to synthesize 2,3-oxidosqualene, which is the precursor of triterpenoid saponin; step three is the post-modification of compounds involved in the cyclization, hydroxylation, and glycosidation of 2,3-oxidosqualene [14,15,16,17] (Supplementary Material Figure S1).
Many enzymes, such as SQE, DDS and CAS, involved in biosynthesis of ginsenosides are shown to be CYP450, which determine the diversity of ginsenoside monomers [18,19,20]. NADPH-cytochrome P450 reductases (CPR) is the redox molecular companion of the cytochrome P450 enzyme system, belonging to the flavin protein family. It is responsible for transferring the electrons in reductive NADPH via flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) cofactors to the ferrous erythrocyte group of cytochrome P450 enzyme, providing electrons for a series of oxidation reactions involving cytochrome P450 enzyme [21]. Plant cytochrome P450 reductase-cytochrome P450 enzyme electron transfer system plays an important role in plant growth, development, biosynthesis of natural products and resistance to external stimuli [22].
In the plant genome, CYP450 is a superfamily of heme mercaptan protein, but CPRs are only encoded by a small number of genes. The amino acid sequences of CPR in different higher plants are highly homologous, so CPRs from different plants can be used as a molecular chaperone of P450 in yeast or other plants. However, different CPRs have different enzyme activities for the same substrate or P450 monooxygenase [23,24]. Sangkyu et al. showed that OsCPR2 showed the highest catalytic activity among the three rice CPRs [25]. In recent years, there are more and more CPRs of non-model plants. Two kinds of CPR were isolated and examined in cotton. The results showed that mechanical injury and inducer could induce the expression of GhCPR2, not GhCPR1 [26]. The sequence and function of ginseng cytochrome P450 have been analyzed one after another, but up to now, there has been no systematic study on ginseng CPR. Arabidopsis CPR (AtCPR1) is still used as a reductase in the study of CYP450 function and ginsenoside synthesis biology.
In this study, two CPR gene sequences PgCPR1(GenBank: MT160178) and PgCPR2(GenBank: MT160179) were selected from the MeJA-induced ginseng root transcription database and bioinformatics analysis was performed. The catalytic activities of PgCPR1, PgCPR2 and AtCPR1 were compared through an enzyme activity test. The expression patterns of PgCPR1 and PgCPR2 under different MeJA treatments were detected by qPCR. The systematic analysis of ginseng CPR is worthwhile in order to deepen our understanding of the CPR family and provide a reference for the study of ginseng secondary metabolites, and a more suitable CPR gene for the study of ginsenoside biosynthesis. This is of great significance for the construction of a suitable microbial expression system and the large-scale industrial production of ginsenosides.

2. Results

2.1. Molecular Cloning and Bioinformatics Analysis of PgCPR1 and PgCPR2 in Panax ginseng

In order to study the biosynthesis of ginsenosides, we performed transcriptome sequencing on Panax ginseng hairy roots and obtained the transcriptome information of ginseng. Two candidate CPR genes were screened from the transcriptome database and were named PgCPR1 and PgCPR2. According to the information of the gene sequence, we designed specific primers to clone PgCPRs (Supplementary Material Table S1). The relative molecular weight of PgCPR1 is 78.56 kDa, and that of PgCPR2 is 75.47 kDa. Blast results show that PgCPR1 and PgCPR2 are 99% similar to ginseng CPR1 (AHA50098.1) and ginseng CPR2 (AIC73829.1), which have been uploaded into the NCBI database. The NCBI conserved domains finder results show that both CPRs contain all functional domains that bind to CYP450 and interact with coenzyme NADPH, FMN, and FAD (Supplementary Material Figure S2). TMHMM prediction results show that both PgCPR1 and PgCPR2 contain transmembrane anchored domains (Supplementary Material Figure S3).
Comparing the amino acid sequence of ginseng PgCPRs with the known CPRs amino acid sequence of model plants, PgCPRs have the same functional characteristics as CPRs in other plants (Figure 1). According to their N-terminal sequence, plant CPRs are classified into class I and class II [27]. Class I includes members of Arabidopsis AtCPR1 and Andrographis APCPR1, and Class II includes AtCPR2, APCPR2 and APCPR3. Here, multiple alignments show that, in the red box, class I is under a longer sequence, while class II has a shorter sequence (Figure 1). MEGA 7.0 program and the neighbor-joining method were used to analyze the phylogeny of the CPR amino acid sequence of ginseng and other plants. The results showed that PgCPR1 and PgCPR2 were both grouped in class II (Figure 2A). The CPR of class II can be induced by stress or induction and is more likely to participate in secondary metabolism.
The three dimensional structure of the SWISS-MODEL showed (Figure 2B,C) that the similarity between PgCPR1 and ATR2 was 76.34%, and GMQE was 0.74 The similarity between PgCPR1 and ATR2 was 76.34%, and GMQE was 0.75 (Figure S4). Through bioinformatics analysis, it is found that the two CPRs of ginseng are very close to those of Arabidopsis CPR2. It is speculated that they may have the same redox function.

2.2. Heterologous Expression and Catalytic Activity Analysisof Recombinant PgCPRs and AtCPR1

The Coomassie blue staining results showed that the truncated PgCPR1, PgCPR2, and AtCPR1 were soluble after IPTG induction at low temperatures, and the purified CPR protein was obtained and used for subsequent enzyme activity experiments (Figure 3). Purified recombinant CPR protein was used for enzyme activity determination. The reaction of CPR protein was expected to be completed in a short time (3–5 min) [28,29]. In order to determine the demand of PgCPRs for electronic donors, NADH or NADPH dependent cytochrome c reduction activities of the recombinant proteins were determined respectively. For PgCPR1, PgCPR2 and AtCPR1, NADPH is effectively used as an electron donor to reduce cytochrome c. On the contrary, NADH cannot be recognized as an electron donor by them, which indicates that three CPRs especially use NADPH as an electron donor (Table 1). The kinetic properties of PgCPR1, PgCPR2 and AtCPR1 on cytochrome c, NADPH and ferricyanide were determined. Three CPRs activities of three different substrates all follow the typical Michaelis Menten curve (Figure 4). The catalytic efficiency of PgCPR1 to cytochrome c is almost the same as that of PgCPR2, which is higher than that of AtCPR1. For NADPH, relative catalytic activity (Kcat/Km value) of PgCPR1 is the highest, and that of AtCPR1 is the lowest. For K3Fe(CN)6, the relative catalytic activity of PgCPR1 is five times that of PgCPR2 and AtCPR1(Table 2). The purified PgCPR1 and PgCPR2 proteins have redox ability, which is higher than that of AtCPR1. In the future study of ginsenoside synthesis biology, using ginseng CPR instead of AtCPR1 may obtain a higher yield.

2.3. Expression Analysis of PgCPRs

Quantitative real-time PCR (qPCR) analysis was utilized to detect gene expression levels in different tissues. Transcripts of PgCPR1 and PgCPR2 exist in all examined tissues (Figure 5A), while PgCPR1 and PgCPR2 were preferentially expressed in the main root and stem. The expression of PgCPR2 was generally higher than that of PgCPR1 in all tissues.
As a plant hormone, MeJA [30] can induce ginsenoside synthesis [31]. The content of PgCPR 1 and PgCPR 2 in the callus was detected by qPCR after being induced in a concentration-dependent manner for 24 h. The results showed that the expression of PgCPR1 and PgCPR2 increased significantly with the increase of MeJA concentration and reached the maximum value in 100 μM MeJA treatments, and then decreased significantly in the 150 μM MeJA treatment. (Figure 5B). Therefore, we added 100 μM MeJA to the ginseng callus medium to study the expression level of PgCPR1 and PgCPR2 at different time points (0,4,8,12,24,48). The qPCRanalysis showed that PgCPR1 and PgCPR2 were slightly expressed before treatment, but are markedly increased to the highest point in 24 or 12 h and then decrease. The expression of PgCPR2 was higher than that of PgCPR1 at different induction times (Figure 5C). The results showed that PgCPR1 and PgCPR2 were induced by MeJA. It is speculated that PgCPR1 and PgCPR2 might be involved in ginsenoside synthesis.

3. Discussion

The current research shows that the cytochrome P450 is the key rate-limiting enzyme of ginsenoside synthesis. Previous reports have shown that CPR is the sole electron donor of all microsomes, and it is also the rate-limiting factor of P450 monooxygenase reaction [32,33]. In this study, we isolated two Panax ginseng cytochrome P450 reductase genes, PgCPR1 and PgCPR2, and identified their protein structure, catalytic activity in vitro for cytochrome c and K3Fe(CN)6, and gene expression level.
In this study, the comparison of two PgCPRs by bioinformatics analysis revealed an 85% similarity in their amino acid sequence. CPR contains several functional domains that are relatively conservative, but N-terminal sequences vary greatly [34]. Using the CPR amino acid sequence of the truncated N-terminal sequence, we can construct phylogenetic trees that are similar to using the full-length sequence [27,35]. The truncation of the N-terminal transmembrane region can increase the solubility of CPR protein but will remain unchanged in its catalytic activity. We conjectured that the N-terminal transmembrane region might be responsible for P450s coupling rather than the reductase activity. Conserved domains in other regions were related to NADPH dependent cytochrome P450 reductase activity.
Based on sequence comparison, PgCPR1 and PgCPR2 belong to class II. Ro et al. believe that CPR class I in general has a short N-terminal, but CPR class II invariably contains a longer N-terminal [27]. Nevertheless, PgCPR2 with N-terminal shorter than PgCPR1 is classified as class II. Amusingly, CPR2 in hybrid poplar, classified as class II, also has a short N-terminal [27]. Andrographis paniculata CPR3 has a very short N-terminal, which is also classified as class II [29]. In their research, Ro et al. think hybrid poplar CPR2 is an exception. According to our results, this appears to be no exception, although rarely seen. Cloning more CPR in other species may lead to more class II CPR with short N-terminal. So, we think that the length of N-terminal is an essential standard, but it is not the unique criterion divided into two classes. Additional sequence information is equally significant for the classification of the two classes.
MeJA is a plant-specific signal regulator, which has been proven to be a valid activator for the biosynthesis of secondary metabolites of several medicinal plants [36]. In ginseng, the up-regulation of critical enzyme genes, for instance, PgCYP450 after MeJA induction led to the accumulation of ginsenoside [37,38,39]. Both CPRs from Panax ginseng were induced after MeJA induction, but PgCPR2 was significantly higher than PgCPR1. The role of MeJA and its correlation with ginsenoside production connote that PgCPRs may play an important role in ginsenoside synthesis. In future experiments, we plan to use transgenic technology to over express and silence PgCPR1 and PgCPR2 genes, which will further clarify their functions in ginsenoside synthesis. The identification of CPRs of Panax ginseng is helpful for illustrating the function of P450 monooxygenase and the biosynthesis pathway of ginsenoside.
The endogenous CPR of yeast is not compatible with plant P450s [40]. Therefore, the WAT11 strain containing AtCPR1 is usually used to validation the function of P450s in different plants [41,42]. The functional plant CPRs can provide electrons for many CYP450 of many different plants, but the matching of P450 with suitable CPRs can significantly improve the production of a secondary metabolite.While the PgCPR1, PgCPR2 and AtCPR2-1 were introduced into engineering yeast respectively, the PPD produced by PgCPR1 and PgCPR2 yeast is greater than twice that of AtCPR2-1 yeast [43]. When CPRs from various plants were paired with CYP450: Uni25647, the titers of 11-oxo-β-amyrin were increased by 1422 times, and the glycyrrhetinic acid and were increased by 945.5 times [44,45]. Although the amino acid sequence of CPRs in plants is highly homologous, their enzyme activities in the same substrate or CYP450 appear to be different. In the production of ginsenosides by yeast metabolic engineering, different CYP450 paired with an appropriate CPR can increase ginsenoside production significantly.
Acting in accordance with the above results, ginseng CPRs can be used in yeast engineering to produce secondary metabolites. Through the establishment of yeast cell factories, more optimization and feasible strategies are offered to improve the supply of ginsenoside. PgCPRs can be used as an ideal biological module in synthetic biology. It can be combined with the existing biological system to redesign or develop completely different ginsenoside production methods. In the future, we will combine PgCPR1 and PgCPR2 with PgDDS, PgCYPs and PgUGTs to determine the best match. This provides a basis for the construction and synthetic biology of ginsenoside yeast engineering bacteria.

4. Materials and Methods

4.1. Plant Materials and Culture Conditions

The Arabidopsis cDNA used in this study was provided by Hunan University. The ginseng used was four-year-old ginseng from the Changbai Mountain Forest of Jilin Province. The callus and adventitious roots of ginseng were induced and subcultured in our laboratory, which laid a foundation for the subsequent experiments.

4.2. Total RNA Extraction and cDNA Synthesis

In this experiment, the total RNA of all plant tissues was extracted with the Tiangen RNA prep Pure Plant kit (Polysaccharides&Polyphenolics-rich) (Tiangen, Beijing, China). Then the concentration and purity of RNA were measured by a spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). The first cDNA was synthesized with the HiScript® II Reverse Transcriptase kit (Vazyme, Nanjing, China) and was stored at −20 °C for future use.

4.3. Isolation and Cloning of PgCPRs Coding Sequences

The primers used in this study are listed in the table (Supplementary Material Table S1). The specific primers were synthesized by the Invitrogen company. According to the PgCPRs gene sequence obtained from the transcriptome database of the laboratory, the target fragment was amplified by PrimeSTAR® Max DNA Polymerase (Takara, Beijing, China). The PCR product was cloned into pCE2 TA/Blunt-Zero Vector by 5minTM TA/Blunt-Zero Cloning Kit (Vazyme, Nanjing, China) and transformed into E. coli DH5 α. The recombinant plasmid was isolated and sequenced by Tsingke Biotechnology.

4.4. Bioinformatic Analysis of PgCPR1 and PgCPR2

Use ORF finder (https://www.ncbi.nlm.nih.gov/orffinder/ (accessed on 7 October 2021)) to find the Open Reading Frame of PgCPR2. The NCBI conserved domains finder (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi (accessed on 7 October 2021)) was used to find the conserved domain of DNA. Calculate assumed molecular weight using ExPASytool (http://ca.expasy.org/tools (accessed on 7 October 2021)). The protein transmembrane helices were predicted by TMHMM Server v. 2.0 (http://www.cbs.dtu.dk/services/TMHMM/ (accessed on 7 October 2021)). Amino acid residues of CPRs were compared with the ClustalW program. Then, DNAMAN software was employed for multiple sequence alignment analysis. Phylogenetic trees are connected by the neighbor-joining method using MEGA 7.0 software. Bootstrap analysis of 1000 replicates was also performed to obtain the level of confidence of each branch. Based on the Arabidopsis CPR2 (ATR2) protein as a three-dimensional homologous model, the third-order structure of ginseng CPR protein was predicted by using the SWISS-MODEL online tool (https://swissmodel.expasy.org/ (accessed on 7 October 2021)).

4.5. Prokaryotic Expression and Protein Purification of Recombinant PgCPR Proteins

All plant CPRs have nearly 60 amino acid transmembrane structures, which are located at the N-terminal. Considering the optimization of the solubility of prokaryotic expression protein, the coding sequence of PgCPRs without the N-terminal membrane-anchored domain was amplified. Using ClonExpress® II one-step cloning Kit (Vazyme, Nanjing, China), the target sequence was inserted into the middle of the Hind III site and EcoRI site of expression vector pET28a, and the expression, purification and kinetic analysis of the binding protein were as described above [31]. After sequencing and validation, the recombinant vector was transformed into E. coli BL21 (DE3) active cells (Qingke, Beijing, China). Positive colony in LB medium (containing 100 μg/mL Kana) at 37 °C and 200 rpm. When OD600 = 1.0, add IPTG to the final concentration of 1 mM and induce at 16 °C for 12 h. The cells were then collected and resuspended in a binding buffer (50mM NaH2PO4, 300 mM NaCl, 10mM imidazole, pH 8.0) and ultrasonically treated on ice for 5 s to 20 times. After centrifugation at 4 °C for 30 min at 12,000 rpm (13,201× g), the Ni-NTA resins were added to the supernatant for affinity purification. The protein samples were separated by SDS-PAGE and stained by Coomassie brilliant blue R250. The concentration of recombinant protein was determined by Pierce™ BCA Protein Assay Kit (Bio-Rad, Hercules, CA, USA).

4.6. Enzyme Activity Assays In Vitro

The enzyme activities of AtCPR1 and PgCPRs in vitro were determined. The experimental methods refer to the reported articles [27]. All incubation and reaction were carried out at 25 °C in 50 mmol/L Tris HCl buffers (pH7.4, containing 0.1 mmol/L EDTA). Cytochrome reduction was monitored by an increase in absorbance at 550 nm. The molar absorption coefficient of cytochrome c was 21 mM−1cm−1. The reduction of K3Fe(CN)6 was monitored at 424 nm (1.02 mM−1cm−1). To determine the kinetic parameters of cytochrome c, 100 μmol/L, NADPH was added to the reaction mixture containing different concentrations of cytochrome c (0–800 μmol/L). The kinetic parameters of NADPH were determined with cytochrome c solution of 100 μmol/L, the concentration range of NADPH was 0–400 μmol/L. The kinetic parameters of K3Fe(CN)6 were determined and 100 μmol/L NADPH was added to the reaction mixture containing different concentrations of K3Fe(CN)6 (0–660 μmol/L). The absorption changes with time were recorded on the epoch spectrophotometer (BioTek, Winooski, VT, USA). GraphPad Prism 5 software was used to carry out nonlinear regression analysis and calculation. At least three independent copies of each sample were tested.

4.7. Quantitative Real-Time PCR

In this paper, the main root, lateral root, stem, leaf, and fruit of Panax ginseng were collected, and the callus of Panax ginseng was induced by MeJA. Expression patterns of PgCPRs were detected by using the cDNA of different tissues or processing stages as the template of qPCR amplification. The real-time PCR was carried out on a Mastercycler ep realplex Real-time fluorescence quantitative PCR instrument (Eppendorf, Germany) by using a Bioteke miRNA qPCR detection kit (SYBR Green) miRNA (Bioteke, Beijing, China). According to the requirements of the real-time PCR system, the primers of the qPCR specific reference gene, PgCPR1 and PgCPR2 are listed in Table S1. The calibration curve was drawn for each primer pair and was quantified by GraphPad prism 6 analysis software. Each sample was tested for at least three independent biological repeats.

Supplementary Materials

The following are available online. Figure S1: Main characteristic steps of the ginsenoside biosynthesis pathway; Figure S2: The conserved domains of PgCPR1 and PgCPR2; Figure S3: The transmembrane domain of PgCPR1 and PgCPR2; Figure S4: The protein tertiary structure analysis of PgCPR1 and PgCPR2; Table S1: Specific primers used in the study.

Author Contributions

Conceptualization, Y.Z.; Data curation, X.Z. (Xian Zou); Formal analysis, T.L. and G.L.; Funding acquisition, X.Z. (Xu Zeng); Software, Y.D.; Supervision, Z.L.; Writing—original draft, X.Z. (Xian Zou); Writing—review & editing, Y.Z., X.Z. (Xu Zeng), T.L., Y.X. and Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

The work was funded by the National Natural Science Foundation of China (Grant No.81673544, 81874332 and 81071821), Central South University Innovation and Entrepreneurship Project by Teacher and Students (Grant No.2017gczd028) and Independent Exploration and Innovation Project for Postgraduates of Central South University (Grant No.2018zzts395).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We sincerely thank Xiaoying Zhao (Hunan University) for providing Arabidopsis cDNA.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yun, T.K. Brief introduction of Panax ginseng C.A. Meyer. J. Korean Med. Sci. 2001, 16, S3–S5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Wen, J.; Zimmer, E.A. Phylogeny and biogeography of Panax, L. (the ginseng genus, araliaceae): Inferences from ITS sequences of nuclear ribosomal DNA. Mol. Phylogenet. Evol. 1996, 6, 167–177. [Google Scholar] [CrossRef]
  3. Qi, L.; Wang, C.; Yuan, C. Isolation and analysis of ginseng: Advances and challenges. Nat. Prod. Rep. 2011, 28, 467–495. [Google Scholar] [CrossRef]
  4. Zheng, S.; Wu, H.; Wu, D. Roles and mechanisms of ginseng in protecting the heart. Chin. J. Integr. Med. 2012, 18, 548–555. [Google Scholar] [CrossRef] [PubMed]
  5. Lee, K.J.; Ji, G.E. The effect of fermented red ginseng on depression is mediated by lipids. Nutr. Neurosci. 2014, 17, 7–15. [Google Scholar] [CrossRef] [PubMed]
  6. Wang, Y.; Chen, J.; Chu, S.; Wang, Y.; Wang, X.; Chen, N.; Zhang, J. Improvement of memory in mice and increase of hippocampal excitability in rats byginsenoside Rg1’s metabolites ginsenoside Rh1 and protopanaxatriol. J. Pharm. Sci. 2009, 109, 504–510. [Google Scholar] [CrossRef] [Green Version]
  7. Wang, Y.; Liu, J.; Zhang, Z.; Bi, P.; Qi, Z.; Zhang, C. Anti-neuroinflammation effect of ginsenoside Rbl in a rat model of Alzheimerdisease. Neurosci. Lett. 2011, 487, 70–72. [Google Scholar] [CrossRef] [PubMed]
  8. Wei, S. Potential therapeutic action of natural products from traditional Chinesemedicine on Alzheimer’s disease animal models targeting neurotrophic factors. Fund. Clin. Pharmacol. 2016, 30, 490–501. [Google Scholar] [CrossRef] [PubMed]
  9. Wang, C.; Anderson, S.; Du, W.; He, T.; Yuan, C. Red ginseng and cancer treatment. Chin. J. Nat. Med. 2016, 14, 7–16. [Google Scholar]
  10. Lee, S.Y.; Kim, G.T.; Roh, S.H.; Song, J.; Kim, H.; Hong, S.; Kwon, S.W.; Park, J.H. Proteome changes related to the anti-cancer activity of HT29 cells by thetreatment of ginsenoside Rd. Die Pharmazie 2009, 64, 242–247. [Google Scholar]
  11. Son, K.; Choi, K.R.; Lee, S.J.; Lee, H. Immunogenic cell death Induced by ginsenoside Rg3: Significance in dendritic cell-based anti-tumor immunotherapy. Immune Netw. 2016, 16, 75–84. [Google Scholar] [CrossRef] [Green Version]
  12. Quan, K.; Liu, Q.; Wan, J.; Zhao, Y.; Guo, R.; Alolga, R.N.; Li, P.; Qi, L. Rapid preparation of rare ginsenosides by acid transformation and theirstructure-activity relationships against cancer cells. Sci. Rep. 2015, 5, 8598. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Shibata, S. Chemistry and cancer preventing activities of ginseng saponins and some relatedtriterpenoid compounds. J. Korean Med. Sci. 2001, 16, S28–S37. [Google Scholar] [CrossRef] [Green Version]
  14. Han, J.Y.; Hwang, H.S.; Choi, S.W.; Kim, H.J.; Choi, Y.E. Cytochrome P450 CYP716A53v2 catalyzes the formation of protopanaxatriol from protopanaxadiol during ginsenoside biosynthesis in Panax ginseng. Plant Cell Physiol. 2012, 53, 1535–1545. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Han, J.; Kim, H.; Kwon, Y.; Choi, Y. The Cyt P450 enzyme CYP716A47 catalyzes the formation of protopanaxadiol from dammarenediol-II during ginsenoside biosynthesis in Panax ginseng. Plant Cell Physiol. 2011, 52, 2062–2073. [Google Scholar] [CrossRef]
  16. Han, J.Y.; Kim, M.J.; Ban, Y.W.; Hwang, H.S.; Choi, Y.E. The involvement of beta-amyrin 28-oxidase (CYP716A52v2) in oleanane-type ginsenoside biosynthesis in Panax ginseng. Plant Cell Physiol. 2013, 54, 2034–2046. [Google Scholar] [CrossRef] [PubMed]
  17. Han, J.Y.; Kwon, Y.S.; Yang, D.C.; Jung, Y.R.; Choi, Y.E. Expression and RNA interference-induced silencing of the dammarenediol synthase gene in Panax ginseng. Plant Cell Physiol. 2006, 47, 1653–1662. [Google Scholar] [CrossRef]
  18. Schuhegger, R.; Nafisi, M.; Mansourova, M.; Petersen, B.L.; Olsen, C.E.; Svatos, A.; Halkier, B.A.; Glawischnig, E. CYP71B15 (PAD3) catalyzes the final step in camalexin biosynthesis. Plant Physiol. 2006, 141, 1248–1254. [Google Scholar] [CrossRef] [Green Version]
  19. Isin, E.M.; Guengerich, F.P. Substrate binding to cytochromes P450. Anal. Bioanal. Chem. 2008, 392, 1019–1030. [Google Scholar] [CrossRef] [Green Version]
  20. Renault, H.; Bassard, J.; Hamberger, B.; Werck-Reichhart, D. Cytochrome P450-mediated metabolic engineering: Current progress and futurechallenges. Curr. Opin. Plant Biol. 2014, 19, 27–34. [Google Scholar] [CrossRef] [Green Version]
  21. Nebert, D.W.; Gonzalez, F.J. P450 genes: Structure, evolution, and regulation. Annu. Rev. Biochem. 1987, 56, 945–993. [Google Scholar] [CrossRef] [PubMed]
  22. Bhat, W.W.; Rana, S.; Dhar, N.; Razdan, S.; Pandith, S.A.; Vishwakarma, R.; Lattoo, S.K. An inducible NADPH-cytochrome P450 reductase from Picrorhiza kurrooa—animperative redox partner of cytochrome P450 enzymes. Funct. Integr. Genomic. 2014, 14, 381–399. [Google Scholar] [CrossRef]
  23. Mizutani, M.; Ohta, D. Diversification of P450 genes during land plant evolution. Annu. Rev. Plant Biol. 2010, 61, 291–315. [Google Scholar] [CrossRef] [PubMed]
  24. Bak, S.; Beisson, F.; Bishop, G.; Hamberger, B.; Hofer, R.; Paquette, S.; Werck-Reichhart, D. Cytochromes p450. Arab. Book 2011, 9, e144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Park, S.; Kim, Y.; Rupasinghe, S.G.; Schuler, M.A.; Back, K. Rice P450 reductases differentially affect P450-mediated metabolism in bacterial expression systems. Bioproc. Biosyst. Eng. 2013, 36, 325–331. [Google Scholar] [CrossRef]
  26. Yang, C.; Lu, S.; Mao, Y.; Wang, L.; Chen, X. Characterization of two NADPH: Cytochrome P450 reductases from cotton (Gossypium hirsutum). Phytochemistry 2010, 71, 27–35. [Google Scholar] [CrossRef]
  27. Ro, D.; Ehlting, J.; Douglas, C.J. Cloning, functional expression, and subcellular localization of multiple NADPH-Cytochrome P450 reductases from Hybrid Poplar. Plant Physiol. 2002, 130, 1837–1851. [Google Scholar] [CrossRef] [Green Version]
  28. Tsou, C.; Matsunaga, S.; Okada, S. Molecular cloning and functional characterization of NADPH-dependent cytochrome P450 reductase from the green microalga Botryococcus braunii, B race. J. Biosci. Bioeng. 2018, 125, 30–37. [Google Scholar] [CrossRef]
  29. Lin, H.; Wang, J.; Qi, M.; Guo, J.; Rong, Q.; Tang, J.; Wu, Y.; Ma, X.; Huang, L. Molecular cloning and functional characterization of multiple NADPH-cytochrome P450 reductases from Andrographis paniculata. Int. J. Biol. Macromol. 2017, 102, 208–217. [Google Scholar] [CrossRef]
  30. Surguchov, A.; Bernal, L.; Surguchev, A.A. Phytochemicals as regulators of genes involved in synucleinopathies. Biomolecules 2021, 11, 624. [Google Scholar] [CrossRef]
  31. Cheong, J.; Choi, Y.D. Methyl jasmonate as a vital substance in plants. Trends Genet. 2003, 19, 409–413. [Google Scholar] [CrossRef]
  32. Jensen, K.; Moller, B.L. Plant NADPH-cytochrome P450 oxidoreductases. Phytochemistry 2010, 71, 132–141. [Google Scholar] [CrossRef]
  33. Riddick, D.S.; Ding, X.; Wolf, C.R.; Porter, T.D.; Pandey, A.V.; Zhang, Q.; Gu, J.; Finn, R.D.; Ronseaux, S.; McLaughlin, L.A.; et al. NADPH-cytochrome P450 oxidoreductase: Roles in physiology, pharmacology, and toxicology. In Drug Metabolism and Disposition: The Biological Fate of Chemicals, 1st ed.; ASPET: Maryland, MD, USA, 2013; Volume 41, pp. 12–23. [Google Scholar]
  34. Ming, W.; David, L.R.; Rosemary, P.; Thomas, M.S. Three-dimensional structure of NADPH-cytochrome P450 reductase. Proc. Natl. Acad. Sci. USA 1997, 94, 8411. [Google Scholar]
  35. Zhao, H.; Wang, J.; Tang, Q.; Mo, C.; Guo, J.; Chen, T.; Lin, H.; Tang, J.; Guo, L.; Huang, L.; et al. Functional expression of two NADPH-cytochrome P450 reductases from Siraitia grosvenorii. Int. J. Biol. Macromol. 2018, 120, 1515–1524. [Google Scholar] [CrossRef] [PubMed]
  36. Wasternack, C. Jasmonates: An update on biosynthesis, signal transduction and action in plant stress response, growth and development. Ann. Bot. 2007, 100, 681–697. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Lee, M.; Jeong, J.; Seo, J.; Shin, C.; Kim, Y.; In, J.; Yang, D.; Yi, J.; Choi, Y. Enhanced triterpene and phytosterol biosynthesis in Panax ginseng overexpressing the squalene synthase gene. Plant Cell Physiol. 2004, 45, 976–984. [Google Scholar] [CrossRef] [Green Version]
  38. Kim, Y.S.; Hahn, E.J.; Murthy, H.N.; Paek, K.Y. Adventitious root growth and ginsenoside accumulation in Panax ginseng cultures as affected by methyl jasmonate. Biotechnol. Lett. 2004, 26, 1619–1622. [Google Scholar] [CrossRef]
  39. Kim, O.T.; Bang, K.H.; Kim, Y.C.; Hyun, D.Y.; Kim, M.Y.; Cha, S.W. Upregulation of ginsenoside and gene expression related to triterpene biosynthesis in ginseng hairy root cultures elicited by methyl jasmonate. Plant Cell Tiss. Org. Cult. 2009, 98, 25–33. [Google Scholar] [CrossRef]
  40. Guo, J.; Zhou, Y.J.; Hillwig, M.L.; Shen, Y.; Yang, L.; Wang, Y.; Zhang, X.; Liu, W.; Peters, R.J.; Chen, X.; et al. CYP76AH1 catalyze the turnover of miltiradiene in tanshinones biosynthesis and enables heterologous production of ferruginol in yeasts. Proc. Natl. Acad. Sci. USA 2013, 110, 12108–12113. [Google Scholar] [CrossRef] [Green Version]
  41. Guo, J.; Ma, X.; Cai, Y.; Ma, Y.; Zhan, Z.; Zhou, Y.J.; Liu, W.; Guan, M.; Yang, J.; Cui, G.; et al. Cytochrome P450 promiscuity leads to a bifurcating biosynthetic pathway for tanshinones. New Phytol. 2016, 210, 525–534. [Google Scholar] [CrossRef] [Green Version]
  42. Zhou, C.; Li, J.; Li, C.; Zhang, Y. Improvement of betulinic acid biosynthesis in yeast employing multiple strategies. BMC Biotechnol. 2016, 16, 59. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Yan, X.; Fan, Y.; Wei, W.; Wang, P.; Liu, Q.; Wei, Y.; Zhang, L.; Zhao, G.; Yue, J.; Zhou, Z. Production of bioactive ginsenoside compound K in metabolically engineered yeast. Cell Res. 2014, 24, 770–773. [Google Scholar] [CrossRef] [PubMed]
  44. Zhu, M.; Wang, C.; Sun, W.; Zhou, A.; Wang, Y.; Zhang, G.; Zhou, X.; Huo, Y.; Li, C. Boosting 11-oxo-β-amyrin and glycyrrhetinic acid synthesis in Saccharomyces cerevisiae via pairing novel oxidation and reduction system from legume plants. Metab. Eng. 2018, 45, 43–50. [Google Scholar] [CrossRef]
  45. Seki, H.; Sawai, S.; Ohyama, K.; Mizutani, M.; Ohnishi, T.; Sudo, H.; Fukushima, E.O.; Akashi, T.; Aoki, T.; Saito, K.; et al. Triterpene functional genomics in licorice for identification of CYP72A154 involved in the biosynthesis of glycyrrhizin. Plant Cell. 2011, 23, 4112–4123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Molecules 26 06654 g001 550
Figure 1. The amino acid sequence of PgCPRs was compared with that of other CPR homologs. The CPR for sequence alignment is from Panax ginseng (PgCPR1: MT160178, PgCPR2: MT160179), Artemisia annua (AaCPR: ABM88789), Siraitia grosvenorii (SgCPR: AYE89265), Andrographis paniculate (ApCPR1: AQT38168, ApCPR2: AQT38169) and Arabidopsis thaliana (AtCPR1: NP194183, AtCPR2: NP194750) by the program Clustal W. The conservative domains are underlined as follows: Membrane anchor (blue), FMN-binding domain (purple), P450-binding domain (red), connecting domain (aurantium), FAD-binding domain (green) and NADPH binding domain (brown).
Figure 1. The amino acid sequence of PgCPRs was compared with that of other CPR homologs. The CPR for sequence alignment is from Panax ginseng (PgCPR1: MT160178, PgCPR2: MT160179), Artemisia annua (AaCPR: ABM88789), Siraitia grosvenorii (SgCPR: AYE89265), Andrographis paniculate (ApCPR1: AQT38168, ApCPR2: AQT38169) and Arabidopsis thaliana (AtCPR1: NP194183, AtCPR2: NP194750) by the program Clustal W. The conservative domains are underlined as follows: Membrane anchor (blue), FMN-binding domain (purple), P450-binding domain (red), connecting domain (aurantium), FAD-binding domain (green) and NADPH binding domain (brown).
Molecules 26 06654 g001
Molecules 26 06654 g002 550
Figure 2. Phylogenetic analysis of PgCPRs. (A) The phylogenetic tree was constructed by the neighbor-joining method of MEGA 7.0, and 27 amino acid sequences of CPR protein were analyzed, which come from the following plant species: Arachis hypogaea (AhCPR: XP025668510), Artemisia annua (AaCPR: ABM88789), Actinidia Chinensis (AcCPR: PSS36215), Andrographis paniculate (ApCPR1: AQT38168, ApCPR2: AQT38169), Arabidopsis thaliana (AtCPR1: NP194183, AtCPR2: NP194750), Camptotheca acuminate (CaCPR: AJW67229), Dendrobium catenatum (DcCPR: PKU78870), Durio zibethinus (DzCPR: XP022729926), Gossypium hirsutum (GhCPR2: NP001314398), Glycine soja (GsCPR: KHN17598), Lotus japonicus (LjCPR: BAG68945 ), Malus domestica (MdCPR: XP028963051), Medicago truncatula (MtCPR: RHN70186), Nicotiana attenuate (NaCPR: OIT38777), Ocimum basilicum (ObCPR: ANW46529), Panax ginseng (PgCPR1: MT160178, PgCPR2: MT160179), Populus trichocarpa (PtCPR: XP002324469), Rosa chinensis (RcCPR: XP024171228), Santalum album (SaCPR: ANQ46483), Siraitia grosvenorii (SgCPR: AYE89265), Solanum Lycopersicum (SlCPR: XP004242931), Salvia miltiorrhiza (SmCPR: CBX24555), Stevia rebaudiana (SrCPR: ABB88839), Vigna radiata (VrCPR: NP001304239). The scale represents the phylogenetic distance calculated from the number of differences. (B,C), The three-dimensional structure diagrams of the PgCPR1 and PgCPR2, constructed by Swiss-Model software.
Figure 2. Phylogenetic analysis of PgCPRs. (A) The phylogenetic tree was constructed by the neighbor-joining method of MEGA 7.0, and 27 amino acid sequences of CPR protein were analyzed, which come from the following plant species: Arachis hypogaea (AhCPR: XP025668510), Artemisia annua (AaCPR: ABM88789), Actinidia Chinensis (AcCPR: PSS36215), Andrographis paniculate (ApCPR1: AQT38168, ApCPR2: AQT38169), Arabidopsis thaliana (AtCPR1: NP194183, AtCPR2: NP194750), Camptotheca acuminate (CaCPR: AJW67229), Dendrobium catenatum (DcCPR: PKU78870), Durio zibethinus (DzCPR: XP022729926), Gossypium hirsutum (GhCPR2: NP001314398), Glycine soja (GsCPR: KHN17598), Lotus japonicus (LjCPR: BAG68945 ), Malus domestica (MdCPR: XP028963051), Medicago truncatula (MtCPR: RHN70186), Nicotiana attenuate (NaCPR: OIT38777), Ocimum basilicum (ObCPR: ANW46529), Panax ginseng (PgCPR1: MT160178, PgCPR2: MT160179), Populus trichocarpa (PtCPR: XP002324469), Rosa chinensis (RcCPR: XP024171228), Santalum album (SaCPR: ANQ46483), Siraitia grosvenorii (SgCPR: AYE89265), Solanum Lycopersicum (SlCPR: XP004242931), Salvia miltiorrhiza (SmCPR: CBX24555), Stevia rebaudiana (SrCPR: ABB88839), Vigna radiata (VrCPR: NP001304239). The scale represents the phylogenetic distance calculated from the number of differences. (B,C), The three-dimensional structure diagrams of the PgCPR1 and PgCPR2, constructed by Swiss-Model software.
Molecules 26 06654 g002
Molecules 26 06654 g003 550
Figure 3. Expression and purification of recombinant PgCPRs and AtCPR1. The SDS-PAGE analysis of PgCPRs and AtCPR1 stained with Coomassie blue R-250. Lane M: protein marker; Lane 1: The crude enzyme in BL21 without IPTG. Lane 2: IPTG induced crude enzyme; Lane 3: The purified protein with his-tag.
Figure 3. Expression and purification of recombinant PgCPRs and AtCPR1. The SDS-PAGE analysis of PgCPRs and AtCPR1 stained with Coomassie blue R-250. Lane M: protein marker; Lane 1: The crude enzyme in BL21 without IPTG. Lane 2: IPTG induced crude enzyme; Lane 3: The purified protein with his-tag.
Molecules 26 06654 g003
Molecules 26 06654 g004 550
Figure 4. Kinetic parameters of the recombinant protein. (A–C), the kinetic parameters of recombinant protein to cytochrome c were measured by 550 nm absorbance increase method; (D–F), the kinetic parameters of NADPH were measured by cytochrome reduction method; (G–I) the kinetic parameters of K3Fe(CN)6 were studied by 424 nm absorbance change method.
Figure 4. Kinetic parameters of the recombinant protein. (A–C), the kinetic parameters of recombinant protein to cytochrome c were measured by 550 nm absorbance increase method; (D–F), the kinetic parameters of NADPH were measured by cytochrome reduction method; (G–I) the kinetic parameters of K3Fe(CN)6 were studied by 424 nm absorbance change method.
Molecules 26 06654 g004
Molecules 26 06654 g005 550
Figure 5. Expression analysis of PgCPRs. (A) The relative expression level of PgCPRs in different tissues of ginseng was calculated by using the expression level of PgCPR1 at the main root as standard 1. (B) The expression level of PgCPR1 at 0 μM MeJA was used as standard 1 to calculate the relative expression level of PgCPRs in different concentrations of MeJA at 24 h. (C) The expression level of PgCPRs at different time points was induced by 100 μM MeJA. The expression level of PgCPR1 at 0 h was used as the standard. All the above expression levels were standardized according to the expression level of β-actin. The standard deviation of three independent repeated tests is expressed by error bars. The asterisk indicates that there is significant difference in t-test, * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 5. Expression analysis of PgCPRs. (A) The relative expression level of PgCPRs in different tissues of ginseng was calculated by using the expression level of PgCPR1 at the main root as standard 1. (B) The expression level of PgCPR1 at 0 μM MeJA was used as standard 1 to calculate the relative expression level of PgCPRs in different concentrations of MeJA at 24 h. (C) The expression level of PgCPRs at different time points was induced by 100 μM MeJA. The expression level of PgCPR1 at 0 h was used as the standard. All the above expression levels were standardized according to the expression level of β-actin. The standard deviation of three independent repeated tests is expressed by error bars. The asterisk indicates that there is significant difference in t-test, * p < 0.05, ** p < 0.01, *** p < 0.001.
Molecules 26 06654 g005
Table
Table 1. Specific activities of recombinant CPRs in reducing cytochrome c (100 μmol/L) with 100 μmol/L of NADPH or NADH. All data represent means ± standard deviations (SD). ND: Not detectable.
Table 1. Specific activities of recombinant CPRs in reducing cytochrome c (100 μmol/L) with 100 μmol/L of NADPH or NADH. All data represent means ± standard deviations (SD). ND: Not detectable.
Specific Activity (μmol/min/mg Protein)
NADPHNADH
PgCPR17.23 ± 0.45ND
PgCPR26.53 ± 0.37ND
ATCPR15.56 ± 0.26ND
Table
Table 2. Kinetic analysis of recombinant CPRs. All data represent the mean ± SDS calculated by nonlinear regression analysis using GraphPad Prism 7.0 software.
Table 2. Kinetic analysis of recombinant CPRs. All data represent the mean ± SDS calculated by nonlinear regression analysis using GraphPad Prism 7.0 software.
Vma (μmol/min/mg)Km (μmol/L)Kcat (min−1)Kcat/Km
Cytochrome cPgCPR119.85 ± 0.5344159.4 ± 12.1397 ± 10.6882.49059
PgCPR218.72 ± 0.4907157.3 ± 11.69374.4 ± 9.8142.380165
ATCPR112.49 ± 0.3651151.9 ± 12.78249.8 ± 7.3021.644503
NADPHPgCPR19.962 ± 0.290562.9 ± 5.807199.24 ± 5.813.167568
PgCPR210.14 ± 0.261185.27 ± 6.07202.8 ± 5.2222.378328
ATCPR18.235 ± 0.299884.61 ± 8.548164.7 ± 5.9961.946578
K3Fe(CN)6PgCPR116.67 ± 0.279164.91 ± 5.639333.4 ± 5.5825.136343
PgCPR24.441 ± 0.0890785.3 ± 7.60188.82 ± 1.78141.041266
ATCPR14.119 ± 0.09382.73 ± 8.44282.38 ± 1.860.995769
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zou, X.; Zhang, Y.; Zeng, X.; Liu, T.; Li, G.; Dai, Y.; Xie, Y.; Luo, Z. Molecular Cloning and Identification of NADPH Cytochrome P450 Reductase from Panax ginseng. Molecules 2021, 26, 6654. https://doi.org/10.3390/molecules26216654

AMA Style

Zou X, Zhang Y, Zeng X, Liu T, Li G, Dai Y, Xie Y, Luo Z. Molecular Cloning and Identification of NADPH Cytochrome P450 Reductase from Panax ginseng. Molecules. 2021; 26(21):6654. https://doi.org/10.3390/molecules26216654

Chicago/Turabian Style

Zou, Xian, Yue Zhang, Xu Zeng, Tuo Liu, Gui Li, Yuxin Dai, Yuanzhu Xie, and Zhiyong Luo. 2021. "Molecular Cloning and Identification of NADPH Cytochrome P450 Reductase from Panax ginseng" Molecules 26, no. 21: 6654. https://doi.org/10.3390/molecules26216654

Article Metrics

Back to TopTop