Burden Imposed by Heterologous Protein Production in Two Major Industrial Yeast Cell Factories: Identifying Sources and Mitigation Strategies

Production of heterologous proteins, especially biopharmaceuticals and industrial enzymes, in living cell factories consumes cellular resources. Such resources are reallocated from normal cellular processes toward production of the heterologous protein that is often of no benefit to the host cell. This competition for resources is a burden to host cells, has a negative impact on cell fitness, and may consequently trigger stress responses. Importantly, this often causes a reduction in final protein titers. Engineering strategies to generate more burden resilient production strains offer sustainable opportunities to increase production and profitability for this growing billion-dollar global industry. We review recently reported impacts of burden derived from resource competition in two commonly used protein-producing yeast cell factories: Saccharomyces cerevisiae and Komagataella phaffii (syn. Pichia pastoris). We dissect possible sources of burden in these organisms, from aspects related to genetic engineering to protein translation and export of soluble protein. We also summarize advances as well as challenges for cell factory design to mitigate burden and increase overall heterologous protein production from metabolic engineering, systems biology, and synthetic biology perspectives. Lastly, future profiling and engineering strategies are highlighted that may lead to constructing robust burden-resistant cell factories. This includes incorporation of systems-level data into mathematical models for rational design and engineering dynamical regulation circuits in production strains.

Introduction

Biomanufacturing of heterologous proteins from genetically engineered cell factories is a growing industry. The first such biopharmaceutical product was approved by the FDA four decades ago (U.S. Food Drug Administration, 1982; Nielsen, 2013), and now over 300 different biopharmaceuticals are available (Walsh, 2018). Beyond delivering many life-saving drugs, biopharmaceutical protein production serves as an economically profitable industry representing an expanding multi-billion-dollar global market (Walsh, 2018; Mordor Intelligence, 2021). Apart from this, commercial protein products for industrial and research applications are also produced in engineered cells and likewise constitute a growing industry (Baghban et al., 2019). However, despite many previous and ongoing successes, continued bioengineering of cell factories is required for more robust and reliable production of heterologous protein.

Today, an assorted palette of living organisms with distinct advantages and disadvantages are used for industrial-scale heterologous protein production, including bacteria, fungi, and mammalian cell lines (Dumont et al., 2016). Depending on the protein product of interest, developers choose the most appropriate host organism based on quality and quantity (Porro et al., 2011). Yeast represents an attractive host for the production of many types of heterologous proteins when compared to other host systems.

The first heterologous protein product released to the market was human insulin produced in an engineered strain of Escherichia coli using recombinant DNA technology and branded as Humulin™ (U.S. Food Drug Administration, 1982). Since then, bacterial expression systems have been a popular choice for producing multiple heterologous proteins since they replicate quickly in cheap growth media and, in general, produce high protein titers (Karbalaei et al., 2020). On the other hand, they lack the capacity to perform necessary post-translational modifications required, especially for many pharmaceutical proteins to be active (Ghaderi et al., 2012). Most bacteria, including E. coli, are also unable to successfully secrete such proteins, requiring added cell harvesting, cell disruption, and product isolation steps before protein purification (Lalor et al., 2019). Thus, protein recovery from lysed bacteria cells often requires time-consuming and costly downstream processing (Vieira Gomes et al., 2018). Alternatively, mammalian cell lines can be used to achieve proper protein modifications, folding, and secretion for simpler downstream purification (Dumont et al., 2016). Although several valuable pharmaceutical proteins are expressed from mammalian cell platforms, including therapeutic monoclonal antibodies (Kunert and Reinhart, 2016; Walsh, 2018), these cells grow slower than microbial cell factories and are much more expensive to cultivate (Karbalaei et al., 2020).

Since yeast are unicellular eukaryotes, these fungi require less expensive culturing conditions and possess many of the post-translational and secretion pathways present in higher eukaryotes. As such, yeast cell factories are commonly used to produce diverse heterologous proteins for pharmaceutical, industrial, and research applications (Table 1, Baghban et al., 2019; Kulagina et al., 2021). The two most favored yeast production hosts are Saccharomyces cerevisiae and Komagataella phaffii (syn. Pichia pastoris) (Tripathi and Shrivastava, 2019) and will be the focus of this review.

TABLE 1

Table 1. Examples of heterologous proteins produced in Saccharomyces cerevisiae (Sc) and Komagataella phaffii (Kp).

Despite many commercial success stories, challenges still remain when engineering heterologous protein production in yeast. Yeast cells have evolved to readily adapt to fluctuating environmental conditions, and may down-regulate protein production to match intracellular demands imposed by various external and internal conditions and stresses (Gasch, 2003). For production strains, this coordination is important to ensure a competitive trade-off between fitness and cost-intensive protein production which consumes cellular resources. Thus, resource intensive protein production can easily disrupt this delicate intracellular balance due to the redirection of cellular resources that are normally distributed among native cellular activities, including biomass formation and growth (Heyland et al., 2011b). Such production burden (Glick, 1995; Heyland et al., 2011b; de Ruijter et al., 2018), typically results in reduced protein titers, limiting process efficiency and production profitability.

Burden is not the only limitation that may negatively impact production of protein in yeast cells. Product toxicity, caused by protein aggregation for example, can also have an impact on the health and productivity of yeast cell factories (Eguchi et al., 2018; Chen et al., 2020; Romero-Suarez et al., 2021). However, not all heterologous proteins have toxic effects on host cells. In contrast, burden represents a universal challenge for heterologous protein production since it incurs resource demands for the host cells at potentially all levels of gene expression.

In this review, we will focus on burden conferred by heterologous soluble protein production in yeast cell factories from recombinant DNA replication through to the final secreted product. Since achieving desirable heterologous protein production is paramount, we will highlight recent strategies to mitigate burden in two of the most common yeast systems used for this purpose and propose new paths for future research. Although many concepts discussed here will be relevant to burden in other production organisms, the unique challenges presented by non-yeast systems remain outside the scope of this review.

Two Favored Budding Yeast Production Strains

Saccharomyces cerevisiae

Several yeasts are used as host organisms for protein production, including S. cerevisiae, Kluyveromyces lactis, Yarrowia lipolytica, Hansenula polymorpha, Ogataea polymorpha, and K. phaffii (Rebello et al., 2018; Vieira Gomes et al., 2018; Kulagina et al., 2021). The budding yeast S. cerevisiae is among the most well-established organisms for heterologous protein production. Often called a “conventional yeast,” S. cerevisiae is highly domesticated and has been used by humans for thousands of years to produce bread and alcoholic beverages (Mattanovich et al., 2014; Patra et al., 2021). The extensive use of S. cerevisiae has facilitated its establishment as an important model organism for life sciences (Chen et al., 2020). Notably, it was the first eukaryotic organism to have its genome sequenced (Goffeau et al., 1996). Subsequent years of research have resulted in a well-annotated genome and an extensive genetic toolbox for S. cerevisiae production strain engineering.

Beyond genetic engineering, S. cerevisiae possesses many additional advantages for heterologous protein production. It possesses eukaryotic pathways to post-translationally process and secrete proteins to the extracellular medium, while simultaneously secreting few endogenous proteins in low amounts. S. cerevisiae can tolerate harsh growth conditions, such as low pH, and finally, it has gained the generally recognized as safe (GRAS) status allowing for easier approval processes (U.S. Food Drug Administration, 2018). Indeed, For these reasons, S. cerevisiae has become an established host for production of heterologous proteins such as human insulin, glucagon, and lignocellulosic enzymes (Table 1, Baghban et al., 2019; Kulagina et al., 2021).

Komagataella phaffii (syn. Pichia pastoris)

Other non-conventional yeasts have been explored for their suitability in heterologous protein production. An attractive alternative budding yeast is K. phaffii. Ecallantide, a drug used to treat hereditary angioedema, was the first pharmaceutic protein produced and approved in K. phaffii in 2009. It has also become established for industrial-scale protein production and manufacturing of insulin, human serum albumin, antibody fragments, phytases, and trypsin (Research Corporation Technologies, 2009, 2019). K. phaffii possesses many of the same characteristic advantages mentioned for S. cerevisiae, such as GRAS status (Vogl et al., 2013). Compared to S. cerevisiae, it can grow to higher cell densities and has a greater secretory capacity (Duman-Özdamar and Binay, 2021). Another important difference involves lower glucose uptake in K. phaffii compared to S. cerevisiae (Peña et al., 2018). K. phaffii is a crab tree-negative yeast, which can be an advantage over the crab tree-positive yeasts, e.g., S. cerevisiae, that can produce toxic levels of ethanol during aerobic cultivations (Cereghino and Cregg, 2000).

Another attractive feature of K. phaffii is its methylotrophic nature. It thrives on simple carbon sources like methanol, which is otherwise toxic to many other microorganisms including S. cerevisiae (Riley et al., 2016). Genetic engineering often exploits this unique feature by placing genes for desired heterologous proteins under the control of a strong methanol-inducible promoter from the Alcohol oxidase 1 gene (AOX1). This allows for a biomass formation phase prior to an induced expression production phase, where the medium carbon source is switched from glycerol, glucose, or sorbitol to methanol (Celik et al., 2009; Noseda et al., 2013). Challenges such as inducer-toxicity, fire hazard, and heat production are related to the use of methanol, present concerns for large-scale production and have caused researchers to explore the engineering of methanol-free expression systems in K. phaffii (Wang et al., 2017). Despite its strengths, K. phaffii is less established and genetically characterized, meaning that metabolic and genome engineering of K. phaffii lags behind that of S. cerevisiae (Kalender and Çalik, 2020; Duman-Özdamar and Binay, 2021). However, genome sequences and annotations are available for the two most relevant K. phaffii strains GS115 and CBS7435 (De Schutter et al., 2009; Küberl et al., 2011; Love et al., 2016; Valli et al., 2016) and substantial progress in genetic engineering has in recent years enhanced heterologous protein production (Yang and Zhang, 2018). For example, employing the constitutive glyceraldehyde-3-phosphate dehydrogenase (GAP) promoter (Qin et al., 2011). Overall, continued advances help overcome important initial drawbacks of using K. phaffii as a production host and underpin the relevance of using it for industrial-scale heterologous protein production.

Common Challenges

Both of these yeast cell factories face similar challenges. For example, the challenge of N-glycan hyper-mannosylation of heterologous protein in S. cerevisiae, although less pronounced in K. phaffii (Mizukami et al., 2018; Vieira Gomes et al., 2018), remains problematic since it can render certain pharmaceutical proteins inactive, such as monoclonal antibodies (Jung and Kim, 2018). Yet, there are still opportunities to use yeast cell factories to produce some glycoprotein products. This requires either additional biochemical processing steps before they can be used as pharmaceutical proteins, or further glycoengineering of production strains with more humanized glycosylation patterns (Liu et al., 2018; Kulagina et al., 2021).

Together, these two budding yeasts represent highly attractive and cost-effective platforms for producing a plethora of different therapeutic and industrial proteins and ongoing research continues to enhance their performance. Despite their advantages, burden conferred by heterologous protein production remains a challenge encountered in all production hosts, and yeasts are no exception (Heyland et al., 2011a,b; Niklas et al., 2013; de Ruijter et al., 2018; Zou et al., 2018). Like the challenge of N-glycan hyper-mannosylation, we expect genome engineering will be critical for building more burden resilient production strains.

Burden—The Competition for Cellular Resources

Heterologous protein production in yeast consumes a variety of cellular resources for expression, processing, and transport. This includes consumption of metabolic precursors, redox co-factors, and energy sources (Mattanovich et al., 2014; Klein et al., 2015; Zahrl et al., 2019). Such metabolic remodeling impacts cell growth and limits normal cellular metabolic processes, often curbing heterologous protein production (Glick, 1995; Kazemi Seresht et al., 2013). This reallocation was first defined as metabolic burden in bacteria (Glick, 1995) but the nomenclature describing the phenomenon in yeast varies from metabolic burden to protein burden or fitness burden (Harrison et al., 2012; Kafri et al., 2016; Deparis et al., 2017; Huang et al., 2017; Saeki et al., 2020; Garrigós-Martínez et al., 2021; Wright et al., 2021). In this review, we will refer to this concept simply as “burden.”

