Identification and Characterization of a Ginsenoside-Transforming β-Glucosidase from Pseudonocardia sp. Gsoil 1536 and Its Application for Enhanced Production of Minor Ginsenoside Rg2(S)

Abstract

The ginsenoside Rg2(S), which is one of the pharmaceutical components of ginseng, is known to have neuroprotective, anti-inflammation, and anti-diabetic effects. However, the usage of ginsenoside Rg2(S) is restricted owing to the small amounts found in white and red ginseng. To enhance the production of ginsenoside Rg2(S) as a 100 gram unit with high specificity, yield, and purity, an enzymatic bioconversion method was developed to adopt the recombinant glycoside hydrolase (BglPC28), which is a ginsenoside-transforming recombinant β-glucosidase from Pseudonocardia sp. strain Gsoil 1536. The gene, termed bglPC28, encoding β-glucosidase (BglPC28) belonging to the glycoside hydrolase family 3 was cloned. bglPC28 consists of 2,232 bp (743 amino acid residues) with a predicted molecular mass of 78,975 Da. This enzyme was overexpressed in Escherichia coli BL21(DE3) using a GST-fused pGEX 4T-1 vector system. The optimum conditions of the recombinant BglPC28 were pH 7.0 and 37°C. BglPC28 can effectively transform the ginsenoside Re to Rg2(S); the Km values of PNPG and Re were 6.36±1.10 and 1.42±0.13 mM, respectively, and the Vmax values were 40.0±2.55 and 5.62±0.21 µmol min−1 mg−1 of protein, respectively. A scaled-up biotransformation reaction was performed in a 10 L jar fermenter at pH 7.0 and 30°C for 12 hours with a concentration of 20 mg/ml of ginsenoside Re from American ginseng roots. Finally, 113 g of Rg2(S) was produced from 150 g of Re with 84.0±1.1% chromatographic purity. These results suggest that this enzymatic method could be usefully exploited in the preparation of ginsenoside Rg2(S) in the cosmetics, functional food, and pharmaceutical industries.

About the Authors:

Juan Du

Contributed equally to this work with: Juan Du, Chang-Hao Cui

Affiliations KAIST Institute for Biocentury, Korea Advanced Institute of Science and Technology, Yuseong-gu, Daejeon, Republic of Korea, College of Biotechnology, Dalian Polytechnic University, Ganjingzi-qu, Dalian, P. R. China

Chang-Hao Cui

Contributed equally to this work with: Juan Du, Chang-Hao Cui

Affiliation: Department of Biological Sciences, Korea Advanced Institute of Science and Technology, Yuseong-gu, Daejeon, Republic of Korea

Sung Chul Park

Affiliation: Intelligent Synthetic Biology Center, Yuseong-gu, Daejeon, Republic of Korea

Jin-Kwang Kim

Affiliation: KAIST Institute for Biocentury, Korea Advanced Institute of Science and Technology, Yuseong-gu, Daejeon, Republic of Korea

Hong-Shan Yu

Affiliation: College of Biotechnology, Dalian Polytechnic University, Ganjingzi-qu, Dalian, P. R. China

Feng-Xie Jin

Affiliation: College of Biotechnology, Dalian Polytechnic University, Ganjingzi-qu, Dalian, P. R. China

Changkai Sun

Affiliation: Institute for Brain Disorders, Dalian Medical University, Dalian, P.R. China

Sun-Chang Kim

Affiliations KAIST Institute for Biocentury, Korea Advanced Institute of Science and Technology, Yuseong-gu, Daejeon, Republic of Korea, Department of Biological Sciences, Korea Advanced Institute of Science and Technology, Yuseong-gu, Daejeon, Republic of Korea, Intelligent Synthetic Biology Center, Yuseong-gu, Daejeon, Republic of Korea

Wan-Taek Im

* E-mail: [email protected]

Affiliation: Department of Biotechnology, Hankyoung National University, Chungang-no Anseong-si, Republic of Korea

Introduction

Ginseng, an important herbal medicine, has been widely used for thousands of years in East Asia and has been popularized in the West during the past decades [1], [2]. Reports show that ginseng has a range of pharmacological and therapeutic uses [3], [4], [5], [6], [7]. The ginseng root consists of ginsenosides, polysaccharides, peptide, polyacetylenic alcohols, and fatty acids [8]. Ginsenosides are the major active components of ginseng, and they appear to be responsible for the principle pharmacological activities of ginseng, including vasorelaxation and anti-neoplastic, anti-diabetic, anti-inflammation, and anti-oxidant effects [9], [10], [11]. Ginsenosides can be categorized as protopanaxadiol (PPD), protopanaxatriol (PPT), and oleanane saponins, based on the structure of the aglycon, with a dammarane skeleton [12]. The PPD- and PPT-type ginsenosides are further classified into subgroups based on the position and number of sugar moieties attached to the aglycon at positions C3 or C6 and C20. For example, one of the largest PPD-type ginsenosides, Rb1, contains 4 glucose moieties, with two attached via glycosidic linkages to the C3 and C20 positions, respectively, of the aglycon (Fig. 1). They may also be classified as major or minor ginsenosides based on the amount found in cultivated ginseng.

[Figure omitted. See PDF.]

Figure 1. Chemical structures of protopanaxadiol and protopanaxatriol ginsenosides.

The ginsenosides represented here are all (S)-type ginsenosides. glc, β-D-glucopyranosyl; arap, α-L-arabinopyranosyl; araf, α-L-arabinofuranosyl; rha, α-L-rhamnopyranosyl; Gyp, gypenoside; C, compound; C-O, compound O; C-Y, compound Y, C-Mc1, compound Mc1; C-Mc, compound Mc; C-K, compound K.

https://doi.org/10.1371/journal.pone.0096914.g001

Ginsenosides Rb1, Rb2, Rc, Rd, Re, and Rg1 are major ginsenosides thatmake up more than 90% of total ginseng ginsenosides [3], [12]. Minor ginsenosides (Rg3, Rh2, F2, C-K, Rg2, Rh1, and F1) that are deglycosylated from the major ginsenosides exist in smaller amounts. As secondary metabolite compounds, the minor ginsenosides have some chemical reactivity that the major ginsenosides do not. Furthermore, there is emerging evidence that the minor ginsenosides have more important pharmaceutical effects, such as anti-cancer, anti-diabetic, anti-oxidative, and anti-aging effects, than the glycosylated major ginsenosides [10], [11], [13], [14], [15].

20(S)-GinsenosideRg2 [hereafter Rg2(S)] is a native ginsenoside found in Korean white and red ginseng (Panax ginseng, C.A. Meyer), and accounts forless than 0.02% of dried mass in red ginseng (a heat-treated ginseng with more minor ginsenosides) [16]. Ginsenoside Rg2(S) has significant neuroprotective pharmaceutical effects. It inhibits mitochondrial permeability transition pores in rat brain as a neuroprotective agent with mechanisms of anti-oxidation and anti-apoptosis [17]. Rg2 could hence potentially be used in a treatment strategy of Alzheimer’s disease [18]. Second, Rg2 improved neurological performance and memory ability of vascular dementia rats [19]. Other investigations suggest that Rg2 can protect cells against UVB-induced genotoxicityby increasing DNA repair and decreasing apoptosis [20], [21]. Rg2 was noted to have a significant impacton improving hemodynamic state and activity SOD (superoxide dismutase) in hemorrhagic shock dogs [22]. It was further reported that Rg2 inhibits leukocyte adhesion on the vascular wall in human cells, providing protection against vascular inflammatory disease [23]. Finally, Rg2 inhibits hepatic glucose production and showed potential for therapeutic treatment for type 2 diabetic patients [24].

Rg2(S) can be easily produced through a weakacid/high temperature process, i.e. the Korean red ginseng process. However, this process creates by products such as Rg2(R), Rg4, and Rg6. Rg2(S) consequently has low concentrations in Korean red ginseng or red ginseng extract [16], [25], [26]. Another method of obtaining Rg2(S) is to use microbes or enzymatic processes. One example is Mucilaginibacter composti TR6-03T, which has β-glucosidase activity and showed an ability to convert ginsenoside Re into Rg2(S) [27]. Although several ginsenoside hydrolyzing recombinant enzymes have been constructed, the majority of these have hydrolysis ability with respect to protopanaxadiol type ginsenosides (i.e. Rb1, Rb2, Rc, Rd) rather than protopanaxatriol type ginsenosides (i.e. Re, Rf, Rg1) [28], [29], [30], [31], [32], [33], [34], [35], [36]. Among these, three enzymes were reported to have hydrolysis activity for the glucose moiety at the C20 position of PPT aglycon [28], [29], [35]. However, Cui et al [28] and Quanet al [35] only conducted a simple enzyme characterization, without further scale-up or process engineering. Cui et al [29] also conducted enzyme characterization and subsequently obtained a gram scale Rg2(S) using a PPT type ginsenoside mixture. However, the yield was low (33.4%) and the final purification step using silica resin is complicated.

