Research article | Open | Published:
Imbalance of heterologous protein folding and disulfide bond formation rates yields runaway oxidative stress
BMC Biologyvolume 10, Article number: 16 (2012)
The protein secretory pathway must process a wide assortment of native proteins for eukaryotic cells to function. As well, recombinant protein secretion is used extensively to produce many biologics and industrial enzymes. Therefore, secretory pathway dysfunction can be highly detrimental to the cell and can drastically inhibit product titers in biochemical production. Because the secretory pathway is a highly-integrated, multi-organelle system, dysfunction can happen at many levels and dissecting the root cause can be challenging. In this study, we apply a systems biology approach to analyze secretory pathway dysfunctions resulting from heterologous production of a small protein (insulin precursor) or a larger protein (α-amylase).
HAC1-dependent and independent dysfunctions and cellular responses were apparent across multiple datasets. In particular, processes involving (a) degradation of protein/recycling amino acids, (b) overall transcription/translation repression, and (c) oxidative stress were broadly associated with secretory stress.
Apparent runaway oxidative stress due to radical production observed here and elsewhere can be explained by a futile cycle of disulfide formation and breaking that consumes reduced glutathione and produces reactive oxygen species. The futile cycle is dominating when protein folding rates are low relative to disulfide bond formation rates. While not strictly conclusive with the present data, this insight does provide a molecular interpretation to an, until now, largely empirical understanding of optimizing heterologous protein secretion. This molecular insight has direct implications on engineering a broad range of recombinant proteins for secretion and provides potential hypotheses for the root causes of several secretory-associated diseases.
The protein secretory pathway is an extensive process in eukaryal cells, as it is responsible for processing approximately one-third of all proteins. Substantial cellular resources are therefore utilized to maintain this pathway's functions, and stressed conditions in the secretory pathway have consequences for the whole cell . Distress in secretory pathway organelles has been implicated as the molecular basis for several diseases, for example, β cell apoptosis in diabetes, cystic fibrosis, and prion-related disease, among others . In biotechnology, efficient secretion of useful recombinant proteins in yeast and fungi is a key industrial objective with applications in enzyme production required for the production of biofuels, detergents, fabrics, food, and biologics, such as imunoglobulins, hormones, and vaccines. Significant effort has gone into engineering yeast for increasing protein secretion . Strategies, such as changing environmental parameters (for example, temperature, media composition)  or altering genetics, can increase secretion for some proteins, but they rarely represent generic solutions for improving protein secretion [5, 6]. The lack of a single engineering strategy that improves protein secretion across the board implies that there are several possible bottlenecks in the secretory pathway, and different proteins may be constrained in different ways. There is therefore a requirement for more fundamental insight into this complex pathway that involves a very large number of components.
In yeast, the secretory pathway is a multi-organelle system that is responsible for trafficking proteins to the extracellular space, cell membrane, or vacuole . During this transit, multiple processes must be coordinated, including folding, specific proteolytic cleavage, glycosylation, and disulfide bond formation, all with a layer of quality control at key check points. The pathway requires substantial cellular resources to perform these tasks, such as glycans, electron acceptors, electron donors, and ATP. In the ER, the nascent peptide is folded into its native structure while disulfide bonds are formed. The rate of protein folding is dependent upon the complexity of the protein to be folded, the availability of chaperones to assist folding, and ATP used by the chaperones . Proteins that are slow to fold or terminally misfolded proteins are removed from the ER via the ER-associated degradation (ERAD) pathway . Disulfide bond formation requires the removal of electrons from cysteine thiols via protein disulfide isomerase (PDI) and Ero1p to the final electron acceptor, typically oxygen [9, 10]. This process produces reactive oxygen species (ROS) in stoichiometric amounts to the number of disulfide bonds formed . Disulfide bond formation is random, and incorrect bond pairs must be exchanged for native bonds via PDI-based processes . In addition, reduced glutathione (GSH) acts as a buffer for the redox state of the ER . A more detailed description of oxidative protein folding can be found in the reviews by Sevier et al. and Chakravarthi et al. [14, 15].
The secretory pathway must adjust the chaperone capacity, oxidizing equivalents, ATP, glycan, and other metabolic requirements, as well as trafficking patterns, based on the portfolio of proteins that need to be expressed at a given time, and the resources required to process that set of proteins. In yeast, the unfolded protein response (UPR) is one transcriptional mechanism that adjusts secretory resources and controls to handle overload of the folding machinery in the ER . In the UPR, accumulation of unfolded proteins in the ER signals a pathway that results in translation of Hac1p, a transcription factor (TF) known to activate or repress over 100 genes, including many ER-associated proteins such as Kar2p, Pdi1p, and Ero1p .
In this study, we identified biological mechanisms which alter the secretory pathway in response to secretion of recombinant proteins with different properties (size, number of disulfide bonds, and glycans) in a Hac1p-dependent and independent manner. The secretory pathway was perturbed by secreting a small protein, human insulin precursor (IP), or a comparatively larger protein, α-amylase, in wild-type (WT) and Δhac1 Saccharomyces cerevisiae. These proteins were chosen because the two proteins elicit different behavior in the secretory pathway. These differences will arise because α-amylase is a relatively larger (and likely more difficult to fold), has an odd number of cysteines (which may complicate disulfide isomerization) and has glycosylation, compared to insulin which is small, has even number of cysteines, and is not glycosylated. As well, α-amylase has one more disulfide bond than IP. To identify biological mechanisms, we characterized changes in physiological properties (specific growth rate, carbon utilization efficiency, and recombinant protein secretion), TF activity (as inferred from transcriptome analysis) and metabolic demand (as inferred by changes in metabolic flux diversion). Through this, we identified the following biological processes: amino acids recycling from degraded proteins, trans-Golgi network (TGN) sorting changes, overall expression repression, and oxidative stress. Motivated by secretory-related oxidative stress observations, we present a model for disulfide bond formation and electron transfer in the ER which takes into account thermodynamic irreversibilities caused by differences in electron affinity. The proposed model explains the non-stoichiometric ROS formation that we observed that results from disulfide bond formation and causes oxidative stress under folding-stress conditions. If proven by genetic and biochemical results, the futile cycle model yields insight into a fundamental problem in secretory stress and reveals new avenues to reduce oxidative stress and increase productivity in industrial protein production.
Protein size and Hac1p activity affect protein secretion quantity and cell growth
Yeast strains were constructed that produce and secrete (a) IP or (b) α-amylase and were compared to yeast strains containing (c) an empty vector in both wild-type and HAC1 deletion backgrounds. IP and α-amylase were chosen because they are very different types of proteins to secrete. IP is 51 amino acids in length, with six cysteines forming three disulfide bonds, and no glycosylation. α-amylase is 478 amino acid in length, with nine cysteines forming only four disulfide bonds and one glycosylation. The odd number of cysteines in α-amylase complicates disulfide pairing, as the random isomerization process may incorporate the cysteine that should not be incorporated into a disulfide bond. Both proteins were targeted for secretion using a YAP3 pre sequence (21 amino acids, cleaved off in the ER) and a rationally designed pro sequence (TA57, 42 amino acids, no glycosylation or disulfides) were cloned behind a TDH3 promoter in a high copy 2 micron plasmid . α-amylase was expressed using the same plasmid, promoter, and leader sequences. These strains are named WN (WT with empty vector), WI (WT secreting IP), WA (WT secreting α-amylase), dN (Δhac1 with empty vector), dI (Δhac1 secreting IP), and dA (Δhac1 secreting α-amylase). Strains were characterized in batch fermentation to understand the effects on cell physiology.
The cellular burden induced by (a) synthesizing and secreting IP and α-amylase and (b) deleting the key TF for the UPR, Hac1p, substantially affected the cells. Protein titers in WT strain were 9 mg/L and 20 mg/L, for IP and α-amylase, respectively (Figure 1a). On a per biomass basis, this is approximately half the insulin produced, and one-third the α-amylase reported for rich media [19, 20]. Rich media appears to be favorable for heterologous protein production, but may present complications in downstream separations. Comparing the small and larger proteins, α-amylase was secreted in higher levels on a mass basis, but six-fold more insulin molecules were secreted (1.52 μM IP in WI compared to 0.26 μM α-amylase in WA). Δhac1 strains secreted significantly less protein than WT, confirming that Hac1p is important for efficient secretion (Figure 1a) .
Reduced specific growth rates imply impairment of cellular processes (Figure 1b). In WT yeast, IP production did not affect growth; however, α-amylase production reduced growth by 25%. This, combined with the differences in protein titers, implies that α-amylase is more challenging to fold and secrete than IP. In the Δhac1 background, recombinant protein strains dI and dA had approximately 20% lower growth rates compared to dN. This growth reduction occurs despite no change in specific glucose uptake rate (Additional file 1, Tables S1 and S2) pointing toward higher energy requirements to maintain homeostasis in Δhac1 while trying to secrete recombinant proteins. Δhac1 strains had overall lower final cell densities. Δhac1 strains produced more glycerol than WT strains implying impaired oxidative processes in the Δhac1 strains (Additional file 2).
Secretory stress shifts metabolism to increase oxygen and ATP requirements
The physiological changes due to the secretory perturbations affect the distribution of resources through the metabolic network. The glucose uptake and range of products produced were altered by the protein production conditions (Table 1). Changes in the underlying metabolic network were estimated by flux balance analysis (FBA) using a yeast central carbon metabolism model, constrained by measured extracellular fluxes (Additional file 1, Tables S1 and S2, Additional files 3 and 4) . Figure 2a shows a metabolic map of central carbon metabolism for each of the six conditions based on the exchange fluxes in Table 1 and the FBA analysis. The shift in metabolic fluxes were correlated with changes in redox requirements. As expected, the catabolic functions of the TCA cycle was predicted to have very low activity due to glucose repression . Figure 2b shows that the oxygen uptake was twice as high in the strains that were growth inhibited (for example, WA, dI, dA) than those that were not. This increased oxygen uptake was not used for oxidative phosphorylation, as the biomass yields on glucose were lower in WA, dI, and dA, and it may therefore be a result of increased oxidation in connection with formation of disulfide bonds.
Figure 2c shows that the maintenance ATP consumption is increased in WA, dI, and dA according to FBA calculations. In WT background, WI did not consume a detectable increase in ATP, likely because IP is short and easily folded, thereby minimally taxing the translation and folding machinery. WA did increase two-fold in ATP consumption, most likely because α-amylase is 10-fold larger and likely more difficult to fold and has more disulfide bond pairing possibilities. In the Δhac1 background, folding efficiency is likely decreased due to ER dysfunction. With native secretion, dN did not require higher ATP maintenance consumption compared to WT. However, even the smaller, easier to fold IP resulted in ER stress that required significant ATP consumption compared to WT. dA, which was already stressed under WT, continued to show high ATP consumption. Despite the increased ATP consumption in dI and dA, little protein was secreted.
Transcription factors controlling oxidative stress, amino acid salvaging, and expression repression are linked to secretory response
Growth phase transcriptomics measurements were carried out to identify cellular processes that were activated under the stresses of HAC1 deletion and recombinant protein production. HAC1 deletion resulted in 339 significantly changed genes in the no recombinant protein case (WN vs. dN). HAC1 deletions in the insulin strain and α-amylase strain resulted in much larger cellular responses of 1628 (WI vs. dI) and 1511 (WA vs. dA) significantly expressed genes, respectively. KAR2 (ER chaperone) expression was significantly reduced upon HAC1 deletion (↓ three-fold dN vs WTN, P = 1 × 10-4) and the four yeast protein disulfide isomerases (PDI1, EUG1, MPD1, MPD2) reduced an average of 2.9-fold (P < 0.05).
The effects of producing IP or α-amylase within a strain background (WT or HAC1) were not as pronounced as the effect of HAC1 deletion, 40 and 194 genes were significantly changed in WI (compared to WN) and WA (compared to WN). Likwsise, 74 and 90 genes were significantly changed for dI (compared to dN) and dA (compared to dN).
To reduce the dimensionality of the data and identify putative TFs involved in protein secretion, the Reporter Transcription Factor algorithm was used . TFs were scored by the modulation in expression level of genes that the TFs bind in the upstream region according to ChIP-chip data . Therefore, the score is not indicative of change in the TF expression level itself, but of the genes under its influence. Reporter TF algorithm is useful, because although the statistical significance of an individual gene may not meet an arbitrary threshold, if several genes linked to the same TF have similar behavior, the likelihood of observing the group of genes is low, making TF identification very sensitive. Figure 3 shows significant secretory process TFs shown to be involved in up- and down-regulating different cellular process under their control. Interestingly, different TFs were identified for the two different proteins. This is likely the combined effect of different protein size and number of disulfide bonds. A complete list of significant transcription factors is provided in Additional files 5 and 6.
In WT (Figure 3a), several TFs were activated by protein secretion. Oxidative and osmotic stress pathway up-regulation was common to both proteins. Oxidative stress is likely caused by ROS that is formed when Ero1p shuttles electrons to oxygen in disulfide bond formation . Osmotic stress response, particular hypo-osmotic stress, strengthens the cell wall to counteract internal turgor pressure by changing the cell wall composition. This change in composition requires remodeling the secretory pathway by changing which components are trafficked to the cell wall . Surprisingly, the Reporter TF algorithm found several Hac1p-influenced genes down-regulated. Genes that Hac1p binds from the ChIP-chip data that are significantly down-regulated are KEG1, MCD4, and ERJ5. KEG1 and MCD4 genes are involved in glycan modifications and ERJ5 is a secondary ER chaperone [27–29]. These genes may be influenced by other TFs not included in the ChIP-chip network. Genes known to be regulated by Hac1p (KAR2 and ERO1) were not significantly changed upon secreting recombinant protein, indicating that there is not an actual Hac1p response in the WT.
Clear differences between large and small protein secretion emerge in WT. IP stimulated modification of the TGN through MCM1 and STE12. Overall expression is reduced by altering mRNA degradation pathways via STO1. α-amylase had a much larger effect on the cell, as compared to IP, as was implied by physiological parameters of Figure 1 and number of altered genes. Additional oxidative and osmotic stress pathways were activated in WA, as well as a down-regulation in some amino acid synthesis pathways and overall reduction of transcription.
In the Δhac1 background (Figure 3b), many of the effects found in WA, have become common to both IP and α-amylase producing strains. HAC1 deletion clearly makes the cell more susceptible to recombinant secretion overload. Both insulin and α-amylase secretion cause considerable oxidative stress response and down-regulation of amino acid synthesis, including the general amino acid synthesis TF, Gcn4p. In dI, translational capacity repression is also employed (via Fhlp/Rap1p) and adjustments in amino acid metabolism. dA shows a mix of up- and down-regulation of genes that are controlled by Hac1p. Other TFs appear to be controlling these genes in the absence of HAC1. Some oxidative and osmotic stress pathways appear independent of HAC1. Skn7p and Cin5p were similarly activated in both WT and Δhac1. Oxidative and hypo-osmotic stress, while important for managing the secretory pathway, appears not to be directly managed through the UPR.
Thermodynamic irreversibilities in redox reactions can explain increased oxidative stress in slow protein folding conditions
The increases in oxidative stress, oxygen consumption, and reduced growth observed in the study can be explained by electron transfer in ER redox pathways. Disulfide bond formation has been established to consume oxygen and produce ROS (and thereby consume cellular resources to protect against the ROS) in stoichiometric quantities with the number of disulfide bonds formed . When non-native disulfide linkages are formed, these linkages must be rearranged to the correct disulfide pairings for the native protein to be folded, a process called disulfide isomerization .
Disulfide isomerization involves (a) breaking the non-native bond by transferring electrons to the non-native bond creating a cysteine linkage with the PDI, and (b) creating a new disulfide linkage in the nascent protein by transferring the electrons to break the PDI-nascent protein linkage. By random pairing, the native disulfide bonds are found.
Directionality in these redox reactions is determined by thermodynamic favorability through electron affinity of the potential disulfide bonds. Disulfide isomerization is redox neutral, not requiring electron donors or acceptors. However, it does require each disulfide pairing to have a lower electron affinity than the next (non-native disulfide in folding protein < PDI-folding protein disulfide < native disulfide in folding protein) to allow the electrons to transfer. Under slow folding conditions, PDI may hold the disulfide bond (oxidized state) for extended time because a native disulfide cannot be found, resulting in PDI being reduced by other moieties, likely GSH.
Given the observations in our experiments, and the thermodynamic reasoning immediately above, we propose a simple thermodynamic model of disulfide bond formation and breaking that explains increased oxidative stress, oxygen consumption, and reduced growth observed in our experiments. This model expands upon the mechanism by Cuozzo and Kaiser . The thermodynamic model assumes there are PDI disulfide bonds that have electron affinities above and below the nascent proteins disulfide bonds (Figure 4a). The disulfide is formed by the typical oxidation pathway (Figure 4a, green) catalyzed by high electron affinity PDI (called PDIA here). Instead of isomerization, the incorrect disulfide is reduced by an electron donor with a low electron affinity (most likely a different PDI paralogue, called PDIB here) (Figure 4a, blue). The difference in electron affinity between the folding protein's cysteines and a specific PDI's cysteines can only allow the electrons to flow in one direction (toward the higher electron affinity cysteines) (Figure 4a). Therefore, a different PDI is required to form and break the incorrect disulfide bond. This futile cycle relies on a strong electron affinity gradient to complete an isomerization-like process. The net result of the futile cycle is GSH consumption and ROS production. This model implies that the ROS produced is not stoichiometrically linked to the number of disulfide bonds formed, but varies by the number of futile cycles before the correct bond is formed.
The metabolic and transcriptional data supports this model. Upon HAC1 deletion, ER chaperones (KAR2) and PDIs (PDI1, FUG1, MPD1, and MPD2) expression is reduced. This downregulation of ER chaperones and PDIs results in suppressed ER folding and disulfide bond formation in the Δhac1 mutants. In the dN case, minimal oxidation stress is seen. However, when there is an increased demand for protein folding and disulfide bond formation, as is the case for dI and dA case, we see high oxygen consumption, ATP requirements, and many oxidative stress pathways being activated transcriptionally. Although both folding and disulfide bond formation is down, an imbalance toward faster disulfide bond formation compared to folding will result in futile cycles. Therefore, this disulfide/folding imbalance acts as a catalyst for drastically increasing ROS production.
Based on this thermodynamic model, the relative rates of protein folding and disulfide bond formation for nascent peptides have important consequences for oxidative stress (Figure 4b). When folding is faster than disulfide bond formation, ROS is produced in near one-to-one amounts with the disulfide bonds formed. Under these conditions, isomerization may be more efficient to resort incorrect disulfide bonds, as native structures with low electron affinity disulfide pairs are favored, and isomerization does not produce ROS. However, when folding is slow compared to disulfide bond formation, as is the case when the protein folding machinery gets overloaded, the nascent peptides cycles through the futile redox cycle producing ROS and consuming GSH in excess to the final number of disulfide bonds formed. The physiological result of a high disulfide bond formation to ER folding rate is oxidative damage to a broad range of cellular proteins and consumption of reducing equivalents that could otherwise be used for anabolism.
In this study, we have identified biological mechanisms related to protein synthesis and secretion by introducing perturbations to the cell, in the form of HAC1 deletion and different recombinant protein expression, and measuring the system level cellular responses, via transcriptomics and metabolic fluxes. These measurements, combined with data analysis algorithms, Reporter TF algorithm and FBA, were able to identify cellular adjustments in (a) overall expression level, (b) post-Golgi sorting, (c) amino acid biosynthesis and savaging, and (d) oxidative stress. These biological effects are a result of the combined influence of protein synthesis and trafficking through the secretory pathway.
Overall transcription and translation were repressed in response to α-amylase expression (a larger protein) and in the Δhac1 strains with any recombinant protein secretion. Repressing overall expression is a broad spectrum response used to adjust the rates of all other cellular processes to match the reduced folding capacity in the ER. Several mechanisms were used to alter overall expression: repressing mRNA synthesis, increasing mRNA degradation rates, and repressing protein translation rates through reducing ribosome numbers. Specifically, mRNA concentrations are lowered by decreasing RNA polymerase accessibility (HIR2), inhibiting transcriptional elongation (THO2), and controlling RNA degradation (STO1) [31, 32]. Ribosome concentration, and thereby translation rates, can be reduced by the TFs Fhl1p and Rap1p which control expression of rRNA and ribosomal proteins . This is seen in IP production in Δhac1 strain, both by the reporter TFs (Figure 3) and by expression of ribosomal proteins (Additional file 7). In this context, extrachromosomal plasmids offer advantages over chromosomal expression. HIR2, whose mechanism is to silence the chromosome, would not affect extrachromosomal plasmids. Increased recombinant protein secretion would be accomplished by silencing native ER genes, while recombinant, plasmid-born gene would not be affected.
Pronounced adjustments to the TGN were observed in the transcriptome in all conditions. TFs involved in pheromone responses (STE12, MCM1, ASH1), invasive/pseudohyphal growth (STE12, MSN1, PHD1, RIM101), and osmotic stress (CIN5, SKN7, SKO1, YAP6, MSN1) were all identified by the Reporter TF algorithm and point to an underlying set of activities that are required to increase the traffic of secretory vesicles to the membrane. Invasive, pseudohyphal, and filamentous growth morphologies have a high surface to volume ratio and inherently require higher Golgi-to-cell membrane trafficking rates to supply cell membrane and cell wall components for growth. These altered morphologies can be activated through the filamentous and invasive response elements (FREs)  bound by STE12 and used to regulate PHD1 . HAC1 deletion has been shown to cause filamentous growth .
Osmotic stress TFs are also responsible for affecting protein secretion, as the external cell wall must be strengthened in response to hypo-osmotic conditions, thereby requiring an efficient secretory pathway to ferry cell wall proteins . MSN1 is known to induce starch degradation, requiring the actions necessary to secrete the appropriate enzymes through filamentous growth activation . SKN7 has a dual role in invasive growth and osmotic stress . Although osmotic stress TFs are commonly associated with the hyper-osmotic glycerol (HOG) pathway, Ypd1p can phosphorylate Skn7p, signaling the hypo-osmotic stress pathway . Because there were no apparent hypo-osmotic conditions in this study, this indicates that these TFs are not directly controlled by osmotic conditions, but possibly through a secondary response to upregulation and increased secretion of cell wall proteins.
TGN TFs and/or the genes they regulate are possible targets for increasing Golgi-to-cell membrane trafficking. In S. cerevisiae, recombinant protein intended for secretion has been found mis-trafficked to the vacuole. This has been shown for insulin and green fluorescent protein secretion in yeast [40, 41]. Proteins involved in vesicle trafficking, namely Sly1 and Munc18 have been found to increase recombinant secretory rates in Chinese hamster ovarian (CHO) and several mammalian cell lines [42, 43]. It is likely that similar proteins are present in yeast and could be exploited for improving protein secretion.
Significant alterations in amino acid metabolism were observed, particularly in the Δhac1 strains. De novo amino acid synthesis (GCN4, BAS1, MET32, ARG81, RTG3) was suppressed. On the surface, this appears contradictory, as increased amino acid requirements should be observed with recombinant protein production. However, this decrease in de novo amino acid synthesis is accompanied by observed increases in scavenging mechanisms for amino acids (SNT2, CUP9, PUT3). High scavenging rates and decrease synthesis imply high protein degradation rates where the degraded proteins result in available amino acids for scavenging; reducing the need for newly synthesized amino acids. This is consistent with either ERAD, a process where proteins that are stalled in the ER are transported back into the cytoplasm for degradation by the proteosome, or vacuolar-localized protein degradation. In either case, the cell is expending energy on synthesizing proteins that are ultimately degraded. These effects appear in the strains that are the slowest growing with the highest ATP requirements (Figures 1b and 2b). In these cases the ER folding capacity is likely saturated, resulting in ER holdup and ERAD.
Oxidative stress TFs were also found in all conditions. Several were dual oxidative/osmotic stress TFs (CIN5, SKN7, SKO1), and others were dedicated to oxidative stress only (AFT2, YAP1). TFs were found in all three of the major oxidative stress signaling pathways, (a) the Hog1 MAPK pathway (where SKO1 is the DNA binding agent), (b) Sln1 pathway (where SKN7 is the DNA binding agent), and (c) YAP1 and CIN5, which directly sense oxidative stress and bind DNA . The cell's control machinery appears to have hard-wired oxidative stress responses to increased secretory demand, as oxidative/hypo-osmotic pathways have a high degree of overlap, which is appropriate because increased secretion of cell wall proteins will result in higher oxidative stress. In particular, Skn7p, which has already been mentioned for its role in managing secretory pathway directly in an osmotic stress pathway, can also activate oxidative stress response genes .
Oxidative stress was pronounced with all secretory perturbations and has been identified in other studies to be associated with secretory stress [1, 17]. Futile cycling may be the dominant disulfide resorting pathway when folding is limited. In previous studies, oxidative stress, induced by tunicamycin, a N-linked glycosylation inhibitor, increased with ER stress, despite no increase in the net disulfide bond formation demand . The futile cycle does predict non-stoichiometric ROS formation, while isomerization does not. ROS can be formed at potentially limitless amounts through multiple rounds of disulfide formation and breaking. This will occur under conditions where the rate of folding is slow, a result of proteins that are specifically difficult to fold, or a result of the overall ER folding capacity being saturated. As well, futile cycling will increase as the number of available cysteine residues available for disulfide bonding increase, as is the case for α-amylase, due to the extended amount of isomerization that may be needed to form the correct disulfide bonds.
One implication of the proposed thermodynamic model is that PDI paralogues, or cysteines within a PDI, must exist at different electron affinities that are above and below the electron affinity of the protein to be folded. Although in vivo redox potentials of PDI cysteine pairs were not measured, from first principles it would appear highly likely that these PDIs would need different redox potentials to carry out isomerization. In Figure 4a, we assume that only PDIs interact with the folding protein. This appears the case, as kinetic rates for direct glutathione oxidation/reduction are too slow to be physiologically relevant . Electron affinity (and therefore redox potential) is broadly determined by the proximity of the two cysteines, with the proximity determined by the current structure of the protein . Cysteines that are in the correct orientation will have a low electron affinity and easily form disulfide bonds, while cysteines that are not in the correct orientation will have a high electron affinity and will have unstable disulfide bonds. Therefore, the electron affinity of a correctly folded/correct disulfide bond would be lower than that of a misfolded or incorrect disulfide bond. This difference in electron affinity may allow PDIs to selectively break disulfides with high electron affinity (incorrect bonds), but not disulfide bonds with low affinity (correct bonds).
The need for different PDIs to form or break disulfide bonds may explain the need for many PDI homologues in the ER, each with different structures, and therefore different electron affinities. These PDIs can only span a finite range of electron affinities, and there may be implications for proteins that have disulfide pairs with electron affinities higher than the highest PDI or lower than the lowest PDI. If no PDI has a lower electron affinity than an incorrect disulfide bond, then the disulfide bond cannot be broken and the protein is terminally misfolded. As well, a protein that has a native disulfide pairing with an electron affinity higher than any PDI cannot form a bond. This may be the case when recombinant proteins are being processed in the ER.
Futile cycling as a large potential ROS source has broad implications on the cell. Tu and Weissman predict Ero1p-produced ROS that is one-to-one with disulfide bond formation could attribute approximately 25% of cellular ROS to the secretory pathway . Therefore, even larger ROS production is likely if the futile cycle is the dominant disulfide resorting pathway under folding stress. This also has implications on GSH and possibly NADPH availability, as it is doubly consumed (a) by the reduction of ROS and (b) directly in the futile cycle. The futile cycle limits reducing equivalents needed for anabolic processes, and may explain the reduced growth rates observed in folding stressed strains (WA, dI, and dA).
In all, Figure 4b highlights that the relative rates of two processes, protein folding and disulfide bond formation, must be kept in balance to avoid significant cellular stress. If disulfide bond formation is fast compared to folding, high futile cycle use will result in high ROS formation, NADPH loss, and high protein degradation as a result of ERAD. This scenario is observed in the Δhac1 strains dI and dA.
The engineering implications for protein secretion become much clearer with this understanding of protein folding to disulfide bond formation ratio. When overexpressing a recombinant protein, an optimal expression must be found, where transcription is as high as possible without overloading the ER folding capacity and sending the cell into an oxidative stressed state. This optimal expression level will be different for different proteins, as protein folding rates will vary according to the protein size and structure. We see this in comparing IP and α-amylase expression. The concept of an optimal expression has been identified heuristically, in the present study we identify the competing molecular effects that could define these phenomena . This optimal expression ratio extends to recombinant proteins that do not have disulfide bonds. For recombinant proteins without disulfide bonds, recombinant protein folding in the ER will consume folding resources, thus slowing down folding rates. Although the recombinant protein has no disulfide bonds, many native proteins still require disulfide bonds. Because of this, the folding to disulfide bond formation ratio will be disturbed, resulting in similar ROS stress.
To maintain an optimal ratio, either protein folding rates must increase or oxidation rates decrease. Overexpression of chaperones that increase folding capacity has successfully been used to increase protein secretion [6, 48]. For particularly large or difficult to fold proteins this may not be adequate. A new approach may be to limit the oxidation rate of Ero1p to slow down the first step of the futile cycle. This would be done in concert with repressing ERAD, as proteins would have long retention times in the ER. In this scenario, recombinant proteins would be slowly folded, albeit without high cellular stresses. This would result in longer overall process times, but may be required for difficult to fold proteins.
In this study, we identified post-Golgi vesicle sorting, high protein degradation rates, repressed overall expression, and oxidative stress in response to +/- UPR strains secreting different sized recombinant protein. These processes were identified through scoring TFs and estimating alteration to the metabolic network. These observations imply our proposed futile cycling is the dominant disulfide resorting pathway in the ER and explains non-stoichiometric ROS formation seen in our study and elsewhere. The futile cycle model, producing ROS and consuming GSH, has a clear thermodynamic driving force compared to disulfide bond isomerization. If correct, futile cycling is likely the dominant mechanism under secretory stress. This interplay between protein folding and futile cycling sheds light on a largely empirical understanding of engineering protein secretion and implies the relative rates of protein folding and disulfide bond formation are critical to maintaining cellular homeostasis. This increased molecular understanding of the secretory pathway should allow for more insightful design of secretory engineering strategies.
Strains and media
All experiments were performed in the background of CEN.PK 113-5D (MAT a SUC2 MAL2-8c ura3-52, P. Kötter, Frankfurt, Germany) . Genomic DNA from Y05650 (BY4741; Mat a; his3D1; leu2D0; met15D0; ura3D0; YFL031w::kanMX4, obtained from EUROSCARF) was used as a template for the HAC1 knockout cassette. Standard molecular biology techniques were used  and all plasmids were maintained in Escherichia coli DH5α in Luria Bertani (LB) broth with 80 mg/L ampicillin. PCR primers are listed in Additional file 8.
Genomic DNA was purified from Y05650 using Fast DNA Spin Kit for Soil (MP Biomedicals Solon, OH, USA). A 2.6 kb DNA fragment containing the genomic replacement of HAC1 with KanMX and flanking regions was amplified by PCR using primers KT007/KT008 (Additional file 8). The HAC1::kanMX4 fragment was integrated at the HAC1 loci of CEN.PK 113-5D by standard yeast transformation  and selected on 200 mg/L G418 to create the Δhac1 strain. Correct integration was confirmed by PCR.
DNA coding for an insulin precursor with a Yap3 pre-leader sequence and the TA57 pro-leader sequence and spacers as described  for correct secretory processing was synthesized with optimal codon usage for yeast and delivered on plasmid pUC57-Yap3Insulin (GenScript Co. Piscataway, NJ, USA) (Additional file 9 for sequence). α-Amylase DNA was amplified from Saccharomyces kluyveri YKM37  using LZH018 and LZH039. The pre-pro-leader was amplified from pUC57-Yap3Insulin using primers LZH015 and LZH016. The pre-pro-leader was connected to the α-amylase by fusion PCR of the two segments together using primers LZH015 and LZH039 . The pre-pro-insulin and pre-pro-amylase were cloned into the SpeI/SalI or SpeI/EcoRI sites of p426GPD, respectively, downstream of the constitutive GAPDH promoter , to create pYapIns and pYapAmy. Plasmids p426GPD, pYapIns, and pYapAmy were transformed into CEN.PK 113-5D and Δhac1 strains by standard methods .
Strains were grown in SD-2xSCAA , containing 20 g/L glucose, 6.7 g/L yeast nitrogen base minus amino acids (Formedium, Norfolk, UK), 2 g/L KH2PO4 (pH = 6 by NaOH), 190 mg/L Arg, 108 mg/L Met, 52 mg/L Tyr, 290 mg/L Ile, 440 mg/L Lys, 200 mg/L Phe, 1260 mg/L Glu, 400 mg/L Asp, 380 mg/L Val, 220 mg/L Thr, 130 mg/L Gly, 400 mg/L Leu, 40 mg/L Trp, 140 mg/L His, 1 g/L bovine serum albumin. Five hundred mL of medium was inoculated in a 1 L bioreactor (DasGip, Jülich, Germany) at 30°C, 600 rpm agitation, 30 standard L/h air flow, pH controlled at 6 by KOH (2 M). Strains were inoculated to an A600 = 0.01 from late exponential phase cultures and A600 was measured throughout the cultivation. Dry cell weight (DCW) was measured by filtering 5 mL of culture broth through a 0.45 μm nitrocellulose filter and measuring the increased weight of the dry filter. Glucose, ethanol, glycerol, and acetate were measured using a Summit HPLC (Dionex, Thermo Scientific, Waltham, MA, USA) with an Aminex HPX-87H column (Bio-Rad, Hercules, CA, USA). Carbon dioxide and oxygen levels were measured in the off-gas and dissolved oxygen was monitored. Transcriptome samples were taken after 5+ doublings at A600 = 1.0-1.4. Triplicate fermentations were carried out for each strain.
Insulin was measured by a modification of the assay by Snell et al. . One mL of cell culture was centrifuged at 4000 × g for 4 min. Eight parts supernatant was added to one part 0.1 N HCl and 5.5 μM sodium azide and stored at 4°C until measurement. Insulin concentration was determined by HPLC using a Luna 5 μ C18(2) (250 mm × 4.6 mm) (Phenomenex, Torrance, CA, USA) column and gradient-based elution. Buffer A contained 68 mM phosphoric acid, 0.2 M sodiumsulphate and 10% (w/v) acetonitrile in water, and Buffer B contained 50% acetonitrile in water. HPLC was run with 25 μL injections at 1 mL/min and 50°C. Gradient protocol: 20% B for 10 min. Linear gradient from 20% B to 60% B over 10 min. Hold at 60% B for 5 min and then to 20% B for 3 min to re-equilibrate for next sample. Insulin standards eluted at 22.6 min and insulin precursor at 20.0 min. HPLC peaks were verified to be the correct protein by SDS-PAGE. Human insulin was used as a standard (Sigma, St. Louis, MO, USA).
α-amylase concentration was calculated from enzyme activity. α-amylase activity was measured using the Ceralpha kit (Megazyme K-CERA, Wicklow, Ireland) using α-amylase from Aspergillus oryzae (Sigma, St. Louis, MO, USA) as a standard. This conversion was calculated using a 1.79 U/mg (weight includes salts and purified protein) standard from Sigma using the Protein 80 chip on the Bioanalyzer (Agilent, Santa Clara, CA, USA). By this, α-amylase was found to be 0.0257 g α-amylase/g total.α-amylase activity was then converted to mass using 70 U/mg α-amylase protein.
Samples for microarray were taken as described previously and stored at -80°C until processing . RNA was isolated using the RNeasy Minikit (Qiagen, Valencia, CA, USA). Cells were lysed in RNeasy RLT buffer using Lysing Matrix C (MP Biomedicals Solon, OH, USA) in a Fast Prep 24 (MP Biomedicals Solon, OH, USA) as follows: 20 s at speed 6, 1 min at 4°C, 20 s at speed 6. RNA was processed to aRNA using the Genechip 3' IVT Express Kit (Affymetrix, Santa Clara, CA, USA) and hybridized/scanned on the Yeast Genome 2.0 Array (Affymetrix, Santa Clara, CA, USA) following commercial protocols to create CEL files.
Images were analyzed using R 2.10.1 statistical software and the 'affy' and 'limma' packages as described previously . Briefly, background normalization was carried out using robust multi-array (RMA) average method with perfect match (PM) probes only. Interchip normalization used the qspline algorithm with median polish summary method. Statistical analysis was carried out by comparison of triplicate bioreactor measurements for each strain. Emperical Bayesian statistics were used to moderate standard errors within each gene and Benjamini-Hochberg's method to adjust for multiple testing. Microarray data was submitted to the GEO database and have accession number GSE27062 (see http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=dpyzfywysoqecbk&acc=GSE27062).
Reporter transcription factor analysis
Transcription factor activity was scored using the Reporter Effector algorithm . Transcription factor-DNA interactions were gathered from ChIP-chip with P < 0.001 . Significant interactions were found for 176 transcription factors regulating 3,796 genes for a total of 10,849 unique interactions. Gene P values from comparing different strains were used to score transcription factors that were known to bind to the upstream DNA. Transcription factors with P < 0.05 of being activated between conditions are reported.
Flux balance analysis
Estimates of intracellular reaction rates were performed using measured exchange fluxes of glucose, ethanol, acetate, glycerol, and carbon dioxide. Model-based error correction was used to close carbon and electron balances . Flux balance analysis was carried out using a 85 reaction model of yeast central carbon metabolism and biomass yield were used . Additional file 4 contains the complete results of the analysis which are used to estimate ATP consumption in the different strains.
Δhac1 secreting α-amylase
Δhac1 secreting IP
Δhac1 with empty vector
flux balance analysis
filamentous and invasive responsive elements
nicotinamide adenine dinucleotide phosphate
protein disulfide isomerase
reactive oxygen species
unfolded protein response
WT secreting α-amylase
WT secreting IP
WT with empty vector
Tu BP, Weissman JS: Oxidative protein folding in eukaryotes: mechanisms and consequences. Journal of Cell Biology. 2004, 164 (3): 341-346. 10.1083/jcb.200311055.
Petranovic D, Tyo K, Vemuri GN, Nielsen J: Prospects of yeast systems biology for human health: integrating lipid, protein and energy metabolism. FEMS Yeast Res. 2010, 10 (8): 1046-1059. 10.1111/j.1567-1364.2010.00689.x.
Idiris A, Tohda H, Kumagai H, Takegawa K: Engineering of protein secretion in yeast: strategies and impact on protein production. Applied Microbiology and Biotechnology. 2010, 86 (2): 403-417. 10.1007/s00253-010-2447-0.
Wedekind A, O'Malley MA, Niebauer RT, Robinson AS: Optimization of the human adenosine A(2)a receptor yields in Saccharomyces cerevisiae. Biotechnology Progress. 2006, 22 (5): 1249-1255. 10.1002/bp050431r.
Valkonen M, Penttila M, Saloheimo M: Effects of inactivation and constitutive expression of the unfolded-protein response pathway on protein production in the yeast Saccharomyces cerevisiae. Applied and Environmental Microbiology. 2003, 69 (4): 2065-2072. 10.1128/AEM.69.4.2065-2072.2003.
Shusta EV, Raines RT, Pluckthun A, Wittrup KD: Increasing the secretory capacity of Saccharomyces cerevisiae for production of single-chain antibody fragments. Nature Biotechnology. 1998, 16 (8): 773-777. 10.1038/nbt0898-773.
Ellgaard L, Molinari M, Helenius A: Setting the standards: Quality control in the secretory pathway. Science. 1999, 286 (5446): 1882-1888. 10.1126/science.286.5446.1882.
Bhamidipati A, Denic V, Quan EM, Weissman JS: Exploration of the topological requirements of ERAD identifies Yos9p as a lectin sensor of misfolded glycoproteins in the ER lumen. Molecular Cell. 2005, 19 (6): 741-751. 10.1016/j.molcel.2005.07.027.
Tu BP, Ho-Schleyer SC, Travers KJ, Weissman JS: Biochemical basis of oxidative protein folding in the endoplasmic reticulum. Science. 2000, 290 (5496): 1571-1574.
Frand AR, Cuozzo JW, Kaiser CA: Pathways for protein disulphide bond formation. Trends in Cell Biology. 2000, 10 (5): 203-210. 10.1016/S0962-8924(00)01745-1.
Tu BP, Weissman JS: The FAD- and O-2-dependent reaction cycle of Ero1-mediated oxidative protein folding in the endoplasmic reticulum. Molecular Cell. 2002, 10 (5): 983-994. 10.1016/S1097-2765(02)00696-2.
Kulp MS, Frickel EM, Ellgaard L, Weissman JS: Domain architecture of protein-disulfide isomerase facilitates its dual role as an oxidase and an isomerase in Ero1p-mediated disulfide formation. Journal of Biological Chemistry. 2006, 281 (2): 876-884.
Cuozzo JW, Kaiser CA: Competition between glutathione and protein thiols for disulphide-bond formation. Nature Cell Biology. 1999, 1 (3): 130-135. 10.1038/11047.
Sevier CS, Kaiser CA: Ero1 and redox homeostasis in the endoplasmic reticulum. Biochimica Et Biophysica Acta-Molecular Cell Research. 2008, 1783 (4): 549-556. 10.1016/j.bbamcr.2007.12.011.
Chakravarthi S, Jessop CE, Bulleid NJ: The role of glutathione in disulphide bond formation and endoplasmic-reticulum-generated oxidative stress. EMBO Rep. 2006, 7 (3): 271-275. 10.1038/sj.embor.7400645.
Patil C, Walter P: Intracellular signaling from the endoplasmic reticulum to the nucleus: the unfolded protein response in yeast and mammals. Current Opinion in Cell Biology. 2001, 13 (3): 349-356. 10.1016/S0955-0674(00)00219-2.
Kimata Y, Ishiwata-Kimata Y, Yamada S, Kohno K: Yeast unfolded protein response pathway regulates expression of genes for anti-oxidative stress and for cell surface proteins. Genes to Cells. 2006, 11 (1): 59-69.
Kjeldsen T, Pettersson AF, Hach M: The role of leaders in intracellular transport and secretion of the insulin precursor in the yeast Saccharomyces cerevisiae. Journal of Biotechnology. 1999, 75 (2-3): 195-208. 10.1016/S0168-1656(99)00159-5.
Kjeldsen T: Yeast secretory expression of insulin precursors. Appl Microbiol Biotechnol. 2000, 54 (3): 277-286. 10.1007/s002530000402.
Nagashima T, Yamamoto Y, Gomi K, Kitamoto K, Kumagai C: A novel culture method for high level production of heterologous protein in Saccharomyces cerevisiae. Biosci Biotechnol Biochem. 1994, 58 (7): 1292-1296. 10.1271/bbb.58.1292.
Forster J, Gombert AK, Nielsen J: A functional genomics approach using metabolomics and in silico pathway analysis. Biotechnology and Bioengineering. 2002, 79 (7): 703-712. 10.1002/bit.10378.
Blank LM, Sauer U: TCA cycle activity in Saccharomyces cerevisiae is a function of the environmentally determined specific growth and glucose uptake rates. Microbiology. 2004, 150 (Pt 4): 1085-1093.
Oliveira A, Patil K, Nielsen J: Architecture of transcriptional regulatory circuits is knitted over the topology of bio-molecular interaction networks. BMC Systems Biology. 2008, 2 (1): 17-10.1186/1752-0509-2-17.
Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW, Hannett NM, Tagne JB, Reynolds DB, Yoo J, et al: Transcriptional regulatory code of a eukaryotic genome. Nature. 2004, 431 (7004): 99-104. 10.1038/nature02800.
Riemer J, Bulleid N, Herrmann JM: Disulfide Formation in the ER and Mitochondria: Two Solutions to a Common Process. Science. 2009, 324 (5932): 1284-1287. 10.1126/science.1170653.
Levin DE: Cell wall integrity signaling in Saccharomyces cerevisiae. Microbiology and Molecular Biology Reviews. 2005, 69 (2): 262-+. 10.1128/MMBR.69.2.262-291.2005.
Carla Fama M, Raden D, Zacchi N, Lemos DR, Robinson AS, Silberstein S: The Saccharomyces cerevisiae YFR041C/ERJ5 gene encoding a type I membrane protein with a J domain is required to preserve the folding capacity of the endoplasmic reticulum. Biochimica Et Biophysica Acta-Molecular Cell Research. 2007, 1773 (2): 232-242. 10.1016/j.bbamcr.2006.10.011.
Nakamata K, Kurita T, Bhuiyan MSA, Sato K, Noda Y, Yoda K: KEG1/YFR042w encodes a novel Kre6-binding endoplasmic reticulum membrane protein responsible for beta-1,6-glucan synthesis in Saccharomyces cerevisiae. Journal of Biological Chemistry. 2007, 282 (47): 34315-34324. 10.1074/jbc.M706486200.
Gaynor EC, Mondesert G, Grimme SJ, Reed SI, Orlean P, Emr SD: MCD4 encodes a conserved endoplasmic reticulum membrane protein essential for glycosylphosphatidylinositol anchor synthesis in yeast. Molecular Biology of the Cell. 1999, 10 (3): 627-648.
Laboissière M, Sturley S, Raines R: The essential function of protein-disulfide isomerase is to unscramble non-native disulfide bonds. Journal of Biological Chemistry. 1995, 270 (47): 28006-28009. 10.1074/jbc.270.47.28006.
Piruat JI, Aguilera A: A novel yeast gene, THO2, is involved in RNA pol II transcription and provides new evidence for transcriptional elongation-associated recombination. Embo Journal. 1998, 17 (16): 4859-4872. 10.1093/emboj/17.16.4859.
Das B, Butler JS, Sherman F: Degradation of normal mRNA in the nucleus of Saccharomyces cerevisiae. Molecular and Cellular Biology. 2003, 23 (16): 5502-5515. 10.1128/MCB.23.16.5502-5515.2003.
Rudra D, Mallick J, Zhao Y, Warner JR: Potential interface between ribosomal protein production and pre-rRNA processing. Molecular and Cellular Biology. 2007, 27 (13): 4815-4824. 10.1128/MCB.02062-06.
Madhani HD, Styles CA, Fink GR: MAP kinases with distinct inhibitory functions impart signaling specificity during yeast differentiation. Cell. 1997, 91 (5): 673-684. 10.1016/S0092-8674(00)80454-7.
Erdman S, Lin L, Malczynski M, Snyder M: Pheromone-regulated genes required for yeast mating differentiation. Journal of Cell Biology. 1998, 140 (3): 461-483. 10.1083/jcb.140.3.461.
Schroder M, Chang JS, Kaufman RJ: The unfolded protein response represses nitrogen-starvation induced developmental differentiation in yeast. Genes & Development. 2000, 14 (23): 2962-2975. 10.1101/gad.852300.
Lambrechts MG, Bauer FF, Marmur J, Pretorius IS: Muc1, a mucin-like protein that is regulated by Mss10, is critical for pseudohyphal differentiation in yeast. Proceedings of the National Academy of Sciences of the United States of America. 1996, 93 (16): 8419-8424. 10.1073/pnas.93.16.8419.
Palecek SP, Parikh AS, Kron SJ: Genetic analysis reveals that FLO11 upregulation and cell polarization independently regulate invasive growth in Saccharomyces cerevisiae. Genetics. 2000, 156 (3): 1005-1023.
Tao W, Deschenes RJ, Fassler JS: Intracellular glycerol levels modulate the activity of Sln1p, a Saccharomyces cerevisiae two-component regulator. Journal of Biological Chemistry. 1999, 274 (1): 360-367. 10.1074/jbc.274.1.360.
Zhang BY, Chang A, Kjeldsen TB, Arvan P: Intracellular retention of newly synthesized insulin in yeast is caused by endoproteolytic processing in the golgi complex. Journal of Cell Biology. 2001, 153 (6): 1187-1197. 10.1083/jcb.153.6.1187.
Huang DG, Gore PR, Shusta EV: Increasing Yeast Secretion of Heterologous Proteins by Regulating Expression Rates and Post-Secretory Loss. Biotechnology and Bioengineering. 2008, 101 (6): 1264-1275. 10.1002/bit.22019.
Peng RW, Guetg C, Tigges M, Fussenegger M: The vesicle-trafficking protein munc18b increases the secretory capacity of mammalian cells. Metabolic Engineering. 2010, 12 (1): 18-25. 10.1016/j.ymben.2009.08.007.
Peng RW, Fussenegger M: Molecular Engineering of Exocytic Vesicle Traffic Enhances the Productivity of Chinese Hamster Ovary Cells. Biotechnology and Bioengineering. 2009, 102 (4): 1170-1181. 10.1002/bit.22141.
Ikner A, Shiozaki K: Yeast signaling pathways in the oxidative stress response. Mutation Research-Fundamental and Molecular Mechanisms of Mutagenesis. 2005, 569 (1-2): 13-27. 10.1016/j.mrfmmm.2004.09.006.
Li S, Dean S, Li ZJ, Horecka J, Deschenes RJ, Fassler JS: The eukaryotic two-component histidine kinase Sln1p regulates OCH1 via the transcription factor, Skn7p. Molecular Biology of the Cell. 2002, 13 (2): 412-424. 10.1091/mbc.01-09-0434.
Hatahet F, Ruddock LW: Protein Disulfide Isomerase: A Critical Evaluation of Its Function in Disulfide Bond Formation. Antioxidants & Redox Signaling. 2009, 11 (11): 2807-2850. 10.1089/ars.2009.2466.
Wittrup K, Robinson A, Parekh R, Forrester K: Existence of an Optimum Expression Level for Secretion of Foreign Proteins in Yeast. Annals of the New York Academy of Sciences. 1994, 745 (1): 321-330.
Robinson AS, Bockhaus JA, Voegler AC, Wittrup KD: Reduction of BiP levels decreases heterologous protein secretion in Saccharomyces cerevisiae. Journal of Biological Chemistry. 1996, 271 (17): 10017-10022. 10.1074/jbc.271.17.10017.
van Dijken JP, Bauer J, Brambilla L, Duboc P, Francois JM, Gancedo C, Giuseppin MLF, Heijnen JJ, Hoare M, Lange HC, et al: An interlaboratory comparison of physiological and genetic properties of four Saccharomyces cerevisiae strains. Enzyme and Microbial Technology. 2000, 26 (9-10): 706-714. 10.1016/S0141-0229(00)00162-9.
Sambrook J, Fritsch EF, Maniatis T: Molecular Cloning: A Laboratory Manual. 1989, Cold Springs, NY, USA: Cold Spring Harbor Laboratory Press, 2
Gietz RD, Woods RA: Transformation of yeast by lithium acetate/single-stranded carrier DNA/polyethylene glycol method. Guide to Yeast Genetics and Molecular and Cell Biology, Pt B. 2002, 350: 87-96.
Moller K, Sharif MZ, Olsson L: Production of fungal alpha-amylase by Saccharomyces kluyveri in glucose-limited cultivations. J Biotechnol. 2004, 111 (3): 311-318. 10.1016/j.jbiotec.2004.04.013.
Pont-Kingdon G: Creation of Chimeric Juctions, Deletions, and insertions by PCR. PCR protocols. Edited by: Bartlett JMS, Stirling D. 2003, Totowa, NJ, USA: Humana Press Inc., 511-515. 2
Mumberg D, Müller R, Funk M: Yeast vectors for the controlled expression of heterologous proteins in different genetic backgrounds. Gene. 1995, 156 (1): 119-122. 10.1016/0378-1119(95)00037-7.
Wittrup KD, Benig V: Optimization of Amino-Acid Supplements for Heterologous Protein Secretion in Saccharomyces-Cerevisiae. Biotechnology Techniques. 1994, 8 (3): 161-166. 10.1007/BF00161582.
Snel L, Damgaard U: Proinsulin heterogeneity in pigs. Horm Metab Res. 1988, 20 (8): 476-480. 10.1055/s-2007-1010862.
Usaite R, Patil KR, Grotkjaer T, Nielsen J, Regenberg B: Global transcriptional and physiological responses of Saccharomyces cerevisiae to ammonium, L-alanine, or L-glutamine limitation. Applied and Environmental Microbiology. 2006, 72 (9): 6194-6203. 10.1128/AEM.00548-06.
Salazar M, Vongsangnak W, Panagiotou G, Andersen MR, Nielsen J: Uncovering transcriptional regulation of glycerol metabolism in Aspergilli through genome-wide gene expression data analysis. Molecular Genetics and Genomics. 2009, 282 (6): 571-586. 10.1007/s00438-009-0486-y.
Stephanopoulos GN, Aristidou AA, Nielsen J: Material Balances and Data Consistency. Metabolic Engineering: Principles and Methodologies. 1998, New York, New York: Academic Press
We thank Dr Intawat Nookaew for assistance with transcriptome analysis and FBA. We thank NIH F32 Kirschstein NRSA fellowship (F32 GM083647), The Knut and Alice Wallenberg Foundation, EU Framework VII project SYSINBIO (Grant no. 212766), European Research Council project INSYSBIO (Grant no. 247013), and the Chalmers Foundation for funding.
The authors declare that they have no competing interests.
KT, DP, and JN designed the experiment. KT and ZL carried out all cloning, fermentations, and analytical measurements. KT did primary calculations in transcriptomics and metabolic flux data. KT, ZL, DP, and JN analyzed data and wrote the manuscript. DP and JN supervised the research. All authors have read and approved of the final manuscript.