To what extent burden affects the cell is case-dependent, since heterologous proteins of differing size and biochemical complexity demand different resources from the host cell (Heyland et al., 2011b; de Ruijter et al., 2018). Therefore, the competition for resources varies depending on several factors, such as amino acid composition, post-translational modifications, metabolism of the host, and carbon availability (Heyland et al., 2011b; Tyo et al., 2012). Differing expression levels for any given heterologous protein can also have a difficult-to-predict impact on production. For example, Mellitzer et al. described different classes of heterologous proteins based on how gene dosage affected expression yield. They found that increasing gene copy number for some heterologous proteins resulted in higher active protein yields, while for others yields remained constant or even decreased (Mellitzer et al., 2012, 2014). Revealing the underlying mechanisms responsible for these differences would provide critical information about how burden is established in different production strains. In line with this, different engineering strategies will likely be needed to mitigate burden depending on the protein of interest (Gu et al., 2015). But first we must be able to assess the ways in which cells are burdened by heterologous protein expression.

Assessing burden imposed by heterologous protein production is crucial for engineering more resilient and higher-producing strains. But how is burden measured? Widely used approaches involve analyzing physiological parameters like cell growth rate, biomass yield, and respiratory capacity (Kazemi Seresht et al., 2013; Liu et al., 2014; de Ruijter et al., 2018). Apart from monitoring cells through these parameters, additional methods and -omics technologies have been applied in recent years to assess burden and its sources (Heyland et al., 2011b; Jordà et al., 2012; Liu et al., 2016; Huang et al., 2017; Wright et al., 2020). These methods provide measurements of the internal cellular state, including, but not limited to, metabolites, transcriptional regulation, translational efficiency, and carbon flux. Together, these approaches help describe the observed burdened phenotype.

Figure 1 illustrates critical steps for the production of heterologous proteins in yeast cell factories and highlights engineering targets for mitigating burden. In the following sections, we will review how burden manifests in production cells from heterologous DNA to secreted protein product.

FIGURE 1

Figure 1. Burden-triggering bottlenecks in yeast cell factories during heterologous protein synthesis and secretion. Schematic representation of yeast cell factory engineered to synthesize and secrete a heterologous protein of interest (POI), encoded by a gene of interest (GOI). The GOI is transcribed into messenger RNA (mRNA), later translated by ribosomes into a peptide that is translocated into endoplasmic reticulum (ER), where it is folded and modified. Via anterograde transport, the POI reaches the golgi complex for further modifications and secretion signal cleavage. Via vesicle exocytosis, the POI is secreted to the extracellular medium. Meanwhile, heterologous protein production is fueled by metabolic pathways such as glycolysis, pentose phosphate pathway (PPP), the tricarboxylic acid cycle (TCA), and the electron transport chain (ETC) which deliver energy (ATP), redox cofactors (NADPH, NADH), nucleotides, and amino acids required for heterologous protein production. Different engineering targets (ETs) that have been highlighted in literature to mitigate the burden that heterologous protein production imposes on the host cell are marked with a target symbol. ET1: Gene dosage, promoter strength, plasmid vs. genome integration. ET2: Codon optimization, co- or post-translational translocation. ET3: Tuning genes involved in unfolded protein response, oxidative stress response, ER-associated protein degradation pathways. ET4: Direct POI for secretion or express intracellularly. ET5: Process engineering by adding certain amino acids to the medium, changing the medium or carbon source. ET6: Redirecting metabolic fluxes by tuning relevant target genes. Synthetic circuit engineering for dynamic regulation may be simultaneously aimed at any combination of ETs 1–4 and/or 6.

The Impacts of Heterologous DNA on Production Strains

Heterologous DNA is maintained in host cells on plasmids or as genome-integrated expression cassettes. Although the energetic cost of replicating extra plasmid DNA is deemed negligible in haploid S. cerevisiae (Krogh et al., 2008; Rugbjerg and Sommer, 2019), researchers studying the instability of plasmid copy number in S. cerevisiae in both haploid and diploid strains producing heterologous aprotinin found that diploid cells adapted to burden through plasmid loss (Krogh et al., 2008). In contrast, despite an observed decrease in aprotinin production, haploid strains maintained constant plasmid copy number while adapting to extended growth in minimal medium (Krogh et al., 2008). We expect this to be similar in K. phaffii but to our knowledge, no research is available in this organism on this topic. However, genomic integration is a more common approach in K. phaffii, which we will get back to. Maintaining plasmids, often requires selection marker expression, such as genes encoding for antibiotic resistance or metabolic enzymes, the expression of such drains even more resources from cells and adds to burden (Karim et al., 2013; He et al., 2020; Shen et al., 2021). Therefore, chromosomal integration of heterologous genes without the need for selection markers can help mitigate the fitness cost associated with carrying high-copy number plasmids (Harrison et al., 2012; Vieira Gomes et al., 2018). Alternatively, a weaker promoter driving expression of the selection marker can be utilized, which can help relieve burden as recently reported in K. phaffii (Shen et al., 2021). In this case, a 3.3-fold higher plasmid copy number and significantly higher expression levels of insulin precursor were achieved when employing the 300 bp promoter region upstream 2-deoxyglucose-6-phosphate phosphatase gene (P_Dog2p300), compared to the popular TEF1 promoter (Shen et al., 2021).

In S. cerevisiae production strains, plasmid-based expression and chromosomal integration strategies are both widely employed for heterologous protein production. For K. phaffii, chromosomal integration of expression cassettes is preferred since plasmids are much less stable in this species. However, high-copy integration strains can be genetically unstable, likely as a consequence of increased burden from heterologous protein production (Zhu et al., 2009; Yu et al., 2020). In rare instances, loss of gene copies through recombination events in genetically unstable cells can occur, providing a growth advantage over burdened cells with higher heterologous protein expressions levels (Zhu et al., 2009). Thus, chromosomal integrations provide many advantages for heterologous protein expression in both S. cerevisiae and K. phaffii, but they also present potential challenges if copies of identical sequences or similar genetic elements are used.

The primary biosynthetic cost from high-copy plasmids or chromosomal integrations encoding heterologous genes is protein overexpression rather than the cost of replication (Kim et al., 2001; Harrison et al., 2012; Eguchi et al., 2018; Huang et al., 2018; He et al., 2020; Yu et al., 2020). Consistent with this, a non-linear correlation between protein production and gene dosage is often observed in yeast production strains (Mellitzer et al., 2014; Cámara et al., 2016; Fang et al., 2017; Huang et al., 2018; He et al., 2020). The optimal threshold for gene dosage varies between studies and depends on the protein of interest (Mellitzer et al., 2012). For example, a two-copy strain of K. phaffii producing Rhizopus oryzae lipase performed better (12.73UA-10¹²/cell) than strains with more than two copies (8.19-12.09UA-10¹²/cell) (Cámara et al., 2016). In another study, a 10-copy strain of K. phaffii producing 1.2 g/L Thermomyces lanuginosus xylanase A performed better than 6- and 18-copy strains producing ~0.85 and 0.25 g/L, respectively (Mellitzer et al., 2012). In a third study, K. phaffii production of the porcine insulin precursor reached a maximum of 0.9 g/L at 6 copies, showing approximately the same production titers for 6, 12, and 18 copies (Chen et al., 2017). These findings indicate that cells possess mechanisms for downregulating heterologous genes when reaching a certain gene dosage plateau, a Goldilocks range of maximum expression, above which burden is triggered (Mellitzer et al., 2014).

Strategies for delivering and maintaining heterologous DNA in production strains, as discussed above, correspond to engineering target 1 in Figure 1. Overall, plasmids or integrated genes in high copy number present potential sources of strain instability. Heterologous DNA maintenance also clearly requires cellular resources from the host. However, the biosynthetic cost of protein production from such foreign DNA has a significantly larger impact on triggering a burden state in the production host. This can be divided into the different steps linked to metabolic pathways, which we will discuss in the following sections.

Transcription of Heterologous Genes: A Source of Potential Burden

The biosynthetic process of transcription has been reported to be a limiting process for growing S. cerevisiae and K. phaffii under certain conditions, including in phosphate limiting conditions or when grown in standard medium with increasing gene-copy number (Kafri et al., 2016; Cámara et al., 2017; Farkas et al., 2018). Engineering target 1 in Figure 1 also highlights the potential strategies for controlling optimal heterologous gene expression. It is conceivable that burden could be induced when employing strong promoters due to competition from the heterologous gene promoter for essential transcription factors against native promoters, such as the K. phaffii transcription factor Mxr1p responsible for stimulating expression from the AOX1 promoter (Görgens et al., 2001; Cámara et al., 2017). It is equally possible that high expression levels of the heterologous gene drain other transcriptional resources, such as the transcriptional machinery or ribonucleotides, thereby reducing growth rates and biomass yield (Görgens et al., 2001; Cámara et al., 2017; Rugbjerg and Sommer, 2019). Having said that, eukaryotic genomes are widely transcribed, including many rapidly degraded non-coding RNAs from intergenic and anti-sense regions (Jensen et al., 2013), implying additional transcription from a heterologous gene present in one or more copies should have a negligible impact (cellular protein content is 5–10-fold higher than nucleic acid content). Instead, burden from heterologous genes and their strong transcription is more likely indirect since the mRNAs produced might overload costly downstream processes involved in protein synthesis (Saeki et al., 2020). Nevertheless, increased transcriptional activity of recombinant genes likely contributes, however modestly, to burden (Farkas et al., 2018).

Protein Synthesis and tRNA Availability Place a Burden on Production Strains

Protein translation is a costly process consuming energy, nutrients, and other resources. Overall, this step is reported to be an important limiting process for protein production in yeast (Kafri et al., 2016; Cámara et al., 2017). Although different resources are critical during translation, the most obvious components for this process include ribosomes and amino acyl tRNAs. Intuitively, high expression levels from heterologous genes leaves fewer ribosomes available to translate native proteins, eventually leading to growth defects (Saeki et al., 2020). However, multiple studies in yeast suggest that the major limiting resources during translation are in fact metabolic products like energy and metabolic precursors, with the occupation of the translational machinery representing a minor source of burden (Heyland et al., 2011b; Kafri et al., 2016). Thus, identifying and enhancing the availability of limiting resources for protein synthesis is expected to reduce burden and enhance overall protein production. This critical step corresponds to engineering target 2 in Figure 1.

The availability of certain tRNAs can also be limited during translation. Aminoacyl-tRNAs formation consumes considerable amounts of cellular energy through hydrolyzation of ATP molecules, while GTP molecules are also subsequently required for amino acid polymerization (Schimmel, 1993). We will discuss causative factors of energy limitations later in this review. The copy number for genes encoding different tRNAs controls to some extent the availability of tRNAs for translation. tRNA gene copy numbers encoding specific anti-codons varies almost 20-fold in S. cerevisiae (Chan and Lowe, 2016) and 10-fold in K. phaffii (De Schutter et al., 2009). Thus, codon optimizing a heterologous gene is expected to be important to avoid the depletion of specific tRNAs, especially rare ones. Based on differing tRNA gene copy numbers in S. cerevisiae and K. phaffii, it is however plausible that cognate tRNA availability is a more important consideration for burden mitigation in S. cerevisiae production strains. Indeed, tRNA availability and recharging has been suggested to be a non-limiting factor in burdened K. phaffii (Mellitzer et al., 2014). Mellitzer et al. applied different codon optimization strategies and observed that certain genes resulted in higher titers independent of the promoter used to drive transcription and carbon source used, conditions under which the aminacyl-tRNA pool was otherwise expected to change and thus have an influence on what optimized genes would result in the highest protein titer or burden the translation machinery. Changes in tRNA concentrations can impact cellular mRNA expression profiles by altering mRNA stability (Presnyak et al., 2015). With this parameter also comes translational speed, which can indirectly contribute to burden since fast translation kinetics can cause incorrect protein folding triggering cellular stress responses (Drummond and Wilke, 2008; Yang et al., 2014). For example, slowing translation through the use of less “optimal” codons can help correct protein folding for some difficult to express heterologous proteins (Mellitzer et al., 2014; Yang et al., 2014). On the other hand, Mellitzer et al. observed a 23-fold difference in expression levels for cellobiohydrolase 2 with different synonymous codon changes in K. phaffii (Mellitzer et al., 2014).

