Subscriber access provided by University of Texas Libraries
Article
Regeneration of NADPH coupled with HMG-CoA reductase activity increases squalene synthesis in Saccharomyces cerevisiae Kalaivani Paramasivan, and Sarma Mutturi J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.7b02945 • Publication Date (Web): 28 Aug 2017 Downloaded from http://pubs.acs.org on August 29, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Journal of Agricultural and Food Chemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 39
Journal of Agricultural and Food Chemistry
1
Regeneration of NADPH coupled with HMG-CoA reductase activity increases squalene
2
synthesis in Saccharomyces cerevisiae
3
Kalaivani Paramasivan†,‡ and Sarma Mutturi†,‡,*
4 5
†
6
Research Institute, Mysore, India.
7
‡
Microbiology and Fermentation Technology Department, CSIR-Central Food Technological
Academy of Scientific and Innovative Research, Mysore, New Delhi, India.
8 9 10
11
12
13
14
15
16
17
18
19
20
1 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
21
Page 2 of 39
ABSTRACT
22
23
Although overexpression of tHMG1 gene is a well-known strategy for terpene synthesis in
24
Saccharomyces cerevisiae, the optimal level for tHMG1p has not been established. In the present
25
study, it was observed that two copies of tHMG1 gene on a dual gene expression cassette
26
improved squalene synthesis in laboratory strain by 16.8-fold in comparison to single copy
27
expression. It was also observed that tHMG1p is limited by its co-factor (NADPH), as the
28
NADPH regenerating genes’ viz., ZWF1 and POS5 (full length and without mitochondrial
29
presequence) has led to its increased enzyme activity. Further, it was demonstrated that
30
overexpression of full length POS5 has improved squalene synthesis in cytosol. Finally, when
31
tHMG1 and full length POS5 were co-overexpressed there was a net 27.5-fold increase in
32
squalene when compared to control strain. These results suggest novel strategies to increase
33
squalene accumulation in S. cerevisiae.
34 35
KEYWORDS: S. cerevisiae, Squalene, tHMG1, ZWF1, POS5
36
37 38 39 40 41 42 43 2 ACS Paragon Plus Environment
Page 3 of 39
44
Journal of Agricultural and Food Chemistry
INTRODUCTION:
45 46
S. cerevisiae is a well-studied eukaryotic model organism for production of terpenes due to the
47
presence of specific precursor molecules in its sterol pathway.1-3 S. cerevisiae is advantageous
48
over its prokaryotic counterpart, E. coli, for production of terpenoids owing to its ability to
49
express cytochrome P450 oxidase and presence of cognate reductase enzymes.4 S. cerevisiae has
50
been a cell factory for biosynthesis of diverse terpene molecules.5 Squalene which is an
51
intermediate of ergosterol pathway, is an unsaturated triterpene hydrocarbon (C30H50) and is a
52
precursor for sterol molecules in S. cerevisiae. Squalene has a demand of >1200 tonnes per year6
53
and has an estimated market size of 211 million USD by 2021.7 Squalene is used in cosmetic
54
industry as a skin moisturizer and also acts as antioxidant. In pharmaceutical industry, it is used
55
as adjuvant for vaccines and also used for drug delivery. Squalene has shown to exhibit tumor-
56
suppressive activity by preventing breast cancer and acts as a cardio-protective agent by reducing
57
serum cholesterol level8. Further, squalene has also been used as a dietary supplement (Naziri et
58
al. 2011).6
59
In most metabolic engineering studies of S. cerevisiae towards a terpene production,
60
overexpression of HMG-CoA reductase has become imperative as it is a key regulatory enzyme
61
in the ergosterol pathway.9 HMG-CoA reductase is present as Hmg1p and Hmg2p isomers in S.
62
cerevisiae, among which Hmg1p contributes to at least 83% of net enzyme activity in the wild-
63
type strain.10 Overexpression of a membrane bound Hmg1p leads to endoplasmic reticulum
64
membrane stacks otherwise called as karmellae formation.11 To overcome the regulation caused
65
due to the feedback inhibition by ergosterol, a truncated form of HMG-CoA reductase enzyme
66
has been overexpressed.12 Overexpression of a truncated form of Hmg1p coding only for the
3 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 4 of 39
67
catalytic domain of the protein to bypass its regulation has not increased ergosterol content as
68
expected rather it improved squalene synthesis over 30 times in yeast cells.13 During production
69
of a terpene molecule, different researchers have overexpressed tHMG1 either on an episomal
70
plasmid or integrated single or multiple copies.14,15 However, the copy number of this gene
71
during overexpression was not optimized.
72
Reactions involving HMG-CoA reductase (HMG1) and squalene synthase (ERG9) in the sterol
73
pathway utilizes NADPH as a cofactor for squalene synthesis. NADPH is one of the universal
74
electron carrier in yeast cells.16 S. cerevisiae has three different NAD kinase enzymes namely
75
Utr1p, Yef1p and Pos5p which phosphorylates NAD+ to form NADP+ in cytoplasm and
76
mitochondria.17 NADP+ synthesized by NAD kinases is further reduced to NADPH by NADP+-
77
dependent dehydrogenase.18 NADPH synthesis and regeneration in yeast depends on activities of
78
NADH kinases and NADP+-dependent dehydrogenases. NADH kinase phosphorylates NADH to
79
form NADPH while the dehydrogenases reduces NADP+ to form NADPH. Pos5p has the
80
activity of both NAD kinase as well as NADH kinase, however the latter activity is observed
81
more.19 NADPH regeneration in the oxidative part of pentose phosphate pathway is carried out
82
by two key enzymes such as glucose-6-phosphate dehydrogenase (G6PD coded by ZWF1) and 6-
83
phospho gluconate dehydrogenase. ZWF1 is a key target for increasing regeneration of cytosolic
84
NADPH not only in S. cerevisiae but also in Corynebacterium glutamicum for amino acid
85
production.20 Although ZWF1 and POS5 have been targeted in some studies for increased
86
regeneration of NADPH, neither of these genes were studied exclusively for squalene
87
improvement. Furthermore, overexpression of POS5 was carried by eliminating the
88
mitochondrial presequence, assuming the full length expression of this gene would not enhance
89
cytosolic NADPH.21 Hence the aim of the present study is to address the following questions:
4 ACS Paragon Plus Environment
Page 5 of 39
Journal of Agricultural and Food Chemistry
90
a) Is single copy tHMG1 overexpression optimal for improving squalene?
91
b) Role of regeneration of cytosolic NADPH in squalene synthesis
92
c) Will the Pos5p present in mitochondria affect the cytosolic squalene synthesis?
93 94
MATERIALS AND METHODS:
95
Chemicals: Squalene and ergosterol standards were obtained from Sigma-Aldrich (Bangalore,
96
India). Other HPLC grade chemicals were procured from SRL (Mumbai, India) and Qualigens
97
(Mumbai, India). Gene amplifications were performed using Phusion® high-fidelity DNA
98
polymerase (New England Biolabs, Ipswich, MA, USA) and Taq polymerase (Sigma-Aldrich,
99
Bangalore, India). All restriction and ligase enzymes were from obtained from New England
100
BioLabs (New England Biolabs, Ipswich, MA, USA). GenElute™ Plasmid Miniprep Kit
101
(#PLN70, Sigma-Aldrich) was used for plasmid DNA purification. PureLink® Quick Gel
102
Extraction and PCR Purification Combo Kit (#K220001, Invitrogen) was used for gel extraction
103
and DNA purification, respectively.
104 105
Strain and vectors: The information on strains, vectors and primers used in the present study
106
are shown in Table 1. S. cerevisiae laboratory strain BY4741 (ATCC 201388) (Euroscarf,
107
Germany) was kindly gifted by Prof. Ram Rajasekharan (Lipid Science Dept., CSIR-CFTRI).
108
BY4741 was maintained as glycerol stocks stored at -80oC and was subcultured on YPD (yeast
109
extract, peptone, dextrose) agar plates.
110
medium22 lacking uracil or with G418 wherever appropriate. Ampicillin and G418 were used at
111
100 µg/mL and 200 µg/mL, respectively. Nitrogen source in the SD medium was either
112
(NH4)2SO4 (5 g/L) or glutamic acid (1 g/L). E. coli DH5α was used for propagating the plasmids
Engineered yeast strains were grown either in SD
5 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 6 of 39
113
before transformation into yeast strains. The following expression cassettes were used for
114
transformation in yeast cells: pYES2 (URA3, 2µ, PGal1; Invitrogen, Bangalore); pCEV-G1-Km
115
(kanMX, 2µ, PTEF1, PPGK1) was obtained from Lars Nielsen and Claudia Vickers (Addgene
116
plasmid #46813).23
117 118
Instrumental analysis: An isocratic RP-HPLC was performed using a semi preparative
119
reversed-phase C-18 column (Phenomenox Kinetex, Hyderabad, India) (particle size 5 µm, 250
120
X 4.6 mm i.d.). The column was mounted to a HPLC unit (Shimadzu Scientific instruments,
121
Kyoto, Japan) and was maintained at 35°C. The pump flow rate was 1.5 mL/min with 100%
122
acetonitrile as the mobile phase. The injection volume was fixed to 20 µL. The chromatograph
123
was coupled to a UV-VIS diode array multiple wavelength detector. Detection and quantification
124
of squalene was carried out at 195 nm. Peak identification was achieved by comparing the
125
retention time with the authentic standards and confirmation of spectral data. The
126
chromatograms were processed using LC solution software. Quantification was accomplished
127
with the aid of standard curves calculated using five point linear regression analysis.
128
LC-MS was performed using Waters Q-Tof Ultima system equipped with an electrospray
129
ionization source (Waters, Milford, MA, USA). All the samples and standards were diluted in
130
chloroform prior sampling. The samples were injected using the autosampler (Waters 2695
131
separation module) into a semi preparative reversed-phase C-18 column (Phenomenex Kinetex,
132
Hyderabad, India) (particle size 5 µm, 250 X 4.6 mm i.d.) using isocratic elution with 100%
133
acetonitrile as a mobile phase with flow rate of 1.5 mL/min. The separated components were
134
detected using photo diode array detector (Waters 2996 PDA detector). Electro spray ionization
135
in positive mode was used for the metabolite fragmentation. The source temperature was held at
6 ACS Paragon Plus Environment
Page 7 of 39
Journal of Agricultural and Food Chemistry
136
120°C. High-resolution data was obtained in the mass range of 200-1000 Da with a scan time of
137
1 s and interscan of 0.1 s using 4 GHz TDC detector. The mass spectrum was obtained at a cone
138
voltage of 100V. For positive mode, the capillary voltage was set to 3.5 kV. The high-resolution
139
data was analyzed by MassLynx 4.0 software (Waters, Milford, MA, USA). Nitrogen was used
140
as collision gas in the collision cell. The analysis was carried out at a desolvation temperature of
141
350 oC and a cone gas flow of 50 L/h and a desolvation gas flow of 500 L/h. The Q-tof was
142
tuned in positive mode allowing the passage of [M (410.7) + H (1) + CH3CN (41)]+. The results
143
from LC-MS are provided in supporting information (Figure S1).
144 145
Gene cloning and construction of strains: Template genomic DNA from BY4741 for PCR
146
reactions was isolated using the standard glass beads method.22 PCR amplification of the tHMG1
147
gene from template was performed using oligonucleotides tHmg1-Gal-F and tHmg1-Gal-R with
148
flanking restriction sites BamHI and EcoRI (Table 1). This fragment was restriction digested and
149
ligated into pYES2/NT-C to obtain pYGH. All DNA sub cloning steps were performed with E.
150
coli DH5α using standard methods as described24 and nucleotide sequences of cloned genes were
151
verified by sequencing (ABI-310 DNA sequencer, Applied Biosystems, CA, USA). pYGH was
152
then introduced into the parent strains BY4741 to obtain SK3. PCR amplification of the tHMG1
153
gene from template was performed using oligonucleotides tHmg1-pTEF-F and tHmg1-pTEF-R
154
with flanking restriction sites NotI and SpeI (Table 1). The fragment was restriction digested and
155
ligated with pCEV-G1-km to obtain pCTH. PCR amplification of the tHMG1 gene from
156
template was performed using oligonucleotides tHmg1-pPGK-F and tHmg1-pPGK-R with
157
flanking restriction sites BamHI and SalI (Table 1). The fragment was restriction digested and
158
ligated with pCEV-G1-km to obtain pCTH. PCR amplification of the ZWF1, POS5(wM: with
7 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 8 of 39
159
Mitochondrial presequence), POS5(w/oM: without mitochondrial presequence), and tHMG1
160
gene from the genomic DNA template of BY4741 were performed using oligonucleotides pairs
161
Zwf1 – F, Zwf1 – R; Pos5 – F (wM), Pos5 – R; Pos5-F (w/oM), Pos5 – R and tHmg1-pPGK-F,
162
tHmg1-pPGK-R, respectively with flanking restriction sites BamHI and SalI (Table 1). The
163
fragments were restriction digested and ligated with pCEV-G1-km to obtain pCPZ, pCPP(wM),
164
pCPP(w/oM) and pCPH. These fragments were also inserted into pCTH to obtain pCTHPZ,
165
pCTHPP(wM), pCTHPP(w/oM), and pCTHPH, respectively. The developed plasmids were
166
introduced into BY4741 strain using Sc easy Com transformation kit (Invitrogen, Bangalore,
167
India). The plasmid maps of vectors used in the present study are provided in supplementary
168
figure (Figure S2).
169 170
Squalene extraction from S. cerevisiae: Positive transformants were picked and pre-cultured in
171
50 mL of yeast synthetic medium at 30°C, 180 rpm for 16 h. The inoculum from pre-culture was
172
used to cultivate yeast cells with an initial OD of 0.1 in production culture at 30°C, 180 rpm for
173
24 h. Cells were harvested by centrifugation to obtain cell pellets and the supernatants were
174
discarded. The cells were frozen at -20°C overnight and later subjected to lyophilisation at -50°C
175
under vacuum for 4 h. The lyophilized cells were dispersed in chloroform/ methanol (2:1, v/v)
176
mixture (~50mL). The dispersed cells were then subjected to sonication for a 3 s pulse for 10
177
cycles at 70% amplitude using sonicator (Qsonica, Newtown, CT, USA). The lysed cells were
178
incubated under shaking condition for extraction at 30˚C, 180 rpm. Subsequently, the extract was
179
filtered and subjected to flash evaporation at 50°C. The samples were dispersed in chloroform
180
solvent and subjected to high speed centrifugation to remove debris prior to RP-HPLC analysis.
181
To obtain the dry cell weight (DCW), cells harvested by centrifugation from the 50 mL culture
8 ACS Paragon Plus Environment
Page 9 of 39
Journal of Agricultural and Food Chemistry
182
were oven-dried. We obtained the following relationship between OD600 and DCW per
183
millilitre (g-DCW/mL): DCW=(OD*0.53429)-0.22216. All optical densities at 600nm (OD600)
184
measurements were taken using a Shimadzu UV-2550 spectrophotometer. All the experiments
185
were carried out in triplicates and the experimental data represents the averages of biological
186
triplicates. Standard deviations of the data were less than 7% among the replicates.
187 188
HMG-CoA reductase assay: The cells were broken by vortexing with glass beads in 100 mM
189
phosphate buffer (pH – 7.0), 1% Triton X-100 and 1 mM phenyl methyl sulphonyl fluoride
190
(PMSF) following the method described in Polakowski et al. (1998).12 The supernatant was
191
collected and the cytosolic fraction was further used for determining protein concentration and
192
HMG-CoA reductase activity. The enzyme activity was carried out using HMG-CoA reductase
193
assay kit from Sigma-Aldrich (St. Louis, MO, USA). Reactions were set using 20 µL of
194
reconstituted NADPH (to obtain a final concentration of 400 µM) and 60 µL of HMG-CoA
195
substrate solution (to obtain a final concentration of 0.3 mg/mL). The reaction was initiated
196
either by addition of HMG-CoA reductase (concentration of 0.5 – 0.7 mg/mL) provided with the
197
kit or by addition of protein extracted from the sample. The samples were mixed thoroughly
198
prior to spectrophotometric analysis. The activity was calculated according to the manufacturer’s
199
protocol and expressed as µmol of NADPH oxidized/min/mg protein. The protein content of the
200
cell extracts was calculated using Bradford assay with bovine serum albumin (BSA) as a
201
standard. Absorbance was measured at 340 nm against the blank sample using Shimadzu UV-
202
2550 UV/Vis spectrophotometer connected to a computer. The absorbance was measured every
203
15 s up to 5 min using a kinetic program. The A340 value gradually decreased over the time-
204
course due to the oxidation of NADPH.
9 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
205
Page 10 of 39
RESULTS AND DISCUSSION
206 207
Effect of promoters and copy number on HMG-CoA overexpression during squalene
208
synthesis
209
In the present study, effect of different promoters (both inducible and constitutive) on tHMG1
210
gene overexpression towards squalene synthesis has been evaluated. Towards this, three of the
211
widely used promoters PGAL1, PTEF1 and PPGK1 were selected using a 2µ episomal plasmid vector.
212
The expression level of tHMG1 gene from these plasmids has been indirectly measured by
213
quantifying squalene. Quantification was carried out with HPLC and the peak was verified using
214
LC-MS. The LC-MS results are provided in supporting information (Figure S1). The strains
215
SK3, SK10 and SK12 harboring tHMG1 gene with PGAL1, PTEF1 and PPGK1 promoters,
216
respectively, have significantly improved both squalene concentration and squalene yield (cf.
217
Figure 1A) when compared to that of base strain, BY4741. A maximum squalene concentration
218
of 6.1, 6.2 and 4.8 mg/L was observed in SK3, SK10 and SK12, respectively. A fold increase of
219
3.0, 3.1 and 2.3 in squalene concentration was observed in these strains when compared to that of
220
the base strain, substantiating the efficacy of HMG-CoA reductase on sterol pathway flux (Table
221
2). The maximum squalene yield of 6.5 (mg/g DW) was observed in SK3, where PGAL1 was used
222
as the promoter. The relative promoter efficiency ranking in terms of squalene titer observed to
223
follow PGAL1=PTEF1>PPGK1. This result is in accordance with Vickers et al. (2013)23 where they
224
have observed that expression strength of PPGK1 was poor compared to PTEF1. Partow et al.
225
(2010)25 used lacZ as reporter gene to evaluate efficiency of seven different promoters and have
226
observed that PTEF1 and PPGK1 performed similarly in glucose consumption phase and the
227
performance of PTEF1> PPGK1 during ethanol consumption phase. These results also indicates that
10 ACS Paragon Plus Environment
Page 11 of 39
Journal of Agricultural and Food Chemistry
228
the performance of constitutive promoter (PTEF1) is as good as strong inducible promoter (PGAL1).
229
The main advantage of having constitutive promoter such as PTEF1, is that the carbon source
230
could be a simple sugar such as glucose instead of expensive galactose for large scale
231
cultivations.
232 233
Polakowski et al. (1998)12 and several other researchers have established that truncated HMG-
234
CoA (tHMG1) overexpression improved the flux towards the ergosterol synthesis in S.
235
cerevisiae. Since then overexpression of tHMG1 has more or less become a ubiquitous target for
236
producing terpene compounds from S. cerevisiae. In order to assay whether the overexpression
237
of single copy of tHMG1 (in an episomal plasmid) is optimal or sub-optimal, two copies of this
238
gene was cloned into pCEV-G1-Km to generate pCTHPH for expression studies. Therefore,
239
strain SK16 has been generated which harbors tHMG1 under PTEF1 and PPGK1 promoters (cf.
240
Figure 1). From Figure (1A) it can be observed that there is a significant increase in the squalene
241
concentration as well as squalene yield in SK16 strain, suggesting enhanced activity of Hmg1p.
242
A maximum squalene titer of 28.1 mg/L and yield of 35.7 (mg/g DW) has been achieved with
243
SK16 strain. There was 4.5- and 14-fold increase in squalene concentration in SK16 compared to
244
that of SK10 and BY4741 strains, respectively (Table 2). Prior investigations have overexpressed
245
tHMG1 gene either using episomal plasmids26 or by integrating this gene at a specific loci in the
246
chromosome.27,28 Generally, S. cerevisiae harboring an episomal plasmid will usually have 10-40
247
plasmid copies per cell.29 In case of commercial cultivations, expression of genes on plasmid is
248
not preferred due to instability issues. Hence, key target genes are integrated into chromosome.
249
For instance, tHMG1 gene has been integrated at three different loci to enhance the rate-limiting
250
step in the terpene biosynthesis pathway27. In the current study, the strain harboring two copies
11 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 12 of 39
251
of tHMG1 gene has led to a 4.5-fold increase in squalene yield over the strains harboring a single
252
copy of the same gene. This suggests that even an episomal plasmid carrying a single copy of
253
tHMG1 is sub-optimal for increasing the flux towards squalene or any other terpene precursor.
254
Hence there is a scope for further improvement of terpene synthesis in S. cerevisiae if copy
255
number of tHMG1 gene is targeted either by multiple integrations or by expressing multiple
256
copies on an episomal plasmid to achieve results closer to optimum. This result also shows that
257
the overexpression of HMG-CoA reductase using two copies scavenges the NADPH in cytosol,
258
and might lower the fluxes of other anabolic reactions which require NADPH as cofactor. It has
259
been proved that increasing copy number of the episomal plasmid will increase the protein levels
260
translated from it.30 Likewise, it can be concluded that expression of the same gene as multiple
261
copies in a single cassette could also increase the protein levels.
262 263
Effect of ZWF1 and POS5 overexpression for cofactor improvement on squalene synthesis
264
The iMM904 genome-scale metabolic model of S. cerevisiae has been screened for NADPH
265
synthesizing reactions which led to the identification of 18 enzymes carrying irreversible
266
reactions for NADPH regeneration and 7 reactions involved in reversible reactions (Table S1).31
267
Among the 18 enzymes, one is present in peroxisome, six in mitochondria and eleven enzymes
268
are present in cytoplasm. The model iMM904 has 1575 reactions among which 91 reactions
269
involve NADPH. Out of 91 reactions, 18 reactions (irreversible) produce NADPH while 75
270
reactions consume NADPH which implies that NADPH is one of the important cofactor in the
271
yeast cells (cf. Table S1). Enzymes such as cytosolic Zwf1p, Ald6p and Idp2p contribute to the
272
regeneration of NADPH from NADP+ in the cytosol of S. cerevisiae (cf. Table S1).32 The
273
NADP+ required for these dehydrogenases is supplied by two key NAD kinases, Yef1p and
12 ACS Paragon Plus Environment
Page 13 of 39
Journal of Agricultural and Food Chemistry
274
Utr1p in the cytosol. In case of mitochondria, the NADPH is regenerated using Pos5p which
275
executes both NADH kinase and NAD kinase activities, however the former activity was found
276
to be higher.32,33 Hence, in the present study ZWF1 and POS5 which regenerate NADPH pool in
277
cytoplasm and mitochondria, respectively, have been overexpressed to see the effect on squalene
278
synthesis. ZWF1 and POS5 (with mitochondrial presequence) and POS5 (without mitochondrial
279
presequence) were overexpressed in BY4741 to generate strains SK13, SK14 and SK15,
280
respectively. Significant increase in squalene concentration and squalene yield was observed
281
when compared to BY4741 (Figure 2A). The squalene yield obtained from SK13, SK14 and
282
SK15 were found to be 14.9, 9.0 and 7.5 (mg/g DW), respectively (Table 2). Owing to the stable
283
expression of the HMG-CoA reductase under PTEF1 promoter, the PTEF1-tHMG1 cassette was
284
combined with PPGK1- ZWF1 and PPGK1-POS5 gene cassettes to form the strains, SK17, SK18
285
and SK19 respectively. When these strains were analysed for squalene improvement, squalene
286
yield of 47.4, 58.6 and 33 mg/g DCW was observed in SK17, SK18 and SK19 strains,
287
respectively (Figure 2B). There was 2 and 3.7-folds increase in squalene concentration in SK17
288
and SK18 strains compared to that of SK13 and SK14, respectively. The increase in
289
concentration and yield with respect to BY4741 was 12.2 and 22.2 folds in case of SK17 strain,
290
and 14 and 27.4 folds in case of SK18 strain (Table 2). Interestingly, the dry cell weight (DCW)
291
of SK16, SK17 and SK18 were observed to be lower than all other strains (cf. Table 2). This
292
suggests that the flux towards biomass formation is lowered whenever squalene flux was
293
increased. Moreover, Donald et al. (1997)34 also suggests that increased squalene or
294
farnesyldiphosphate (a presqualene intermediate) due to overexpression of tHMG1 could be
295
cytotoxic to the cell and hence the reduced biomass.
13 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 14 of 39
296
Naziri et al. (2011)35 observed maximum squalene yield of 10.02 mg/g DCW using terbinafine
297
and methyl jasmonate as selection pressure during cultivation of wild-type S. cerevisiae.
298
Whereas, Mantzouridou and Tsimidou (2010)36 achieved maximum yield of 18.5 mg/g DCW
299
when HMG2 was overexpressed along with K6R mutation in Hmg2p.
300
resulted in lesser squalene yield than the present study. Zhuang and Chappell (2015)37 reported a
301
maximum squalene titer of 250 mg/L in ZXB yeast strain at end of 12 days, which is higher than
302
maximum titer obtained in present study (46 mg/L in SK19 at the end of 24 h). ZXB was derived
303
from parental BY4741 strain having mutations for sterol uptake enhancement (EMS treatment
304
and selection), overexpression of hamster tHMG1, deletion of ERG1 and ERG9, and finally
305
heterologous expression of squalene synthase from Botryococcus braunii.
306
NADPH is tightly regulated in the cell in comparison to NADH.38 Studies have been conducted
307
to manipulate the levels of NADP/NADPH ratio in the cytosol.39 Reducing equivalents in form
308
of NADPH are essential for several metabolic reactions involved in the biosynthesis of cellular
309
macromolecules.40 NADPH is a major cofactor involved in ergosterol biosynthesis and is
310
required by HMG-CoA reductase (HMG1) and squalene synthase (ERG9) for squalene
311
synthesis. Increased squalene content in SK-13, SK-14 and SK-15 can be correlated to the
312
increased NADPH regeneration by overexpressing ZWF1 and POS5 genes. When
313
overexpression of these NADPH regenerating genes was coupled to tHMG1 overexpression, as
314
seen in strains SK-17, SK-18 and SK-19 there was a synergistic effect on squalene synthesis.
315
This indicates that either of Hmg1p or Erg9p, or both may be are limited by NADPH supply and
316
hence improving regeneration of NADPH significantly improved squalene yield. Overexpression
317
of POS5 gene was more effective during such coupled co-expression and has shown 1.2 fold
318
higher squalene yield when compared to that of ZWF1 overexpression. This could be due to the
Both these studies
14 ACS Paragon Plus Environment
Page 15 of 39
Journal of Agricultural and Food Chemistry
319
fact that NADPH regeneration is higher in POS5 overexpressed strain owing to the dual
320
functionality of Pos5p as a NAD+ and NADH kinase. The mechanism through which it improves
321
squalene flux has been discussed subsequently.
322
Overexpression of ZWF1 was proven to improve lysine production in C. glutamicum. 41 In case
323
of S. cerevisiae, ZWF1 overexpression has improved tolerance to furfural, which could be
324
exploited for lignocellulose- based bioethanol production.42 There are very few reports on its
325
overexpression for terpene production in S. cerevisiae. Brown et al. (2015)43 has overexpressed
326
ZWF1 for strictosidine production and Zhao et al. (2015)21 studied its effect on carotene
327
production. Kwon et al. (2006)44 observed there was 6.0 fold increase in Zwf1p activity in S.
328
cerevisiae strain when ZWF1 was overexpressed. However NADPH levels were not elevated
329
indicating rapid consumption of this co-factor. In their study the regenerated NADPH was
330
utilized by NADPH- dependent xylose reductase for converting xylose to xylitol as there was
331
significant improvement in xylitol (86 g/L) when compared to control (71 g/L). In contrast to
332
Zwf1p, Pos5p is localized in mitochondrial matrix and plays a key role in protection of
333
mitochondria from oxidative stress.45 Hou et al. (2009)46 overexpressed both full length Pos5 and
334
Pos5∆17 (lacking the mitochondrial target signal) in S. cerevisiae to generate cNDK (equivalent
335
to SK15 and SK19 in present study) and mNDK (equivalent to SK14 and SK18 in present study)
336
strains. They observed improved NADH kinase activity in mNDK and cNDK strains when
337
compared to reference strain and no significant changes in NAD kinase activity when cultivated
338
aerobically using glucose as carbon source.
339
The rationale behind improved squalene due to POS5 overexpression is not straightforward as in
340
case of ZWF1. The mitochondrial membrane is impermeable to NADH and NADPH. Hence the
341
cofactors are restricted to their respective organelles.33 The utr1yef1 double mutant S. cerevisiae
15 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 16 of 39
342
strain showed viability implying that cytosolic NADP+ can also be supplied by Pos5p.32 Thus
343
supplied NADP+ can be reduced to NADPH by cytosolic dehydrogenases (Zwf1p, Ald6p and
344
Idp2p). However, the mechanism of NADP+ transport between cytosol and mitochondria has not
345
been well characterized and the function of several mitochondrial carrier proteins are yet to be
346
uncovered42. Furthermore, the mitochondrial transport protein Yhm2p which is involved in
347
citrate transport to cytosol has also enhanced the NADPH in cytosol via citrate-oxoglutarate
348
NADPH redox shuttle between mitochondria and cytosol.45 These two findings suggest that the
349
regeneration of NADPH (via NADH kinase) and NADP+ (via NAD kinase) inside the
350
mitochondria could possibly enhance the NADPH pool in the cytosol. Thus improved NADPH
351
regeneration in the cytosol in turn could have possibly increased the rates of HMG-CoA
352
reductase (Hmg1p) and squalene synthase (erg9p) thereby increasing the squalene flux. This
353
could be the plausible rationale for improved squalene during overexpression of mitochondrial
354
POS5 either in isolation or when coupled with tHMG1. Zhao et al. (2015)21 overexpressed ZWF1
355
and POS5 to observe improvement in carotenoid biosynthesis, however, the POS5 gene has been
356
overexpressed without the mitochondrial presequence. Hence the POS5 gene product was
357
restricted only to cytoplasm. Moreover the fold increase in the β-carotene yield was only 1.8 and
358
in lycopene it was 1.6 when compared to that of the control strain, however significant. NAD
359
kinase activity was found to be significantly improved in the strain overexpressing POS5 in their
360
study. However, NADH kinase activity was not assayed which is supposedly seen more when
361
compared to NAD kinase activity of Pos5p. Shi et al. (2013)47 overexpressed POS5 and ZWF1
362
from S. cerevisiae in C. glutamicum spp. lactofermentum to observe increased isoleucine
363
production. They have extended this study to proteomic analysis in Shi et al. (2015)48 to observe
364
24 differentially expressed genes in case of POS5 overexpression and 7 differentially expressed
16 ACS Paragon Plus Environment
Page 17 of 39
Journal of Agricultural and Food Chemistry
365
in case of ZWF-ppnk overexpression. Prior investigations were focused on the expression of
366
truncated POS5 gene and its influence in biochemical production.21,48 Strand et al. (2003)19 has
367
reported that even the full length POS5 overexpression has yielded much higher protein when
368
compared to truncated POS5 and the protein from full length expression has accumulated only in
369
mitochondria. To the best of our knowledge this is the first study where we report full length
370
POS5 overexpression has significantly improved squalene in S. cerevisiae. If we carefully
371
observe, the squalene levels were higher in the case of strains harboring plasmid containing both
372
HMG1 and Zwf1/Pos5, i.e. SK17, SK18 & SK19 in comparison to only Zwf1/Pos5, i.e. SK13,
373
SK14 & SK15. This indicates that overexpression of HMG1 has served as an additional sink for
374
the regenerated NADPH. Similar result was observed by Kwon et al. (2006)44 where expression
375
of NADPH- requiring xylose reductase has improved xylitol production when NADPH was
376
regenerated using overexpression of Zwf1. Hou et al. (2009)46 observes that NADPH
377
regeneration using oxidative pentose phosphate pathway is coupled with CO2 production
378
whereas NADPH regeneration via POS5 is a decoupled process. Probably this could be the
379
reason for higher concentration of squalene in SK18 and SK19 than SK17 (Table 2). From the
380
biomass data it observed that the strains SK14 and SK18 had lower dry cell values (DCW) in
381
comparison to strains SK15 and SK19. Especially, SK18 had the lowest biomass (0.48 g/L)
382
where maximum squalene yield was observed. Here we hypothesize that in strains SK15 and
383
SK19, the regenerated NADPH was involved in the biosynthesis of cellular macromolecules
384
thereby contributing to biomass. Moreover as mentioned in Hou et al. (2009)46, as the redox
385
cofactors are impermeable across mitochondrial membrane, the phenotypes could vary when
386
cytosolic and mitochondrial redox are perturbed.
387
17 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 18 of 39
388
In-vitro HMG-CoA reductase assay
389
HMG1 contributes about 83% of the HMG-CoA reductase activity in wild-type cells while its
390
isomer, HMG2, contributes the remaining activity thus indicating the significance of HMG1.10
391
HMG-CoA is transported to endoplasmic reticulum where it is reduced to mevalonate by
392
membrane bound glycoprotein HMG-CoA reductase. The reductive cleavage of HMG-CoA to
393
mevalonate and CoA requires two molecules of NADPH. HMG-CoA reductase activity was
394
determined based on the reduction in absorbance at 340 nm which in turn is due to the NADPH
395
oxidation by the catalytic subunit of HMG-CoA reductase in the presence of the substrate HMG-
396
CoA. In the present study, the strains conferring HMG-CoA reductase activity was ascertained
397
based on the estimation of NADPH changes as described earlier in methods section. When
398
compared to BY4741, the crude cell extracts of all the transformants displayed significantly
399
higher enzyme activity (cf. Figure 3). Recombinant strains, SK10 and SK12 attributed similar
400
activity indicating that TEF1 and PGK1 promoter performance was almost same. Whereas SK3
401
strain, which harbors tHMG1 under PGAL1 showed slightly higher activity when compared to that
402
of SK10 and SK12. Maximum increase in the reductase activity was observed in the strains
403
SK17 and SK18 with fold increase of 10 and 11, respectively, when compared to that of BY4741
404
(Table 2). Strains SK17 and SK18 harbors tADH1-tHMG1-PTEF1-PPGK1-ZWF1-tCYC1 and
405
tADH1-tHMG1-PTEF1-PPGK1-POS5-tCYC1 cassettes, respectively, indicating the tHMG1
406
overexpression and NADPH regenerating ZWF1/POS5 improves the HMG-CoA reductase
407
activity significantly in comparison to other recombinant strains. Interestingly, the SK16 strain
408
harboring tADH1-tHMG1-PTEF1-PPGK1-tHMG1-tCYC1 increased the reductase activity by 1.4-
409
fold in comparison to SK10/SK12 and 1.2- fold in comparison to SK3 indicating that presence of
410
two tHMG1 gene copies on a single cassette will improve its activity when compared to single
18 ACS Paragon Plus Environment
Page 19 of 39
Journal of Agricultural and Food Chemistry
411
copy. Hence the strain SK16 produced more squalene in comparison to SK3, SK10 and SK12.
412
Interestingly, ZWF1 (strain SK13) and POS5 (strain SK14) overexpression improved the HMG-
413
CoA reductase activity by 6.5- and 5.5- folds in comparison to BY4741 suggesting the HMG1
414
up-regulation. Moreover, this fold increase is alike to that of strains SK10 and SK12 indicating
415
same effect as that of tHMG1 overexpression. Thus the enzyme activity results corroborates the
416
squalene estimation results found in the earlier section. Polakowski et al. (1998)12 has estimated
417
the HMG-CoA reductase activity in the tHMG1 overexpressed strains. In their study it was
418
observed that the episomal plasmid- based overexpression had significantly higher activity in
419
comparison to the single copy integration of tHMG1. In case of episomal expression the specific
420
activity obtained in their study was 0.18 U/mg protein whereas in the present study it ranged
421
between 0.13-0.15 U/mg protein. The specific activity values obtained in the present study is in
422
agreement with Polakowski et al. (1998)12 for both wild-type and tHMG1 overexpressed strains.
423
A maximum specific activity of 0.218 U/mg protein was obtained in the strain SK18 harboring
424
tADH1-tHMG1-PTEF1-PPGK1-POS5-tCYC1 cassette.
425 426
Effect of nitrogen sources on G418 – based vector for tHMG1 overexpression
427
Aminoglycoside 3′-phosphotransferase coded by aphA1 gene from E. coli confers
428
aminoglycoside, G418 (Geneticin®) resistance in S. cerevisiae and is a commonly implemented
429
selection marker during transformation in S. cerevisiae.23,49,50 The kanMX module on a G418
430
selection cassette contains the aphA1 gene under a specific promoter and terminator. The
431
plasmid pCEV-G1-Km used in the present study harbors kanMX module for selection of
432
transformants using G418 antibiotic. It has been established that if ammonium sulphate is used
433
as nitrogen source in the growth medium for cultivation of S. cerevisiae strain harboring kanMX
19 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 20 of 39
434
cassette, the resistance towards G418 is lost.23,50 The exact reason for such loss is not established,
435
however, it was suggested that the ammonium sulphate assimilation might downregulate the
436
reactions involved in import and activity of G418.50 So when selective pressure is lost during
437
successive generations of cell cycle, the percentage of cells harboring plasmid will be lowered as
438
the cells which are free of plasmids will have higher growth rate due to selective advantage. This
439
in turn lowers the expression of genes present on the plasmid. The very role of selective pressure
440
(addition of G418 to the medium for present study) is to maintain the plasmid in cell population.
441
It was also established that if ammonium sulphate is substituted by glutamate (1 g/L) as the
442
nitrogen the resistance to G418 was retained and therefore the selection pressure was
443
maintained.23,50 To corroborate this reasoning experimentally, it was observed that in the strain
444
SK10 squalene content was meagre 1.42 mg/L when (NH4)2SO4 was used as nitrogen source
445
whereas substituting with glutamate the squalene synthesis was improved to 6.3 mg/L (cf. Figure
446
4A). Squalene yield and titers improved significantly indicating the positive effect of glutamate
447
as nitrogen source on G418- based gene expression (cf. Figure 4A). Hence, the SD medium used
448
for cultivation has been modified by replacing (NH4)2SO4 with glutamic acid (1 g/L) for the
449
strains harboring plasmids with geneticin as a selection marker for all shake-flask experiments.
450
Vickers et al. (2013)23 provided comprehensive information on factors affecting the G418
451
selection, and medium pH and its composition were suggested to be key factors.
452
Similarly the carbon sources, glucose and galactose, have been evaluated for squalene production
453
when constitutive promoter PTEF1 was utilized. When strain SK10 was cultivated on 2% glucose
454
or 2% galactose using glutamic acid (1 g/L) as the nitrogen source, the squalene production
455
levels were slightly higher in case of 2% glucose when compared to 2% galactose (Figure 4B).
456
Partow et al. (2010) reported that TEF1 is a truly constitutive promoter is not affected by change
20 ACS Paragon Plus Environment
Page 21 of 39
Journal of Agricultural and Food Chemistry
457
in the carbon source. Therefore, the results obtained are also in accordance with Partow et al.
458
(2010).25
459 460
In conclusion, we demonstrate that single copy episomal overexpression of tHMG1 is suboptimal
461
for producing squalene in S. cerevisiae. We also show that NADPH regenerating enzymes can
462
improve squalene synthesis in S. cerevisiae when overexpressed either in isolation or in tandem
463
with tHMG1. In case of latter, the effect was highly synergistic. We also observed that full length
464
POS5 overexpression improved squalene in cytosol even though its protein, pos5p, localizes in
465
mitochondria. Although further studies are needed to elucidate the exact reason for this
466
observation. Finally we conclude that the strategies employed in present study could be
467
potentially applied for engineering S. cerevisiae for squalene improvement and also for other
468
industrially relevant terpenes.
469
470
471
AUTHOR INFORMATION
472
*Corresponding author
473
Microbiology and Fermentation Technology Department, CSIR-Central Food Technological
474
Research Institute, Mysore, India.
475
E-mail:
[email protected];
[email protected] 476
Telephone: +91-821-2517539
477
ORCID: 0000-0002-3251-8643
21 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 22 of 39
478
Author contributions:
479
Funding: The financial support was provided by the Department of Biotechnology, India (award
480
of research fellowship) and by the Science Engineering and Research Board (SERB), India
481
(YSS/2014/000565).
482
Notes:
483
The authors declare no competing financial interest
484
485
ACKNOWLEDGEMENTS
486
The authors thank Mr. Sandeep Kumar, RGNF-JRF for his useful comments and suggestions on
487
cloning experiments and Director Prof. Ram Rajasekharan (Lipid Science Dept., CSIR-CFTRI)
488
for providing us the BY4741 strain and the plasmid, pYES2/NTC. The authors also acknowledge
489
Mr. P Mukund Lakman, Central Instrumentation Facilities, CSIR-CFTRI for and HPLC and LC-
490
MS studies.
491
492
Supporting information
493
Figure S1: Shows the LC-MS profile of squalene against a standard
494
Figure S2: Shows the plasmid maps constructed for the current study
495
Table S1: Provides the list of reactions in Saccharomyces cerevisiae involved in NADPH
496
synthesis based on iMM904 in-silico model
497
22 ACS Paragon Plus Environment
Page 23 of 39
Journal of Agricultural and Food Chemistry
498
REFERENCES
499
(1)
500
Horning, T.; Tsuruta, H.; Melis, D. J.; Owens, A.; Fickes, S.; Diola, D.; Benjamin, K. R.;
501
Keasling, J. D.; Leavell, M. D.; McPhee, D. J.; Renninger, N. S.; Newman, J. D.; Paddon, C. J.,
502
Production of amorphadiene in yeast, and its conversion to dihydroartemisinic acid, precursor to
503
the antimalarial agent artemisinin. Proc. Natl. Acad. Sci. U S A. 2012, 109, E111-8.
504
(2)
505
terpenoids. Comput. Struct. Biotechnol. J. 2012, 3, e201210006.
506
(3)
507
Makris, A. M.; Roussis, V.; Kampranis, S. C., Reconstructing the chemical diversity of labdane-
508
type diterpene biosynthesis in yeast. Metab. Eng. 2015, 28, 91-103.
509
(4)
510
Ho, K. a.; Eachus, R. a.; Ham, T. S.; Kirby, J.; Chang, M. C. Y.; Withers, S. T.; Shiba, Y.;
511
Sarpong, R.; Keasling, J. D., Production of the antimalarial drug precursor artemisinic acid in
512
engineered yeast. Nature. 2006, 440, 940-3.
513
(5)
514
of Saccharomyces cerevisiae. Crit. Rev. Biotechnol., 2017,
515
DOI:10.1080/07388551.2017.1299679
516
(6)
517
potential of biotechnology. Lipid. Technol. 2011, 23, 270-273.
518
(7)
519
(Accessed on 24th August 2017)
Westfall, P. J.; Pitera, D. J.; Lenihan, J. R.; Eng, D.; Woolard, F. X.; Regentin, R.;
Kampranis, S. C.; Makris, A. M., Developing a yeast cell factory for the production of
Ignea, C.; Ioannou, E.; Georgantea, P.; Loupassaki, S.; Trikka, F. A.; Kanellis, A. K.;
Ro, D.-K.; Paradise, E. M.; Ouellet, M.; Fisher, K. J.; Newman, K. L.; Ndungu, J. M.;
Paramasivan, K.; Mutturi, S., Progress in terpene synthesis strategies through engineering
Naziri, E.; Mantzouridou, F.; Tsimidou, M. Z., Squalene resources and uses point to the
http://www.marketsandmarkets.com/Market-Reports/squalene-market-542345.html
23 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 24 of 39
520
(8)
Spanova, M.; Daum, G., Squalene - biochemistry, molecular biology, process
521
biotechnology, and applications. Eur. J. Lipid. Sci. Tech. 2011, 113, 1299-1320.
522
(9) Rodwell, V. W.; Nordstrom, J. L.; Mitschelen, J. J., Regulation of HMG-CoA reductase. Adv.
523
Lipid Res. 1976, 14, 1-74.
524
(10) Basson, M. E.; Thorsness, M.; Rine, J., Saccharomyces cerevisiae contains two functional
525
genes encoding 3-hydroxy-3-methylglutaryl-coenzyme A reductase. Proc. Natl. Acad.
526
Sci. U.S.A. 1986, 83, 5563-67.
527
(11)
528
A proliferation of stacked membrane pairs surrounding the yeast nucleus. J. Cell Biol. 1988,
529
107, 101-114.
530
(12)
531
glutaryl-CoA reductase leads to squalene accumulation in yeast. Appl. Environ. Microbiol.1998,
532
49, 66-71.
533
(13)
534
cerevisiae. i-ACES 2014, 1, 57-63.
535
(14)
536
gene encoding HMG-CoA reductase in Saccharomyces cerevisiae for production of prenyl
537
alcohols. Appl. Microbiol. Biotechnol. 2009, 82, 837–845.
538
(15)
539
precursor supply in Saccharomyces cerevisiae. J. Biotechnol. 2013, 168, 446–451.
540
(16)
541
of NADP. Redox Rep. 2010, 15, 2-10.
Wright, R.; Basson, M., Increased amounts of HMG-CoA reductase induce "Karmellae":
Polakowski, T.; Stahl, U.; Lang, C., Overexpression of a cytosolic hydroxyl methyl
Thompson, A.; Kwak, S.; Jin, Y.-s., Squalene Production using Saccharomyces
Ohto, C.; Muramatsu, M.; Obata, S.; Sakuradani, E.; Shimizu, S. Overexpression of the
Liu, J.; Zhang, W.; Du, G.; Chen, J.; Zhou, J. Overproduction of geraniol by enhanced
Agledal, L.; Niere, M.; Ziegler, M., The phosphate makes a difference: cellular functions
24 ACS Paragon Plus Environment
Page 25 of 39
Journal of Agricultural and Food Chemistry
542
(17)
Shi, F.; Kawai, S.; Mori, S.; Kono, E.; Murata, K., Identification of ATP-NADH kinase
543
isozymes and their contribution to supply of NADP(H) in Saccharomyces cerevisiae. FEBS. J.
544
2005, 272, 3337-3349.
545
(18)
546
in the respiratory chain activity of Saccharomyces cerevisiae. Acta. Biochim. Biophys. Sin. 2011,
547
43, 989-995.
548
(19)
549
Copeland, W. C., POS5 gene of Saccharomyces cerevisiae encodes a mitochondrial NADH
550
kinase required for stability of mitochondrial DNA. Eukaryot. Cell. 2003, 2, 809-820.
551
(20)
552
flux engineering of L-lysine production in Corynebacterium glutamicum: over expression and
553
modification of G6P dehydrogenase. J Biotechnol. 2007, 132, 99–109.
554
(21)
555
biosynthesis in recombinant Saccharomyces cerevisiae. Lett. Appl. Microbiol. 2015, 61, 354-
556
360.
557
(22)
558
Cold Spring Harbor Laboratory Course Manual. Cold Spring Harbor Laboratory Press, Cold
559
Spring Harbor, NY 2005.
560
(23)
561
vectors with antibiotic selection markers for engineering in Saccharomyces cerevisiae. Microb.
562
Cell. Fact. 2013, 12, 96.
563
(24)
564
Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY 2001.
Shi, F.; Li, Z.; Sun, M.; Li, Y., Role of mitochondrial NADH kinase and NADPH supply
Strand, M. K.; Stuart, G. R.; Longley, M. J.; Graziewicz, M. A.; Dominick, O. C.;
Judith, B.; Corinna, K.; Andrea, H.; Oskar, Z.; Christoph, J.B.; Christoph, W., Metabolic
Zhao, X.; Shi, F.; Zhan, W., Overexpression of ZWF1 and POS5 improves carotenoid
Amberg, D. C.; Burke, D. J.; Dawson, D.; Stearns, T., Methods in Yeast Genetics: A
Vickers, C. E.; Bydder, S. F.; Zhou, Y.; Nielsen, L. K., Dual gene expression cassette
Sambrook, J., and D. E. Russell., Molecular cloning: a laboratory manual, No. Ed. 2.
25 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 26 of 39
565
(25)
Partow, S.; Siewers, V.; Bjørn, S.; Nielsen, J.; Maury, J., Characterization of different
566
promoters for designing a new expression vector in Saccharomyces cerevisiae. Yeast 2010, 27,
567
955-64.
568
(26)
569
cerevisiae by metabolic engineering. FEMS. Microbiol. Lett. 2013, 345, 94-101.
570
(27)
571
diphosphate pool as direct precursor of sesquiterpenes through metabolic engineering of the
572
mevalonate pathway in Saccharomyces cerevisiae. Biotechnol. Bioeng. 2010, 106, 86-96.
573
(28)
574
engineering applications in Saccharomyces cerevisiae. FEMS. Yeast. Res. 2012, 12, 197-214.
575
(29)
576
Yeast. 1992, 8, 423–488.
577
(30)
578
number of episomal plasmids in Saccharomyces cerevisiae for improved protein production.
579
FEMS Yeast Res. 2012, 12, 598-607.
580
(31) Mo, M. L.; Palsson, B. Ø.; Herrgård, M. J. Connecting extracellular metabolomic
581
measurements to intracellular flux states in yeast. BMC Syst. Biol. 2009, 3, 37.
582
(32)
583
Saccharomyces cerevisiae. J. Biol. Chem. 2009, 284, 7553-7560.
584
(33)
585
NADPH in Saccharomyces cerevisiae. EMBO. J. 2003, 22, 2015-2024.
Li, Q.; Sun, Z.; Li, J.; Zhang, Y., Enhancing beta-carotene production in Saccharomyces
Asadollahi, M. A.; Maury, J.; Schalk, M.; Clark, A.; Nielsen, J., Enhancement of farnesyl
Da Silva, N. A., and Srikrishnan, S., Introduction and expression of genes for metabolic
Romanos, M.A.; Scorer, C.A.; Clare, J.J., Foreign gene expression in yeast – a review.
Chen, Y., Partow, S., Scalcinati, G., Siewers, V., & Nielsen, J. Enhancing the copy
Miyagi, H.; Kawai, S.; Murata, K., Two sources of mitochondrial NADPH in the yeast
Outten, C. E.; Culostta, V. C., A novel NADH kinase is the mitochondrial source of
26 ACS Paragon Plus Environment
Page 27 of 39
Journal of Agricultural and Food Chemistry
586
(34) Donald, K. A. G.; Hampton, R. Y.; Fritz, I. B. Effects of overproduction of the catalytic
587
domain of 3-hydroxy-3- methylglutaryl coenzyme A reductase on squalene synthesis in
588
Saccharomyces cerevisiae. Appl. Environ. Microbiol. 1997, 63, 3341–3344.
589
(35) Naziri, E.; Mantzouridou, F.; Tsimidou, M. Z. Enhanced squalene production by wild-type
590
Saccharomyces cerevisiae strains using safe chemical means. J. Agric. Food Chem. 2011, 59,
591
9980–9989.
592
(36) Mantzouridou, F.; Tsimidou, M. Z. Observations on squalene accumulation in
593
Saccharomyces cerevisiae due to the manipulation of HMG2 and ERG6. FEMS Yeast Res. 2010,
594
10, 699–707.
595
(37) Zhuang, X.; Chappell, J. Building terpene production platforms in yeast. Biotechnol.
596
Bioeng. 2015, 112, 1854–1864.
597
(38) Moreira Dos Santos, M.; Raghevendran, V.; Kötter, P.; Olsson, L.; Nielsen, J. Manipulation
598
of malic enzyme in Saccharomyces cerevisiae for increasing NADPH production capacity
599
aerobically in different cellular compartments. Metab. Eng. 2004, 6, 352–363.
600
(39) Asadollahi, M. A.; Maury, J.; Patil, K. R.; Schalk, M.; Clark, A.; Nielsen, J., Enhancing
601
sesquiterpene production in Saccharomyces cerevisiae through in silico driven metabolic
602
engineering. Metab. Eng. 2009, 11, 328-34.
603
(40) Grabowska, D.; Chelstowska, A. The ALD6 gene product is indispensable for providing
604
NADPH in yeast cells lacking glucose-6-phosphate dehydrogenase activity. J. Biol. Chem. 2003,
605
278, 13984–13988.
606
(41) Becker, J.; Klopprogge, C.; Herold, A.; Zelder, O.; Bolten, C. J.; Wittmann, C., Metabolic
607
flux engineering of L-lysine production in Corynebacterium glutamicum—over expression and
608
modification of G6P dehydrogenase. J. Biotechnol. 2007, 132, 99-109.
27 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 28 of 39
609
(42) Gorsich, S.; Dien, B.; Nichols, N.; Slininger, P.; Liu, Z.; Skory, C., Tolerance to furfural-
610
induced stress is associated with pentose phosphate pathway genes ZWF1, GND1, RPE1, and
611
TKL1 in Saccharomyces cerevisiae. Appl. Environ. Microbiol. 2006, 71, 339-349.
612
(43) Brown, S.; Clastre, M.; Courdavault, V.; O’Connor, S. E., De novo production of the plant-
613
derived alkaloid strictosidine in yeast. Proc. Natl. Acad. Sci. U.S.A. 2015, 3205-10.
614
(44) Kwon, D.-H.; Kim, M.-D.; Lee, T.-H.; Oh, Y.-J.; Ryu, Y.-W.; Seo, J.-H. Elevation of
615
glucose 6-phosphate dehydrogenase activity increases xylitol production in recombinant
616
Saccharomyces cerevisiae. J. Mol. Catal. B Enzym. 2006, 43, 86–89.
617
(45) Castegna, A.; Scarcia, P.; Agrimi, G.; Palmieri, L.; Rottensteiner, H.; Spera, I.; Germinario,
618
L.; Palmieri, F., Identification and functional characterization of a novel mitochondrial carrier for
619
citrate and oxoglutarate in Saccharomyces cerevisiae. J. Biol. Chem. 2010, 285, 17359-17370.
620
(46) Hou, J.; Vemuri, G. N.; Bao, X.; Olsson, L. Impact of overexpressing NADH kinase on
621
glucose and xylose metabolism in recombinant xylose-utilizing Saccharomyces cerevisiae. Appl.
622
Microbiol. Biotechnol. 2009, 909–919.
623
(47) Shi, F.; Li, K.; Huan, X.; Wang, X., Expression of NAD(H) kinase and glucose-6-phosphate
624
dehydrogenase improve NADPH supply and L-isoleucine biosynthesis in Corynebacterium
625
glutamicum ssp. lactofermentum. Appl. Biochem. Biotechnol. 2013, 171, 504-521.
626
(48) Shi, F.; Li, K.; Li, Y., Comparative proteome analysis of global effect of POS5 and zwf-
627
ppnK
628
lactofermentum. Biotechnol. Lett. 2015, 37, 1063-1071.
629
(49) Webster, T. D.; Dickson, R. C., Direct selection of Saccharomyces cerevisiae resistant to
630
the antibiotic G418 following transformation with a DNA vector carrying the kanamycin-
631
resistance gene of Tn903. Gene. 1983, 26, 243-252.
overexpression
in
l-isoleucine
producing
Corynebacterium
glutamicum
ssp.
28 ACS Paragon Plus Environment
Page 29 of 39
Journal of Agricultural and Food Chemistry
632
(50) Cheng, T.H.; Chang, C.-R.; Joy, P.; Yablok, S.; Gartenberg, M. R., Controlling gene
633
expression in yeast by inducible site-specific recombination. Nucleic. Acids. Res. 2000, 28, e108-
634
e108.
635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654
29 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
655
Page 30 of 39
Figure Legends
656 657
Figure 1: Squalene synthesis in tHMG1 overexpressed S. cerevisiae strains. (A) The comparison
658
of BY4741 with respect to SK13, SK10 and SK12 represents the effect of single copy
659
overexpression of tHMG1 using promoters PGAL1, PTEF1 and PPGK1, respectively. Whereas, the
660
comparison of SK16 to other engineered strains (SK10, SK12 and SK13) reveals the effect of
661
plasmid harboring two copies tHMG1 to single copy tHMG1 on squalene. (B) Pictorial
662
representation of gene cassettes used in the strains.
663
Figure 2: Squalene synthesis in S. cerevisiae strains with overexpression of NADPH
664
regenerating genes (A) Comparison among strains SK13, SK14 and SK15 indicates the effect of
665
overexpression of ZWF1, POS5 (with mitochondrial presequence) and POS5 (without
666
mitochondrial presequence), respectively. (B) Comparison between the strains SK17, SK18 and
667
SK19 indicates the synergistic effect of simultaneous expression of tHMG1 and NADPH
668
regenerating genes. Here SK17, SK18 and SK19 represents SK13, SK14 and SK15 strains,
669
respectively, with additional tHMG1 coexpression.
670
Figure 3: Comparison of HMG-CoA reductase activity in various S. cerevisiae strains used in
671
the present study
672
Figure 4: Effect of medium components on squalene synthesis (A) Effect of nitrogen source on
673
squalene synthesis using the strain SK10 harboring kanMX cassette for conferring G418
674
resistance when cultivated in (NH4)2SO4 and glutamate. (B) Effect of carbon source on squalene
675
synthesis using the SK10 having tHMG1 gene under constitutive promoter PTEF1 when cultivated
676
in glucose and galactose.
677
30 ACS Paragon Plus Environment
Page 31 of 39
Journal of Agricultural and Food Chemistry
Tables Table 1: Primers, plasmids and strains used in this study Name
Description
Reference
tHmg1-Gal-F
CGGGATCCGACCAATTGGTGAAAACTGAAG
This study
tHmg1-Gal-R
CGGAATTCGTAACACATGGTGCTGTTGT
This study
Primers
ATAAGAATAGGCGGCCGCATGGACCAATTGG tHmg1-pTEF - F
TGAAAACTGAAG
This study
CGGTCGACTAGTGTAACACATGGTGCTGTTG tHmg1 -pTEF- R
T
This study
CGGTCGGGATCCATGGACCAATTGGTGAAAA tHmg1-pPGK - F
CTGAAG
This study
ACGCGCGTCGACGTAACACATGGTGCTGTTG tHmg1 -pPGK- R
T
This study
Zwf1 – F
CGGTCGGGATCCATGAGTGAAGGCCCCGT
This study
ACGCGCGTCGACGCCGATAAATGAATGTGCT Zwf1 – R
TGC
This study
CGGTCGGGATCCACTTCCACAGTTCTCAACTC Pos5 – F(wM)
TTC
This study
CGGTCGGGATCCATGCACCATGAGTACGTTG Pos5 – F (w/oM)
GATTCACA
This study
ACGCGCGTCGACAGAGAATCTCATTGAATCT Pos5 – R
TTGCAT
This study
Yeast expression vector with URA selection,
Invitrogen
Plasmids pYES2-NTC
(Vickers et pCEV-G1-Km
Yeast expression vector with G418 selection
al. 2013)
pYGH
pYES2NT/C- PGAL1-tHmg1
This study
pCTH
pCEV-G1-Km-PTEF1-tHmg1
This study
pCPZ
pCEV-G1-Km-PPGK1-Zwf1
This study 31
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 32 of 39
pCPP(wM)
pCEV-G1-Km-PPGK1-Pos5(wM)
This study
pCPH
pCEV-G1-Km-PPGK1-tHmg1
This study
pCPP(w/oM)
pCEV-G1-Km-PPGK1-Pos5 (w/oM)
This study
pCTHPZ
pCEV-G1-Km-PTEF1-tHmg1-PPGK1-Zwf1
This study
pCTHPP(wM)
pCEV-G1-Km-PTEF1-tHmg1-PPGK1-Pos5(wM)
This study
pCTHPH
pCEV-G1-Km-PTEF1-tHmg1-PPGK1-tHmg1
This study
pCTHPP(w/oM)
pCEV-G1-Km-PTEF1-tHmg1-PPGK1-Pos5(w/oM)
This study
MATa; his3∆ 1; leu2∆ 0; met15∆ 0; ura3∆ 0
Euroscarf,
Strains BY4741
SRD GmbH SK1
Yeast BY4741 harboring pYES2/NTC
This study
SK3
Yeast BY4741 harboring pYGH
This study
SK8
Yeast BY4741 harboring pCEV-G1-Km
This study
SK10
Yeast BY4741 harboring pCTH
This study
SK12
Yeast BY4741 harboring pCPH
This study
SK13
Yeast BY4741 harboring pCPZ
This study
SK14
Yeast BY4741 harboring pCPP(wM)
This study
SK15
Yeast BY4741 harboring pCPP(w/oM)
This study
SK16
Yeast BY4741 harboring pCTHPH
This study
SK17
Yeast BY4741 harboring pCTHPZ
This study
SK18
Yeast BY4741 harboring pCTHPP(wM)
This study
SK19
Yeast BY4741 harboring pCTHPP(w/oM)
This study
32 ACS Paragon Plus Environment
Page 33 of 39
Journal of Agricultural and Food Chemistry
Table 2: Fold improvement in squalene production for different S. cerevisiae strains Strain
Overexpressed protein
Yield (mg/g
Yield fold
Conc.
Titer fold
DW)
increase a
(mg/L)
increase a
BY4741
NA
2.1±0.36
NA
2±0.21
NA
SK3
tHmg1p (1X)
6.5±2.6
3
6.1±0.05
3
SK10
tHmg1p (1X)
3.5±0.13
1.6
6.3±0.03
3
SK12
tHmg1p (1X)
4.4±2.07
2
4.8±2.39
2.3
SK13
Zwf1p
14.9±0.32
6.9
12.6±1.72
6.1
SK14
Pos5p (wM)
9.05±0.47
4.3
7.6±1.05
3.8
SK15
Pos5p (w/oM)
7.5±1.48
3.6
10±3.48
5
SK16
tHmg1p (2X)
35.7±0.92
16.8
28.1±3
13.7
SK17
tHmg1p & Zwf1p
47.4±3.2
22.2
24.9±0.32
12.2
tHmg1p & Pos5p
58.6±1.43
27.5
28.4±1.08
13.9
33±2.96
16.5
46±4.08
23
SK18
(wM) tHmg1p & Pos5p
SK19
a
(w/oM)
With respect to BY4741
NA: Not Applicable
33 ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 34 of 39
34 ACS Paragon Plus Environment
Page 35 of 39
Journal of Agricultural and Food Chemistry
Figure 1. Squalene synthesis in tHMG1 overexpressed S. cerevisiae strains. (A) The comparison of BY4741 with respect to SK13, SK10 and SK12 represents the effect of single copy overexpression of tHMG1 using promoters PGAL1, PTEF1 and PPGK1, respectively. Whereas, the comparison of SK16 to other engineered strains (SK10, SK12 and SK13) reveals the effect of plasmid harboring two copies tHMG1 to single copy tHMG1 on squalene. (B) Pictorial representation of gene cassettes used in the strains. 352x142mm (96 x 96 DPI)
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Figure 2. Squalene synthesis in S. cerevisiae strains with overexpression of NADPH regenerating genes (A) Comparison among strains SK13, SK14 and SK15 indicates the effect of overexpression of ZWF1, POS5 (with mitochondrial presequence) and POS5 (without mitochondrial presequence), respectively. (B) Comparison between the strains SK17, SK18 and SK19 indicates the synergistic effect of simultaneous expression of tHMG1 and NADPH regenerating genes. Here SK17, SK18 and SK19 represents SK13, SK14 and SK15 strains, respectively, with additional tHMG1 co-expression. 245x282mm (300 x 300 DPI)
ACS Paragon Plus Environment
Page 36 of 39
Page 37 of 39
Journal of Agricultural and Food Chemistry
Figure 3. Comparison of HMG-CoA reductase activity in various S. cerevisiae strains used in the present study 272x208mm (300 x 300 DPI)
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Figure 4. Effect of medium components on squalene synthesis (A) Effect of nitrogen source on squalene synthesis using the strain SK10 harboring kanMX cassette for conferring G418 resistance when cultivated in (NH4)2SO4 and glutamate. (B) Effect of carbon source on squalene synthesis using the SK10 having tHMG1 gene under constitutive promoter PTEF1 when cultivated in glucose and galactose. 278x282mm (300 x 300 DPI)
ACS Paragon Plus Environment
Page 38 of 39
Page 39 of 39
Journal of Agricultural and Food Chemistry
TOC 324x144mm (96 x 96 DPI)
ACS Paragon Plus Environment