Codon Optimization Significantly Improves the Expression Level of α -Amylase Gene from Bacillus licheniformis in Pichia pastoris

A Systematic Review of the Potential of Pichia pastoris (Komagataella phaffii) as an Alternative Host for Biologics Production

The methylotrophic yeast Pichia pastoris is garnering interest as a chassis cell factory for the manufacture of recombinant proteins because it effectively satisfies the requirements of both laboratory and industrial set up. The optimisation of P. pastoris cultivation is still necessary due to strain- and product-specific problems such as promoter strength, methanol utilisation type, and culturing conditions to realize the high yields of heterologous protein(s) of interest. Techniques integrating genetic and process engineering have been used to overcome these problems. Insight into the Pichia as an expression system utilizing MUT pathway and the development of methanol free systems are highlighted in this systematic review. Recent developments in the improved production of proteins in P. pastoris by (i) diverse genetic engineering such as codon optimization and gene dosage; (ii) cultivating tactics including co-expression of chaperones; (iii) advances in the use of the 2A peptide system, and (iv) CRISPR/Cas technologies are widely discussed. We believe that by combining these strategies, P. pastoris will become a formidable platform for the production of high value therapeutic proteins.

Komagataella phaffii as a Platform for Heterologous Expression of Enzymes Used for Industry

In the 1980s, Escherichia coli was the preferred host for heterologous protein expression owing to its capacity for rapid growth in complex media; well-studied genetics; rapid and direct transformation with foreign DNA; and easily scalable fermentation. Despite the relative ease of use of E. coli for achieving the high expression of many recombinant proteins, for some proteins, e.g., membrane proteins or proteins of eukaryotic origin, this approach can be rather ineffective. Another microorganism long-used and popular as an expression system is baker's yeast, Saccharomyces cerevisiae. In spite of a number of obvious advantages of these yeasts as host cells, there are some limitations on their use as expression systems, for example, inefficient secretion, misfolding, hyperglycosylation, and aberrant proteolytic processing of proteins. Over the past decade, nontraditional yeast species have been adapted to the role of alternative hosts for the production of recombinant proteins, e.g., Komagataella phaffii, Yarrowia lipolytica, and Schizosaccharomyces pombe. These yeast species' several physiological characteristics (that are different from those of S. cerevisiae), such as faster growth on cheap carbon sources and higher secretion capacity, make them practical alternative hosts for biotechnological purposes. Currently, the K. phaffii-based expression system is one of the most popular for the production of heterologous proteins. Along with the low secretion of endogenous proteins, K. phaffii efficiently produces and secretes heterologous proteins in high yields, thereby reducing the cost of purifying the latter. This review will discuss practical approaches and technological solutions for the efficient expression of recombinant proteins in K. phaffii, mainly based on the example of enzymes used for the feed industry.

An exploratory in silico comparison of open-source codon harmonization tools

BACKGROUND

Not changing the native constitution of genes prior to their expression by a heterologous host can affect the amount of proteins synthesized as well as their folding, hampering their activity and even cell viability. Over the past decades, several strategies have been developed to optimize the translation of heterologous genes by accommodating the difference in codon usage between species. While there have been a handful of studies assessing various codon optimization strategies, to the best of our knowledge, no research has been performed towards the evaluation and comparison of codon harmonization algorithms. To highlight their importance and encourage meaningful discussion, we compared different open-source codon harmonization tools pertaining to their in silico performance, and we investigated the influence of different gene-specific factors.

RESULTS

In total, 27 genes were harmonized with four tools toward two different heterologous hosts. The difference in %MinMax values between the harmonized and the original sequences was calculated (ΔMinMax), and statistical analysis of the obtained results was carried out. It became clear that not all tools perform similarly, and the choice of tool should depend on the intended application. Almost all biological factors under investigation (GC content, RNA secondary structures and choice of heterologous host) had a significant influence on the harmonization results and thus must be taken into account. These findings were substantiated using a validation dataset consisting of 8 strategically chosen genes.

CONCLUSIONS

Due to the size of the dataset, no complex models could be developed. However, this initial study showcases significant differences between the results of various codon harmonization tools. Although more elaborate investigation is needed, it is clear that biological factors such as GC content, RNA secondary structures and heterologous hosts must be taken into account when selecting the codon harmonization tool.

Zika Virus Envelope Protein Domain III Produced in K. phaffii Has the Potential for Diagnostic Applications