Notably, if codon optimality is around the average tRNA adaptation index score [tAI score developed by dos Reis et al. (2004)] of a gene in the specific host organism−0.37 in S. cerevisiae (Eguchi et al., 2018); 0.47 in K. phaffii (Xu et al., 2021)—the translation elongation rate is often not high enough for an overexpressed protein to reach the burden-limit triggering growth defects. Authors using S. cerevisiae as a model organism estimate this to be around 15% of total cellular protein (Eguchi et al., 2018). Similar estimations for K. phaffii have not yet been experimentally determined or modeled. Whether tRNA abundance is a factor that should be considered in burdened yeast cell factories and the extent to which codon optimization is a useful strategy to significantly mitigate burden remains unanswered in most cases.

Amino Acid Supply as a Limiting Factor During Translation

As highlighted above, heterologous gene expression triggers an increased demand for nucleotides, tRNAs, and amino acids, removing resources from the pool of free intracellular building blocks required for mRNA and protein synthesis (Gonzalez et al., 2003; de Ruijter et al., 2018). By comparing S. cerevisiae metabolite profiles or flux distribution in production strains to wild type, an increased rate of consumed amino acids has been observed (Gonzalez et al., 2003; de Ruijter et al., 2018). This leaves fewer metabolic precursors for downstream pathways like the tricarboxylic acid cycle (TCA) cycle and pentose phosphate pathway (PPP), leading to decreased flux into energy and reducing power formation. This explains decreased growth rates observed and why supplementing cultivation medium with certain amino acids improves heterologous protein production (Görgens et al., 2005; Heyland et al., 2011b; Van Rensburg et al., 2012; Liu et al., 2016; Huang et al., 2017).

Glutamate and glutamine are important building blocks for nascent polypeptides and amino donors for the biosynthesis of other amino acids (Huang et al., 2017). By supplementing these two amino acids in growth medium, their increased availability helps relieve metabolic limitations by contributing also as additional carbon sources (Heyland et al., 2011b; Nie et al., 2014). Addition of glutamine alone to standard defined medium was reported to cause an increase in β-aminopeptidase 3-2W4 BapA from 0.37 to 0.42 mg/g (Heyland et al., 2011b), whereas in another experiment addition of glutamate caused an increase in β-galactosidase (in this case measured as activity) from 14,570 to 20,460 U/g_DCW (Nie et al., 2014). In line with this, RNA-seq revealed that genes involved in de novo biosynthesis of glutamate and glutamine were upregulated despite supplementation, as well as amino group transfer genes to help convert them into other amino acids and genes encoding membrane transporters (Huang et al., 2017). Thus, glutamine and glutamate especially might play a central role in adapting to the changing environment, but also amino acid composition of the heterologous protein produced might change the need and limitations in the burdened host cells. Huang and coworkers also found a significant upregulation of YCT1 and ERC1 gene transcription encoding cysteine transporters. In their study they produced a-amylase which amino acid composition has a 9.3-fold higher requirement for cysteine compared to that of the average yeast cell protein (Huang et al., 2017) indicating amino acid composition can affect limitations related to amino acids. On the other hand, the authors of another study, producing insulin precursor from S. cerevisiae, found that heterologous protein production had an impact on the cellular free amino acid pool, but that this had no correlation with the relative amino acid abundance in native yeast protein and insulin precursor protein (Kazemi Seresht et al., 2013). Thus, cells appear to adapt to changing demands from heterologous protein production by enhancing both amino acid uptake via membrane transporters as well as de novo synthesis.

Glutamate addition changes flux distributions of central carbon metabolism and increases the pool of intermediates in the TCA cycle to improve energy generation in K. phaffii (Nie et al., 2014; Liu et al., 2016). Notably, TCA cycle activity appears to react differently in K. phaffii and S. cerevisiae. For K. phaffii, an upper limit of absolute TCA cycle activity to catabolize carbon was reported (Heyland et al., 2011b). Thus, the capacity of the TCA cycle is not always robust enough to catabolize sufficient carbon to meet the increased demands for energy, which will limit growth rates and protein production in K. phaffi (Heyland et al., 2011b). On the contrary, the relative TCA cycle flux for S. cerevisiae correlates with glucose uptake (Blank and Sauer, 2004). Therefore, S. cerevisiae might not have such an upper limit as has been reported for K. phaffii, highlighting a potentially important metabolic difference to consider when engineering more resilient strains of each organism.

In summary, amino acids play an essential role in protein production as both building blocks as well as additional carbon sources for central carbon metabolism. Evidence supports supplementation of certain amino acids in growth medium to alleviate burden in yeast cell factories, as indicated for engineering strategy 5 in Figure 1. Particularly, glutamate and glutamine enhance protein production by relieving metabolic limitation as an added nitrogen and carbon source. However, due to the economic cost of using such supplementations in large-scale bioreactors (Heyland et al., 2011b), this process engineering strategy can be suboptimal at industry-level. Instead, genetic engineering strategies to overcome or bypass such metabolic limitations likely offer more cost-effective alternatives to mitigate burden at this stage in the protein production process.

NADPH Demands During Translation Contribute to Burden

High demand for amino acids can also contribute to establishing a burden state in production cells due to limitations in the redox cofactor nicotinamide adenine dinucleotide phosphate (NADPH). NADPH is required in many anabolic processes, including de novo synthesis of amino acids, and demand for this redox cofactor predictably rises in production strains (Zahrl et al., 2019).

NADPH can be generated through two main pathways: the oxidative phase of PPP or the acetate formation pathway (Grabowska and Chelstowska, 2003; Kwolek-Mirek et al., 2019). Primary generation of NADPH via PPP was reported for burdened K. phaffii production strains where acetate did not accumulate despite increased demands for NADPH (Nie et al., 2014). Therefore, upregulating genes in the oxidative phase of PPP might be expected to help cells meet the growing demand for NADPH. Such approaches are highlighted as engineering strategy 6 in Figure 1 and have been confirmed in K. phaffii strains producing human superoxide dismutase, whereby enhancing the first two steps in the oxidative phase of PPP proved to increase production by 3.8-fold (Nocon et al., 2016). However, it did not alleviate the observed decreases in specific growth rate and biomass yield. This potential solution might not work in burdened S. cerevisiae strains since RNA-seq and proteomic analyses show that genes involved in PPP are not upregulated, but rather downregulated (Huang et al., 2017; Wright et al., 2020). Instead, increased acetate formation has been reported from an insulin-producing S. cerevisiae strain, indicating preferential NADPH production via the acetate formation pathway to meet increased demand resulting from augmented amino acids needs (Wright et al., 2020). Others however argue that S. cerevisiae production strains meet NAPDH demands by simply reducing biomass yield in general (Huang et al., 2017).

In contrast to S. cerevisiae, researchers studying protein production in K. phaffii observed an increase in biomass yield (from 0.39 to 0.49 g_CDW${\text{g}}_{\text{glucose}}^{-1}$), and a decrease in growth rate (from 0.3 to 0.18h⁻¹) and glucose uptake rate (from 4.2 to 2.2 mmol g⁻¹h⁻¹) for BapA production strains compared to the reference strain (Heyland et al., 2011b). The PPP flux correlated with biomass yield, whereas correlations between PPP and protein production were low. This indicates a higher demand for NADPH due to an increase in biomass synthesis rather than protein production. In this context, it is also important to mention that flux through PPP also depends on substrates used for cultivation. In general, methylotrophic yeasts like K. phaffii have retained high PPP activity during evolution (Riley et al., 2016). The assimilation of methanol requires a higher carbon flux through PPP (Jordà et al., 2012) since the ketose sugar PPP intermediate, xylulose 5-phosphate, is needed to drive this reaction. Thus, increased PPP flux does not only reflect higher NADPH demand in cases where methanol is used during cultivation (Nie et al., 2014). Meanwhile, S. cerevisiae is unable to grow on methanol, using simple sugars such as glucose as carbon sources.

In summary, these two organisms rely on different pathways for NADPH production in burdened production strains: the acetate formation pathway is favored in S. cerevisiae, while the oxidative phase of PPP in K. phaffii.

General Challenges With Redox Imbalance in Burdened Production Strains

Another important redox cofactor in cells is NAD+/NADH. During oxidative phosphorylation in mitochondria, NADH is required for ATP energy formation through the respiratory chain reaction. ATP is needed for almost all growth-related activities, cellular maintenance, and heterologous protein production. Hence, the burden from heterologous protein production increases the demand for carbon allocated to the TCA cycle and oxidative phosphorylation for ATP production (Heyland et al., 2010). For this reason, ATP is likely to decrease in high producing strains as a response to burdensome protein production, compensated by slower growth and carbon uptake (Tyo et al., 2012; Nie et al., 2014). This further highlights the central role of carbon metabolism on burden.

High NADH demand in production strains requires greater redox balance. Under aerobic conditions, S. cerevisiae and K. phaffii can use glycerol as an electron acceptor to oxidize NADH. This was used as an explanation for glycerol production in K. phaffii to help maintain redox balance and avoid cytosolic NAD+ depletion (Nie et al., 2014). In the same study, the authors also found that higher producing strains reduced by-product formation, such as glycerol production, and concluded this is a strategy for the cell to compensate for a higher demand of redox-cofactors and energy (Nie et al., 2014). In line with this, an increase in glycerol production can also lead to lower biomass yield as carbon is removed from central carbon metabolism, contributing to burden (Krogh et al., 2008). Of note for K. phaffii, problems associated with NADH limitations will likely be worse in glucose cultivation as opposed to growth on methanol in co-fed cultures because of the energetic regulation of enzymes involved in methanol oxidation and assimilation pathway (Jordà et al., 2012). While redox requirements for translation are higher in production strains, post-translational processes also require redox factors and other resources. The impact on burden from post-translational processes is discussed next.

Burdensome Protein Folding, Export, and Secretion Triggers Cellular Stress Responses

Secretion of post-translationally modified proteins by yeast cell factories is advantageous as it facilitates simple downstream purification and processing steps (Nielsen, 2013). The downside is that it increases demand for heterologous protein translocation into the endoplasmic reticulum (ER), folding in the ER, and export via the Golgi apparatus. Together this can introduce bottlenecks that place additive burden on yeast cell factories and trigger cellular stress responses (Jordà et al., 2012; Tyo et al., 2012; Huang et al., 2017; de Ruijter et al., 2018). Indeed, the impact of post-translational processes on the burden state is thought to be a major limitation in yeast, especially considering relatively lower expression levels for these organisms compared to bacteria (Mattanovich et al., 2004).

Proteins are modified, folded, and transported in the ER. Yeast cells must compensate for the burden conferred by the energy-consuming secretory pathway. This is achieved by allocating more resources to cellular maintenance for energy production and away from growth (Jordà et al., 2012; Chen et al., 2020). Even when directing expressed proteins to mitochondria, in contrast to cellular secretion, resources for transportation are a limiting factor causing plasmid expression repression and growth defects (Kintaka et al., 2016; Eguchi et al., 2018). Similarly, anchoring heterologous β-glucosidase to the yeast cell membrane induces burden (Ding et al., 2018). Thus, in order to fully understand burden in yeast, we must better elucidate the burdensome steps during protein transportation and secretion.

Translocation Into the ER