The experiment in this study was designed to overcome the above disadvantagesand to meet the industrial demand for mass production of ginsenoside Rg2(S) and fulfill the original purpose of its application as are combinant enzyme. From a dozen bacterial strains capable of transforming ginsenosides, we selected Pseudonocardia sp. Gsoil 1536 on the basis of its high activity for transforming Re into Rg2(S). A genome fosmid library was constructed, the key ginsenoside-transforming β-glucosidase (BglPC28) was identified, its novel gene was cloned for production as a recombinant enzyme, and the enzymatic properties and substrate specificities of the recombinant enzyme were thoroughly investigated. Treatment of purified ginsenoside Re with BglPC28 followed by purification yielded 100 gram-scale Rg2(S) with high purity. This is the first report of the 100 gram-scale production of high-purity Rg2(S) by the application of a recombinant glycoside hydrolase.

Materials and Methods

2.1. Materials

Theginsenoside Re (purity: 87.6%) from the root of Panax quinquefolius from FusongBiotech Co. Ltd. (China) was used as the substrate in the currentinvestigation. Ginsenosides standards which are over 98% purity such as Rb1, Rb2, Rc, Rd, Rg2(S), Rh2(S), F2, compound K (C-K), protopanaxadiol (PPD), Rg1, Re, Rg2(S), Rh1(S), and protopanaxatriol (PPT) were purchased from Nanjing Zelang Medical Technology Co. Ltd. (China). Methanol and acetonitrile with HPLC grade were obtained from Merck (Darmstadt, Germany). The other chemicals used in this study were a minimum of analytical reagent grade, and the sources are noted individually in the Methods section. Pseudonocardia sp. Gsoil 1536, which has ginsenoside-hydrolyzing activity, was isolated from the soil sample in ginseng filed, Pocheon Province, South Korea, and cultivated on R2A agar (BD, USA) under aerobic conditions at 30°C and used for the gene cloning experiment. Escherichia coli BL21 (DE3) and pGEX 4T-1 plasmid (GE Healthcare, USA) were used as host, and expression vectorsources, respectively. The recombinant E. coli for protein expression was cultivated in a Luria-Bertani (LB) medium supplemented with ampicillin (100 mg/l).

2.2. Fosmid Library Construction and Fosmid Sequencing

A CopyControl™ Fosmid Production kit (Epicentre Technologies, WI) was used to clone the ginsenoside hydrolyzing glycosidase gene from Pseudonocardia sp. Gsoil 1536. A fosmid library was constructed according to the manufacturer’s protocol. Infected E. coli was transferred onto LB plates supplemented with 40 µg/ml X-Glc (5-bromo-4-chloro-3-indolyl β-D-glucopyranoside) and 12.5 µg/ml chloramphenicol and then incubated at 37°C for 16 h. The blue color clones were selected as putative ginsenoside hydrolyzing clones. After confirmation of the ginsenoside-hydrolyzing activity by a TLC assay, one clone was selected for fosmid sequencing. Fosmid DNA was purified according to the manufacturer’s protocols (Fosmid MAX DNA purification kit, Epicentre, WI) and was sequenced by Macrogen Co. Ltd. (Korea). The final sequences assembly procedure was conducted by the SeqMan program in the DNASTAR package (DNASTAR, WI), which yieldedtwo contigs (14.3 - and 4.7 kb).

2.3. Phylogenetic Analysis of BglPC28