Zika virus (ZIKV) represents a global human health threat and it is related to severe diseases such as congenital Zika syndrome (CZS) and Guillain-Barré syndrome (GBS). There is no vaccine available nor specific antiviral treatment, so developing sensitive, specific, and low-cost diagnostic tests is necessary. Thus, the objective of this work was to produce the Zika virus envelope protein domain III (ZIKV-EDIII) in Komagataella phaffii KM71H and evaluate its potential for diagnostic applications. After the K. phaffii had been transformed with the pPICZαA-ZIKV-EDIII vector, an SDS-PAGE and Western Blot were performed to characterize the recombinant protein and an ELISA to evaluate the antigenic potential. The results show that ZIKV-EDIII was produced in the expected size, with a good purity grade and yield of 2.58 mg/L. The receiver operating characteristic (ROC) curve showed 90% sensitivity and 87.5% specificity for IgM, and 93.33% sensitivity and 82.76% specificity for IgG. The ZIKV-EDIII protein was efficiently produced in K. phaffi, and it has the potential for diagnostic applications.

Optimization of Pichia pastoris Expression System for High-Level Production of Margatoxin

Margatoxin (MgTx) is a high-affinity blocker of voltage-gated potassium (Kv) channels. It inhibits Kv1.1–Kv1.3 ion channels in picomolar concentrations. This toxin is widely used to study physiological function of Kv ion channels in various cell types, including immune cells. Isolation of native MgTx in large quantities from scorpion venom is not affordable. Chemical synthesis and recombinant production in Escherichia coli need in vitro oxidative refolding for proper disulfide bond formation, resulting in a very low yield of peptide production. The Pichia pastoris expression system offers an economical approach to overcome all these limitations and gives a higher yield of correctly refolded recombinant peptides. In this study, improved heterologous expression of recombinant MgTx (rMgTx) in P. pastoris was obtained by using preferential codons, selecting the hyper-resistant clone against Zeocin, and optimizing the culturing conditions. About 36 ± 4 mg/L of >98% pure His-tagged rMgTx (TrMgTx) was produced, which is a threefold higher yield than has been previously reported. Proteolytic digestion of TrMgTx with factor Xa generated untagged rMgTx (UrMgTx). Both TrMgTx and UrMgTx blocked the Kv1.2 and Kv1.3 currents (patch-clamp) (K_d for Kv1.2 were 64 and 14 pM, and for Kv1.3, 86 and 50 pM, respectively) with comparable potency to the native MgTx. The analysis of the binding kinetics showed that TrMgTx had a lower association rate than UrMgTx for both Kv1.2 and Kv1.3. The dissociation rate of both the analogues was the same for Kv1.3. However, in the case of Kv1.2, TrMgTx showed a much higher dissociation rate with full recovery of the block than UrMgTx. Moreover, in a biological functional assay, both peptides significantly downregulated the expression of early activation markers IL2R and CD40L in activated CD4⁺ T_EM lymphocytes whose activation was Kv1.3 dependent. In conclusion, the authors report that the Pichia expression system is a powerful method to produce disulfide-rich peptides, the overexpression of which could be enhanced noticeably through optimization strategies, making it more cost-effective. Since the presence of the His-tag on rMgTx only mildly altered the block equilibrium and binding kinetics, recombinant toxins could be used in ion channel research without removing the tag and could thus reduce the cost and time demand for toxin production.

Bioproduced Proteins On Demand (Bio-POD) in hydrogels using Pichia pastoris

Traditional production of industrial and therapeutic proteins by eukaryotic cells typically requires large-scale fermentation capacity. As a result, these systems are not easily portable or reusable for on-demand protein production applications. In this study, we employ Bioproduced Proteins On Demand (Bio-POD), a F127-bisurethane methacrylate hydrogel-based technique that immobilizes engineered Pichia pastoris for preservable, on-demand production and secretion of medium- and high-molecular weight proteins (in this case, SEAP, α-amylase, and anti-HER2). The gel samples containing encapsulated-yeast demonstrated sustained protein production and exhibited productivity immediately after lyophilization and rehydration. The hydrogel platform described here is the first hydrogel immobilization using a P. pastoris system to produce recombinant proteins of this breadth. These results highlight the potential of this formulation to establish a cost-effective bioprocessing strategy for on-demand protein production.

Genetic and metabolic engineering challenges of C1-gas fermenting acetogenic chassis organisms