Once a heterologous protein has been effectively translated in yeast, the nascent peptide is translocated into the ER. The timing of nascent peptide translocation, e.g., co- or post-translational, can affect burden on yeast cells (Tang et al., 2015; Barrero et al., 2021) and is directed by the pre-sequence of the leader (Hou et al., 2012). For example, the K. phaffii pre-alpha-mating factor (MF), which is commonly fused to heterologous proteins to drive their post-translational translocation into the ER, augmented burden during Rhizopus oryzae lipase (ROL) expression. This was observed as a reduction in growth rate, biomass yield, and final protein titer (Barrero et al., 2021). Notably, burden was less pronounced when using an alternative translocation signal called pre-Ost which stimulates co-translational translocation of the nascent peptide into ER. In fact, compared to ROL fused to pre-alpha-MF signal, the final ROL titer was improved from 192.4 to 291.9 U/ml (Barrero et al., 2021). Similarly, co-translational translocation enhanced secretion from S. cerevisiae production stains expressing β-glucosidase, endoglucanase, and α-amylase (Tang et al., 2015). In some cases, native leader sequences will be advantageous to use instead of yeast-specific leaders as for example reported for human serum albumin (Sleep et al., 1990). However, predicting which leader sequence most efficiently facilitates translocation into ER and further secretion is difficult (Hou et al., 2012). The use of non-efficient leader sequences for directing a heterologous protein through the secretory pathway can contribute to burden. Secretion via co-translational translocation of heterologous proteins into the ER, in contrast to post-translational translocation, was reported to alleviate some of the burden introduced at this step and can be a strategy to mitigate burden. Such strategy is highlighted as engineering target 2 in Figure 1.

Protein Folding in the ER

Major resource consuming steps in the secretory pathway include protein folding and disulfide bond formation, a common stabilizing feature in many secreted proteins. Protein chaperones consume significant amounts of ATP to facilitate protein folding and thereby prevent their aggregation and/or degradation (Umebayashi et al., 1997; Walter and Buchner, 2002; Hartl et al., 2011). Despite added energy demands, co-expressing certain chaperones in yeast production strains has been shown to enhance heterologous protein by helping to overcome burden at this crucial folding step (Zhang et al., 2006; Gu et al., 2015; Yu et al., 2017). Addition of affinity tags, usually added for detection and purification, can also help facilitate folding. These fusion partners might however induce burden at different levels, such as widely used His-6 tag is proposed to induce low burden and glutathione S-transferase tag to induce high burden (Waugh, 2005). Therefore, it is advised to consider this when adding affinity tags.

The formation of disulfide bonds in the folding process is redox-driven, which can lead to severe oxidative stress in production cells if imbalanced (Margittai and Sitia, 2011). De Ruijter et al. recently reported that production of heterologous antibody fragments in S. cerevisiae, which require disulfide bonds, clearly burdened ER redox balance and increased flux of glutathione metabolism (de Ruijter et al., 2018). This is because protein folding becomes slower than disulfide bond formation as the protein folding machinery becomes overloaded and this will consume glutathione for disulfide bond breaking, producing reactive oxygen species (ROS) that induce oxidative stress (Tyo et al., 2012). This has also been proposed as an explanation for proteome changes over time in S. cerevisiae cells burdened by insulin production, which also requires disulfide bond formation (Wright et al., 2020). In line with this, overexpressing glutathione peroxidase 1, the enzyme responsible for detoxifying ROS from reduced glutathione, was observed to improve disulfide bond-containing AppA phytase production by 1.3-fold (from ~350 mU/ml) in K. phaffii (Navone et al., 2021). Alternatively, in another study burden was mitigated by over-expressing HAP1 gene encoding a transcription factor activating oxidative stress response genes and helped to increase growth rate and respiration. This alpha-amylase producing S. cerevisiae strain performed better (45 U${\text{g}}_{\text{DCW}}^{-1}$h⁻¹) than the reference strain (34 U${\text{g}}_{\text{DCW}}^{-1}$h⁻¹) in a chemostat with dilution rate near the maximum of the reference strain but still far from maximum for the optimized strain (Martínez et al., 2016). Such approaches correspond to engineering target 3 in Figure 1.

Engineering strategies focused on co-expressing protein chaperones, peroxidases or transcription factors to resolve secretory pathway bottlenecks might seem counter-intuitive, since expressing additional proteins might be expected to have an additive effect on burden by competing for some of the same translational and post-translational resources as the heterologous protein (Liu et al., 2014; Huang et al., 2018; Yu et al., 2020). This can explain why such strategies do not always work or sometimes only modestly enhance expression (Yu et al., 2017). Tightly regulating and fine-tuning the overall amount and timing of expression for such proteins is expected to have better results (Navone et al., 2021) than utilizing high copy number plasmids for over-expression (Huang et al., 2018).

The cofactor NADPH is also used as a reducing agent of glutathione disulfide to regenerate reduced glutathione (Ayer et al., 2014). The need for NADPH was suggested to cause similar levels of flux through the oxidative PPP branch in lipase-producing K. phaffii compared to the reference strain, despite lower biomass yield in the production strains. The high PPP flux ensures adequate levels of the glutathione electron donor for protein folding (Jordà et al., 2012). Moreover, for an S. cerevisiae strain producing a human insulin precursor, burden-induced lower PPP flux was considered to have a negative impact on folding in ER, because another enzyme involved in disulfide bond formation, termed protein disulfide isomerase (PDI1), also requires NADPH (Kazemi Seresht et al., 2013). Therefore, NADPH limitations become more problematic in production strains engineered to secrete proteins, as opposed to intracellularly expressed proteins, due to the additional demand for NADPH during later export processes after protein folding in ER (Liu et al., 2016). While tuning the oxidative phase of PPP in burdened K. phaffii enhances heterologous protein production (Nocon et al., 2016), it is currently unclear whether this PPP engineering strategy would improve production in S. cerevisiae.

Quality Control Systems in the ER

Burden associated with increased oxidative stress response and the accumulation of unfolded protein in ER can trigger the activation of the unfolded protein response (UPR). UPR is a quality control system that helps re-establish cellular proteome homeostasis and is therefore activated as a response to burden in the ER (Huang et al., 2017; Chen et al., 2020). While chaperones that act on misfolded proteins may buffer against burden (Farkas et al., 2018), UPR acts as a defense system against burden. However, in some instances UPR activation might contribute additional burden on cells (Lamour et al., 2019), due to the upregulation in many UPR genes (de Ruijter et al., 2018). In contrast to this view, a recent study in K. phaffii found that proteins involved in UPR decreased when switching from glycerol batch pre-cultivation to methanol fed-batch (Vanz et al., 2014). Glycerol cultivation is associated with high UPR induction, and although cultivation in methanol does not induce the same level of UPR activation, the authors suggest that heterologous protein production may benefit from initially high UPR activity. In line with this, the basal UPR level in K. phaffii appears to be higher compared to other yeasts (Vanz et al., 2014). Therefore, enhanced UPR activity or the co-expression of UPR factors might mitigate burden upon the strong induction of heterologous protein synthesis following the switch to methanol as a carbon source.

The main regulator of UPR in yeast is the transcription factor Homologous to Atf/Creb1 (HAC1). Prolonged overexpressing of this transcription factor improved production in burdened K. phaffii strains, increasing heterologous Xylanase from ~140 to 200 U/ml and glucose oxidase from ~55 to 140 U/ml (Lin et al., 2013; Yu et al., 2020). However, it also reduced cell growth. The negative correlation between protein secretion capacity and cellular growth was ascribed to stress but meets our definition of burden. And despite efforts to reduce burden with this strategy, the co-expression of a heterologous protein and the UPR regulator was able to enhance expression but contributed its own burden (Yu et al., 2020). This might be explained by the fact that increased ROS levels were associated with HAC1 overexpression, suggesting that prolonged UPR activation triggered oxidative stress (Yu et al., 2020). Therefore, while this strategy successfully enhanced heterologous protein expression, a new source of burden is elicited.

In addition to UPR, the ER-associated-degradation (ERAD) pathway is another important quality control system in yeast. ERAD targets misfolded proteins for cytosolic degradation and ensures that ER homeostasis is maintained (Hwang and Qi, 2018). Co-expressing ERAD-related genes has been explored as a strategy to mitigate burden in yeast. For example, co-expressing the ERAD ubiquitin ligase called HMG-coA Reductase Degradation (HRD1) together with a heterologous glucose oxidase helped limit burden in K. phaffii (Gu et al., 2015). While co-expressing single genes such as HRD1 or HAC1 alleviates some of the cellular stress induced by protein folding burden, yeast can adapt to prolonged ER stress induced by folding burden through chromosomal duplication (Beaupere et al., 2018; Beaupere and Labunskyy, 2019). Such genomic instability might present additional problems for production strains, indicating that more tightly controlled regulatory circuits that selectively turn specific genes in the UPR and ERAD pathways on or off might be more desirable than constitutively over-expressing them.

To summarize, cellular stress is tightly coupled to burden. As heterologous proteins are directed to the secretory pathway, burden is induced by overload and subsequent disruptions in protein folding and redox balance in the ER which trigger the onset of different cellular stress responses. Engineering target 3 in Figure 1 highlights this important step. Evidence supports increased heterologous protein titers when co-expressing one or a few specific stress-response proteins, indicating that stress-response pathways provide attractive targets for engineering. However, more robust, and tightly regulated circuits can potentially be more efficient to mitigate burden accumulation and allow cells to maintain a healthier homeostasis.

Future Profiling and Engineering Strategies

Several process and genetic engineering strategies have been proposed to improve production from burdened host cells (Deparis et al., 2017). We summarize these in Figure 1 and Table 2. These strategies encompass: enhancing the availability of certain co-factors; improving protein folding and/or folding capacity in the ER; optimizing gene copy number and codon usage; engineering improved gene promoters; and finally, optimizing cultivation conditions to support increased resource demand. Moving forward, engineering metabolic networks represents a promising approach to improve heterologous protein titers by helping alleviate redox cofactor and energy limitations in burdened cells, such as leading flux toward necessary precursors (Mattanovich et al., 2014; Nie et al., 2014).

TABLE 2

Table 2. Strategies to improve production from burdened yeast cell factories.

Systems biology offers a promising interdisciplinary approach to discover new targets for metabolic engineering (Hou et al., 2012). Apart from analyzing systems level data, it can also be integrated into mechanistic models. Combining in silico modeling with -omics techniques and bioinformatic analyses allows for the identification of factors that are critical but might not necessarily stand out when the data is first observed. Thus, modeling represents a useful tool for quantifying the bioproduction limitations of living cells (Mattanovich et al., 2014; Volkova et al., 2020; Patra et al., 2021). Constraint-based stoichiometric models, in particular, can be used to calculate the metabolic differences between strains and allow for quantification of burden on metabolism (Gonzalez et al., 2003). For example, carbon flux distribution in central carbon metabolic pathways can be traced by applying C¹³-metabolic flux analysis and additional in silico data into stoichiometric models (Jordà et al., 2012; Nie et al., 2014). To broaden the applicability, burden should be addressed and incorporated into such models like large genome-scale models (Saitua et al., 2017).

Indeed, diverse systems-based efforts have been introduced to alleviate burden in production hosts. Notably, such efforts have inspired the construction of dynamic control systems that respond to changing metabolic and energy demands. These have proven highly effective in E. coli production strains and show strengths over the more classic genetic engineering approach of simply deleting or over-/under-expressing candidate genes. As heterologous protein production requires resources from shared pools of metabolites with native cellular processes, static overexpression of specific genes leads to a less robust systems since it does not help minimize the changing levels of burden on production host during cultivation (Boo et al., 2019). Instead, dynamic control can stabilize heterologous protein production over time. For example, Ceroni et al. engineered a burden-driven feedback system in E. coli using CRISPR technology which dynamically regulated heterologous protein production in response to burden (Ceroni et al., 2018). In this system, transcription of a single guide RNA (sgRNA) is controlled by a burden-responsive promoter identified by RNA-seq analysis. In response to burden, the sgRNA is transcribed and directs a catalytically-dead Cas9 protein (dCas9) to the heterologous gene promoter, suppressing its expression until burden is alleviated (Ceroni et al., 2018). By temporarily reducing heterologous protein production, burdened E. coli cells have time to recover and achieve greater production levels over the total course of the cultivation. This study highlights the potential for engineering similar burden-driven negative feedback systems for heterologous protein production in yeast.

