ABSTRACT
Plant polyphenols are of great interest for drug discovery and drug development since many of these compounds have health-promoting activities as treatments against various diseases, such as diabetes, cancer, or heart diseases. However, the limited availability of polyphenols represents a major obstacle to clinical applications that must be overcome. In comparison to the quantities of these compounds obtained by isolation from natural sources or costly chemical synthesis, the microbial production of these compounds could provide sufficient quantities from inexpensive substrates. In this work, we describe the development of an Escherichia coli platform strain for the production of pinosylvin, a stilbene found in the heartwood of pine trees which could aid in the treatment of various cancers and cardiovascular diseases. Initially, several configurations of the three-step biosynthetic pathway to pinosylvin were constructed from a set of two different enzymes for each enzymatic step. After optimization of gene expression and evaluation of different construct environments, low pinosylvin concentrations up to 3 mg/liter could be detected. Analysis of the precursor supply and a comparative analysis of the intracellular pools of pathway intermediates and product identified the limited malonyl coenzyme A (malonyl-CoA) availability and low stilbene synthase activity in the heterologous host to be the main bottlenecks during pinosylvin production. Addition of cerulenin for increasing intracellular malonyl-CoA pools and the in vivo evolution of the stilbene synthase from Pinus strobus for improved activity in E. coli proved to be the keys to elevated product titers. These measures allowed product titers of 70 mg/liter pinosylvin from glucose, which could be further increased to 91 mg/liter by the addition of l-phenylalanine.
INTRODUCTION
Plant polyphenols are plant secondary metabolites which have important functions in plant defense and signaling or serve as pigments. Furthermore, many polyphenols have antioxidant, anticancer, and anti-inflammatory properties, rendering these metabolites an invaluable source of bioactive compounds for drug discovery and drug development in the pharmaceutical industry (1).
The polyphenol pinosylvin (trans-3,5-dihydroxystilbene) is a stilbene which can be predominantly found in the heartwood of conifer trees of the genus Pinus (2). Pinosylvin protects the plants against microbial and fungal decay (3), but it has also attracted attention due to its health benefits for humans. Several recent studies indicate a positive effect of pinosylvin in the treatment of various cancers (4, 5), cardiovascular inflammatory diseases (6), and adjuvant arthritis (7). Pinosylvin is synthesized from the aromatic amino acid l-phenylalanine in three enzymatic steps (Fig. 1) (8). l-Phenylalanine is first deaminated to trans-cinnamic acid by a phenylalanine ammonia lyase (PAL; EC 4.3.1.24). The resulting acid is subsequently coenzyme A (CoA) activated by a 4-coumarate-CoA ligase (4CL; EC 6.2.1.12), yielding the CoA thioester trans-cinnamoyl-CoA. A stilbene synthase (STS; EC 2.3.1.146), which is a type III polyketide synthase (PKS), catalyzes the successive condensation of three malonyl-CoA molecules with trans-cinnamoyl-CoA, forming a linear tetraketide intermediate. Finally, the STS also cyclizes the tetraketide intermediate via an intramolecular C-2 → C-7 aldol condensation to form the stilbene pinosylvin (8).
Biosynthetic pathway from l-phenylalanine to pinosylvin. Abbreviations: PAL, phenylalanine ammonia lyase; 4CL, 4-coumarate-CoA ligase; STS, stilbene synthase.
The concentration of pinosylvin in pine heartwood ranges from 1 to 40 mg/g (dry weight) (including its monomethyl ether) (3), rendering access to pinosylvin through extraction difficult. Furthermore, purification of pinosylvin from these extracts would require separation from a multitude of often structurally very similar compounds, such as other stilbenes or flavonoids (9). Various catalytic and noncatalytic chemical synthetic routes for the synthesis of stilbenes have been devised. However, these strategies suffer from the demand for expensive substrates, catalyst instability, or degradation or the formation of by-products and various isomers (10). In contrast, microbial production of pinosylvin from glucose would have the advantage of being much more environmentally friendly than chemical synthesis, since it avoids the use of organic solvents, heavy metals, and strong acids or bases. To date, several microorganisms have been engineered for microbial stilbene production, especially for the synthesis of resveratrol (11–13). However, stilbene production with these strains relies on the supplementation of phenylpropanoids (e.g., coumaric acid or cinnamic acid) as stilbene precursors. This strategy has led to some considerable success with regard to production of the stilbene resveratrol, where feeding of 15 mM p-coumaric acid resulted in maximum product titers of 1.4 g/liter in Escherichia coli (11), which could be further improved to 1.6 g/liter by increasing the carbon flux into malonyl-CoA (14). More recently, E. coli has also been engineered to produce 35 mg/liter resveratrol from 3 mM its amino acid precursor l-tyrosine (15). This was achieved by following a modular metabolic engineering strategy for balancing the pathway with the matB-matC system from Rhizobium trifolii, which synthesizes the precursor malonyl-CoA from supplemented malonate. In comparison, product titers reported for pinosylvin are modest, with 0.6 mg/liter being achieved with Streptomyces venezuelae (when supplementing 1.2 mM trans-cinnamic acid) (16) and 155 mg/liter being achieved during biotransformation with E. coli (when providing 1 mM trans-cinnamic acid and augmenting intracellular malonyl-CoA levels) (17).
In this study, we examined multiple genes and pathway configurations for microbial pinosylvin production by E. coli from glucose. For this purpose, we optimized heterologous gene expression, examined precursor availability, and performed laboratory protein evolution to enhance the activity of a plant-derived STS in the microbial host.
MATERIALS AND METHODS
Bacterial strains, plasmids, and growth conditions.The bacterial strains and plasmids used or constructed in the course of this work are listed in Table 1. E. coli DH5α, solely used for cloning purposes, was routinely cultivated on a rotary shaker (170 rpm) at 37°C in LB medium or on LB plates (LB medium with 1.5% [wt/vol] agar) (18). If appropriate, carbenicillin (50 μg/ml), kanamycin (30 μg/ml), or spectinomycin (50 μg/ml) was added to the medium. Growth was determined by measuring the optical density at 600 nm (OD600).
Strains and plasmids used in this study
Recombinant DNA techniques.The enzymes for recombinant DNA work were obtained from Thermo Scientific (Waltham, MA, USA) and Merck Millipore (Billerica, MA, USA). Routine methods like PCR, restriction, or ligation were carried out according to standard protocols (18). Oligonucleotides used for cloning, error-prone PCR (epPCR), and site-directed mutagenesis (19) were obtained from Eurofins Genomics (Ebersberg, Germany) and are listed in Table 2. E. coli was transformed by the RbCl method (20). The PAL, 4CL, and STS genes were chemically synthesized by GeneArt Gene Synthesis Services from Life Technologies, a Thermo Fisher Scientific company (Waltham, MA, USA). Genes were codon optimized for expression in E. coli by employing the proprietary GeneOptimizer software. For subcloning purposes, all PAL and STS genes were synthesized with a 5′ NcoI restriction site and a 3′ BamHI restriction site, whereas the synthetic 4CL genes incorporate a 5′ NdeI site and a 3′ XhoI site. The 5′ NdeI and 5′ NcoI restriction sites already include the start codon ATG as part of the recognition sites. However, in the case of the sts2 gene from Pinus strobus, subcloning using the NcoI site required the addition of the glycine codon GGA after the start codon since the NcoI recognition site (5′-CCATGG-3′) comprises a G following the start ATG codon.
Oligonucleotides used in this study
Heterologous gene expression.Recombinant E. coli BL21(DE3) strains expressing genes for the pinosylvin synthetic pathway were cultivated in 25 ml LB medium in 100-ml baffled shake flasks at 37°C and 120 rpm until an OD600 of 0.6 was reached. After induction with 1 mM isopropyl-β-d-thiogalactopyranoside (IPTG), the cultivation continued at 30°C for 24 h. Cells were then harvested by centrifugation (3,900 × g, 15 min, 4°C) and washed with 5 ml 0.9% (wt/vol) NaCl solution. After centrifugation (3,900 × g, 15 min, 4°C), the cells were resuspended in 1 ml phosphate-buffered saline (1.9 mM NaH2PO4, 8.1 mM Na2HPO4, 154 mM NaCl, pH 7.4) and disrupted by sonication using a Branson sonifier 250 (intensity, 2; duty cycle, 20%; 4 min; Branson, Danbury, CT, USA). After removal of the cellular debris by centrifugation (16,000 × g, 30 min, 4°C), the cell extract (CE) was subjected to SDS-PAGE analysis using NuPAGE bis-Tris gels from Life Technologies (Waltham, MA, USA).
For purification of the His6-tagged STS proteins, E. coli BL21(DE3) cultures were cultivated in 100 ml LB medium in 500-ml baffled shake flasks on a rotary shaker (120 rpm) at 37°C until an OD600 of 0.2 was reached. The cultures were shifted to 30°C, and gene expression was induced with 0.1 mM IPTG as soon as an OD600 of 0.6 was reached. After 4 h, the cells were harvested by centrifugation (6,200 × g, 10 min, 4°C) and washed with an ice-cold 0.9% (wt/vol) NaCl solution. After centrifugation (6,200 × g, 10 min, 4°C), the cells were resuspended in 1 ml NNI10 lysis buffer (50 mM NaH2PO4, 300 mM NaCl, 10 mM imidazole, pH 8) and subsequently disrupted by sonication (intensity, 2; duty cycle, 20%; 4 min). The resulting CE was clarified by centrifugation (16,000 × g, 30 min, 4°C). A column volume (CV) of 1 ml nickel-nitrilotriacetic acid resin (Qiagen, Hilden, Germany) was used for protein purification. After equilibration of the column by applying 5 CVs of NNI10 buffer, the CE was loaded on the column. Subsequently, the column was washed three times with 10 CVs of NNI70 wash buffer (50 mM NaH2PO4, 300 mM NaCl, 70 mM imidazole, pH 8). The elution of the STS proteins was achieved by applying 1 CV of NNI250 elution buffer (50 mM NaH2PO4, 300 mM NaCl, 250 mM imidazole, pH 8) six times. The eluted fractions were analyzed by SDS-PAGE.
Furthermore, the successful expression of all PAL, 4CL, and His-STS proteins was verified by matrix-assisted laser desorption ionization–time of flight (MALDI-TOF) mass spectrometry (MS) as previously described (21).
Pinosylvin production with E. coli.Recombinant E. coli BL21(DE3) strains were cultivated either in 50 ml LB broth or in 50 ml yeast nitrogen base (YNB) defined medium in 500-ml baffled shake flasks on a rotary shaker (120 rpm) at 26°C. For the preparation of 1 liter of YNB medium, 100 ml of 10-times-concentrated YNB was added to 900 ml base medium [6 g K2HPO4, 3 g KH2PO4, 10 g 3-(N-morpholino)propanesulfonic acid (MOPS), pH 7], and 5 g/liter glucose was added as a carbon source. Precultures were cultivated overnight at 37°C (170 rpm) in LB or YNB medium and then diluted to an OD600 of 0.1 in 50 ml fresh LB or YNB medium, respectively. When the OD600 of these main cultures reached 0.6, gene expression was induced with 1 mM IPTG (when necessary), and the cultivation was continued for 36 h at 26°C. Samples were taken at regular time intervals for liquid chromatography (LC)-MS analysis to follow the formation of pinosylvin.
A two-stage cultivation procedure was also employed for pinosylvin production. In the first stage, cultivation and gene expression were performed in LB as described above. After 5 h of heterologous gene expression, cells were collected by centrifugation (4,300 × g, 15 min, 4°C) and subsequently resuspended in 10 ml YNB medium–5 g/liter glucose, which was supplemented with 1 mM IPTG. The cultivation of this second phase was continued for 36 h at 26°C. Again, samples were taken at regular time intervals for LC-MS analysis.
Metabolite extraction.Metabolite extracts from supernatants and cells were prepared for LC-MS analysis by mixing 1 ml culture with 1 ml ethyl acetate and vigorous shaking (1,400 rpm, 10 min, 20°C) in an Eppendorf thermomixer (Hamburg, Germany). This suspension (800 μl) was centrifuged for 5 min at 16,000 × g, and the ethyl acetate layer was transferred to an organic solvent-resistant deep-well plate (Eppendorf, Hamburg, Germany). After evaporation of the ethyl acetate, 100 µl acetonitrile was added prior to LC-MS analysis, and this concentrated the extract 8-fold.
Determination of cytoplasmic metabolite concentrations.E. coli cells were separated from the medium and inactivated by silicone oil centrifugation for the determination of cytoplasmic concentrations of trans-cinnamic acid, trans-cinnamoyl CoA, malonyl-CoA, and pinosylvin as described previously (22). The resulting upper phase was used for further analysis by LC-MS/MS.
Chemical synthesis of trans-cinnamoyl-CoA.trans-Cinnamic acid was esterified in the presence of CoA and dicyclohexyl carbodiimide as previously described to prepare a standard for the LC-MS/MS analysis of trans-cinnamoyl-CoA (23). All required chemicals were purchased from Sigma-Aldrich (St. Louis, MO, USA).
Determination of intermediates and pinosylvin by LC-MS and LC-MS/MS.The trans-cinnamic acid and pinosylvin in the cell-free culture medium as well as in the ethyl acetate-extracted samples were quantified by LC-MS using an ultra-high-performance LC (uHPLC) 1290 Infinity system coupled to a 6130 quadrupole LC-MS system (Agilent, Santa Clara, CA, USA). LC separation was carried out with a Kinetex 1.7u C18 100-Å-pore-size column (50 mm by 2.1 mm [internal diameter]; Phenomenex, Torrance, CA, USA) at 50°C. For elution, 0.1% acetic acid (solvent A) and acetonitrile supplemented with 0.1% acetic acid (solvent B) were applied as the mobile phases at a flow rate of 0.5 ml/min. A gradient was used, where solvent B was increased stepwise (minute 0 to 1, 15% to 22%; minute 1 to 2, 22% to 40%; minute 2 to 2.5, 40% to 50%; minute 2.5 to 3, 50% to 100%). The mass spectrometer was operated in the negative electrospray ionization (ESI) mode, and data acquisition was performed in selected-ion-monitoring mode. LC-MS/MS was used for the quantification of trans-cinnamic acid, trans-cinnamoyl-CoA, malonyl-CoA, and pinosylvin in the cytoplasmic extracts following a previously described method (24). trans-Cinnamic acid, malonyl-CoA, and pinosylvin standards were purchased from Sigma-Aldrich (St. Louis, MO, USA).
Construction of Pstrsts2 epPCR libraries and microtiter plate screening for increased pinosylvin production.All epPCRs of the sts2 gene from P. strobes (Pstrsts2), with mutation rates varying from 2 to 4.6 mutations/kb, were performed using a Diversify PCR random mutagenesis kit from Clontech (Mountain View, CA, USA). The mutated Pstrsts2 genes were subcloned into the pRSFDuet1 backbone carrying the genes for PAL isozyme 1 from Petroselinum crispum PAL1 (PcPAL1) and 4CL from Streptomyces coelicolor (Sc4CL), resulting in the vector pR-HisPstrsts2*-Sc4cl-Pcpal1, and this vector was transformed into E. coli BL21(DE3). At least 90 clones for each mutation rate, always accompanied by a negative control [E. coli BL21(DE3)/pRSFDuet1] and a strain expressing the wild-type sts2 gene [E. coli BL21(DE3)/pR-HisPstrsts2-Sc4cl-Pcpal1], were screened for improved pinosylvin titers. For this purpose, cultivation and screening were performed in the 96-well microtiter plate format. The libraries were cultivated in 600 μl YNB medium supplemented with 5 g/liter glucose, 3 mM l-phenylalanine, 0.5 mM IPTG, and 30 μg/ml kanamycin in deep-well plates from Eppendorf (Hamburg, Germany) and a Microtron shaker from Infors (Bottmingen, Switzerland) at 30°C, 900 rpm, and 75% humidity for 24 h. Subsequently, 100 μl culture from each well was transferred to a UV-transparent 96-well flat-bottom plate from Corning (Corning, NY, USA). An Infinite M200 Pro microplate reader (Tecan, Männedorf, Switzerland) was used to compare relative pinosylvin titers on the basis of the fluorescent properties of pinosylvin (excitation λ = 316 nm, emission λ = 388 nm).
RESULTS
Design of a synthetic pathway for pinosylvin.For the construction of a synthetic pinosylvin pathway, we first consulted the recent literature and the BRaunschweig ENzyme DAtabase (BRENDA) to identify suitable phenylalanine ammonia lyases (PALs), 4-coumarate-CoA ligases (4CLs), and stilbene synthases (STSs) originating from either plants or microorganisms (25). Enzymes were considered suitable if they had already been functionally expressed in E. coli and if the determined basic catalytic parameters (e.g., kcat, Km) appeared to be promising for the construction of an efficient route toward the biosynthesis of pinosylvin. In the course of that survey, we selected PAL isozyme 1 from Petroselinum crispum (PcPAL1) and PAL isozyme 2 from Arabidopsis thaliana (AtPAL2), since both enzymes are characterized by high levels of activity when expressed in E. coli. Furthermore, we decided to express a variant of a 4CL-like enzyme from the Gram-positive bacterium Streptomyces coelicolor (Sc4CL), which favors trans-cinnamic acid over p-coumaric acid as the substrate, and a variant of 4CL isozyme 2 of Arabidopsis thaliana (At4CL2). As candidates for stilbene synthases, we selected STS isozyme 3 of Pinus densiflora (PdSTS3) and STS isozyme 2 of Pinus strobus (PstrSTS2), as these isozymes were reported to have high specificity and activity toward cinnamoyl-CoA.
Optimization of heterologous gene expression.For initial gene expression experiments, all PAL genes were cloned into the pETDuet1 vector, whereas the STS and 4CL genes were cloned into pCDFDuet1 vector under the individual control of a T7 promoter. The successful expression of soluble protein could be verified for all PALs and 4CLs by SDS-PAGE analysis. However, no visible bands indicating successful expression of PdSTS3 or PstrSTS2 could be observed in the soluble or the insoluble protein fractions (data not shown). Neither an incremental reduction of the expression temperature from 30°C to 15°C nor variation of the IPTG concentration (0.4 mM to 1 mM) improved STS expression. Furthermore, a first analysis of the supernatant of the expression cultures showed that up to 1.7 mM trans-cinnamic acid but no pinosylvin was produced by the recombinant E. coli strains. In order to improve the expression of PdSTS3 and PstrSTS2, the genes for both were individually cloned into the pRSFDuet1 vector, another member of the pDuet vector family bearing the origin of replication from the plasmid RSF1030, which allows for a higher number of plasmid copies per cell (>100 copies per cell). In addition, both genes were cloned with an N-terminal His6 tag, since such a tag improved the soluble expression of PdSTS3 without affecting the activity of STS in in vitro enzyme assays (26). Despite these alterations, no soluble STS expression was observable in CEs, but both His6-tagged enzymes (HisPstrSTS2 and HisPdSTS3) could be identified in the whole-cell fraction. Purification of the enzymes from the CE by affinity chromatography finally revealed a low level of expression of soluble HisPstrSTS2 protein but not HisPdSTS3 in the CE (see Fig. S1 in the supplemental material). Possibly, native E. coli proteins concealed the band corresponding to HisPstrSTS2 on SDS-PAGE analysis of the CE. Finally, the identity of all PALs, 4CLs, and HisSTS proteins was confirmed by MALDI-TOF MS (data not shown). The accompanying chemical analysis of the cultures during the optimization of heterologous gene expression revealed that the modification of STS expression already resulted in the accumulation of low concentrations of pinosylvin (0.068 ± 0.004 mg/liter), but only in extracts of strains expressing HisSTS of P. strobus.
Identification of the most suitable pathway configuration.After the identification of His6-tagged PstrSTS2 as the only functional STS in E. coli, four different synthetic pathways with all combinations of the available PALs and 4CLs were constructed. Subsequently, the pinosylvin production of all four strains was compared using the two-phase cultivation conditions described in Materials and Methods. As it turned out, the combination of the PAL from P. crispum, the 4CL from S. coelicolor, and the STS from P. strobus proved to be the most productive pathway variant, but only a low product concentration of up to 0.64 mg/liter pinosylvin could be detected (Table 3). Interestingly, the combination of PcPAL1 and At4CL2 did not lead to any detectable product formation, even though both enzymes functioned in other pathway configurations. The positive effect of the His6 tag on PstrSTS2 expression was also reflected by increased pinosylvin titers, as the same strain expressing PstrSTS2 without the His6 tag accumulated only up to 0.22 mg/liter pinosylvin (Table 3).
Pinosylvin production with E. coli strains harboring different variants of the pinosylvin pathway in mono- or tricistronic transcriptional units
At this point, all genes of the pinosylvin pathway were expressed under the control of their own T7 promoter to rule out any polar expression effects. With the aim of evaluating the expression of the pathway from a single three-gene operon, two additional plasmids were constructed in which the operon was expressed under the control of either the IPTG-inducible T7 promoter (using the pRSFDuet1 vector backbone) or the constitutive gap promoter (Pgap) of the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from E. coli (using the pUC18 vector backbone) (27). The gap promoter was selected, as it induces gene expression during cultivation on glucose as the carbon source, thus allowing the continuous transcription of the three heterologous genes. Surprisingly, no pinosylvin formation could be detected when the synthetic operon was constitutively expressed during single-phase cultivation in YNB supplemented with 5 g/liter glucose. Only the accumulation of up to 10 mg/liter trans-cinnamic acid during the cultivation could be measured. In contrast, IPTG-induced expression from the T7 promoter of the tricistronic pRSFDuet1 construct led to a 6-fold increased pinosylvin titer with E. coli BL21(DE3)/pR-HisPstrsts2-Sc4cl-Pcpal1 (up to 3.37 mg/liter) compared to that achieved with the monocistronic organization of the same genes in E. coli BL21(DE3)/pE-Pcpal1-Sc4cl pR-HisPstrsts2 (Table 3).
Availability of precursors and intermediates during pinosylvin synthesis.At this stage, insufficient levels of the precursors l-phenylalanine and malonyl-CoA could be limiting for the overall product titers. In order to find out whether the availability of l-phenylalanine or malonyl-CoA was limiting, cultivations were performed in which 3 mM l-phenylalanine and/or cerulenin (up to 200 μM) was supplemented during the production phase (Fig. 2). Cerulenin is an antifungal antibiotic produced by Cephalosporium caerulens, which blocks fatty acid biosynthesis by inhibiting the β-ketoacyl-acyl carrier protein (ACP) synthases FabB and FabF, thereby preventing channeling of malonyl-CoA into the pathway for fatty acid synthesis (28). The exclusive supplementation of l-phenylalanine had no positive effect on pinosylvin production (3.29 ± 0.11 mg/liter without supplementation versus 3.49 ± 0.42 mg/liter with supplementation) but resulted in the accumulation of the intermediate trans-cinnamic acid in the supernatant (data not shown). In contrast, addition of cerulenin drastically increased product titers up to 18-fold, reaching 59 mg/liter pinosylvin at a concentration of 200 μM cerulenin (and no addition of l-phenylalanine).
Pinosylvin production of E. coli BL21(DE3)/pR-HisPstrsts2-Sc4cl-Pcpal1 with and without supplementation of l-phenylalanine and/or addition of cerulenin. The experiments were performed in triplicate.
Subsequently, the cytoplasmic levels of malonyl-CoA in E. coli BL21(DE3) were determined to verify that the addition of cerulenin had an effect on the intracellular malonyl-CoA pool. These experiments revealed that the intracellular malonyl-CoA level increased from 2 pmol/mg cell dry weight (CDW) (no cerulenin) to 105 pmol/mg CDW at 1.5 h after the addition of 200 μM cerulenin. In contrast, E. coli BL21(DE3)/pR-HisPstrsts2-Sc4cl-Pcpal1 accumulated malonyl-CoA only up to 59 pmol/mg CDW in the presence of 200 μM cerulenin and during pinosylvin synthesis (see Fig. S2 in the supplemental material), indicating malonyl-CoA consumption during pinosylvin formation. l-Phenylalanine supplementation had a positive impact on product formation during cerulenin-mediated inhibition of fatty acid synthesis, indicating that the availability of this precursor becomes limiting when more malonyl-CoA is available (Fig. 2). At this stage, supplementation of 3 mM l-phenylalanine and addition of 200 μM cerulenin led to the highest pinosylvin titer of 70 mg/liter after 36 h of cultivation.
Two additional strategies for improving the intracellular malonyl-CoA availability in E. coli to circumvent cerulenin addition were pursued: overexpression of the intrinsic β-ketoacyl-ACP synthase II (FabF) for reducing the drain of malonyl-CoA into fatty acid synthesis, as well as heterologous expression of the acetyl-CoA carboxylase (ACC) from Corynebacterium glutamicum for increasing the intracellular generation of malonyl-CoA from acetyl-CoA. Unfortunately, heterologous expression of either enzyme led only to reduced pinosylvin titers (data not shown).
The availability of the intermediate trans-cinnamoyl-CoA during pinosylvin synthesis could also be limiting due to low 4CL activity. Therefore, the relative intracellular levels of trans-cinnamic acid, trans-cinnamoyl-CoA, and pinosylvin were determined during cultivation with 3 mM l-phenylalanine and with and without addition of cerulenin in the production phase (Fig. 3). During these experiments, l-phenylalanine was always supplemented to ensure that trans-cinnamic acid, as the direct precursor of trans-cinnamoyl-CoA, was not limiting in the presence of cerulenin. All measurements were conducted 2 h after addition of cerulenin, but only the relative levels of intermediates and product could be determined, since no trans-cinnamoyl-CoA standard of sufficient purity could be chemically synthesized. As a result, addition of 200 μM cerulenin led to an 18-fold increase of pinosylvin, accompanied by a drop of the trans-cinnamoyl-CoA level of only 30%. This observation hints that the levels of trans-cinnamoyl-CoA are sufficient to achieve higher product titers, as the pool of this intermediate is not depleted or even significantly reduced.
Relative intracellular levels of trans-cinnamic acid, trans-cinnamoyl-CoA, and pinosylvin determined by LC-MS/MS. The experiments with E. coli BL21(DE3)/pR-HisPstrsts2-Sc4cl-Pcpal1 were performed in duplicate with supplementation of 3 mM l-phenylalanine and 0 μM or 200 μM cerulenin. The relative levels are given as the area (percent), with the areas for the intermediates determined with 0 μM cerulenin being set equal to 100%.
Directed protein evolution of HisPstrSTS2.The PAL and 4CL activity in E. coli appeared to be sufficient to support higher product titers, as the intracellular pools of trans-cinnamic acid and trans-cinnamoyl-CoA were not depleted during pinosylvin production. In contrast, expression of the STS was still low, even though the sequences of all STS genes were optimized for expression in E. coli prior to gene synthesis. Hence, the poor soluble expression of HisPstrSTS2 in the microbial host might be limiting the performance of the pinosylvin pathway. With the aim to adapt the STS2 of P. strobus for optimal expression in E. coli, we performed directed evolution of the HisPstrSTS2 in vivo. For this purpose, we randomly mutated the codon-optimized open reading frame of HisPstrsts2 by epPCR and subcloned the resulting STS library into pR-Sc4cl-Pcpal to complete the pinosylvin pathway. Cultivation and screening of the library were performed in the microtiter plate (MTP) format, taking advantage of the fluorescence of pinosylvin for identifying clones with increased product titers in comparison to the titer of the starting strain, E. coli BL21(DE3)/pR-HisPstrsts2-Sc4cl-Pcpal1. Screening of 450 clones yielded 2 clones exhibiting a higher fluorescence under screening conditions. DNA sequencing revealed that each STS gene carries a single point mutation (A742G and A1082G), both of which lead to amino acid substitutions (T248A and Q361R, respectively). After recloning of the STS variants and retransformation to exclude any effect of the plasmid or strain background, the product titers of both variants were compared to the titer of the starting strain using the optimized cultivation conditions with and without the addition of cerulenin or the supplementation of l-phenylalanine.
Under all cultivation conditions tested, both amino substitutions in HisPstrSTS2 individually improved the overall product titers (Fig. 4). In a direct comparison, the positive effect of T248A for pinosylvin formation was more pronounced, as more pinosylvin could be detected in the extracts. In the presence of cerulenin, but without supplementation of l-phenylalanine, 70 mg/liter pinosylvin accumulated in the E. coli BL21(DE3)/pR-HisPstrsts2T248A-Sc4cl-Pcpal1 cultures, whereas an additional supplementation of 3 mM l-phenylalanine increased the pinosylvin titer to 91 mg/liter. However, use of the combination of both amino acid substitutions in an attempt to further boost product formation failed and most likely resulted in an inactive STS variant, as no pinosylvin synthesis of E. coli BL21(DE3)/pR-HisPstrsts2T248A/Q361R-Sc4cl-Pcpal1 was detected.
Comparison of the pinosylvin titers of E. coli BL21(DE3)/pR-HisPstrsts2-Sc4cl-Pcpal1 expressing either wild-type HisPstrsts2 (light gray bars), HisPstrsts2 with the Q361R substitution (dark gray bars), or HisPstrsts2 with the T248A substitution (black bars) without any supplementation, with addition of 200 μM cerulenin, or with addition of 200 μM cerulenin and supplementation of 3 mM l-phenylalanine. The experiments were performed in duplicate.
DISCUSSION
The development of efficient microbial platform organisms for the production of (plant) natural products requires the identification of suitable enzymes, the design of stable genetic constructs for gene expression, and the optimal adaptation of the synthetic pathway to the host cell metabolism (or vice versa). In this work, we systematically explored several strategies to engineer E. coli for the microbial production of the stilbene pinosylvin. The efforts included the optimization of gene expression conditions, the comparison of different expression constructs, and protein engineering to adapt the expression of the plant-derived STS2 from P. strobus to the microbial host system.
Since construction of stilbene-producing microorganisms has mostly been limited to strains requiring the supplementation of phenylpropanoid intermediates, we decided to design and implement the complete three-step pathway to pinosylvin starting from the amino acid l-phenylalanine. We set off with the identification of the most suitable enzymes from various microbial or plant sources in the literature and databases and selected two enzymes for each of the three required enzymatic steps based on their kinetic properties. PALs, catalyzing the first committed step toward phenylpropanoid synthesis, can be found in all higher plants but are also present in yeasts, fungi, and some prokaryotic species (29). PAL isozyme 1 from Petroselinum crispum (PcPAL1) and PAL isozyme 2 from Arabidopsis thaliana (AtPAL2) were selected for the construction of the pinosylvin pathway since both enzymes are characterized by high levels of activity when expressed in E. coli (30, 31). In comparison to PALs, 4CLs are specific to the secondary metabolism of plants. Most 4CLs described are substrate specific, since all described natural substrates are characterized by a 4′-hydroxyl group on the phenyl ring (32). However, a synthetic pathway leading to pinosylvin requires a CoA-ligase that accepts trans-cinnamic acid lacking this para-hydroxy moiety. Interestingly, Kaneko and coworkers discovered a bacterial 4CL-like enzyme from the filamentous, soil-dwelling, Gram-positive bacterium Streptomyces coelicolor (33). The enzyme showed a distinct 4CL activity, favoring trans-cinnamic acid over p-coumaric acid as the substrate. Mutational analysis of the substrate binding pocket indicated that the amino acid substitution A294G conferred even higher catalytic activity toward trans-cinnamic acid than the native enzyme. This particular enzyme variant (Sc4CL) was selected for evaluation in the synthetic pinosylvin pathway. The second 4CL candidate was a variant of 4CL isozyme 2 of Arabidopsis thaliana (At4CL2) bearing three amino acid substitutions (N256A, M293P, K320L). In in vitro enzyme assays, this mutein demonstrated a 30-fold higher level of conversion of trans-cinnamic acid than the native enzyme (34). Pinosylvin-forming stilbene synthases (EC 2.3.1.146) are typical for gymnosperms, and to date several such enzymes have been identified in three pine species, Pinus sylvestris (35), Pinus densiflora (26), and Pinus strobus (36). STS isozyme 3 of P. densiflora (PdSTS3) and STS isozyme 2 of P. strobus (PstrSTS2) were selected as STS candidates, as these isozymes were reported to have higher specificity and activity toward cinnamoyl-CoA than the other isozymes in in vitro enzyme assays (26, 36). Interestingly, compared to the other characterized STSs, PdSTS3 has a truncated C terminus, which is believed to be responsible for the observed release from feedback inhibition by pinosylvin (26).
The combination of PcPAL1 and Sc4CL in the synthetic pinosylvin pathway turned out to be more beneficial than the combination of AtPAL2 and Sc4CL to achieve high product titers. This is interesting, because At4CL2 performed better with AtPAL2 than with the PAL from P. crispum. This could simply be due to the common origin of both enzymes from A. thaliana but also supports the notion that the performance of an entire synthetic pathway instead of the performance of a pathway just assembled from individually characterized best parts should be evaluated. Pinosylvin production could be improved by organizing the three genes of the final pinosylvin pathway configuration in an operon under the control of a single promoter. The same effect was also observed during the development of a microbial production strain for resveratrol, where cotranscription of the genes for 4CLs and STSs was believed to lead to more balanced gene expression (11).
The availability of malonyl-CoA turned out to be crucial for the overall product titers, as the addition of the fatty acid production inhibitor cerulenin boosted the pinosylvin titers 18- to 20-fold, whereas the supplementation of l-phenylalanine alone had no effect. In order to avoid the addition of cerulenin, two strategies for increasing the availability of intracellular malonyl-CoA were followed: overexpression of the E. coli FabF (37) as well as heterologous expression of the ACC from C. glutamicum, as was successfully demonstrated during microbial polyketide synthesis and flavanone synthesis, respectively (17, 38, 39). However, none of these approaches had the desired effect, as pinosylvin production significantly decreased or ceased entirely. No improvement of microbial flavonoid production was also observed during microbial flavonoid production with E. coli upon expression of the ACC from C. glutamicum, suggesting that this strategy is not generally applicable (40). A comparative analysis of intracellular intermediate availability and product concentration in the presence or absence of 200 μM cerulenin was conducted, but the results obtained did not reveal a limiting enzymatic step when more malonyl-CoA was available. Additional experiments will be necessary to shed more light on this subject.
Finally, we performed directed protein evolution to adapt the expression of the pine-derived HisPstrSTS2 to the microbial host system. This was directly done in the genetic context of the pinosylvin pathway. Only one round of diversity generation by epPCR and MTP-based screening of 450 clones for increased pinosylvin-related fluorescence identified two clones, each bearing a single amino acid substitution, with improved product formation. The combination of both mutations had a detrimental effect on product formation, but the T248A substitution alone increased the final pinosylvin titer of the best strain, E. coli BL21(DE3)/pR-HisPstrsts2-Sc4cl-Pcpal1, from 59 mg/liter to 70 mg/liter in the presence of 200 μM cerulenin and by an additional 30% to 91 mg/liter when l-phenylalanine was also supplemented. With the aim to understand the structural consequences of both amino acid substitutions, we generated a structure model for a PstrSTS2 dimer using the SWISS-MODEL server (41). A structure model for the pinosylvin synthase of P. sylvestris (PDB accession number 1U0U) (42), which was not among the pinosylvin synthases tested by us, served as the template for the calculated model since this enzyme shares 87% sequence identity to PstrSTS2 on the protein level. A structural alignment with structure models for the closely related chalcone synthase from Medicago sativa with coenzyme A bound (PDB accession number 1BQ6) (43) and the stilbene synthase from Arachis hypogaea with resveratrol bound (PDB accession number 1Z1F) (44) helped us to pinpoint the active site in the calculated PstrSTS2 model. Neither substitution is in close proximity to the assumed active sites of the PstrSTS2 dimer (Fig. 5A). T248 is located in the middle of a long β strand, and its side chain faces the interior of the protein at the dimer interface. The replacement by an l-alanine at this position could simply improve the flexibility of the protein and might ultimately lead to the observed improvement of STS activity in E. coli. In contrast, Q361 is positioned at the end of an α helix on the protein surface and located next to the amino acid sequence N362-G363-C364. This pattern follows the typical N-glycosylation sequence N-X-S/T/C of l-asparagine (N)-linked protein glycosylation, where X can be any amino acid except l-proline (45, 46). Analysis of the PstrSTS2 amino acid sequence with the web server GlycoEP also identified N362 to be a potential site for N-glycosylation in this particular STS (47). E. coli cannot perform l-asparagine N-linked protein glycosylation. Hence, replacement of l-glutamine with the more hydrophilic l-arginine at position 361 might simply lead to an improved stability or solubility of PstrSTS2 in the absence of any glycosylation in the heterologous host. When both amino acid substitutions were combined, no pinosylvin synthesis could be observed, suggesting that this enzyme variant is inactive. The β strand, in which T248 is located, is adjacent to the C-terminal β strand of PstrSTS2, and the α helix of Q361 interacts indirectly with the same C-terminal β strand via a loop (Fig. 5B). Possibly, the combination of both mutations simply destabilizes the architecture of PstrSTS2, rendering this enzyme variant inactive.
(A) Cartoon representation of a close-up of the pinosylvin synthase STS2 from P. strobus. The structure model was calculated using the model of the pinosylvin synthase of P. sylvestris as the template. The monomers are shown in green and purple, and the resveratrol (Res) and coenzyme A (CoA) highlighted in one monomer are shown in the stick mode in yellow and cyan/orange, respectively. T248 is located at the dimer interface, whereas Q361 is located on the protein surface. Both amino acids are shown in the stick mode. (B) Both T248 and Q361, shown in orange, interact with the C-terminal β strand. The β strand in which T248 is located is adjacent to this C-terminal β strand, and the α helix of Q361 interacts indirectly with the same C-terminal β strand via a loop.
The key to success for the microbial production of pinosylvin from glucose without supplementation of any phenylpropanoid intermediate was the use of a combination of various metabolic engineering approaches and protein engineering of a heterologously expressed key enzyme. However, at this point the addition of cerulenin is still required for achieving significant product titers. This could be circumvented by carrying out multiple genetic alterations of the central carbon metabolism to redirect the carbon flux to malonyl-CoA, as has been done for the microbial synthesis of the flavanone naringenin with E. coli (48). Here, the successful modifications suggested by the Optforce model included the overexpression of the genes for the pyruvate dehydrogenase multienzyme complex, phosphoglycerate kinase and acetyl-CoA carboxylase, and deletion of the genes encoding the fumarase and succinyl-CoA synthetase. Following this strategy, 5.6 times higher naringenin titers could be attained.
Nonetheless, protein engineering by directed evolution in the genetic context of the synthetic pathway, as presented here, could also be used to harmonize the entire pinosylvin pathway with E. coli for achieving higher product titers. This strategy might also be useful during the development of other microbial strains for the production of small metabolites of high value.
ACKNOWLEDGMENTS
This work was supported by the BOOST fund PNP-EXPRESS of the NRW Strategieprojekt BioSC and the Helmholtz Initiative Synthetic Biology of the Helmholtz Association.
We thank Petra Geilenkirchen (Forschungszentrum Jülich) for help with the LC-MS/MS analyses, Katharina Neufeld (University of Düsseldorf) for assisting with the chemical synthesis of trans-cinnamoyl-CoA, and Marco Bocola (RWTH Aachen University) for useful advice with the structural interpretation of the amino acid substitutions obtained in the HisPstrSTS2 variants.
FOOTNOTES
- Received 12 September 2014.
- Accepted 12 November 2014.
- Accepted manuscript posted online 14 November 2014.
Supplemental material for this article may be found at http://dx.doi.org/10.1128/AEM.02966-14.
- Copyright © 2015, American Society for Microbiology. All Rights Reserved.