Unabated mining and utilisation of petroleum and petroleum resources and their conversion to essential fuels and chemicals have drastic environmental consequences, contributing to global warming and climate change. In addition, fossil fuels are finite resources, with a fast-approaching shortage. Accordingly, research efforts are increasingly focusing on developing sustainable alternatives for chemicals and fuels production. In this context, bioprocesses, relying on microorganisms, have gained particular interest. For example, acetogens use the Wood-Ljungdahl pathway to grow on single carbon C1-gases (CO2 and CO) as their sole carbon source and produce valuable products such as acetate or ethanol. These autotrophs can, therefore, be exploited for large-scale fermentation processes to produce industrially relevant chemicals from abundant greenhouse gases. In addition, genetic tools have recently been developed to improve these chassis organisms through synthetic biology approaches. This review will focus on the challenges of genetically and metabolically modifying acetogens. It will first discuss the physical and biochemical obstacles complicating successful DNA transfer in these organisms. Current genetic tools developed for several acetogens, crucial for strain engineering to consolidate and expand their catalogue of products, will then be described. Recent tool applications for metabolic engineering purposes to allow redirection of metabolic fluxes or production of non-native compounds will lastly be covered.

High-efficiency expression and secretion of human FGF21 in Bacillus subtilis by intercalation of a mini-cistron cassette and combinatorial optimization of cell regulatory components

Background

Recombinant human Fibroblast growth factor 21 (rhFGF21) is an endocrine hormone that has profound effects on treatment of metabolic diseases. However, rhFGF21 is prone to form inclusion body when expressed in bacteria, which results in, the downstream process of purification of bioactive rhFGF21 is time-consuming and labor intensive. The aim of this work is to explore a new method for improving the soluble expression and secretion level of rhFGF21 in B. subtilis.

Results

A codon optimized rhFGF21 gene was expressed under the control of a strong inducible promoter P_malA in B. subtilis. A mini-cistron cassette (from gsiB) was located upstream of rhFGF21 in expression vector (pMATEFc5), which could reduce the locally stabilized mRNA secondary structure of transcripts and enhance the efficiency of translation initiation. Then various chaperones were further overexpressed to improve the expression efficiency of rhFGF21. Results showed that overexpression of the chaperone DnaK contributed to the increase of solubility of rhFGF21. Moreover, an extracellular proteases deficient strain B. subtilis Kno6cf was used to accumulate the secreted rhFGF21 solidly. In addition, eleven signal peptides from B. subtilis were evaluated and the SP_dacB appeared the highest secretion yield of rhFGF21 in B. subtilis. Finally, the combinatorial optimized strain achieved an about ninefold increase of the soluble rhFGF21 production after 24 h of flask fermentation in comparison with the initial production strain.

Conclusion

This work provided a comprehensive strategy for secretory expressing the heterologous protein rhFGF21 in B. subtilis. To our knowledge, this is the first report of the highly efficient production of rhFGF21 in B. subtilis and this approach may provide some suggestions for heterologous proteins production in B. subtilis.

Electronic supplementary material

The online version of this article (10.1186/s12934-019-1066-4) contains supplementary material, which is available to authorized users.

Integrated analysis of individual codon contribution to protein biosynthesis reveals a new approach to improving the basis of rational gene design

Gene codon optimization may be impaired by the misinterpretation of frequency and optimality of codons. Although recent studies have revealed the effects of codon usage bias (CUB) on protein biosynthesis, an integrated perspective of the biological role of individual codons remains unknown. Unlike other previous studies, we show, through an integrated framework that attributes of codons such as frequency, optimality and positional dependency should be combined to unveil individual codon contribution for protein biosynthesis. We designed a codon quantification method for assessing CUB as a function of position within genes with a novel constraint: the relativity of position-dependent codon usage shaped by coding sequence length. Thus, we propose a new way of identifying the enrichment, depletion and non-uniform positional distribution of codons in different regions of yeast genes. We clustered codons that shared attributes of frequency and optimality. The cluster of non-optimal codons with rare occurrence displayed two remarkable characteristics: higher codon decoding time than frequent-non-optimal cluster and enrichment at the 5'-end region, where optimal codons with the highest frequency are depleted. Interestingly, frequent codons with non-optimal adaptation to tRNAs are uniformly distributed in the Saccharomyces cerevisiae genes, suggesting their determinant role as a speed regulator in protein elongation.

Increased dosage of AOX1 promoter-regulated expression cassettes leads to transcription attenuation of the methanol metabolism in Pichia pastoris