Similar and other kinds of regulatory feedback circuits have already proven useful to mitigate or circumvent burden in S. cerevisiae strains engineered to produce valuable compounds such as bioethanol or vanillin-β-glucoside (D'Ambrosio et al., 2020; Qin et al., 2020). For bioethanol production in S. cerevisiae, a more robust strain was achieved by placing target genes in the ethanol pathway under control of stress-regulated promoters to ensure dynamic feedback regulation (Qin et al., 2020). Alternatively, a stable S. cerevisiae cell factory was developed to produce vanillin-β-glucoside by coupling pathway intermediate production to biosensor-controlled expression of an essential gene (D'Ambrosio et al., 2020). The absence of the biosensor-coupled burden pathway intermediate renders host cells inviable. In other words, the cells become addicted to production. Although these two examples highlight alternative engineering strategies to produce valuable chemical compounds in yeast, they also involve tightly regulated changes in protein expression (e.g., enzymes required to synthesize valuable metabolites) in response to burden imposed by production.

The above mentioned strategies represent novel and promising approaches which could be applied for burden-driven regulation of heterologous proteins expression. It might also be useful to place secondary “helper genes” under the control of such dynamic regulation, including genes encoding some combination of factors involved in restoring redox balance, controlling UPR and/or ERAD pathways, or ensuring adequate protein folding. Engineering dynamic systems that tailor heterologous protein expression to ever-changing host cell needs over the course of the cultivation represents a novel approach that aims at multiple engineering targets highlighted in Figure 1 simultaneously. Such an approach is anticipated to generate burden-resistant yeast cell factories to boost biopharmaceutical and industrial enzyme production ensuring more stable production over time.

Discussion

Burden derived from resource competition is often encountered in bioproduction processes triggering a phenotype with compromised production during long-term cultivations. Therefore, engineering production strains that are either more tolerant to burden or able to adjust heterologous protein expression to burden is paramount, and efforts have already been made to engineer such yeast strains as listed in Table 2. To do so requires a better understanding of the different molecular sources of this phenomenon. A reduction in heterologous protein titer, as a consequence of burden, might happen for different reasons. For example, burden itself can impose selective pressure on cells evolving toward lower levels of production through genetic or non-genetic variation (Rugbjerg and Olsson, 2020). This was recently observed in an isogenic population of insulin-producing S. cerevisiae cell factories, where burdensome production and cultivation caused non-genetic heterogeneity in the population to favor the growth of low-producing cells in the population (Wright et al., 2021).

Commonly, the impact of burden is measured as a reduction in physiological parameters like specific growth rate, biomass yield, and respiratory capacity (Kazemi Seresht et al., 2013; Liu et al., 2014; de Ruijter et al., 2018). It is worth mentioning that burden might manifest in different ways depending on scale and mode of fermentation. Experimental setup and fermentation equipment will undoubtedly have an impact on these physiological parameters. While reduced protein titers can be a measure of burden during fermentation more broadly, measures like growth rate and biomass are more relevant in batch-type fermentations (Heyland et al., 2011b; Van Rensburg et al., 2012), often representing a screening scenario for strain development. Transferability of such results to industrially relevant cultivation modes is however difficult. It is questionable whether optimized strains from screening will perform well in fed-batch or continuous fermentation environments (Looser et al., 2014). Proxies for cell maintenance, such as amount of energy that the host cell uses to readjust to protein production, is also relevant to assess burden in controlled and batch fermentations. Further, it will be more relevant to use measures like respiratory quotient, that cannot be accounted for during biomass yield, as a measure for burden in, for example, chemostat cultivation mode (Jordà et al., 2012). To do so, underlying biological processes, such as metabolic fluxes and transcription factor activity, are relevant to quantify.

Sources of host cell burden can be more specifically identified from high-throughput, omics-based approaches since they provide a systems-level overview of the burden response. With tools that profile the metabolome, transcriptome, proteome, and/or fluxome of any specific production strain over time, these complex networks can be dissected to understand the burdened phenotype at a molecular level, identify causal factors responsible for pathway bottleneck, and allow rational engineering strategies to alleviate burden and enhance heterologous protein production (Heyland et al., 2011b; Jordà et al., 2012; Liu et al., 2016; Huang et al., 2017; Wright et al., 2020). While the resource drain originating from all processes during the biosynthetic production of a heterologous protein (transcription, translation, folding, and secretion) can have an impact on overall burden, the main culprit, as we have seen, is often tied back to constraints in central carbon metabolism. This central metabolic network is highly affected by the increased need for energy, redox-cofactors, and free metabolic precursors required for the costly production of an additional non-essential protein (Mattanovich et al., 2014). Therefore, cellular metabolism appears key to understanding the underlying mechanisms for burden state induction in yeast, and therefore suitably defined as metabolic burden.

In line with this, the gene regulatory network of underlying metabolic pathways might play an important role in inducing burden and could be an interesting direction for future studies. As reviewed here, cellular stress responses are intertwined with the burden state. Likely, these stress responses and other pathways have additional functions such as facilitating crosstalk with the metabolic pathways affected by burden (Kalender and Çalik, 2020). For example, UPR is reported to behave as a regulatory network activating cytosolic pathways. In a burdened insulin-producing strain of S. cerevisiae, activation of amino acid biosynthesis was suggested to be elicited following the up-regulation of transcription factor Gcn4, which is induced by UPR (Kazemi Seresht et al., 2013). It is expected that Gcn4 has the same properties in K. phaffii (Gu et al., 2015). Therefore, elucidating such system-wide interactions and regulatory networks will be critical to fully understand the cellular response to burden.

There are still gaps in our understanding of the concept of burden in yeast cell factories. By emphasizing this topic, we hope to draw the attention of future research to elucidate the molecular basis and underlying mechanisms of burden in yeast to mitigate this common production limitation. Despite these challenges, recent advances in bioengineering are already paving the way for improved heterologous protein titers from yeast cell factories by developing strategies that alleviate burden, such as dynamic control system engineering. We expect such approaches to be adopted by the biotech industry and significantly enhance heterologous protein production in the years to come.

Author Contributions

CW and MJ conceived the topic reviewed in this manuscript. CW, MJ, RA, and LK contributed content and ideas for the manuscript structure. LK wrote the manuscript. CW, MJ, and RA contributed with corrections and revised drafts at all stages. All authors approved the submitted version.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We would like to thank Novo Nordisk Foundation for funding this project in the frame of the Fermentation-Based Biomanufacturing education initiative (grant number NNF17SA0031362) and Novo Nordisk Foundation Center for Biosustainability grant number NNF20CC0035580 to MJ. We also would like to thank members of the Krogh and Workman labs for helpful discussions. Additional thanks to the reviewers for valuable feedback that improved this manuscript. The figure was created with www.BioRender.com.

References

Ayer, A., Gourlay, C. W., and Dawes, I. W. (2014). Cellular redox homeostasis, reactive oxygen species and replicative ageing in Saccharomyces cerevisiae. FEMS Yeast Res. 14, 60–72. doi: 10.1111/1567-1364.12114

PubMed Abstract | CrossRef Full Text | Google Scholar

Baghban, R., Farajnia, S., Rajabibazl, M., Ghasemi, Y., Mafi, A. A., Hoseinpoor, R., et al. (2019). Yeast expression systems: overview and recent advances. Mol. Biotechnol. 61, 365–384. doi: 10.1007/s12033-019-00164-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Barrero, J. J., Pagazartaundua, A., Glick, B. S., Valero, F., and Ferrer, P. (2021). Bioreactor-scale cell performance and protein production can be substantially increased by using a secretion signal that drives co-translational translocation in Pichia pastoris. New Biotechnol. 60, 85–95. doi: 10.1016/j.nbt.2020.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Beaupere, C., Dinatto, L., Wasko, B. M., Chen, R. B., VanValkenburg, L., Kiflezghi, M. G., et al. (2018). Genetic screen identifies adaptive aneuploidy as a key mediator of ER stress resistance in yeast. Proc. Natl. Acad. Sci. U. S. A. 115, 9586–9591. doi: 10.1073/pnas.1804264115

PubMed Abstract | CrossRef Full Text | Google Scholar

Beaupere, C., and Labunskyy, V. M. (2019). (Un)folding mechanisms of adaptation to ER stress: lessons from aneuploidy. Curr. Genet. 65, 467–471. doi: 10.1007/s00294-018-0914-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Blank, L. M., and Sauer, U. (2004). TCA cycle activity in Saccharomyces cerevisiae is a function of the environmentally determined specific growth and glucose uptake rates. Microbiology 150, 1085–1093. doi: 10.1099/mic.0.26845-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Boo, A., Ellis, T., and Stan, G. B. (2019). Host-aware synthetic biology. Curr. Opin. Syst. Biol. 14, 66–72. doi: 10.1016/j.coisb.2019.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Cámara, E., Albiol, J., and Ferrer, P. (2016). Droplet digital PCR-aided screening and characterization of Pichia pastoris multiple gene copy strains. Biotechnol. Bioeng. 113, 1542–1551. doi: 10.1002/bit.25916

PubMed Abstract | CrossRef Full Text | Google Scholar

Cámara, E., Landes, N., Albiol, J., Gasser, B., Mattanovich, D., and Ferrer, P. (2017). Increased dosage of AOX1 promoter-regulated expression cassettes leads to transcription attenuation of the methanol metabolism in Pichia pastoris. Sci. Rep. 7, 1–16. doi: 10.1038/srep44302

PubMed Abstract | CrossRef Full Text | Google Scholar

Celik, E., Calik, P., and Oliver, S. G. (2009). Fed-batch methanol feeding strategy for recombinant protein production by Pichia pastoris in the presence of co-substrate sorbitol. Yeast 26, 473–484. doi: 10.1002/yea.1679

PubMed Abstract | CrossRef Full Text | Google Scholar

Cereghino, J. L., and Cregg, J. M. (2000). Heterologous protein expression in the methylotrophic yeast Pichia pastoris. FEMS Microbiol. Rev. 24, 45–66. doi: 10.1111/j.1574-6976.2000.tb00532.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ceroni, F., Boo, A., Furini, S., Gorochowski, T. E., Borkowski, O., Ladak, Y. N., et al. (2018). Burden-driven feedback control of gene expression. Nat. Methods 15, 387–393. doi: 10.1038/nmeth.4635

PubMed Abstract | CrossRef Full Text | Google Scholar

Chan, P. P., and Lowe, T. M. (2016). GtRNAdb 2.0: an expanded database of transfer RNA genes identified in complete and draft genomes. Nucleic Acids Res. 44, D184–D189. doi: 10.1093/nar/gkv1309

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Mohsin, A., Chu, J., Zhuang, Y., Liu, Y., and Guo, M. (2017). Enhanced protein production by sorbitol co-feeding with methanol in recombinant Pichia pastoris strains. Biotechnol. Bioprocess Eng. 22, 767–773. doi: 10.1007/s12257-017-0011-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Ji, B., Hao, X., Li, X., Eisele, F., Nyström, T., et al. (2020). FMN reduces Amyloid-β toxicity in yeast by regulating redox status and cellular metabolism. Nat. Commun. 11, 1–16. doi: 10.1038/s41467-020-14525-4

PubMed Abstract | CrossRef Full Text | Google Scholar

D'Ambrosio, V., Dore, E., Di Blasi, R., van den Broek, M., Sudarsan, S., Horst, J., ter Ambri, F., et al. (2020). Regulatory control circuits for stabilizing long-term anabolic product formation in yeast. Metab. Eng. 61, 369–380. doi: 10.1016/j.ymben.2020.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

