Hybrid metabolic flux analysis and recombinant protein prediction in Pichia pastoris X-33 cultures expressing a single-chain antibody fragment

Precise metabolic modeling in post-omics era: accomplishments and perspectives

Microbes have been extensively utilized for their sustainable and scalable properties in synthesizing desired bio-products. However, insufficient knowledge about intracellular metabolism has impeded further microbial applications. The genome-scale metabolic models (GEMs) play a pivotal role in facilitating a global understanding of cellular metabolic mechanisms. These models enable rational modification by exploring metabolic pathways and predicting potential targets in microorganisms, enabling precise cell regulation without experimental costs. Nonetheless, simplified GEM only considers genome information and network stoichiometry while neglecting other important bio-information, such as enzyme functions, thermodynamic properties, and kinetic parameters. Consequently, uncertainties persist particularly when predicting microbial behaviors in complex and fluctuant systems. The advent of the omics era with its massive quantification of genes, proteins, and metabolites under various conditions has led to the flourishing of multi-constrained models and updated algorithms with improved predicting power and broadened dimension. Meanwhile, machine learning (ML) has demonstrated exceptional analytical and predictive capacities when applied to training sets of biological big data. Incorporating the discriminant strength of ML with GEM facilitates mechanistic modeling efficiency and improves predictive accuracy. This paper provides an overview of research innovations in the GEM, including multi-constrained modeling, analytical approaches, and the latest applications of ML, which may contribute comprehensive knowledge toward genetic refinement, strain development, and yield enhancement for a broad range of biomolecules.

A General Hybrid Modeling Framework for Systems Biology Applications: Combining Mechanistic Knowledge with Deep Neural Networks under the SBML Standard

In this paper, a computational framework is proposed that merges mechanistic modeling with deep neural networks obeying the Systems Biology Markup Language (SBML) standard. Over the last 20 years, the systems biology community has developed a large number of mechanistic models that are currently stored in public databases in SBML. With the proposed framework, existing SBML models may be redesigned into hybrid systems through the incorporation of deep neural networks into the model core, using a freely available python tool. The so-formed hybrid mechanistic/neural network models are trained with a deep learning algorithm based on the adaptive moment estimation method (ADAM), stochastic regularization and semidirect sensitivity equations. The trained hybrid models are encoded in SBML and uploaded in model databases, where they may be further analyzed as regular SBML models. This approach is illustrated with three well-known case studies: the Escherichia coli threonine synthesis model, the P58IPK signal transduction model, and the Yeast glycolytic oscillations model. The proposed framework is expected to greatly facilitate the widespread use of hybrid modeling techniques for systems biology applications.

SBML2HYB: a Python interface for SBML compatible hybrid modeling

SUMMARY

Here, we present sbml2hyb, an easy-to-use standalone Python tool that facilitates the conversion of existing mechanistic models of biological systems in Systems Biology Markup Language (SBML) into hybrid semiparametric models that combine mechanistic functions with machine learning (ML). The so-formed hybrid models can be trained and stored back in databases in SBML format. The tool supports a user-friendly export interface with an internal format validator. Two case studies illustrate the use of the sbml2hyb tool. Additionally, we describe HMOD, a new model format designed to support and facilitate hybrid models building. It aggregates the mechanistic model information with the ML information and follows as close as possible the SBML rules. We expect the sbml2hyb tool and HMOD to greatly facilitate the widespread usage of hybrid modeling techniques for biological systems analysis.

AVAILABILITY AND IMPLEMENTATION

The Python interface, source code and the example models used for the case studies are accessible at: https://github.com/r-costa/sbml2hyb.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-scale modeling of Chinese hamster ovary cells by hybrid semi-parametric flux balance analysis

Flux balance analysis (FBA) is currently the standard method to compute metabolic fluxes in genome-scale networks. Several FBA extensions employing diverse objective functions and/or constraints have been published. Here we propose a hybrid semi-parametric FBA extension that combines mechanistic-level constraints (parametric) with empirical constraints (non-parametric) in the same linear program. A CHO dataset with 27 measured exchange fluxes obtained from 21 reactor experiments served to evaluate the method. The mechanistic constraints were deduced from a reduced CHO-K1 genome-scale network with 686 metabolites, 788 reactions and 210 degrees of freedom. The non-parametric constraints were obtained by principal component analysis of the flux dataset. The two types of constraints were integrated in the same linear program showing comparable computational cost to standard FBA. The hybrid FBA is shown to significantly improve the specific growth rate prediction under different constraints scenarios. A metabolically efficient cell growth feed targeting minimal byproducts accumulation was designed by hybrid FBA. It is concluded that integrating parametric and nonparametric constraints in the same linear program may be an efficient approach to reduce the solution space and to improve the predictive power of FBA methods when critical mechanistic information is missing.

Recombinant protein production in Pichia pastoris: From transcriptionally redesigned strains to bioprocess optimization and metabolic modelling