The methanol-regulated alcohol oxidase promoter (P_AOX1) of Pichia pastoris is one of the strongest promoters for heterologous gene expression in this methylotrophic yeast. Although increasing gene dosage is one of the most common strategies to increase recombinant protein productivities, the increase of gene dosage of Rhizopus oryzae lipase (ROL) in P. pastoris has been previously shown to reduce cell growth, lipase production and substrate consumption in high-copy strains. To better assess that physiological response, transcriptomics analysis was performed of a subset of strains with 1 to 15 ROL copies. The macroscopic physiological parameters confirm that growth yield and carbon uptake rate are gene dosage dependent, and were supported by the transcriptomic data, showing the impact of increased dosage of AOX1 promoter-regulated expression cassettes on P. pastoris physiology under steady methanolic growth conditions. Remarkably, increased number of cassettes led to transcription attenuation of the methanol metabolism and peroxisome biogenesis in P. pastoris, concomitant with reduced secretion levels of the heterologous product. Moreover, our data also point to a block in ROL mRNA translation in the higher ROL-copies constructs, while the low productivities of multi-copy strains under steady growth conditions do not appear to be directly related to UPR and ERAD induction.

Dynamic genome-scale metabolic modeling of the yeast Pichia pastoris

BACKGROUND

Pichia pastoris shows physiological advantages in producing recombinant proteins, compared to other commonly used cell factories. This yeast is mostly grown in dynamic cultivation systems, where the cell's environment is continuously changing and many variables influence process productivity. In this context, a model capable of explaining and predicting cell behavior for the rational design of bioprocesses is highly desirable. Currently, there are five genome-scale metabolic reconstructions of P. pastoris which have been used to predict extracellular cell behavior in stationary conditions.

RESULTS

In this work, we assembled a dynamic genome-scale metabolic model for glucose-limited, aerobic cultivations of Pichia pastoris. Starting from an initial model structure for batch and fed-batch cultures, we performed pre/post regression diagnostics to ensure that model parameters were identifiable, significant and sensitive. Once identified, the non-relevant ones were iteratively fixed until a priori robust modeling structures were found for each type of cultivation. Next, the robustness of these reduced structures was confirmed by calibrating the model with new datasets, where no sensitivity, identifiability or significance problems appeared in their parameters. Afterwards, the model was validated for the prediction of batch and fed-batch dynamics in the studied conditions. Lastly, the model was employed as a case study to analyze the metabolic flux distribution of a fed-batch culture and to unravel genetic and process engineering strategies to improve the production of recombinant Human Serum Albumin (HSA). Simulation of single knock-outs indicated that deviation of carbon towards cysteine and tryptophan formation improves HSA production. The deletion of methylene tetrahydrofolate dehydrogenase could increase the HSA volumetric productivity by 630%. Moreover, given specific bioprocess limitations and strain characteristics, the model suggests that implementation of a decreasing specific growth rate during the feed phase of a fed-batch culture results in a 25% increase of the volumetric productivity of the protein.

CONCLUSION

In this work, we formulated a dynamic genome scale metabolic model of Pichia pastoris that yields realistic metabolic flux distributions throughout dynamic cultivations. The model can be calibrated with experimental data to rationally propose genetic and process engineering strategies to improve the performance of a P. pastoris strain of interest.

Improved Production of Aspergillus usamii endo-β-1,4-Xylanase in Pichia pastoris via Combined Strategies

A series of strategies were applied to improve expression level of recombinant endo-β-1,4-xylanase from Aspergillus usamii (A. usamii) in Pichia pastoris (P. pastoris). Firstly, the endo-β-1,4-xylanase (xynB) gene from A. usamii was optimized for P. pastoris and expressed in P. pastoris. The maximum xylanase activity of optimized (xynB-opt) gene was 33500 U/mL after methanol induction for 144 h in 50 L bioreactor, which was 59% higher than that by wild-type (xynB) gene. To further increase the expression of xynB-opt, the Vitreoscilla hemoglobin (VHb) gene was transformed to the recombinant strain containing xynB-opt. The results showed that recombinant strain harboring the xynB-opt and VHb (named X33/xynB-opt-VHb) displayed higher biomass, cell viability, and xylanase activity. The maximum xylanase activity of X33/xynB-opt-VHb in 50 L bioreactor was 45225 U/mL, which was 35% and 115% higher than that by optimized (xynB-opt) gene and wild-type (xynB) gene. Finally, the induction temperature of X33/xynB-opt-VHb was optimized in 50 L bioreactor. The maximum xylanase activity of X33/xynB-opt-VHb reached 58792 U/mL when the induction temperature was 22°C. The results presented here will greatly contribute to improving the production of recombinant proteins in P. pastoris.