de Ruijter, J. C., Koskela, E. V., Nandania, J., Frey, A. D., and Velagapudi, V. (2018). Understanding the metabolic burden of recombinant antibody production in Saccharomyces cerevisiae using a quantitative metabolomics approach. Yeast 35, 331–341. doi: 10.1002/yea.3298

PubMed Abstract | CrossRef Full Text | Google Scholar

De Schutter, K., Lin, Y. C., Tiels, P., Van Hecke, A., Glinka, S., Weber-Lehmann, J., et al. (2009). Genome sequence of the recombinant protein production host Pichia pastoris. Nat. Biotechnol. 27, 561–566. doi: 10.1038/nbt.1544

PubMed Abstract | CrossRef Full Text | Google Scholar

Deparis, Q., Claes, A., Foulquié-Moreno, M. R., and Thevelein, J. M. (2017). Engineering tolerance to industrially relevant stress factors in yeast cell factories. FEMS Yeast Res. 17, 1–17. doi: 10.1093/femsyr/fox036

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, J., Liang, G., Zhang, K., Hong, J., Zou, S., Lu, H., et al. (2018). Extra metabolic burden by displaying over secreting: Growth, fermentation and enzymatic activity in cellobiose of recombinant yeast expressing B-glucosidase. Bioresour. Technol. 254, 107–114. doi: 10.1016/j.biortech.2017.12.030

PubMed Abstract | CrossRef Full Text | Google Scholar

dos Reis, M., Savva, R., and Wernisch, L. (2004). Solving the riddle of codon usage preferences: a test for translational selection. Nucleic Acids Res. 32, 5036–5044. doi: 10.1093/nar/gkh834

PubMed Abstract | CrossRef Full Text | Google Scholar

Drummond, D. A., and Wilke, C. O. (2008). Mistranslation-induced protein misfolding as a dominant constraint on coding-sequence evolution. Cell 134, 341–352. doi: 10.1016/j.cell.2008.05.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Duman-Özdamar, Z. E., and Binay, B. (2021). Production of industrial enzymes via Pichia pastoris as a cell factory in bioreactor: current status and future aspects. Protein J. 40, 367–376. doi: 10.1007/s10930-021-09968-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Dumont, J., Euwart, D., Mei, B., Estes, S., and Kshirsagar, R. (2016). Human cell lines for biopharmaceutical manufacturing: history, status, and future perspectives. Crit. Rev. Biotechnol. 36, 1110–1122. doi: 10.3109/07388551.2015.1084266

PubMed Abstract | CrossRef Full Text | Google Scholar

Eguchi, Y., Makanae, K., Hasunuma, T., Ishibashi, Y., Kito, K., and Moriya, H. (2018). Estimating the protein burden limit of yeast cells by measuring the expression limits of glycolytic proteins. ELife 7, 1–23. doi: 10.7554/eLife.34595

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, C., Wang, Q., Selvaraj, J. N., Zhou, Y., Ma, L., Zhang, G., et al. (2017). High copy and stable expression of the xylanase XynHB in Saccharomyces cerevisiae by rDNA-mediated integration. Sci. Rep. 7, 1–9. doi: 10.1038/s41598-017-08647-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Farkas, Z., Kalapis, D., Bódi, Z., Szamecz, B., Daraba, A., Almási, K., et al. (2018). Hsp70-associated chaperones have a critical role in buffering protein production costs. ELife 7, 1–23. doi: 10.7554/eLife.29845

PubMed Abstract | CrossRef Full Text | Google Scholar

Garrigós-Martínez, J., Vuoristo, K., Nieto-Taype, M. A., Tähtiharju, J., Uusitalo, J., Tukiainen, P., et al. (2021). Bioprocess performance analysis of novel methanol-independent promoters for recombinant protein production with Pichia pastoris. Microbial Cell Fact. 20, 1–12. doi: 10.1186/s12934-021-01564-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Gasch, A. P. (2003). “The environmental stress response: a common yeast response to diverse environmental stresses,” in Yeast Stress Responses, eds S. Hohmann and W. H. Mager (Berlin, Heidelberg: Springer), 11–70. doi: 10.1007/3-540-45611-2_2

CrossRef Full Text | Google Scholar

Ghaderi, D., Zhang, M., Hurtado-Ziola, N., and Varki, A. (2012). Production platforms for biotherapeutic glycoproteins. Occurrence, impact, and challenges of non-human sialylation. Biotechnol. Genet. Eng. Rev. 28, 147–176. doi: 10.5661/bger-28-147

PubMed Abstract | CrossRef Full Text | Google Scholar

Glick, B. R. (1995). Metabolic load and heterologous gene expression. Biotechnol. Adv. 13, 247–261. doi: 10.1016/0734-9750(95)00004-A

PubMed Abstract | CrossRef Full Text | Google Scholar

Goffeau, A. A., Barrell, B. G., Bussey, H., Davis, R. W., Dujon, B., Feldmann, H., et al. (1996). Life with 6000 genes. Science 274, 4–10. doi: 10.1126/science.274.5287.546

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonzalez, R., Andrews, B. A., Molitor, J., and Asenjo, J. A. (2003). Metabolic analysis of the synthesis of high levels of intracellular human SOD in Saccharomyces cerevisiae rhSOD 2060 411 SGA122. Biotechnol. Bioeng. 82, 152–169. doi: 10.1002/bit.10556

PubMed Abstract | CrossRef Full Text | Google Scholar

Görgens, J. F., Van Zyl, W. H., Knoetze, J. H., and Hahn-Hägerdal, B. (2001). The metabolic burden of the PGK1 and ADH2 promoter systems for heterologous xylanase production by saccharomyces cerevisiae in defined medium. Biotechnol. Bioeng. 73, 238–245. doi: 10.1002/bit.1056

PubMed Abstract | CrossRef Full Text | Google Scholar

Görgens, J. F., Van Zyl, W. H., Knoetze, J. H., and Hahn-Hägerdal, B. (2005). Amino acid supplementation improves heterologous protein production by Saccharomyces cerevisiae in defined medium. Appl. Microbiol. Biotechnol. 67, 684–691. doi: 10.1007/s00253-004-1803-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Grabowska, D., and Chelstowska, A. (2003). The ALD6 gene product is indispensable for providing NADPH in yeast cells lacking glucose-6-phosphate dehydrogenase activity. J. Biol. Chem. 278, 13984–13988. doi: 10.1074/jbc.M210076200

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, L., Zhang, J., Du, G., and Chen, J. (2015). Multivariate modular engineering of the protein secretory pathway for production of heterologous glucose oxidase in Pichia pastoris. Enzyme Microbial Technol. 68, 33–42. doi: 10.1016/j.enzmictec.2014.10.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Harrison, E., Koufopanou, V., Burt, A., and MacLean, R. C. (2012). The cost of copy number in a selfish genetic element: the 2-μm plasmid of Saccharomyces cerevisiae. J. Evol. Biol. 25, 2348–2356. doi: 10.1111/j.1420-9101.2012.02610.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hartl, F. U., Bracher, A., and Hayer-Hartl, M. (2011). Molecular chaperones in protein folding and proteostasis. Nature 475(7356), 324–332. doi: 10.1038/nature10317

PubMed Abstract | CrossRef Full Text | Google Scholar

He, H., Wu, S., Mei, M., Ning, J., Li, C., Ma, L., et al. (2020). A combinational strategy for effective heterologous production of functional human lysozyme in Pichia pastoris. Front. Bioeng. Biotechnol. 8, 118. doi: 10.3389/fbioe.2020.00118

PubMed Abstract | CrossRef Full Text | Google Scholar

Heyland, J., Blank, L. M., and Schmid, A. (2011a). Quantification of metabolic limitations during recombinant protein production in Escherichia coli. J. Biotechnol. 155, 178–184. doi: 10.1016/j.jbiotec.2011.06.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Heyland, J., Fu, J., Blank, L. M., and Schmid, A. (2010). Quantitative physiology of Pichia pastoris during glucose-limited high-cell density fed-batch cultivation for recombinant protein production. Biotechnol. Bioeng. 107, 357–368. doi: 10.1002/bit.22836

PubMed Abstract | CrossRef Full Text | Google Scholar

Heyland, J., Fu, J., Blank, L. M., and Schmid, A. (2011b). Carbon metabolism limits recombinant protein production in Pichia pastoris. Biotechnol. Bioeng. 108, 1942–1953. doi: 10.1002/bit.23114

PubMed Abstract | CrossRef Full Text | Google Scholar

Hou, J., Tyo, K. E. J., Liu, Z., Petranovic, D., and Nielsen, J. (2012). Metabolic engineering of recombinant protein secretion by Saccharomyces cerevisiae. FEMS Yeast Res. 12, 491–510. doi: 10.1111/j.1567-1364.2012.00810.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, M., Bao, J., Hallström, B. M., Petranovic, D., and Nielsen, J. (2017). Efficient protein production by yeast requires global tuning of metabolism. Nat. Commun. 8, 1131. doi: 10.1038/s41467-017-00999-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, M., Wang, G., Qin, J., Petranovic, D., and Nielsen, J. (2018). Engineering the protein secretory pathway of Saccharomyces cerevisiae enables improved protein production. Proc. Natl. Acad. Sci. U. S. A., 115, E11025–E11032. doi: 10.1073/pnas.1809921115

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, J., and Qi, L. (2018). Quality control in the endoplasmic reticulum: crosstalk between ERAD and UPR pathways. Trends Biochem. Sci. 43, 593–605. doi: 10.1016/j.tibs.2018.06.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Jensen, T. H., Jacquier, A., and Libri, D. (2013). Dealing with pervasive transcription. Mol. Cell 52, 473–484. doi: 10.1016/j.molcel.2013.10.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Jordà, J., Jouhten, P., Cámara, E., Maaheimo, H., Albiol, J., and Ferrer, P. (2012). Metabolic flux profiling of recombinant protein secreting Pichia pastoris growing on glucose:methanol mixtures. Microbial Cell Fact. 11, 1–14. doi: 10.1186/1475-2859-11-57

PubMed Abstract | CrossRef Full Text | Google Scholar

Jung, S. T., and Kim, D.-I. (2018). “Mass production of full-length IgG monoclonal antibodies from mammalian, yeast, and bacterial hosts.” in Emerging Areas in Bioengineering, ed H. N. Chang (Weinheim: Wiley-VCH), 679–695. doi: 10.1002/9783527803293.ch39

PubMed Abstract | CrossRef Full Text | Google Scholar

Kafri, M., Metzl-Raz, E., Jona, G., and Barkai, N. (2016). The cost of protein production. Cell Rep. 14, 22–31. doi: 10.1016/j.celrep.2015.12.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalender, Ö., and Çalik, P. (2020). Transcriptional regulatory proteins in central carbon metabolism of Pichia pastoris and Saccharomyces cerevisiae. Appl. Microbiol. Biotechnol. 104, 7273–7311. doi: 10.1007/s00253-020-10680-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Karbalaei, M., Rezaee, S. A., and Farsiani, H. (2020). Pichia pastoris: a highly successful expression system for optimal synthesis of heterologous proteins. J. Cell. Physiol. 235, 5867–5881. doi: 10.1002/jcp.29583

PubMed Abstract | CrossRef Full Text | Google Scholar

Karim, A. S., Curran, K. A., and Alper, H. S. (2013). Characterization of plasmid burden and copy number in Saccharomyces cerevisiae for optimization of metabolic engineering applications. FEMS Yeast Res. 13, 107–116. doi: 10.1111/1567-1364.12016

PubMed Abstract | CrossRef Full Text | Google Scholar