Database homology search was performed with BLAST programprovided by NCBI. Sequences of the characterized glycosyl hydrolases were obtained from the CAZY database [Carbohydrate-Active enZymes database (http://www.cazy.org)], and multiple alignments were performed using the CLUSTAL_X program [37]. Gaps were edited in the BioEdit program [38], and evolutionary distances were calculated using the Kimura two-parameter model [39]. A phylogenetic tree was constructed using the neighbor-joining method [40] in the MEGA5 Program [41], with bootstrap values based on 1000 replicates [42]. Furthermore, the multiple amino acid sequence alignment and the conserved patterns of discrete amino acid sequences of BglPC28 and known the most homologousβ-glucosidases were performed by using ClustalW2 program (http://www.ebi.ac.uk/Tools/msa/clustalw2/).

2.4. Molecular Cloning, Expression, and Purification of Recombinant BglPC28

The assembled DNA sequence was analyzed using the ORF Finder program on the NCBI website (www.ncbi.nlm.nih.gov/gorf). Predicted ORFs were subjected to a similarity search using BLASTP, which identified two putative open reading frames of a β-glucosidase belonging to glycosyl hydrolase family 3. The sequence of the oligonucleotide primers used for gene cloning was based on the DNA sequence of bglPC28 (GenBank accession no. JX960416). Forward (5′- GGTTCCGCGTGGATCCATCGACCCCGTTGATCTCACCCTC-3′) and reverse primers (5′- GATGCGGCCGCTCGAGCTAGTTTTGTCCGACGTGTTGGGG-3′) were designed as the primers to introduce the BamHI and XhoI restriction sites (underline), respectively, and were synthesized by BioneerCo. Ltd. (Daejeon, Korea). The amplified DNA fragment obtained from the PCR was purified and insertedinto the pGEX 4T-1GST fusion vector digested with BamHI and XhoI using an EzCloning Kit (Enzynomics Co. Ltd., Korea). The resulting recombinant pGEX-bglPC28 was transformed into E. coli BL21(DE3). The E. coli BL21(DE3) harboring the recombinant plasmid was grown in an LB-ampicillin medium at 37°C until the culture reached an OD600 of 0.6, at which point the protein expression was induced through the addition of 0.1 mM isopropyl-β-D-thiogalactopyranoside (IPTG). The bacterial cells were incubated for a further 24 h at 22°C and were then harvested via centrifugation at 13,000 rpm for 15 min at 4°C. The cells were washed twice with a solution consisting of 100 mM sodium phosphate and 1% Triton X-100 (pH 7.0); then, they were resuspended in 100 mM sodium phosphate (pH 7.0). The cells were disrupted via ultrasonication (Vibra-cell, Sonics & Materials, CT, USA). The intact cells and debris were removed via centrifugation at 13,000 rpmfor 15 min at 4°C in order to obtain the crude cell extract. The GST tag was purified using the GST bind agarose resin (Elpisbiotech Co. Ltd, Korea). The homogeneity of the protein was assessed using 10% SDS-PAGE and an EZ-Gel staining solution (Daeillab Co. Ltd., Korea).

2.5. Effect of pH, Temperature, Metal Ions and Chemical Reagent on Enzyme Activity

The specific activity of purified BglPC28 was determined using p-nitrophenyl-β-D-glucopyranoside (pNPG) as a surrogate substrate in 50 mM sodium phosphate buffer, pH 7.0 at 37°C. Reactions were stopped after 10 minutes (min) by the addition of Na2CO3 at a final concentration of 0.5 M, and the release of p-nitrophenol was measured immediately using a microplate reader at 405 nm (Bio-Rad model 680; Bio-Rad, Hercules, CA). One unit of activity was defined as the amount of protein required to generate 1 µmol of p-nitrophenol per min. Specific activity was expressed as units per milligram of protein. Protein concentrations were determined using the bicinchoninic acid (BCA) protein assay (Pierce, Rockford, IL), with bovine serum albumin (Sigma) as the standard. All assays were performed in triplicate.

The effect of pH on enzymatic activity was determined using 1.0 Mm pNPG as a substrate in the following buffers (each at 50 mM): KCl-HCl (pH 2.0), glycine-HCl (pH 3.0), sodium acetate (pH 4.0 and 5.0), sodium phosphate (pH 6.0, 7.0 and 7.5), Tris-HCl (pH 8.0, and 9.0) and glycine-sodium hydroxide (pH 10). The pH stability of recombinant BglPC28 was determined by measuring enzymatic activity after incubation in each buffer (containing 2.0 mM pNPG in 50 mM potassium buffer as a substrate) for 12 h at 4°C. The results are expressed as a percentage of the activity obtained at the optimum pH. The effect of temperature on enzymatic activity was tested by incubating the enzyme at various temperatures ranging from 4 to 65°C (4, 10, 18, 25, 30, 37, 45, 55, 65°C) at optimum pH for 5 min in 50 mM potassium phosphate buffer containing 2.0 mM pNPG. The thermo-stability of the enzyme was examined by incubating the enzyme in 50 mM potassium phosphate buffer for 30 min at different temperatures. After cooling the sample on ice for 10 min, activity was determined using pNPG as the substrate.

The effects of metals and other chemicals on BglPC28 activity were also determined. BglPC28 activity was tested in the presence of 1 or 10 mM (final concentration) of HgCl2, MnCl2, CaCl2, CoCl2, MgCl2, EDTA, NaCl, KCl, CuCl2, SDS, dithiothreitol (DTT), or β-mercaptoethanol for 30 min at 37°C. The remaining activity was determined using pNPG as a substrate, and activities are expressed as a percentage of the activity obtained in the absence of the compound.

Substrate preference was examined using 2.0 mM chromogenic o-nitrophenyl (ONP) and p-nitrophenyl (PNP) as substrates at 37°C for 5 min, with one activity unit being defined as the release of 1 µmol o-nitrophenol or p-nitrophenol per min. The following substrates were tested: PNP-β-D-glucopyranoside, PNP-β-D-galactopyranoside, PNP-β-D-fucopyranoside, PNP-N-acetyl-β-D-glucosaminide, PNP-β-L-arabinopyranoside, PNP-β-D-mannopyranoside, PNP-β-D-xylopyranoside, PNP-α-D-glucopyranoside, PNP-α-L-arabinofuranoside, PNP-α-L-arabinopyranoside, PNP-α-L-rhamnopyranoside, PNP-α-D-mannopyranoside, PNP-α-D-xylopyranoside, ONP-β-D-glucopyranoside, ONP-β-D-galactopyranoside, ONP-β-D-fucopyranoside and ONP-α-D-galactopyranoside (all from Sigma).

2.6. Determination of Kinetic Parameters

Kinetic studies were performed with freshly purified enzymes using pNPG at 1–20 mM and Re at concentrations from 0.2 mM to 5.0 mM. One unit of activity was defined as the amount of protein required to generate 1 µmol of p-nitrophenol or to convert 1 µmol of Re per minute. All assays were performed in triplicate. The parameters, Km and Vmax, were determined using the enzyme kinetics program described by Cleland [43].

2.7. Biotransformation Activity of Ginsenosides Using BglPC28

The initial biotransformation experiments using the major ginsenosides Rb1 and Re as substrates revealed that the GST-fused enzyme does not affect the activities of BglPC28. Therefore, the fusion protein (GST-BglPC28) was used to determine the specificity and selectivity of the enzymes for the hydrolysis of the glucose moieties attached at the C3 and C20 positions in the five PPD ginsenosides. The enzyme solutions at a concentration of 0.1 mg/ml in 100 mM of sodium phosphate buffer (pH 7.0) were reacted with an equal volume of Rb1, Rb2, Rb3, Rc, Rd, Re, and Rg1 solution at a concentration of 0.2% (w/v) in 100 mM of sodium phosphate buffer (pH 7.0) at 37°C. The samples were taken at regular intervalsand analyzed via TLC or HPLC after pretreatment (see analytic methods).

2.8. Optimization of Concentration of the Enzyme and Substrate

In order to determine the optimal condition for the biotransformation of Re to Rg2(S), the substrate concentration of Re in the reaction was optimized. The final crude BglPC28 concentration was fixed to 10, 20, and 50 mg/ml and reacted with ginsenoside Redissolved in DMSO (100 mg/ml) in order to have 20 and 30 mg/ml as the final substrate concentration. These 6 types of optimization reactions were performed in a 2 ml eppendorf tube with a 1 ml working volume at 200 rpm for 24 h at 37°C. The samples were taken at regular intervals and analyzed via TLC and HPLC.

2.9. Preparation of the Recombinant Enzymes BglPC28 using High Cell Density Culture

For the production of the recombinant BglPC28, the LB medium supplemented with ampicillin (100 µg/ml) was used to cultivate the E. coli harboring pGEX-bglPC28 in a 10 L stirred-tank reactor (BiotronGX, Hanil science Co., Korea) with a 5 L working volume at 500 rpm. The pH value of the medium was adjusted to 7.0 using 100 mM of sodium phosphate buffer. The culture was incubated at 37°C until the culture reached an OD of 3.0 at 600 nm. The protein expression was induced through the addition of isopropyl-β-D-thiogalactopyranoside (IPTG) with a final concentration of 0.1 mM with feeding 2% (w/v) glucose. The bacterial cells were incubated for a further 18 h at 22°C and were then harvested via centrifugation at 5,000 rpm for 20 min (Component R, Hanil science Co Ltd., Korea) at 4°C.

The cells suspended in 100 mM of phosphate buffer (pH 7.0) were disrupted via sonication (Vibra-cell, Sonics & Materials, CT, USA), and then the intact cells and debris were removed via centrifugation at 5,000 rpm for 20 min (Component R, Hanil science Co Ltd., Korea) at 4°Cin order to obtain the supernatants of the crude enzymes. The crude recombinant BglPC28 was diluted to the desired concentration with 100 mM of sodium phosphate buffer (pH 7.0) and was used to biotransform the Re.

2.10. Scaled-up Biotransformation of Ginsenoside Re via Crude BglPC28

The scaled-up biotransformation was performed in a 10 L stirred-tank reactor (BiotronGX, Hanil science Co., Korea) with a 7.5 L working volumecontaining 10% DMSO (Dimethyl sulfoxide) at 100 rpm for 24 h. The reaction was performed under conditions in pH 7.0 at 30°C. The reaction started with a composition of 20 mg/ml of substrate ginsenoside (Re; total 150 g) as final concentration and 1.5 L of crude recombinant BglPC28 (20 mg/ml) in 100 mM of phosphate buffer (pH 7.0). Samples were collected at regular intervalsand were analyzed by HPLC in order to determine the production of the ginsenoside Rg2(S) from Re.

2.11. Purification of Rg2(S)

Following the 7.5 L reaction of Re with BglPC28, the mixture was cooled at 4°C and centrifuged at 5,000 rpm for 15 min (Component R, Hanil science Co Ltd., Korea). The biotransformed ginsenoside Rg2(S) in the supernatants and precipitates was processed separatelyin order to purify the samples. The precipitate was also dissolved in 5.0 L of 70% ethanol solution twice and filtered through a filter paper (Advantec, Japan). The ethanol extracts were combinedand adjusted to be a 40% ethanol solution. The column chromatography [3,170(L)×128(D) mm; Doointech, Korea] packed with HP20 resin (Mitsubishi, Japan) was adopted in order to remove the impurities, except the ginsenosides. The supernatants and 40% ethanol solution were loaded on to the column together. The free sugar molecules and unwanted hydrophilic compounds from the HP-20 that were adsorbed in beads were washed with 6 bed volumes (BV) of water, and finally the adsorbed ginsenosides were eluted using 6 BV of 95% ethanol (extra pure grade; SK Chemicals, Korea). The ethanol eluent was evaporated in vacuo. The resulting powder was dissolved in 100% methanol and analyzed via HPLC.

2.12. Analytic Methods

2.12.1. Thin Layer Chromatography (TLC) Analysis.

The thin layer chromatography (TLC) was performed using 60f254 silica gel plates (Merck, Germany) with Chcl3-CH3OH-H2O (65∶35∶10, lower phase) as the solvent. The spots on the TLC plates were identified through comparisons with standard ginsenoside after visualization was made by spraying 10% (vol/vol) H2SO4, followed by heating at 110°C for 5 min.

2.12.2. High Performance Liquid Chromatography (HPLC) analysis.

HPLC analysis of the ginsenosides [Rb1, Rd, Rg2(S), Rh2(S), F2, C-K, Rg1, Re and Rg2(S)] was performed using an HPLC system (Younglin Co., Ltd., Korea), with a quaternary pump, automatic injector, single-wavelength UV detector (model 730 D), and Younglin’s AutoChro 3000 software for peak identification and integration. The separation was carried out on a Prodigy ODS (2) C18 column (5 µm, 150×4.6-mm i.d.; Phenomenex, USA) with a guard column (5 µm, 12.5×4.6-mm i.d.; Eclipse XDB C18). The mobile phases were acetonitrile (A) and water (B). Gradient elution started with 17% solvent A and 83% solvent B changed to: A from 17 to 25%, 12–20 min; A from 25 to 32%, 20–30 min; A from 32 to 55%, 30–35 min; A from 55 to 60%, 35–40 min; A from 60 to 80%, 40–45 min; A from 80 to 100%, 45–50 min; A 100%, 50–54 min; A from 100 to 17%, 54.0–54.1 min; and A 17%, 54.1–65 min. The flow rate was 1.0 ml/min; detection was performed by monitoring the absorbance at 203 nm and an injection volume of 25 µl.

2.12.3. LC/MS/MS analysis.

Electrospray ionization mass spectrum (ESI-MS) was measured for the biotransformed Rg2(S) after purification steps on a triple-quadrupole tandemmass spectrometer (API-2000, Applied Biosystems, Foster City, USA) with negative ion mode. The ESI parameters were as follows: ionspray voltage, −4,200 V; ion source gas 1 (GS1), 20; curtain gas (CUR), 20; collision gas (CAD), 2. The declustering potential (DP), focusing potential (FP), entrance potential (EP), collision cell exit potential (CXP) and collision energy (CE) were variant with regard to measured ginsenosides. For full-scan MS analysis, the spectra were recorded in the m/z range from 400 to 1,000.

Results

3.1. Fosmid Library Construction and Cloning of BglPC28

β-Glucosidase activity was detected from a ginsenoside-hydrolyzing bacterium, Pseudonocardiasp. strain Gsoil 1536, isolated from the soil of a ginseng cultivating field. Since strain Gsoil 1536 converted ginsenosides Rb1 and Rd into C-K via F2 as well as converting Re and Rg1 into Rg2(S) and Rg1(S), respectively (data not shown), it could have two types of β-D-glucosidase activities toward β-1,2- and β-1,6-glucosidic linkages. To isolate and clone the respective genes for the ginsenoside-hydrolyzing enzymes from strain Gsoil 1536, we created a fosmid library, screened a clone containing the gene encoding the enzymatic activity, and determined the two contigs (14.3 - and 4.7 kb) partial sequence of the fosmid vector. An ORF FINDER (www.ncbi.nlm.nih.gov/gorf) analysis of the large contig of the 14.3 kb revealed that it encoded a total of 11 putative ORFs of longer than 300 codons (data not shown). Four putative ORFs of longer than 300 codons (data not shown) were revealed in the small size contig of 4.7 kb. Among three ORFs, one was homologousto glycoside hydrolase genes in glycoside hydrolase family 3(GH3). This ORF, termed bglPC28, consisting of 2,232 bp encoding 743 amino acids, was amplified via PCR and then inserted into the pGEX 4T-1 vector.

3.2. Phylogenetic Analysis of BglPC28 Sequences

Analysis of the amino acid sequences of BglPC28 indicated that it was 61.7% identical to the glycoside hydrolase (Bgp1) of Microbacterium esteraromaticum KACC 16318 (GenBank number AEX88466), which belongs to GH3. The enzymatic activity of Bgp1 has been characterized to have ginsenoside conversion activity [35]. BglPC28 has homology to the protein domain of GH3. The Carbohydrate-Active enZymes database (http://www.cazy.org) describes more than 5,000 uncharacterized and 222 characterized GH3 members that are widespread across numerous organisms. Glycosyl hydrolases are classified according to their amino acid sequence similarities, which reflect the structural features and substrate specificities of the enzymes (http://www.cazy.org/fam/acc_GH.html). The GH3 is subdivided into six subfamilies [44]. In order to determine the evolutionary position of BglPC28 within the characterized enzymes in glycoside hydrolasesfamily 3, a phylogenetic analysis was conducted using the neighbor-joining method in the MEGA5 Program with bootstrap values based on 1,000 replications. The resulting consensus tree is presented in Fig. 2. BglPC28 clustered within subfamily 5 and formed a separate, well-supported clade (bootstrap of 100) with β-glucosidase (Bgp1) derived from Microbacterium esteraromaticum KACC 16318, which was characterized. Several ginsenoside-hydrolyzing β-glucosidases in GH3 have previously been cloned, including a β-glucosidase (rApy-H11) from Bifidobacterium longum H-1 [45], a β-glucosidase (BgpA) from Terrabacterginsenosidimutans [46], a β-glucosidase (BGL1) from Aspergillusniger [47], a β-glucosidase (BglAm) from Actinosynnemamirum [28], a β-glucosidase (BglQM) from Mucilaginibacter sp. QM49 [29], a β-glucosidase (BglSk) from Sanguibacter keddieii [31], and a β-glucosidase (BglBX10) from Flavobacterium johnsoniaeUW101T [48]. The relationship between BglPC28 and these ginsenoside-hydrolyzing β-glucosidases is presented in Fig. 2.

[Figure omitted. See PDF.]

Figure 2. Phylogenetic analysis of characterized glycoside hydrolases family 3 (GH3).

Amino acid sequences were obtained from the NCBI/EMBL database and CAZy database (accession numbers are indicated on the tree). This tree was made using the neighbor-joining method with a poisson model and pairwise deletion. Bootstrap values expressed as percentages of 1,000 replications greater than 65% are shown at the branch points. The bar represents 20 amino acid residues substitutions per 100 amino acid residues.

https://doi.org/10.1371/journal.pone.0096914.g002

Comparison of the BglPC28 sequence with ginsenoside-transforming enzyme (Bgp1) from Microbacterium esteraromaticum [35] and known-structure β-glucosidases from Streptomyces venezuelae [49], Thermotoganeapolitana [50] and Kluyveromycesmarxianus [51] are shows in Fig. 3. The amino acids - D242 and E408 of BglPC28 which were marked by asterisk were thought to serve as nucleophile and the acid–base during the hydrolysis reaction through sequence alignment with the structure-determined above mentioned three β-glucosidases. The PA14 domain which is existed in many proteins such as glycosidases, glycosyl-transferases, proteases, amidases, toxins, adhesins and signaling molecules were found in some GH3 enzymes [52]. PA14 domain was found 230 of analyzed 820 sequences of GH3, and may have a carbohydrate-binding or fundamental role in biological events [51]. However, no PA14 domain was found in BglPC28 and Bgp1.

[Figure omitted. See PDF.]

Figure 3. Sequence alignment of ginsenoside-transforming BglPC28 and Bgp1 or structure-determined glycoside hydrolase family 3 enzymes.

The sequence alignment was created using ClustalW at the EBI-server using default settings and visualized with Jalview (http://www.ebi.ac.uk/Vmichele/jalview/). Regions of identity or high similarity among sequences are shown as black or gray columns, respectively. The conserved active site residues (general acid/base and nucleophile residues) were marked by asterisk. The PA14 domains of S.ven and K.mar are boxed. Genbank IDs of the glycoside hydrolases family 3 are as follows, Pseudonocardia sp. Gsoil 1536 β-glucosidase [BglPC28 (This study)], JX960416; Microbacterium esteraromaticum KACC 16318 β-glucosidase (Bgp1), AEX88466; Streptomyces venezuelae β-glycosidase (S.ven), AAC68679; ThermotoganeapolitanaDSM 4359T β-glycosidase (T.nea), ABI29899; Kluyveromycesmarxianus β-glucosidase (K.mar), ACY95404.

https://doi.org/10.1371/journal.pone.0096914.g003

3.3. Expression and Purification of Recombinant BglPC28

In order to maximize the yield of the fusion protein in a soluble form, different induction conditions were tested and it was found that induction with 0.15 mM IPTG at 22°C for 18 h cultivation after initial induction produced the maximum level of soluble active fusion enzyme (data not shown). The GST-BglPC28 fusion protein was purified using the GST•bind agarose resin, and the GST tag was removed using thrombin at room temperature for 12 h incubation. The recombinant enzyme was purified by GST•bind agarose resin, and then supernatant from cell lysates as well as purified protein was applied to SDS-PAGE (Fig. 4). The molecular mass of the GST-BglPC28 calculated via an amino acid sequence was 105 KDa, which is similar mass detected in SDS-PAGE. In addition, the recombinant GST•BglPC28 contains 36.5±1.4% of total soluble protein in E. coli lysate. This high expression level in soluble form makes it more feasible for industrial application.

[Figure omitted. See PDF.]

Figure 4. SDS-PAGE analysis of the recombinant BglPC28 after purification using the GST•bind agarose resin.

Lane: 1, uninduced crude extract; 2, soluble fraction of the crude extract of the induced recombinant BL21 (DE3) cells; 3, precipitated fraction of the crude extract of the induced recombinant BL21 (DE3) cells; 4, GST-BglPC28 after purification with the GST•bind agarose resin; 5, purified recombinant BglPC28 after cleavage via thrombin.

https://doi.org/10.1371/journal.pone.0096914.g004

3.4. Characterization of Recombinant BglPm

The optimum temperature and pH were determined using purified GST-BglPC28. The optimal temperature activity of BglPC28 was 37°C; at 30°C and 25°C, the enzyme showed 90.0% and 67.2% relative activity, respectively, while thermostability decreased quickly beyond 37°C (Fig. 5a). The enzyme was stable at temperatures lower than 30°C, and about 70.6% of the activity was lost after incubation at 37°C for 30 min (Fig. 5b). BglPC28 had optimal pH activity and stability at pH 7.0 in a sodium phosphate buffer at 30°C; from pH 8.0, the enzyme stability decreased swiftly, while at pH 5.0 the enzyme activity decreased to 58.1% (Fig. 5b). The enzyme is probably mesophilic and stable at a neutral pH range. These optimal conditions are consistent with the soil environment from which Pseudonocardia sp. Gsoil 1536 was isolated. In addition, the near-neutral optimal pH and mild optimal temperature of BglPC28 are similar to those of other ginsenoside-hydrolyzing GH3 from bacteria [28], [29], [35], [46], [47]. Although the optimum temperature of BglPC28 for pNPG is 37°C, the ginsenoside-conversion reaction occurred at 30°C for extension of stable transformation activity.

[Figure omitted. See PDF.]

Figure 5. Effects of temperature (a) and pH (b) on the activity and stability of recombinant BglPC28.

https://doi.org/10.1371/journal.pone.0096914.g005

The effects of metal ions, EDTA, β-mercaptoethanol, and SDS, on BglPC28 activity were also investigated (Table 1). The enzyme did not require Mg2+ for activity. Furthermore, it lost its activity in the presence of 50 mM Cu2+, and was significantly inhibited by 1 mM Hg2+ or 10 mM SDS. BglPC28 activity was not affected by DTT and β-mercaptoethanol, which are well-known thiol group inhibitors. These results suggest that sulfhydryl groups may not be involved in the catalytic center of the enzyme. The chelating agent EDTA did not inhibit BglPC28 activity, which indicated that divalent cations are not required for enzymatic activity. Thus, no dramatic positive effects on the activity of the BglPC28 were found for the tested ions (Table 1).

[Figure omitted. See PDF.]

Table 1. Effects of metal ions and reagents on the activity of purified recombinant BglPC28.

https://doi.org/10.1371/journal.pone.0096914.t001

The substrate specificity of BglPC28 was tested using 2.0 mM pNP- and oNP-glycosides with α and β configurations. The results, which are summarized in Table 2, showed that BglPC28 was maximally active against pNP-β-D-glucopyranoside, followed by pNP-β-L-arabinopyranoside, but had little effect on various other PNP- and ONP-glycosides, with the exception of pNP-α-D-glucopyranoside, pNP-β-D-glucosaminide, and pNP-β-D-xylopyranoside. The specific activity of BglPC28 against pNP-β-D-xylopyranoside is not seen in ginsenoside-transforming GH3 members from Terrabacter ginsenosidimutans [46], Microbacterium esteraromaticum [35], and Mucilaginibacter sp. QM49 [29].

[Figure omitted. See PDF.]

Table 2. Relative activity of purified recombinant BglPC28 towards various chromogenic substrates as measured by PNP or PNP release at 37°C. ND: not determined.

https://doi.org/10.1371/journal.pone.0096914.t002

3.5. Determination of Kinetic Parameters

The kinetic parameters of Vmax and Km of BglPC28 were determined by plotting the substrate concentration vs. the initial velocity of each reaction and subjecting the data to a linear regression analysis. The Km, Vmax, kcat, and kcat/Km values for 4 substrates are presented in Table 3. The Km values of BglPC28 for pNPG and Re were 6.36±1.10 and 1.42±0.13 mM, respectively, and the Vmax values were 40.0±2.55 and 5.62±0.21 µmol min−1 mg−1 of protein, respectively. The kcat values were 52.7±3.4, and 7.40±0.28 S−1, respectively measured using the method described by Cleland [43].

[Figure omitted. See PDF.]

Table 3. Kinetic parameters of recombinant BglPC28 on pNPG and ginsenoside Re.

https://doi.org/10.1371/journal.pone.0096914.t003

The catalytic efficiencies (kcat/Km) for pNPG and Re decreased in this order: pNPG (8.63±2.02 mM−1 S−1) >Re (5.28±0.69 mM−1 S−1). The catalytic efficiencies for ginsenosides are higher than those of β-glycosidase from Sulfolobus acidocaldarius [53] for ginsenosides Rd (4.8 mM−1 min−1) and Rb1 (4.8 mM−1 min−1), and β-glucosidase from Mucilaginibacter sp. QM49 [29] for ginsenosides Rg1 (1.64±0.22 mM) and Re (0.51±0.07 mM).

3.6. Biotransformation Characteristics of BglPC28

For verification of the bioconversion pathway of the five PPD type ginsenosides (Rb1, Rb2, Rb3, Rc, and Rd) and two PPT type ginsenosides (Re and Rg1) by GST-BglPC28, TLC and HPLC analyses were performed at regular intervals. GST-BglPC28 could clearly transform four major ginsenosides (Re, Rg1, Rb1, and Rd), as shown by the Rf values of the TLC analysis (Fig. 6) and the retention times in HPLC (data not shown). BglPC28 could efficiently hydrolyze the glucose moiety at the C20 position of Re and Rg1, transforming them into Rg2(S) and Rh1(S), respectively (Fig. 7a, 7b). In addition to these PPT-type ginsenosides, BglPC28 could also transform the PPD-type ginsenosides, Rb1, Rb3, and Rd. The proposed biotransformation pathways for the actions of BglPC28 on the PPD-type ginsenosides are Rb1 or Rd → Rg3(S) via hydrolysis of the outer glucoses directly at the C20 position. Furthermore, BglPC28 transformed Rb3 to Rd weakly, by hydrolyzing the D-xyloside moiety attached at the C20 site, as indicated by pNP-β-D-xylopyranoside activity in BglPC28 (Table 2). Overall, BglPC28 hydrolyzed Re, Rg1, Rb1, Rb3, and Rd as follows: Re → Rg2(S), Rg1 → Rh1(S), Rb1 → Rg3(S), Rd → Rg3(S), and Rb3 → Rg3(S) (Fig. 8).

[Figure omitted. See PDF.]

Figure 6. Thin layer chromatography (TLC) analyses of biotransformation of Rb1, Rb2, Rb3, Rc and Rd by recombinant BglPC28.

The sampling times were 10: CHCl3-CH3OH-H2O (65∶35∶10, v/v, lower phase). S, ginsenoside standards (PPD type ginsenoside mixtures).

https://doi.org/10.1371/journal.pone.0096914.g006

[Figure omitted. See PDF.]

Figure 7. TLC analyses of time course of ginsenosides bioconversion by BglPC28.

(a), transformation of ginsenoside Re; (b), transformation of ginsenoside Rg1. Developing solvent: CHCl3-CH3OH-H2O (65∶35∶10, v/v, lower phase). Lanes S, ginsenosidestandards (PPT type ginsenoside mixtures).

https://doi.org/10.1371/journal.pone.0096914.g007

[Figure omitted. See PDF.]

Figure 8. Transformation pathways of ginsenosides Rb1, Rb3, Rd, Re, and Rg1 by recombinant BglPC28, respectively to become biotransformed ginsenosides.

https://doi.org/10.1371/journal.pone.0096914.g008

The metabolic pathways of ginsenoside hydrolysis by BglPC28 are very similar to those of a β-glucosidase (Bgp1) from Microbacterium esteraromaticum [33], [34], and this is attributed to their high amino acid sequence similarity (Fig. 2). BglPC28 also has similar metabolic pathways of ginsenosides hydrolysis with a β-glucosidase (BglQM) from Mucilaginibacter sp. strain QM49 except that BglQM hydrolyzed Rb2 and Rc into C-Y and C-Mc, respectively, whereas BglPC28 did not hydrolyze them. In addition, the amino acid sequence similarity of both is very low and located in different subgroups (Fig. 2). Re and Rg1 were deglucosylated when exposed to a culture of Pseudonocardia sp. Gsoil 1536, and recombinant BglPC28 has the same effect, presumably via the hydrolyzing pathways. This suggests that hydrolyzing activity for ginsenosides Re and Rg1 of Pseudonocardia sp. Gsoil 1536 arises from the expression of BglPC28.

When Re, Rg1, Rb1, and Rd (1.0 mg/ml) were used as substrates, they were biotransformed within 10 min by 10 mg/ml of crude recombinant GST-BglPC28. Several other ginsenoside-hydrolyzing recombinant enzymes have been reported, but most were only able to hydrolyze the outer or inner glucose moiety at the C3 position of the aglycon [28], [31], [35], [36], [46], [53], [54], [55]. Three previously described enzymes, Bgp1 [35], BglAm [28] and BglQM [29], had the same metabolic pathways of ginsenoside hydrolysis for ginsenosides Re and Rg1 as BglPC28, but the authors conducted only a simple enzymatic characterization without further scale-up or process engineering and in some case low yield and a complicated purification step impeded practical application of these enzymes. Thus, we sought to optimize our BglPC28 reaction for industrially relevant 100 gram-scale production of Rg2(S).

3.7. Optimization of Re and Enzymes Concentration

Two substrate concentrations (20 mg/ml and 30 mg/ml) and three crude enzyme concentrations (10, 20, and 50 mg/ml) as final concentrations were tested in order to determine the appropriate substrate concentration for decreasing the reactor volume and economical enzyme concentration to reduce production costs. The time course of the ginsenoside Re and product Rg2(S) was determined via HPLC analyses in six test conditions (Fig. 9). In the test conditions of a low substrate concentration (20 mg/ml) and a high crude enzyme concentration (50 mg/ml), the ginsenoside Re was completely converted to ginsenoside Rg2(S) within 4 hours; thisreaction speed is three times faster than that of the 20 mg/ml Re with 20 mg/ml crude enzyme concentration and six times faster than that of 30 mg/ml Re with 50 mg/ml crude enzyme concentration. Under the other three reaction conditions the conversion was not completed within 40 hours. Thus, these three reaction conditions were excluded in the next step. Providing the advantages of smaller usage of enzyme and complete conversion of Re, the conditions of 20 mg/ml substrate concentration and 20 mg/ml crude enzyme concentration were adopted for the next scaled-up biotransformation step.

[Figure omitted. See PDF.]

Figure 9. Effect of the concentration of Re and crude BglPC28 on the production of ginsenoside Rg2(S).

https://doi.org/10.1371/journal.pone.0096914.g009

3.8. Preparation of Crude pGEX- BglPC28 and Scaled-up Production of Ginsenosides Rg2(S)

The bacterial cells that harbor pGEX-bglPC28 and were incubated further for 18 h at 22°C after induction were harvested via centrifugation when the culture reached an OD of 38 at 600 nm. 154 g of wet cells was harvested and the pellets were resuspended in 10 volumes (w/v) of 100 mM sodium phosphate buffer (pH 7.0). The cells were disrupted via ultrasonication and the supernatant was used as crude enzymes for the biotransformation of the ginsenosides. As almost half of the expressed pGEX•BglPC28 was in soluble form, the crude recombinant BglPC28 was applied to the biotransformation reactor. The enzyme reaction occurred using the crude recombinant BglPC28 with RG-Re as the substrate with a concentration of 20 mg/ml, and crude recombinant BglPC28 was adjusted to 20 mg/ml as the final concentration in 7.5 L in order to produce Rg2(S). The ginsenoside Re was completely converted to Rg2(S) within 24 hours after the crude GST-BglPC28 was applied to the ginsenoside Re (Fig. 10). The research team behind this paper searched the ginsenoside hydrolyzing bacteria [27], [56], [57] and constructed several ginsenoside-hydrolyzing recombinant enzymes [28], [29], [30], [31], [46], [48], [58], [59] by surveying the ginsenoside-hydrolyzing enzymes. However, these enzymes are not applicable or are not active in transforming the Re. At last, the finding of exceptional BglPC28 that can efficiently convert Re is a key factor in creating Rg2(S).

[Figure omitted. See PDF.]

Figure 10. HPLC analysis of the transformation of Re by GST-BglPC28; S, standard ginsenosides; a, Re as substrate for production for Rg2(S); b, produced Rg2(S) after treating Re and GST-BglPC28 for 24 hours.

https://doi.org/10.1371/journal.pone.0096914.g010

3.9. Purification of Biotransformed Rg2(S)

In order to remove the enzymes, salt, and free sugars from the reaction mixture of the 7.5 L reaction of Re with BglPC28, the mixture was centrifuged at 5,000 rpm for 15 min. Most of the ginsenoside Rg2(S) was precipitated to form a solid, with a small quantity remaining dissolved in the supernatant (data not shown). After a purification step using a column chromatograph packed with HP20 resin, approximately 24 L of the 95% ethanol eluent was evaporated in vacuo in order to create 113 g of ginsenoside Rg2(S). Its chromatographic purity was 84.0±1.1%, as determined via HPLC. Furthermore, electrospray negative ion mass spectrum of produced ginsenoside Rg2(S) was obtained on API 2000 to confirm its identity. Produced Rg2(S) was a white powder with a pseudomolecularion peak [M−H]− at m/z 783.8 in ESIMS, corresponding to C42H72O13 (calculated molecular weight, 784.5). Thus, produced Rg2(S) was confirmed (Fig. 11). Among Re, the total molar amount of Re that could be biotransformed into Rg2(S) using BglPC28 was 138.8 mmol, which corresponds to 131.4 g of 150 g. The residue (18.6 g) was composed of other types of ginsenosides, moisture and unknown impurities. The molar amount of the produced ginsenoside Rg2(S) was 120.9 mmol. This indicates that the recovery ratio through the biotransformation process using ginsenoside Re to Rg2(S) reached 87.1% after bioprocess engineering.

[Figure omitted. See PDF.]

Figure 11. Electrospray negative ion mass spectrum of produced ginsenoside Rg2(S) of which m/z was 783.8 [M−H]−. Glc, glucose moiety; Rha, rhamnose moiety.

https://doi.org/10.1371/journal.pone.0096914.g011

Discussion

Although Rg2(S) has neuroprotective, anti-inflammation, and anti-diabetic effects, the lack of a selective mass-production technology has hampered its commercial use. Red ginseng, which is a very popular health-promoting food in oriental herbal medicine, contains approximately less than 0.02% Rg2(S) based on dry weight. In total, red ginseng products account for 52.8% of the total health food market in Korea, which was estimated to be worth 625 million US dollars in 2011 [60]. Furthermore, one Korean cosmetic company (Somang Co. Ltd.) recently launched a brand name RGII product that contains the ginsenoside Rg2 as the major active component. However, the high cost of purified Rg2 from red ginseng obstructs its widespread use.

To achieve Rg2(S), a number of researchers have sought to biotransform major ginsenosides into larger quantities of Rg2(S) or Rh1(S) using microorganisms [61], [62] and enzymes [29], [35], [63], [64] in laboratory settings. Ko et al [63] characterized the abilities of crude β-galactosidase from Aspergillus oryzae and crude lactase from Penicillium sp. to biotransform PPT-type ginsenosides to Rg2(S) and Rh1(S), but only on a 10-mg scale. The usage of cloned ginsenoside-hydrolyzing GHs in glycoside hydrolase families 1, 3, and 51 are the other method. As shown in Table 4, family 3 GHs has various types of ginsenoside-hydrolyzing pathways. As such, they have strong potential for the production of various types of minor ginsenosides from major ginsenosides for practical application. Recently, three recombinant β-glucosidases belonging to glycoside hydrolase family 3 were shown to transform Re and Rg1 into Rg2(S) and Rh1(S), respectively. BgpA from Microbacterium esteraromaticum and BglAm from Actinosynnemamirum, but the level of research remains only at the level of simple enzyme characterization, and sufficient efficiency for gram-scale applications has not been obtained [28], [35]. The other (BglQM), from Mucilaginibacter sp. QM49, was applied to obtain a gram scale Rg2(S) and Rh1(S) using a PPT type ginsenoside mixture. However, the yield was low (33.4%) and the final purification step using silica resin was complicated. As a family 3GH, BglPC28 exhibited a similar ginsenoside-transforming pathway to those of BgpA and BglQM, but showed high conversion activity using Re, which was not easily dissolved in water. Re was readily dissolved in DMSO up to 200 mg/ml. Here, we report for the first time that BglPC28 can transform up to 20 mg/ml of Re into 100 gram-scale Rg2(S) within 24 h. Our novel identification of BglPC28 as an enzyme capable of converting Re and Rg1 into Rg2(S) and Rh1(S) is expected to facilitate mass production of Rg2(S) from protopanaxatriol-type ginsenoside Re derived from Panax quinquefolius (American ginseng) or Panax ginseng C. A. Meyer (Korean ginseng).

[Figure omitted. See PDF.]

Table 4. Major ginsenosides transformations by the cloned glycoside hydrolases family 3.

https://doi.org/10.1371/journal.pone.0096914.t004

In summary, we herein describe the isolation of a ginsenoside-hydrolyzing β-glucosidase (BglPC28) belonging to glycoside hydrolase family 3 from Pseudonocardia sp. Gsoil 1536, and the production of a recombinant enzyme for the biotransformation of the major ginsenoside Re into the pharmacologically active rare ginsenoside Rg2(S). This enzyme was expressed in E. coli BL21(DE3) in a soluble form. Characterization revealed that its optimum reaction conditions were 37°C and pH 7.0, and that it could be used for 100 gram-scale Rg2(S) production. In terms of yield, 113 g of Rg2(S) with 84.0±1.1% chromatographic purity was obtained via the biotransformation of 150 g of a purified ginsenoside Re followed by purification with HP20 resin. Bioconversion process was taken place in a 10 L jar fermenter in 100 mM sodium phosphate buffer (pH 7.0) containing 10% DMSO at 30°C for 24 h, with an initial substrate concentration of 20 mg/ml. This is the first report of a cloned enzyme that is capable 100 gram-scale production of Rg2(S) through biotransformation of the major ginsenoside Re.

Acknowledgments

This work was supported by the Intelligent Synthetic Biology Center of Global Frontier Project funded by the Ministry of Education, Science and Technology (2011-0031967), Republic of Korea and by the China Major Projects No.2012ZX09503001-003.

Author Contributions

Conceived and designed the experiments: WI. Performed the experiments: JD CC SP JK WI. Analyzed the data: JD CC HY WI. Contributed reagents/materials/analysis tools: FJ CS SK. Wrote the paper: JD CC WI.

Citation: Du J, Cui C-H, Park SC, Kim J-K, Yu H-S, Jin F-X, et al. (2014) Identification and Characterization of a Ginsenoside-Transforming β-Glucosidase from Pseudonocardia sp. Gsoil 1536 and Its Application for Enhanced Production of Minor Ginsenoside Rg2(S). PLoS ONE 9(6): e96914. https://doi.org/10.1371/journal.pone.0096914

References

1. Cho IH (2012) Effects of Panax ginseng in neurodegenerative diseases. J Ginseng Res 36: 342–353.

2. Park HJ, Kim DH, Park SJ, Kim JM, Ryu JH (2012) Ginseng in traditional herbal prescriptions. J Ginseng Res 36: 225–241.

3. Attele AS, Wu JA, Yuan CS (1999) Ginseng pharmacology: multiple constituents and multiple actions. BiochemPharmacol 58: 1685–1693.

4. Buettner C, Yeh GY, Phillips RS, Mittleman MA, Kaptchuk TJ (2006) Systematic review of the effects of ginseng on cardiovascular risk factors. Ann PharmacoTher 40: 83–95.

5. Kang S, Min H (2012) Ginseng, the ‘Immunity Boost’: the effects of Panax ginseng on immune system. J Ginseng Res 36: 354–368.

6. Kim SK, Park JH (2011) Trends in Ginseng Research in 2010. J Ginseng Res 35: 389–398.

7. Hofseth LJ, Wargovich MJ (2007) Inflammation, cancer, and targets of ginseng. J Nutr 137: 183–185.

8. Lee F.C. (1992). Facts about Ginseng, the Elixir of Life. Hollyn International Corp., Elizabeth, NJ.

9. Yuan HD, SJ Kim, Chung SH (2011) Beneficial effects of IH-901 on glucose and lipid metabolisms via activating adenosine monophosphate-activated protein kinase and phosphatidylinositol-3 kinase pathways. Metabolism 60: 43–51.

10. Lee JH, Ahn JY, Shin TJ, Choi SH, Lee BH, et al. (2011) Effects of minor ginsenosides, ginsenoside metabolites, and ginsenoside epimers on the growth of Caenorhabditiselegans. J Ginseng Res 35: 375–383.

11. Leung KW, Wong AS (2010) Pharmacology of ginsenosides: a literature review. Chin Med 5: 20–22.

12. Christensen LP (2009) Ginsenosides chemistry, biosynthesis, analysis, and potential health effects. Adv Food Nutr Res 55: 1–99.

13. Choi S, Kim TW, Singh SV (2009) Ginsenoside Rh2-mediated G1 phase cell cycle arrest in human breast cancer cells is caused by p15 Ink4B and p27 Kip1-dependent inhibition of cyclin-dependent kinases. Pharm Res 26: 2280–2288.

14. Choi SH, Shin TJ, Hwang SH, Lee BH, Kang J, et al. (2011) Differential effects of ginsenoside metabolites on HERG k channel currents. J Ginseng Res 35: 191–199.

15. Kim YJ, Yamabe N, Choi P, Lee JW, Ham J, et al. (2013) Efficient thermal deglycosylation of ginsenoside Rd and its contribution to the improved anticancer activity of ginseng. J Agric Food Chem 61: 9185–9191.

16. In G, Ahn NG, Bae BS, Han ST, Noh KB, et al. (2012) Newmethod for simultaneous quantification of 12 ginsenosides in red ginseng powder and extract: in-house method validation. J Ginseng Res 36: 205–210.

17. Tian J, Zhang S, Li G, Liu Z, Xu B (2009) 20(S)-ginsenoside Rg2, a neuroprotective agent, inhibits mitochondrial permeability transition pores in rat brain. Phytotherpy Research 23: 486–491.

18. Li N, Liu B, Dluzen DE, Jin Y (2007) Protective effects of ginsenoside Rg2 against glutamate-induced neurotoxicity in PC12 cells. J Ethnopharmacol 111: 458–463.

19. Zhang G, Liu A, Zhou Y, San X, Jin T, et al. (2008) Panax ginseng ginsenoside-Rg2 protects memory impairment via anti-apoptosis in a rat model with vascular dementia. J Ethnopharmacol 115: 441–448.

20. Jeong SJ, Han SH, Kim DY, Lee JC, Kim HS, et al. (2007) Effects of mRg2, a mixture of ginsenosides containing 60% Rg2, on the ultraviolet B-induced DNA repair synthesis and apoptosis in NIH3T3 cells. Int J Toxicol 26: 151–158.

21. Ha SE, Shin DH, Kim HD, Shim SM, Kim HS, et al. (2010) Effects of ginsenoside Rg2 on the ultraviolet B-induced DNA damage responses in HaCaT cells. N-S Arch Pharmacol 382: 89–101.

22. Lu J, Lu WW, Tian JM, Li LY (2001) Effect of Rg2 on hemodynamics of hemorrhagic shock and its anti-oxidation in dogs. ZhongguoZhong Yao ZaZhi 26: 556–558.

23. Cho YS, Kim CH, Ha TS, Lee SJ, Ahn HY (2013) Ginsenoside Rg2 inhibits lipopolysaccharide-induced adhesion molecule expression in human umbilical vein endothelial cell. Korean J PhysiolPharmacol 17: 133–137.

24. Yuan HD, Kim do Y, Quan HY, Kim SJ, Jung MS, et al. (2012) Ginsenoside Rg2 induces orphan nuclear receptor SHP gene expression and inactivates GSK3beta via AMP-activated protein kinase to inhibit hepatic glucose production in HepG2 cells. ChemBiol Interact 195: 35–42.

25. Kim SN, Ha YW, Shin H, Son SH, Wu SJ, et al. (2007) Simultaneous quantification of 14 ginsenosides in Panax ginseng C.A. Meyer (Korean red ginseng) by HPLC-ELSD and its application to quality control. J Pharm Biomed Anal 45: 164–170.

26. Sun LQ (2004) Information on research and application of ginseng, the king of traditional and herbal medicines. Asian J Drug MetabPharma 4: 261–284.

27. Cui CH, Choi TE, Yu H, Jin F, Lee ST, et al. (2011) Mucilaginibacter composti sp. nov, with ginsenoside converting activity, isolated from compost. J Microbiol 49: 393–398.

28. Cui CH, Kim SC, Im WT (2013) Characterization of the ginsenoside-transforming recombinant β-glucosidase from Actinosynnema mirum and bioconversion of major ginsenosides into minor ginsenosides. ApplMicrobiolBiotech 97: 649–659.

29. Cui CH, Liu QM, Kim JK, Sung BH, Kim SG, et al. (2013) Identification and characterization of Mucilaginibacter sp. QM49 β-glucosidase and its use in producing the pharmaceutically active minor ginsenosides, Rh1(S) and Rg2(S). Appl Environ Microbiol 79: 5788–5798.

30. Hong H, Cui CH, Kim JK, Jin FX, Kim SC, et al. (2012) Enzymatic biotransformation of ginsenoside Rb1 and gypenoside XVII into ginsenosides Rd and F2 by recombinant β-glucosidase from Flavobacterium johnsoniae. J Ginseng Res 36: 418–424.

31. Kim JK, Cui CH, Yoon MH, Kim SC, Im WT (2012) Bioconversion of major ginsenosides Rg1 to minor ginsenoside F1 using novel recombinant ginsenoside hydrolyzing glycosidase cloned from Sanguibacter keddieii and enzyme characterization. J Biotechnol 161: 294–301.

32. Park CS, Yoo MH, Noh KH, Oh DK (2010) Biotransformation of ginsenosides by hydrolyzing the sugar moieties of ginsenosides using microbial glycosidases. Appl Microbiol Biotechnol 87: 9–19.

33. Quan LH, Jin Y, Wang C, Min JW, Kim YJ, Yang DC (2012) Enzymatic transformation of the major ginsenoside Rb2 to minor compound Y and compound K by a ginsenoside-hydrolyzing β-glycosidase from Microbacterium esteraromaticum. J IndMicrobiolBiotech 39: 1557–1562.

34. Quan LH, Min JW, Jin Y, Wang C, Kim YJ, Yang DC (2012) Enzymatic biotransformation of ginsenoside Rb1 to compound K by recombinant β-glucosidase from Microbacterium esteraromaticum. J Agr Food Chem 60: 3776–3781.

35. Quan LH, Min JW, Sathiyamoorthy S, Yang DU, Kim YJ, Yang DC (2012) Biotransformation of ginsenosides Re and Rg1 into ginsenosides Rg2 and Rh1 by recombinant β-glucosidase. Biotechnollett 34: 913–917.

36. Quan LH, Min JW, Yang DU, Kim YJ, Yang DC (2012) Enzymatic biotransformation of ginsenoside Rb1 to 20(S)-Rg3 by recombinant β-glucosidase from Microbacterium esteraromaticum. ApplMicrobiolBiotech 94: 377–384.

37. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25: 4876–4882.

38. Hall TA (1999) BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids SympSer 41: 95–98.

39. Kimura M (1983) The Neutral Theory of Molecular Evolution Cambridge: Cambridge University Press, Cambridge, New York.

40. Saitou N, Nei M (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. MolBiolEvol 4: 406–425.

41. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. MolBiolEvol 28: 2731–2739.

42. Felsenstein J (1985) Confidence limits on phylogenies: An approach using the bootstrap Evolution. 39: 783–791.

43. Cleland WW (1979) Statistical analysis of enzyme kinetic data. Methods Enzymol 63: 103–138.

44. Harvey AJ, Hrmova M, De G, N Varghese J, Fincher GB (2000) Comparative modeling of the three-dimensional structures of family 3 glycoside hydrolases. Proteins 41: 257–269.

45. Lee JH, Hyun YJ, Kim DH (2011) Cloning and characterization of α-L-arabinofuranosidase and bifunctional α-L-arabinopyranosidase/β-D-galactopyranosidase from Bifidobacterium longum H-1. J ApplMicrobiol 111: 1097–1107.

46. An DS, Cui CH, Lee HG, Wang L, Kim SC, et al. (2010) Identification and characterization of a novel Terrabacterginsenosidimutans sp. nov. β-glucosidase that transforms ginsenoside Rb1 into the rare gypenosides XVII and LXXV. Appl Environ Microbiol 76: 5827–5836.

47. Ruan CC, Zhang H, Zhang LX, Liu Z, Sun GZ, et al. (2009) Biotransformation of ginsenoside Rf to Rh1 by recombinant beta-glucosidase. Molecules 14: 2043–2048.

48. Kim JK, Cui CH, Liu QM, Yoon MH, Kim SC, et al. (2013) Mass production of the ginsenoside Rg3(S) through the combinative use of two glycoside hydrolases. Food Chem 141: 1369–1377.

49. Zmudka MW, Thoden JB, Holden HM (2013) The structure of DesR from Streptomyces venezuelae, a β-glucosidase involved in macrolide activation. Protein Sci 22: 883–892.

50. Pozzo T, Pasten JL, Karlsson EN, Logan DT (2010) Structural and functional analyses of β-glucosidase 3B from Thermotoga neapolitana: a thermostable three-domain representative of glycoside hydrolase 3. J MolBiol 397: 724–739.

51. Yoshida E, Hidaka M, Fushinobu S, Koyanagi T, Minami H, et al. (2010) Role of a PA14 domain in determining substrate specificity of a glycoside hydrolase family 3 β-glucosidase from Kluyveromycesmarxianus. Biochem J 431: 39–49.

52. Rigden DJ, Mello LV, Galperin MY (2004) The PA14 domain, a conserved all-βdomain in bacterial toxins, enzymes, adhesins and signaling molecules. Trends BiochemSci 29: 335–339.

53. Noh KH, Son JW, Kim HJ, Oh DK (2009) Ginsenoside compound K production from ginseng root extract by a thermostable β-glycosidase from Sulfolobus solfataricus. BiosciBiotechnolBiochem 73: 316–321.

54. Noh KH, Oh DK (2009) Production of the rare ginsenosides compound K, compound Y, and compound Mc by a thermostableβ-glycosidase from Sulfolobus acidocaldarius. Biol Pharm Bull 32: 1830–1835.

55. Lee GW, Kim KR, Oh DK (2012) Production of rare ginsenosides (compound Mc, compound Y and aglyconprotopanaxadiol) by β-glucosidase from Dictyoglomusturgidum that hydrolyzes beta-linked, but not alpha-linked, sugars in ginsenosides. BiotechnolLett 34: 1679–1686.

56. An DS, Wang L, Kim MS, Bae HM, Lee ST, Im WT (2011) Solirubrobacterginsenosidimutans sp. nov., isolated from soil of a ginseng field. Int J SystEvolMicrobiol 61: 2606–2609.

57. Wang L, An DS, Kim SG, Jin FX, Kim SC, Lee ST, Im WT (2012) Ramlibacterginsenosidimutans sp. nov., with ginsenoside-converting activity. J MicrobiolBiotech 22: 311–315.

58. An DS, Cui CH, Sung BH, Yang HC, Kim SC, et al. (2012) Characterization of a novel ginsenoside-hydrolyzing α-L-arabinofuranosidase, AbfA, from Rhodanobacterginsenosidimutans Gsoil 3054T. ApplMicrobiolBiotech 94: 673–682.

59. Wang L, Liu QM, Sung BH, An DS, Lee HG, et al. (2011) Bioconversion of ginsenosides Rb1, Rb2, Rc and Rd by novel beta-glucosidase hydrolyzing outer 3-O glycoside from Sphingomonas sp. 2F2: cloning, expression, and enzyme characterization. J Biotechnol 156: 125–133.

60. Korea Food & Drug Administration (KFDA). Production record of health functional food and analysis result of that in 2011 [Seoul, Korea].

61. Bae EA, Shin JE, Kim DH (2005) Metabolism of ginsenoside Re by human intestinal microflora and its estrogenic effect. Biol Pharm Bull 28: 1903–1908.

62. Chi H, Ji GE (2005) Transformation of ginsenosides Rb1 and Re from Panax ginseng by food microorganisms. BiotechnolLett 27: 765–771.

63. Ko SR, Choi KJ, Uchida K, Suzuki Y (2003) Enzymatic preparation of ginsenosides Rg2, Rh1, and F1 from protopanaxatriol-type ginseng saponin mixture. Planta Med 69: 285–286.

64. Liu QM, Jung HM, Cui CH, Sung BH, Kim JK, et al. (2013) Bioconversion of ginsenoside Rc into Rd by a novel α-L-arabinofuranosidase, Abf22-3 from Leuconostoc sp. 22-3: cloning, expression, and enzyme characterization. AntonLeeuwInt J G 103: 747–754.

65. Jin XF, Yu HS, Wang DM, Liu TQ, Liu CY, et al. (2012) Kinetics of a cloned special ginsenosidase hydrolyzing 3-O-glucoside of multi-protopanaxadiol-type ginsenosides, named ginsenosidase type III. J Microbiol Biotechnol 22: 343–351.