Regeneration of NADPH Coupled with HMG-CoA Reductase Activity

Aug 28, 2017 - (12) Overexpression of a truncated form of Hmg1p coding only for the catalytic domain of the protein to bypass its regulation has not i...
0 downloads 11 Views 1MB Size
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