Kazemi Seresht, A., Cruz, A. L., De Hulster, E., Hebly, M., Palmqvist, E. A., Van Gulik, W., et al. (2013). Long-term adaptation of Saccharomyces cerevisiae to the burden of recombinant insulin production. Biotechnol. Bioeng. 110, 2749–2763. doi: 10.1002/bit.24927

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, M. D., Rhee, S. K., and Seo, J. H. (2001). Enhanced production of anticoagulant hirudin in recombinant Saccharomyces cerevisiae by chromosomal δ-integration. J. Biotechnol. 85, 41–48. doi: 10.1016/S0168-1656(00)00376-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Kintaka, R., Makanae, K., and Moriya, H. (2016). Cellular growth defects triggered by an overload of protein localization processes. Sci. Rep. 6, 1–11. doi: 10.1038/srep31774

PubMed Abstract | CrossRef Full Text | Google Scholar

Klein, T., Niklas, J., and Heinzle, E. (2015). Engineering the supply chain for protein production/secretion in yeasts and mammalian cells. J. Indust. Microbiol. Biotechnol. 42, 453–464. doi: 10.1007/s10295-014-1569-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Krogh, A. M., Beck, V., Christensen, L. H., Henriksen, C. M., Møller, K., and Olsson, L. (2008). Adaptation of Saccharomyces cerevisiae expressing a heterologous protein. J. Biotechnol. 137, 28–33. doi: 10.1016/j.jbiotec.2008.07.1787

PubMed Abstract | CrossRef Full Text | Google Scholar

Küberl, A., Schneider, J., Thallinger, G. G., Anderl, I., Wibberg, D., Hajek, T., et al. (2011). High-quality genome sequence of Pichia pastoris CBS7435. J. Biotechnol. 154, 312–320. doi: 10.1016/j.jbiotec.2011.04.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Kulagina, N., Besseau, S., Godon, C., Goldman, G. H., Papon, N., and Courdavault, V. (2021). Yeasts as biopharmaceutical production platforms. Front. Fungal Biol. 2, 733492. doi: 10.3389/ffunb.2021.733492

CrossRef Full Text | Google Scholar

Kunert, R., and Reinhart, D. (2016). Advances in recombinant antibody manufacturing. Appl. Microbiol. Biotechnol. 100, 3451–3461. doi: 10.1007/s00253-016-7388-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwolek-Mirek, M., Maslanka, R., and Molon, M. (2019). Disorders in NADPH generation via pentose phosphate pathway influence the reproductive potential of the Saccharomyces cerevisiae yeast due to changes in redox status. J. Cell. Biochem. 120, 8521–8533. doi: 10.1002/jcb.28140

PubMed Abstract | CrossRef Full Text | Google Scholar

Lalor, F., Fitzpatrick, J., Sage, C., and Byrne, E. (2019). Sustainability in the biopharmaceutical industry: seeking a holistic perspective. Biotechnol. Adv. 37, 698–707. doi: 10.1016/j.biotechadv.2019.03.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Lamour, J., Wan, C., Zhang, M., Zhao, X., and Den Haan, R. (2019). Overexpression of endogenous stress-tolerance related genes in Saccharomyces cerevisiae improved strain robustness and production of heterologous cellobiohydrolase. FEMS Yeast Res. 19, 1–13. doi: 10.1093/femsyr/foz035

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, X. Q., Liang, S. L., Han, S. Y., Zheng, S. P., Ye, Y. R., and Lin, Y. (2013). Quantitative iTRAQ LC-MS/MS proteomics reveals the cellular response to heterologous protein overexpression and the regulation of HAC1 in Pichia pastoris. J. Proteomics 91, 58–72. doi: 10.1016/j.jprot.2013.06.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C. P., Tsai, T. I., Cheng, T., Shivatare, V. S., Wu, C. Y., Wu, C. Y., et al. (2018). Glycoengineering of antibody (Herceptin) through yeast expression and in vitro enzymatic glycosylation. Proc. Natl. Acad. Sci. U. S. A. 115, 720–725. doi: 10.1073/pnas.1718172115

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., Martínez, J. L., Liu, Z., Petranovic, D., and Nielsen, J. (2014). Balanced globin protein expression and heme biosynthesis improve production of human hemoglobin in Saccharomyces cerevisiae. Metab. Eng. 21, 9–16. doi: 10.1016/j.ymben.2013.10.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, P., Huang, M., Guo, M., Qian, J., Lin, W., Chu, J., et al. (2016). Combined 13C-assisted metabolomics and metabolic flux analysis reveals the impacts of glutamate on the central metabolism of high β-galactosidase-producing Pichia pastoris. Bioresour. Bioprocess. 3, 47. doi: 10.1186/s40643-016-0124-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Looser, V., Bruhlmann, B., Bumbak, F., Stenger, C., Costa, M., Camattari, A., et al. (2014). Cultivation strategies to enhance productivity of Pichia pastoris: a review. Biotechnol. Adv. 33, 1177–1193. doi: 10.1016/j.biotechadv.2015.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Love, K. R., Shah, K. A., Whittaker, C. A., Wu, J., Bartlett, M. C., Ma, D., et al. (2016). Comparative genomics and transcriptomics of Pichia pastoris. BMC Genomics 17, 550. doi: 10.1186/s12864-016-2876-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Majeke, B. M., García-Aparicio, M., Biko, O. D., Viljoen-Bloom, M., van Zyl, W. H., and Görgens, J. F. (2020). Synergistic codon optimization and bioreactor cultivation toward enhanced secretion of fungal lignin peroxidase in Pichia pastoris: enzymatic valorization of technical (industrial) lignins. Enzyme Microbial Technol. 139, 109593. doi: 10.1016/j.enzmictec.2020.109593

PubMed Abstract | CrossRef Full Text | Google Scholar

Margittai, É. V., and Sitia, R. (2011). Oxidative protein folding in the secretory pathway and redox signaling across compartments and cells. Traffic 12, 1–8. doi: 10.1111/j.1600-0854.2010.01108.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Martínez, J. L., Meza, E., Petranovic, D., and Nielsen, J. (2016). The impact of respiration and oxidative stress response on recombinant α-amylase production by Saccharomyces cerevisiae. Metab. Eng. Commun. 3, 205–210. doi: 10.1016/j.meteno.2016.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Mattanovich, D., Gasser, B., Hohenblum, H., and Sauer, M. (2004). Stress in recombinant protein producing yeasts. J. Biotechnol. 113, 121–135. doi: 10.1016/j.jbiotec.2004.04.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Mattanovich, D., Sauer, M., and Gasser, B. (2014). Yeast biotechnology: teaching the old dog new tricks. Microbial Cell Fact. 13, 1–5. doi: 10.1186/1475-2859-13-34

PubMed Abstract | CrossRef Full Text | Google Scholar

Matthews, C. B., Kuo, A., Love, K. R., and Love, J. C. (2018). Development of a general defined medium for Pichia pastoris. Biotechnol. Bioeng. 115, 103–113. doi: 10.1002/bit.26440

PubMed Abstract | CrossRef Full Text | Google Scholar

Mellitzer, A., Ruth, C., Gustafsson, C., Welch, M., Birner-Grünberger, R., Weis, R., et al. (2014). Synergistic modular promoter and gene optimization to push cellulase secretion by Pichia pastoris beyond existing benchmarks. J. Biotechnol. 191, 187–195. doi: 10.1016/j.jbiotec.2014.08.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Mellitzer, A., Weis, R., Glieder, A., and Flicker, K. (2012). Expression of lignocellulolytic enzymes in Pichia pastoris. Microbial Cell Fact. 11, 11. doi: 10.1186/1475-2859-11-61

PubMed Abstract | CrossRef Full Text | Google Scholar

Mizukami, A., Caron, A. L., Picanço-Castro, V., and Swiech, K. (2018). “Platforms for recombinant therapeutic glycoprotein production,” in Recombinant Glycoprotein Production: Methods and Protocols, eds V. Picanço-Castro and K. Swiech (New York, NY: Springer), 1–14. doi: 10.1007/978-1-4939-7312-5_1

PubMed Abstract | CrossRef Full Text | Google Scholar

Mordor Intelligence. (2021). Biopharmaceuticals Market - Growth, Trends, COVID-19 Impact, and Forecasts (2021 - 2026). Biopharmaceuticals market - growth, trends, covid-19 impact, and forecasts (2021 - 2026). Available online at: https://www.mordorintelligence.com/industry-reports/global-biopharmaceuticals-market-industry (accessed August 06, 2021).

Google Scholar

Navone, L., Vogl, T., Luangthongkam, P., Blinco, J. A., Luna-Flores, C., Chen, X., et al. (2021). Synergistic optimisation of expression, folding, and secretion improves E. coli AppA phytase production in Pichia pastoris. Microbial Cell Fact. 20, 1–14. doi: 10.1186/s12934-020-01499-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Nie, Y., Huang, M., Lu, J., Qian, J., Lin, W., Chu, J., et al. (2014). Impacts of high β-galactosidase expression on central metabolism of recombinant Pichia pastoris GS115 using glucose as sole carbon source via 13C metabolic flux analysis. J. Biotechnol. 187, 124–134. doi: 10.1016/j.jbiotec.2014.07.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Nielsen, J. (2013). Production of biopharmaceutical proteins by yeast: advances through metabolic engineering. Bioengineered 4, 207–211. doi: 10.4161/bioe.22856

PubMed Abstract | CrossRef Full Text | Google Scholar

Niklas, J., Priesnitz, C., Rose, T., Sandig, V., and Heinzle, E. (2013). Metabolism and metabolic burden by α1-antitrypsin production in human AGE1.HN cells. Metab. Eng. 16, 103–114. doi: 10.1016/j.ymben.2013.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Nocon, J., Steiger, M., Mairinger, T., Hohlweg, J., Rußmayer, H., Hann, S., et al. (2016). Increasing pentose phosphate pathway flux enhances recombinant protein production in Pichia pastoris. Appl. Microbiol. Biotechnol. 100, 5955–5963. doi: 10.1007/s00253-016-7363-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Noseda, D. G., Recúpero, M. N., Blasco, M., Ortiz, G. E., and Galvagno, M. A. (2013). Cloning, expression and optimized production in a bioreactor of bovine chymosin B in Pichia (Komagataella) pastoris under AOX1 promoter. Protein Express. Purif. 92, 235–244. doi: 10.1016/j.pep.2013.08.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Patra, P., Das, M., Kundu, P., and Ghosh, A. (2021). Recent advances in systems and synthetic biology approaches for developing novel cell-factories in non-conventional yeasts. Biotechno. Adv. 47, 107695. doi: 10.1016/j.biotechadv.2021.107695

PubMed Abstract | CrossRef Full Text | Google Scholar

Peña, D. A., Gasser, B., Zanghellini, J., Steiger, M. G., and Mattanovich, D. (2018). Metabolic engineering of Pichia pastoris. Metab. Eng. 50, 2–15. doi: 10.1016/j.ymben.2018.04.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Porro, D., Gasser, B., Fossati, T., Maurer, M., Branduardi, P., Sauer, M., et al. (2011). Production of recombinant proteins and metabolites in yeasts. Appl. Microbiol. Biotechnol. 89, 939–948. doi: 10.1007/s00253-010-3019-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Presnyak, V., Alhusaini, N., Chen, Y. H., Martin, S., Morris, N., Kline, N., et al. (2015). Codon optimality is a major determinant of mRNA stability. Cell 160, 1111–1124. doi: 10.1016/j.cell.2015.02.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Qin, L., Dong, S., Yu, J., Ning, X., Xu, K., Zhang, S. J., et al. (2020). Stress-driven dynamic regulation of multiple tolerance genes improves robustness and productive capacity of Saccharomyces cerevisiae in industrial lignocellulose fermentation. Metab. Eng. 61, 160–170. doi: 10.1016/j.ymben.2020.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Qin, X., Qian, J., Yao, G., Zhuang, Y., Zhang, S., and Chu, J. (2011). GAP promoter library for fine-tuning of gene expression in Pichia pastoris. Appl. Environ. Microbiol. 77, 3600–3608. doi: 10.1128/AEM.02843-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Rebello, S., Abraham, A., Madhavan, A., Sindhu, R., Binod, P., Karthika Bahuleyan, A., et al. (2018). Non-conventional yeast cell factories for sustainable bioprocesses. FEMS Microbiol. Lett. 365, 1–10. doi: 10.1093/femsle/fny222