Pichia pastoris is one of the most widely used host for the production of recombinant proteins. Expression systems that rely mostly on promoters from genes encoding alcohol oxidase 1 or glyceraldehyde-3-phosphate dehydrogenase have been developed together with related bioreactor operation strategies based on carbon sources such as methanol, glycerol, or glucose. Although, these processes are relatively efficient and easy to use, there have been notable improvements over the last twenty years to better control gene expression from these promoters and their engineered variants. Methanol-free and more efficient protein production platforms have been developed by engineering promoters and transcription factors. The production window of P. pastoris has been also extended by using alternative feedstocks including ethanol, lactic acid, mannitol, sorbitol, sucrose, xylose, gluconate, formate, or rhamnose. Herein, the specific aspects that are emerging as key parameters for recombinant protein synthesis are discussed. For this purpose, a holistic approach has been considered to scrutinize protein production processes from strain design to bioprocess optimization, particularly focusing on promoter engineering, transcriptional circuitry redesign. This review also considers the optimization of bioprocess based on alternative carbon sources and derived co-feeding strategies. Optimization strategies for recombinant protein synthesis through metabolic modelling are also discussed.

Model-assisted DoE software: optimization of growth and biocatalysis in Saccharomyces cerevisiae bioprocesses

Bioprocess development and optimization are still cost- and time-intensive due to the enormous number of experiments involved. In this study, the recently introduced model-assisted Design of Experiments (mDoE) concept (Möller et al. in Bioproc Biosyst Eng 42(5):867, https://doi.org/10.1007/s00449-019-02089-7 , 2019) was extended and implemented into a software ("mDoE-toolbox") to significantly reduce the number of required cultivations. The application of the toolbox is exemplary shown in two case studies with Saccharomyces cerevisiae. In the first case study, a fed-batch process was optimized with respect to the pH value and linearly rising feeding rates of glucose and nitrogen source. Using the mDoE-toolbox, the biomass concentration was increased by 30% compared to previously performed experiments. The second case study was the whole-cell biocatalysis of ethyl acetoacetate (EAA) to (S)-ethyl-3-hydroxybutyrate (E3HB), for which the feeding rates of glucose, nitrogen source, and EAA were optimized. An increase of 80% compared to a previously performed experiment with similar initial conditions was achieved for the E3HB concentration.

Towards a widespread adoption of metabolic modeling tools in biopharmaceutical industry: a process systems biology engineering perspective

In biotechnology, the emergence of high-throughput technologies challenges the interpretation of large datasets. One way to identify meaningful outcomes impacting process and product attributes from large datasets is using systems biology tools such as metabolic models. However, these tools are still not fully exploited for this purpose in industrial context due to gaps in our knowledge and technical limitations. In this paper, key aspects restraining the routine implementation of these tools are highlighted in three research fields: monitoring, network science and hybrid modeling. Advances in these fields could expand the current state of systems biology applications in biopharmaceutical industry to address existing challenges in bioprocess development and improvement.

Systems biology approach in the formulation of chemically defined media for recombinant protein overproduction

The cell culture medium is an intricate mixture of components which has a tremendous effect on cell growth and recombinant protein production. Regular cell culture medium includes various components, and the decision about which component should be included in the formulation and its optimum amount is an underlying issue in biotechnology industries. Applying conventional techniques to design an optimal medium for the production of a recombinant protein requires meticulous and immense research. Moreover, since the medium formulation for the production of one protein could not be the best choice for another protein, hence, the most suitable media should be determined for each recombinant cell line. Accordingly, medium formulation becomes a laborious, time-consuming, and costly process in biomanufacturing of recombinant protein, and finding alternative strategies for medium development seems to be crucial. In silico modeling is an attractive concept to be adapted for medium formulation due to its high potential to supersede laboratory examinations. By emerging the high-throughput datasets, scientists can disclose the knowledge about the effect of medium components on cell growth and metabolism, and via applying this information through systems biology approach, medium formulation optimization could be accomplished in silico with no need of significant amount of experimentation. This review demonstrates some of the applications of systems biology as a powerful tool for medium development and illustrates the effect of medium optimization with system-level analysis on the production of recombinant proteins in different host cells.

Methyl-selective isotope labeling using α-ketoisovalerate for the yeast Pichia pastoris recombinant protein expression system

Methyl-detected NMR spectroscopy is a useful tool for investigating the structures and interactions of large macromolecules such as membrane proteins. The procedures for preparation of methyl-specific isotopically-labeled proteins were established for the Escherichia coli (E. coli) expression system, but typically it is not feasible to express eukaryotic proteins using E. coli. The Pichia pastoris (P. pastoris) expression system is the most common yeast expression system, and is known to be superior to the E. coli system for the expression of mammalian proteins, including secretory and membrane proteins. However, this system has not yet been optimized for methyl-specific isotope labeling, especially for Val/Leu-methyl specific isotope incorporation. To overcome this difficulty, we explored various culture conditions for the yeast cells to efficiently uptake Val/Leu precursors. Among the searched conditions, we found that the cultivation pH has a critical effect on Val/Leu precursor uptake. At an acidic cultivation pH, the uptake of the Val/Leu precursor was increased, and methyl groups of Val and Leu in the synthesized recombinant protein yielded intense ¹H-¹³C correlation signals. Based on these results, we present optimized protocols for the Val/Leu-methyl-selective ¹³C incorporation by the P. pastoris expression system.

Hybrid semiparametric systems for quantitative sequence-activity modeling of synthetic biological parts

Predicting the activity of modified biological parts is difficult due to the typically large size of nucleotide sequences, resulting in combinatorial designs that suffer from the "curse of dimensionality" problem. Mechanistic design methods are often limited by knowledge availability. Empirical methods typically require large data sets, which are difficult and/or costly to obtain. In this study, we explore for the first time the combination of both approaches within a formal hybrid semiparametric framework in an attempt to overcome the limitations of the current approaches. Protein translation as a function of the 5' untranslated region sequence in Escherichia coli is taken as case study. Thermodynamic modeling, partial least squares (PLS) and hybrid parallel combinations thereof are compared for different data sets and data partitioning scenarios. The results suggest a significant and systematic reduction of both calibration and prediction errors by the hybrid approach in comparison to standalone thermodynamic or PLS modeling. Although with different magnitudes, improvements are observed irrespective of sample size and partitioning method. All in all the results suggest an increase of predictive power by the hybrid method potentially leading to a more efficient design of biological parts.

Integrating metabolic modeling and population heterogeneity analysis into optimizing recombinant protein production by Komagataella (Pichia) pastoris

The methylotrophic yeast Komagataella (Pichia) pastoris has become one of the most utilized cell factories for the production of recombinant proteins over the last three decades. This success story is linked to its specific physiological traits, i.e., the ability to grow at high cell density in inexpensive culture medium and to secrete proteins at high yield. Exploiting methanol metabolism is at the core of most P. pastoris-based processes but comes with its own challenges. Co-feeding cultures with glycerol/sorbitol and methanol is a promising approach, which can benefit from improved understanding and prediction of metabolic response. The development of profitable processes relies on the construction and selection of efficient producing strains from less efficient ones but also depends on the ability to master the bioreactor process itself. More specifically, how a bioreactor processes could be monitored and controlled to obtain high yield of production. In this review, new perspectives are detailed regarding a multi-faceted approach to recombinant protein production processes by P. pastoris; including gaining improved understanding of the metabolic pathways involved, accounting for variations in transcriptional and translational efficiency at the single cell level and efficient monitoring and control of methanol levels at the bioreactor level.

Engineering strategies for enhanced production of protein and bio-products in Pichia pastoris: A review

Pichia pastoris has been recognized as one of the most industrially important hosts for heterologous protein production. Despite its high protein productivity, the optimization of P. pastoris cultivation is still imperative due to strain- and product-specific challenges such as promoter strength, methanol utilization type and oxygen demand. To address the issues, strategies involving genetic and process engineering have been employed. Optimization of codon usage and gene dosage, as well as engineering of promoters, protein secretion pathways and methanol metabolic pathways have proved beneficial to innate protein expression levels. Large-scale production of proteins via high cell density fermentation additionally relies on the optimization of process parameters including methanol feed rate, induction temperature and specific growth rate. Recent progress related to the enhanced production of proteins in P. pastoris via various genetic engineering and cultivation strategies are reviewed. Insight into the regulation of the P. pastoris alcohol oxidase 1 (AOX1) promoter and the development of methanol-free systems are highlighted. Novel cultivation strategies such as mixed substrate feeding are discussed. Recent advances regarding substrate and product monitoring techniques are also summarized. Application of P. pastoris to the production of biodiesel and other value-added products via metabolic engineering are also reviewed. P. pastoris is becoming an indispensable platform through the use of these combined engineering strategies.

Analysis of culture media screening data by projection to latent pathways: The case of Pichia pastoris X-33

Cell culture media formulations contain hundreds of individual components in water solutions which have complex interactions with metabolic pathways. The currently used statistical design methods are empirical and very limited to explore such a large design space. In a previous work we developed a computational method called projection to latent pathways (PLP), which was conceived to maximize covariance between envirome and fluxome data under the constraint of metabolic network elementary flux modes (EFM). More specifically, PLP identifies a minimal set of EFMs (i.e., pathways) with the highest possible correlation with envirome and fluxome measurements. In this paper we extend the concept for the analysis of culture media screening data to investigate how culture medium components up-regulate or down-regulate key metabolic pathways. A Pichia pastoris X-33 strain was cultivated in 26 shake flask experiments with variations in trace elements concentrations and basal medium dilution, based on the standard BSM+PTM1 medium. PLP identified 3 EFMs (growth, maintenance and by-product formation) describing 98.8% of the variance in observed fluxes. Furthermore, PLP presented an overall predictive power comparable to that of PLS regression. Our results show iron and manganese at concentrations close to the PTM1 standard inhibit overall metabolic activity, while the main salts concentration (BSM) affected mainly energy expenditures for cellular maintenance.