PubMed Abstract | CrossRef Full Text | Google Scholar

Research Corporation Technologies. (2009). FDA Grants First Approval for a Product Derived From Pichia pastoris. Available online at: https://rctech.com/fda-grants-first-approval-for-a-product-derived-from-pichia-pastoris/ (accessed August 06, 2021).

Google Scholar

Research Corporation Technologies. (2019). RTC Pichia Expression Platform. Available online at: https://pichia.com/science-center/commercialized-products/ (accessed August 5, 2021).

Google Scholar

Riley, R., Haridas, S., Wolfe, K. H., Lopes, M. R., Todd, C., and Göker, M. (2016). Comparative genomics of biotechnologically important yeasts. Proc. Natl. Acad. Sci. U. S. A. doi: 10.1073/pnas.1603941113

PubMed Abstract | CrossRef Full Text | Google Scholar

Romero-Suarez, D., Wulff, T., Rong, Y., Jakočiu Nas, T., Yuzawa, S., Keasling, J. D., et al. (2021). A reporter system for cytosolic protein aggregates in yeast. ACS Synth. Biol. 10, 466–477. doi: 10.1021/acssynbio.0c00476

PubMed Abstract | CrossRef Full Text | Google Scholar

Rugbjerg, P., and Olsson, L. (2020). The future of self-selecting and stable fermentations. J. Indus. Microbiol. Biotechnol. 47, 993–1004. doi: 10.1007/s10295-020-02325-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Rugbjerg, P., and Sommer, M. O. A. (2019). Overcoming genetic heterogeneity in industrial fermentations. Nat. Biotechnol. 37, 869–876. doi: 10.1038/s41587-019-0171-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Saeki, N., Eguchi, Y., Kintaka, R., Makanae, K., Shichino, Y., Iwasaki, S., et al. (2020). N-terminal deletion of Swi3 created by the deletion of a dubious ORF YJL175W mitigates protein burden effect in S. cerevisiae. Sci. Rep. 10, 1–11. doi: 10.1038/s41598-020-66307-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Saitua, F., Torres, P., Pérez-Correa, J. R., and Agosin, E. (2017). Dynamic genome-scale metabolic modeling of the yeast Pichia pastoris. BMC Syst. Biol. 11, 27. doi: 10.1186/s12918-017-0408-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Schimmel, P. (1993). GTP hydrolysis in protein synthesis: two for Tu? Science 259, 1264–1265. doi: 10.1126/science.8446896

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Q., Yu, Z., Zhou, X. T., Zhang, S. J., Zou, S. P., Xiong, N., et al. (2021). Identification of a novel promoter for driving antibiotic-resistant genes to reduce the metabolic burden during protein expression and effectively select multiple integrations in Pichia Pastoris. Appl. Microbiol. Biotechnol. 105, 3211–3223. doi: 10.1007/s00253-021-11195-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Sleep, D., Belfield, G. P., and Goodey, A. R. (1990). The secretion of human serum albumin from the yeast Saccharomyces cerevisiae using five different leader sequnces. Biotechnology 8, 42–46. doi: 10.1038/nbt0190-42

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, H., Bao, X., Shen, Y., Song, M., Wang, S., Wang, C., et al. (2015). Engineering protein folding and translocation improves heterologous protein secretion in Saccharomyces cerevisiae. Biotechnol. Bioeng. 112, 1872–1882. doi: 10.1002/bit.25596

PubMed Abstract | CrossRef Full Text | Google Scholar

Tripathi, N. K., and Shrivastava, A. (2019). Recent developments in bioprocessing of recombinant proteins: expression hosts and process development. Front. Bioeng. Biotechnol. 7, 420. doi: 10.3389/fbioe.2019.00420

PubMed Abstract | CrossRef Full Text | Google Scholar

Tyo, K. E. J., Liu, Z., Petranovic, D., and Nielsen, J. (2012). Imbalance of heterologous protein folding and disulfide bond formation rates yields runaway oxidative stress. BMC Biol. 10, 16. doi: 10.1186/1741-7007-10-16

PubMed Abstract | CrossRef Full Text | Google Scholar

U.S. Food Drug Administration. (1982). HUMULIN-R Approval Letter. Available online at: https://www.accessdata.fda.gov/drugsatfda_docs/label/2018/018780s126s134lbl.pdf (accessed September 09, 2021).

Google Scholar

U.S. Food Drug Administration. (2018). Microorganisms and Microbial-Derived Ingredients Used in Food (Partial List). U.S. Food Drug Administration. Available online at: https://www.fda.gov/food/generally-recognized-safe-gras/microorganisms-microbial-derived-ingredients-used-food-partial-list (accessed September 12, 2021).

Google Scholar

Umebayashi, K., Hirata, A., Fukuda, R., Horiuchi, H., Ohta, A., and Takagi, M. (1997). Accumulation of misfolded protein aggregates leads to the formation of russell body-like dilated endoplasmic reticulum in yeast. Yeast 13, 1009–1020.

PubMed Abstract | Google Scholar

Valli, M., Tatto, N. E., Peymann, A., Gruber, C., Landes, N., Ekker, H., et al. (2016). Curation of the genome annotation of Pichia pastoris (Komagataella phaffii) CBS7435 from gene level to protein function. FEMS Yeast Res. 16, 1–12. doi: 10.1093/femsyr/fow051

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Rensburg, E., Den Haan, R., Smith, J., Van Zyl, W. H., and Görgens, J. F. (2012). The metabolic burden of cellulase expression by recombinant Saccharomyces cerevisiae Y294 in aerobic batch culture. Appl. Microbiol. Biotechnol. 96, 197–209. doi: 10.1007/s00253-012-4037-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanz, A. L., Nimtz, M., and Rinas, U. (2014). Decrease of UPR- and ERAD-related proteins in Pichia pastoris during methanol-induced secretory insulin precursor production in controlled fed-batch cultures. Microbial Cell Fact. 13, 1–10. doi: 10.1186/1475-2859-13-23

PubMed Abstract | CrossRef Full Text | Google Scholar

Vieira Gomes, A., Souza Carmo, T., Silva Carvalho, L., Mendonça Bahia, F., and Parachin, N. (2018). Comparison of yeasts as hosts for recombinant protein production. Microorganisms 6, 38. doi: 10.3390/microorganisms6020038

PubMed Abstract | CrossRef Full Text | Google Scholar

Vogl, T., Hartner, F. S., and Glieder, A. (2013). New opportunities by synthetic biology for biopharmaceutical production in Pichia pastoris. Curr. Opin. Biotechnol. 24, 1094–1101. doi: 10.1016/j.copbio.2013.02.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Volkova, S., Matos, M. R. A., Mattanovich, M., and de Mas, I. M. (2020). Metabolic modelling as a framework for metabolomics data integration and analysis. Metabolites 10, 1–27. doi: 10.3390/metabo10080303

PubMed Abstract | CrossRef Full Text | Google Scholar

Walsh, G. (2018). Biopharmaceutical benchmarks 2018. Nat. Biotechnol. 36, 1136–1145. doi: 10.1038/nbt.4305

PubMed Abstract | CrossRef Full Text | Google Scholar

Walter, S., and Buchner, J. (2002). Molecular chaperones - CELLULAR machines for protein folding. Angew. Chem. Int. Ed. 41, 1098–1113. doi: 10.1002/1521-3773(20020402)41:7andlt;1098::AID-ANIE1098andgt;3.0.CO;2-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Wang, X., Shi, L., Qi, F., Zhang, P., Zhang, Y., et al. (2017). Methanol-independent protein expression by AOX1 promoter with trans-acting elements engineering and glucose-glycerol-shift induction in Pichia pastoris. Sci. Rep. 7, 1–12. doi: 10.1038/srep41850

PubMed Abstract | CrossRef Full Text | Google Scholar

Waugh, D. S. (2005). Making the most of affinity tags. Trends Biotechnol. 23, 316–320. doi: 10.1016/j.tibtech.2005.03.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, N. R., Wulff, T., Palmqvist, E. A., Jørgensen, T. R., Workman, C. T., Sonnenschein, N., et al. (2020). Fluctuations in glucose availability prevent global proteome changes and physiological transition during prolonged chemostat cultivations of Saccharomyces cerevisiae. Biotechnol. Bioeng. 117, 2074–2088. doi: 10.1002/bit.27353

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, N. R., Wulff, T., Workman, C. T., Rønnest, N. P., and Sonnenschein, N. (2021). Emergence of phenotypically distinct subpopulations is a factor in adaptation of recombinant andlt;emandgt;Saccharomyces cerevisiaeandlt;/emandgt; under glucose-limited conditions. BioRxiv [Preprint]. bioRxiv 2021.01.18.427050.

Google Scholar

Xu, Y., Liu, K., Han, Y., Xing, Y., Zhang, Y., Yang, Q., et al. (2021). Codon usage bias regulates gene expression and protein conformation in yeast expression system P. pastoris. Microbial Cell Fact. 20, 1–15. doi: 10.1186/s12934-021-01580-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J. R., Chen, X., and Zhang, J. (2014). Codon-by-codon modulation of translational speed and accuracy via mRNA folding. PLoS Biol. 12, e1001910. doi: 10.1371/journal.pbio.1001910

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Z., and Zhang, Z. (2018). Engineering strategies for enhanced production of protein and bio-products in Pichia pastoris: a review. Biotechnol. Adv. 36, 182–195. doi: 10.1016/j.biotechadv.2017.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, S., Miao, L., Huang, H., Li, Y., and Zhu, T. (2020). High-level production of glucose oxidase in Pichia pastoris: effects of Hac1p overexpression on cell physiology and enzyme expression. Enzyme Microb. Technol. 141, 109671. doi: 10.1016/j.enzmictec.2020.109671

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, X. W., Sun, W. H., Wang, Y. Z., and Xu, Y. (2017). Identification of novel factors enhancing recombinant protein production in multi-copy Komagataella phaffii based on transcriptomic analysis of overexpression effects. Sci. Rep. 7, 1–12. doi: 10.1038/s41598-017-16577-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zahrl, R. J., Gasser, B., Mattanovich, D., and Ferrer, P. (2019). “Detection and elimination of cellular bottlenecks in protein-producing yeasts.” in Recombinant Protein Production in Yeast, eds B. Gasser and D. Mattanovich (New York, NY: Springer), 75–96. doi: 10.1007/978-1-4939-9024-5_2

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, W., Zhao, H. L., Xue, C., Xiong, X. H., Yao, X. Q., Li, X. Y., et al. (2006). Enhanced secretion of heterologous proteins in Pichia pastoris following overexpression of Saccharomyces cerevisiae chaperone proteins. Biotechnol. Progress 22, 1090–1095. doi: 10.1021/bp060019r

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, X., Zhang, Y., Zhang, X., Li, C., Liu, X., Lin, Y., et al. (2019). Fhl1p protein, a positive transcription factor in Pichia pastoris, enhances the expression of recombinant proteins. Microbial Cell Fact. 18, 1–12. doi: 10.1186/s12934-019-1256-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, T., Guo, M., Sun, C., Qian, J., Zhuang, Y., Chu, J., et al. (2009). A systematical investigation on the genetic stability of multi-copy Pichia pastoris strains. Biotechnol. Lett. 31, 679–684. doi: 10.1007/s10529-009-9917-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zou, W., Edros, R., and Al-Rubeai, M. (2018). The relationship of metabolic burden to productivity levels in CHO cell lines. Biotechnol. Appl. Biochem. 65, 173–180. doi: 10.1002/bab.1574

PubMed Abstract | CrossRef Full Text | Google Scholar