Kinetic models in industrial biotechnology - Improving cell factory performance

Ogataea polymorpha as a next-generation chassis for industrial biotechnology

Ogataea (Hansenula) polymorpha is a nonconventional yeast with some unique characteristics, including fast growth, thermostability, and broad substrate spectrum. Other than common applications for protein production, O. polymorpha is attracting interest for chemical and protein production from methanol; a promising feedstock for the next-generation biomanufacturing due to its abundant sources and excellent characteristics. Benefiting from the development of synthetic biology, it has been engineered to produce value-added chemicals by extensively rewiring cellular metabolism. This Review discusses recently developed synthetic biology tools of O. polymorpha. The advances of chemicals production and systems biology were reviewed comprehensively. Finally, we look ahead to the developments of biomanufacturing in O. polymorpha to make an overall understanding of this chassis for academia and industry.

Kinetic modeling and optimization of ethanol fermentation by the marine yeast Wickerhamomyces subpelliculosus ZE75

The study aims to enhance ethanol production by Wickerhamomyces subpelliculosus ZE75 isolated from marine sediment. In addition, analyzing the kinetic parameters of ethanol production and optimization of the fermentation conditions was performed. The marine yeast isolate ZE75 was selected as the front runner ethanol-producer, with an ethanol yield of 89.77 gL-1. ZE75 was identified relying on the phenotypic and genotypic characteristics of W. subpelliculosus. The genotypic characterization based on the Internal Transcribed Spacer (ITS) sequence was deposited in the GenBank database with the accession number OP715873. The maximum specific ethanol production rate (vmax) was 0.482 gg-1 h-1 at 175 gL-1 glucose concentration, with a high accuracy of R2 0.95. The maximum growth specific rates (μmax) were 0.141 h-1 obtained at 150 gL-1 glucose concentration with R2 0.91. Optimization of the fermentation parameters such as pH and salinity has been achieved. The highest ethanol yield 0.5637 gg-1 was achieved in a 100% natural seawater-based medium. The maximum ethanol production of 104.04 gL-1 was achieved at pH 4.5 with a specific ethanol rate of 0.1669 gg-1 h-1. The findings of the present study recommend the possibility of ethanol production from a seawater-based medium on a large scale using W. subpelliculosus ZE75.

Pushing the boundaries of phosphorylase cascade reaction for cellobiose production II: Model-based multiobjective optimization

The inherent complexity of coupled biocatalytic reactions presents a major challenge for process development with one-pot multienzyme cascade transformations. Kinetic models are powerful engineering tools to guide the optimization of cascade reactions towards a performance suitable for scale up to an actual production. Here, we report kinetic model-based window of operation analysis for cellobiose production (≥100 g/L) from sucrose and glucose by indirect transglycosylation via glucose 1-phosphate as intermediate. The two-step cascade transformation is catalyzed by sucrose and cellobiose phosphorylase in the presence of substoichiometric amounts of phosphate (≤27 mol% of substrate). Kinetic modeling was instrumental to uncover the hidden effect of bulk microviscosity due to high sugar concentrations on decreasing the rate of cellobiose phosphorylase specifically. The mechanistic-empirical hybrid model thus developed gives a comprehensive description of the cascade reaction at industrially relevant substrate conditions. Model simulations serve to unravel opposed relationships between efficient utilization of the enzymes and maximized concentration (or yield) of the product within a given process time, in dependence of the initial concentrations of substrate and phosphate used. Optimum balance of these competing key metrics of process performance is suggested from the model-calculated window of operation and is verified experimentally. The evidence shown highlights the important use of kinetic modeling for the characterization and optimization of cascade reactions in ways that appear to be inaccessible to purely data-driven approaches.

Pushing the boundaries of phosphorylase cascade reaction for cellobiose production I: Kinetic model development

One-pot cascade reactions of coupled disaccharide phosphorylases enable an efficient transglycosylation via intermediary α-d-glucose 1-phosphate (G1P). Such transformations have promising applications in the production of carbohydrate commodities, including the disaccharide cellobiose for food and feed use. Several studies have shown sucrose and cellobiose phosphorylase for cellobiose synthesis from sucrose, but the boundaries on transformation efficiency that result from kinetic and thermodynamic characteristics of the individual enzyme reactions are not known. Here, we assessed in a step-by-step systematic fashion the practical requirements of a kinetic model to describe cellobiose production at industrially relevant substrate concentrations of up to 600 mM sucrose and glucose each. Mechanistic initial-rate models of the two-substrate reactions of sucrose phosphorylase (sucrose + phosphate → G1P + fructose) and cellobiose phosphorylase (G1P + glucose → cellobiose + phosphate) were needed and additionally required expansion by terms of glucose inhibition, in particular a distinctive two-site glucose substrate inhibition of the cellobiose phosphorylase (from Cellulumonas uda). Combined with mass action terms accounting for the approach to equilibrium, the kinetic model gave an excellent fit and a robust prediction of the full reaction time courses for a wide range of enzyme activities as well as substrate concentrations, including the variable substoichiometric concentration of phosphate. The model thus provides the essential engineering tool to disentangle the highly interrelated factors of conversion efficiency in the coupled enzyme reaction; and it establishes the necessary basis of window of operation calculations for targeted optimizations toward different process tasks.

Rational strain design with minimal phenotype perturbation

Devising genetic interventions for desired cellular phenotypes remains challenging regarding time and resources. Kinetic models can accelerate this task by simulating metabolic responses to genetic perturbations. However, exhaustive design evaluations with kinetic models are computationally impractical, especially when targeting multiple enzymes. Here, we introduce a framework for efficiently scouting the design space while respecting cellular physiological requirements. The framework employs mixed-integer linear programming and nonlinear simulations with large-scale nonlinear kinetic models to devise genetic interventions while accounting for the network effects of these perturbations. Importantly, it ensures the engineered strain's robustness by maintaining its phenotype close to that of the reference strain. The framework, applied to improve the anthranilate production in E. coli, devises designs for experimental implementation, including eight previously experimentally validated targets. We expect this framework to play a crucial role in future design-build-test-learn cycles, significantly expediting the strain design compared to exhaustive design enumeration.

Modelling of microbial interactions in anaerobic digestion: from black to glass box

Anaerobic and microaerophilic environments are pervasive in nature, providing essential contributions to the maintenance of human health, biogeochemical cycles and the Earth's climate. These ecological niches are characterised by low free oxygen and oxidants, or lack thereof. Under these conditions, interactions between species are essential for supporting the growth of syntrophic species and maintaining thermodynamic feasibility of anaerobic fermentation. Kinetic models provide a simplified view of complex metabolic networks, while genome-scale metabolic models and flux-balance analysis (FBA) aim to unravel these systems as a whole. The target of this review is to outline the main similarities, differences and challenges associated with kinetic and metabolic modelling, and describe state-of-the-art modelling practices for studying syntrophies in the anaerobic digestion (AD) case study.

From observational to actionable: rethinking omics in biologics production

As the era of omics continues to expand with increasing ubiquity and success in both academia and industry, omics-based experiments are becoming commonplace in industrial biotechnology, including efforts to develop novel solutions in bioprocess optimization and cell line development. Omic technologies provide particularly valuable 'observational' insights for discovery science, especially in academic research and industrial R&D; however, biomanufacturing requires a different paradigm to unlock 'actionable' insights from omics. Here, we argue the value of omic experiments in biotechnology can be maximized with deliberate selection of omic approaches and forethought about analysis techniques. We describe important considerations when designing and implementing omic-based experiments and discuss how systems biology analysis strategies can enhance efforts to obtain actionable insights in mammalian-based biologics production.

Metabolic control analysis as a strategy to identify therapeutic targets, the case of cancer glycolysis

The use of kinetic modeling and metabolic control analysis (MCA) to identify possible therapeutic targets and to investigate the controlling and regulatory mechanisms in cancer glycolysis is here reviewed. The glycolytic pathway has been considered a target to decrease cancer cell growth; however, its occurrence in normal cells makes it difficult to design therapeutic strategies that target this pathway in pathological cells. Notwithstanding, the over-expression of all enzymes and transporters, as well as the expression of isoenzymes with different kinetic and regulatory properties in cancer cells, suggested a different distribution of the control of glycolytic flux than that observed in normal cells. Kinetic models of glycolysis are constructed with enzyme kinetics experimental data, validated with the steady-state metabolite concentrations and glycolytic fluxes; applying MCA, permitted us to identify the steps with the highest control of glycolysis in cancer cells, but low control in normal cells. The cancer glycolysis main controlling steps under several metabolic conditions were: glucose transport, hexokinase and hexose-6-phosphate isomerase (HPI); whereas in normal cells were: the first two and phosphofructokinase-1. HPI is the best therapeutic target because it exerts high control in cancer glycolytic flux, but not in normal cells. Furthermore, kinetic modeling also contributed to identifying new feed-back and feed-forward regulatory loops in cancer cells glycolysis, and to understanding the mode of metabolic action of glycolytic inhibitors. Thus, MCA and metabolic modeling allowed to propose new strategies for inhibiting glycolysis in cancer cells.

Machine learning in bioprocess development: from promise to practice

Fostered by novel analytical techniques, digitalization, and automation, modern bioprocess development provides large amounts of heterogeneous experimental data, containing valuable process information. In this context, data-driven methods like machine learning (ML) approaches have great potential to rationally explore large design spaces while exploiting experimental facilities most efficiently. Herein we demonstrate how ML methods have been applied so far in bioprocess development, especially in strain engineering and selection, bioprocess optimization, scale-up, monitoring, and control of bioprocesses. For each topic, we will highlight successful application cases, current challenges, and point out domains that can potentially benefit from technology transfer and further progress in the field of ML.

Optimal enzyme utilization suggests that concentrations and thermodynamics determine binding mechanisms and enzyme saturations

Deciphering the metabolic functions of organisms requires understanding the dynamic responses of living cells upon genetic and environmental perturbations, which in turn can be inferred from enzymatic activity. In this work, we investigate the optimal modes of operation for enzymes in terms of the evolutionary pressure driving them toward increased catalytic efficiency. We develop a framework using a mixed-integer formulation to assess the distribution of thermodynamic forces and enzyme states, providing detailed insights into the enzymatic mode of operation. We use this framework to explore Michaelis-Menten and random-ordered multi-substrate mechanisms. We show that optimal enzyme utilization is achieved by unique or alternative operating modes dependent on reactant concentrations. We find that in a bimolecular enzyme reaction, the random mechanism is optimal over any other ordered mechanism under physiological conditions. Our framework can investigate the optimal catalytic properties of complex enzyme mechanisms. It can further guide the directed evolution of enzymes and fill in the knowledge gaps in enzyme kinetics.

Elucidating yeast glycolytic dynamics at steady state growth and glucose pulses through kinetic metabolic modeling

Microbial cell factories face changing environments during industrial fermentations. Kinetic metabolic models enable the simulation of the dynamic metabolic response to these perturbations, but their development is challenging due to model complexity and experimental data requirements. An example of this is the well-established microbial cell factory Saccharomyces cerevisiae, for which no consensus kinetic model of central metabolism has been developed and implemented in industry. Here, we aim to bring the academic and industrial communities closer to this consensus model. We developed a physiology informed kinetic model of yeast glycolysis connected to central carbon metabolism by including the effect of anabolic reactions precursors, mitochondria and the trehalose cycle. To parametrize such a large model, a parameter estimation pipeline was developed, consisting of a divide and conquer approach, supplemented with regularization and global optimization. Additionally, we show how this first mechanistic description of a growing yeast cell captures experimental dynamics at different growth rates and under a strong glucose perturbation, is robust to parametric uncertainty and explains the contribution of the different pathways in the network. Such a comprehensive model could not have been developed without using steady state and glucose perturbation data sets. The resulting metabolic reconstruction and parameter estimation pipeline can be applied in the future to study other industrially-relevant scenarios. We show this by generating a hybrid CFD-metabolic model to explore intracellular glycolytic dynamics for the first time. The model suggests that all intracellular metabolites oscillate within a physiological range, except carbon storage metabolism, which is sensitive to the extracellular environment.

Machine learning in fermentative biohydrogen production: Advantages, challenges, and applications

Hydrogen can be produced in an environmentally friendly manner through biological processes using a variety of organic waste and biomass as feedstock. However, the complexity of biological processes limits their predictability and reliability, which hinders the scale-up and dissemination. This article reviews contemporary research and perspectives on the application of machine learning in biohydrogen production technology. Several machine learning algorithems have recently been implemented for modeling the nonlinear and complex relationships among operational and performance parameters in biohydrogen production as well as predicting the process performance and microbial population dynamics. Reinforced machine learning methods exhibited precise state prediction and retrieved the underlying kinetics effectively. Machine-learning based prediction was also improved by using microbial sequencing data as input parameters. Further research on machine learning could be instrumental in designing a process control tool to maintain reliable hydrogen production performance and identify connection between the process performance and the microbial population.

Comparative Analysis of a Family of Sliding Mode Observers under Real-Time Conditions for the Monitoring in the Bioethanol Production

Online monitoring of fermentation processes is a necessary task to determine concentrations of key biochemical compounds, diagnose faults in process operations, and implement feedback controllers. However, obtaining the signals of all-important variables in a real process is a task that may be difficult and expensive due to the lack of adequate sensors, or simply because some variables cannot be directly measured. From the above, a model-based approach such as state observers may be a viable alternative to solve the estimation problem. This work shows a comparative analysis of the real-time performance of a family of sliding-mode observers for reconstructing key variables in a batch bioreactor for fermentative ethanol production. These observers were selected for their robust performance under model uncertainties and finite-time estimation convergence. The selected sliding-mode observers were the first-order sliding mode observer, the proportional sliding mode observer, and the high-order sliding mode observer. For estimation purposes, a power law kinetic model for ethanol production by Saccharomyces cerevisiae was performed. A hybrid methodology allows the kinetic parameters to be adjusted, and an approach based on inference diagrams allows the observability of the model to be determined. The experimental results reported here show that the observers under analysis were robust to modeling errors and measurement noise. Moreover, the proportional sliding-mode observer was the algorithm that exhibited the best performance.

Advances in fermentative production, purification, characterization and applications of gellan gum

Multiple microbial exopolysaccharides have been reported in recent decade with their structural and functional features. Gellan gum (GG) is among these emerging biopolymers with versatile properties. Low production yield, high downstream cost, and abundant market demand have made GG a high cost material. Hence, an understanding on the various possibilities to develop cost-effective gellan gum bioprocess is desirable. This review focuses on details of upstream and downstream process of GG from an industrial perspective. It emphasizes on GG producing Sphingomonas spp., updates on biosynthesis, strain and media engineering, kinetic modeling, bioreactor design and scale-up considerations. Details of the downstream operations with possible modifications to make it cost-effective and environmentally sustainable have been discussed. The updated regulatory criteria for GG as a food ingredient and analytical tools required to validate the same have been briefly discussed. Derivatives of GG and their applications in various industrial segments have also been highlighted.

Reconstructing Kinetic Models for Dynamical Studies of Metabolism using Generative Adversarial Networks

Kinetic models of metabolism relate metabolic fluxes, metabolite concentrations and enzyme levels through mechanistic relations, rendering them essential for understanding, predicting and optimizing the behaviour of living organisms. However, due to the lack of kinetic data, traditional kinetic modelling often yields only a few or no kinetic models with desirable dynamical properties, making the analysis unreliable and computationally inefficient. We present REKINDLE (Reconstruction of Kinetic Models using Deep Learning), a deep-learning-based framework for efficiently generating kinetic models with dynamic properties matching the ones observed in cells. We showcase REKINDLE's capabilities to navigate through the physiological states of metabolism using small numbers of data with significantly lower computational requirements. The results show that data-driven neural networks assimilate implicit kinetic knowledge and structure of metabolic networks and generate kinetic models with tailored properties and statistical diversity. We anticipate that our framework will advance our understanding of metabolism and accelerate future research in biotechnology and health.

Modeling the Succinic Acid Bioprocess: A Review

Succinic acid has attracted much interest as a key platform chemical that can be obtained in high titers from biomass through sustainable fermentation processes, thus boosting the bioeconomy as a critical production strategy for the future. After several years of development of the production of succinic acid, many studies on lab or pilot scale production have been reported. The relevant experimental data reveal underlying physical and chemical dynamic phenomena. To take advantage of this vast, but disperse, kinetic information, a number of mathematical kinetic models of the unstructured non-segregated type have been proposed in the first place. These relatively simple models feature critical aspects of interest for the design, control, optimization and operation of this key bioprocess. This review includes a detailed description of the phenomena involved in the bioprocesses and how they reflect on the most important and recent models based on macroscopic and metabolic chemical kinetics, and in some cases even coupling mass transport.

Integrating a dynamic central metabolism model of cancer cells with a hybrid 3D multiscale model for vascular hepatocellular carcinoma growth

We develop here a novel modelling approach with the aim of closing the conceptual gap between tumour-level metabolic processes and the metabolic processes occurring in individual cancer cells. In particular, the metabolism in hepatocellular carcinoma derived cell lines (HEPG2 cells) has been well characterized but implementations of multiscale models integrating this known metabolism have not been previously reported. We therefore extend a previously published multiscale model of vascular tumour growth, and integrate it with an experimentally verified network of central metabolism in HEPG2 cells. This resultant combined model links spatially heterogeneous vascular tumour growth with known metabolic networks within tumour cells and accounts for blood flow, angiogenesis, vascular remodelling and nutrient/growth factor transport within a growing tumour, as well as the movement of, and interactions between normal and cancer cells. Model simulations report for the first time, predictions of spatially resolved time courses of core metabolites in HEPG2 cells. These simulations can be performed at a sufficient scale to incorporate clinically relevant features of different tumour systems using reasonable computational resources. Our results predict larger than expected temporal and spatial heterogeneity in the intracellular concentrations of glucose, oxygen, lactate pyruvate, f16bp and Acetyl-CoA. The integrated multiscale model developed here provides an ideal quantitative framework in which to study the relationship between dosage, timing, and scheduling of anti-neoplastic agents and the physiological effects of tumour metabolism at the cellular level. Such models, therefore, have the potential to inform treatment decisions when drug response is dependent on the metabolic state of individual cancer cells.

Impact of Influenza A Virus Infection on Growth and Metabolism of Suspension MDCK Cells Using a Dynamic Model

Cell cultured-based influenza virus production is a viable option for vaccine manufacturing. In order to achieve a high concentration of viable cells, is requirement to have not only optimal process conditions, but also an active metabolism capable of intracellular synthesis of viral components. Experimental metabolic data collected in such processes are complex and difficult to interpret, for which mathematical models are an appropriate way to simulate and analyze the complex and dynamic interaction between the virus and its host cell. A dynamic model with 35 states was developed in this study to describe growth, metabolism, and influenza A virus production in shake flask cultivations of suspension Madin-Darby Canine Kidney (MDCK) cells. It considers cell growth (concentration of viable cells, mean cell diameters, volume of viable cells), concentrations of key metabolites both at the intracellular and extracellular level and virus titers. Using one set of parameters, the model accurately simulates the dynamics of mock-infected cells and correctly predicts the overall dynamics of virus-infected cells for up to 60 h post infection (hpi). The model clearly suggests that most changes observed after infection are related to cessation of cell growth and the subsequent transition to apoptosis and cell death. However, predictions do not cover late phases of infection, particularly for the extracellular concentrations of glutamate and ammonium after about 12 hpi. Results obtained from additional in silico studies performed indicated that amino acid degradation by extracellular enzymes resulting from cell lysis during late infection stages may contribute to this observed discrepancy.

Mathematical Modeling for the Industrial 2-Mercaptobenzothiazole Batch Production Process

As an important chemical intermediate, 2-mercaptobenzothiazole (MBT) is widely used in various processes, especially in the rubber industry. However, there is no first-principles model that describes the synthetic process of MBT. This paper focuses on the formulation of a reliable mathematical model represented by a series of differential and algebraic equations for the industrial batch MBT production process. It is difficult to estimate all of the unknown parameters in the model because of the lack of sufficient industrial/experimental data. Thus, a reduced estimable parameter set is derived by performing estimability analysis on the original estimation problem. A multiple-starting-point strategy is then applied to numerically solve the non-convex parameter estimation problem with the weighted least-squares approach. Subsequently, a cross-validation technique is employed to evaluate the generalizability of the proposed model. Finally, it is confirmed that the proposed model produces encouraging prediction performance with regard to independent test data.

Kinetic and Stoichiometric Modeling-Based Analysis of Docosahexaenoic Acid (DHA) Production Potential by C. cohnii from Glycerol, Glucose and Ethanol

Docosahexaenoic acid (DHA) is one of the most important long-chain polyunsaturated fatty acids (LC-PUFAs), with numerous health benefits. Crypthecodinium cohnii, a marine heterotrophic dinoflagellate, is successfully used for the industrial production of DHA because it can accumulate DHA at high concentrations within the cells. Glycerol is an interesting renewable substrate for DHA production since it is a by-product of biodiesel production and other industries, and is globally generated in large quantities. The DHA production potential from glycerol, ethanol and glucose is compared by combining fermentation experiments with the pathway-scale kinetic modeling and constraint-based stoichiometric modeling of C. cohnii metabolism. Glycerol has the slowest biomass growth rate among the tested substrates. This is partially compensated by the highest PUFAs fraction, where DHA is dominant. Mathematical modeling reveals that glycerol has the best experimentally observed carbon transformation rate into biomass, reaching the closest values to the theoretical upper limit. In addition to our observations, the published experimental evidence indicates that crude glycerol is readily consumed by C. cohnii, making glycerol an attractive substrate for DHA production.

A protocol for dynamic model calibration

Ordinary differential equation models are nowadays widely used for the mechanistic description of biological processes and their temporal evolution. These models typically have many unknown and nonmeasurable parameters, which have to be determined by fitting the model to experimental data. In order to perform this task, known as parameter estimation or model calibration, the modeller faces challenges such as poor parameter identifiability, lack of sufficiently informative experimental data and the existence of local minima in the objective function landscape. These issues tend to worsen with larger model sizes, increasing the computational complexity and the number of unknown parameters. An incorrectly calibrated model is problematic because it may result in inaccurate predictions and misleading conclusions. For nonexpert users, there are a large number of potential pitfalls. Here, we provide a protocol that guides the user through all the steps involved in the calibration of dynamic models. We illustrate the methodology with two models and provide all the code required to reproduce the results and perform the same analysis on new models. Our protocol provides practitioners and researchers in biological modelling with a one-stop guide that is at the same time compact and sufficiently comprehensive to cover all aspects of the problem.

Kinetic Modeling of Saccharomyces cerevisiae Central Carbon Metabolism: Achievements, Limitations, and Opportunities

Central carbon metabolism comprises the metabolic pathways in the cell that process nutrients into energy, building blocks and byproducts. To unravel the regulation of this network upon glucose perturbation, several metabolic models have been developed for the microorganism Saccharomyces cerevisiae. These dynamic representations have focused on glycolysis and answered multiple research questions, but no commonly applicable model has been presented. This review systematically evaluates the literature to describe the current advances, limitations, and opportunities. Different kinetic models have unraveled key kinetic glycolytic mechanisms. Nevertheless, some uncertainties regarding model topology and parameter values still limit the application to specific cases. Progressive improvements in experimental measurement technologies as well as advances in computational tools create new opportunities to further extend the model scale. Notably, models need to be made more complex to consider the multiple layers of glycolytic regulation and external physiological variables regulating the bioprocess, opening new possibilities for extrapolation and validation. Finally, the onset of new data representative of individual cells will cause these models to evolve from depicting an average cell in an industrial fermenter, to characterizing the heterogeneity of the population, opening new and unseen possibilities for industrial fermentation improvement.

Microbial technologies for biotherapeutics production: Key tools for advanced biopharmaceutical process development and control

Current trends in the biopharmaceutical market such as the diversification of therapies as well as the increasing time-to-market pressure will trigger the rethinking of bioprocess development and production approaches. Thereby, the importance of development time and manufacturing costs will increase, especially for microbial production. In the present review, we investigate three technological approaches which, to our opinion, will play a key role in the future of biopharmaceutical production. The first cornerstone of process development is the generation and effective utilization of platform knowledge. Building processes on well understood microbial and technological platforms allows to accelerate early-stage bioprocess development and to better condense this knowledge into multi-purpose technologies and applicable mathematical models. Second, the application of verified scale down systems and in silico models for process design and characterization will reduce the required number of large scale batches before dossier submission. Third, the broader availability of mathematical process models and the improvement of process analytical technologies will increase the applicability and acceptance of advanced control and process automation in the manufacturing scale. This will reduce process failure rates and subsequently cost of goods. Along these three aspects we give an overview of recently developed key tools and their potential integration into bioprocess development strategies.

A genome-scale nutrient minimization forecast algorithm for controlling essential amino acid levels in CHO cell cultures

Mammalian cell culture processes rely heavily on empirical knowledge in which process control remains a challenge due to the limited characterization/understanding of cell metabolism and inability to predict the cell behaviors. This study facilitates control of Chinese hamster ovary (CHO) processes through a forecast-based feeding approach that predicts multiple essential amino acids levels in the culture from easily acquired viable cell density data. Multiple cell growth behavior forecast extrapolation approaches are considered with logistic curve fitting found to be the most effective. Next, the nutrient-minimized CHO genome-scale model is combined with the growth forecast model to generate essential amino acid forecast profiles of multiple CHO batch cultures. Comparison of the forecast with the measurements suggests that this algorithm can accurately predict the concentration of most essential amino acids from cell density measurement with error mitigated by incorporating off-line amino acids concentration measurements. Finally, the forecast algorithm is applied to CHO fed-batch cultures to support amino acid feeding control to control the concentration of essential amino acids below 1-2 mM for lysine, leucine, and valine as a model over a 9-day fed batch culture while maintaining comparable growth behavior to an empirical-based culture. In turn, glycine production was elevated, alanine reduced and lactate production slightly lower in control cultures due to metabolic shifts in branched-chain amino acid degradation. With the advantage of requiring minimal measurement inputs while providing valuable and in-advance information of the system based on growth measurements, this genome model-based amino acid forecast algorithm represent a powerful and cost-effective tool to facilitate enhanced control over CHO and other mammalian cell-based bioprocesses.

Valorisation of fungal hydrolysates of exhausted sugar beet pulp for lactic acid production

BACKGROUND

Exhausted sugar beet pulp pellets (ESBPP) were used as raw material for lactic acid (LA) fermentation. The enzymatic hydrolysis of ESBPP was performed with the solid obtained after the fungal solid-state fermentation of ESBPP as a source of hydrolytic enzymes. Subsequently, a medium rich in glucose and arabinose was obtained, which was used to produce LA by fermentation. For LA production, two Lactobacillus strains were assayed and the effects of the supplementation of the hydrolysate with a nitrogen source and the mode of pH regulation of the fermentation were investigated. Moreover, a kinetic model for LA fermentation by Lactobacillus plantarum of ESBPP hydrolysates was developed.

RESULTS

L. plantarum produced a LA concentration 34% higher than that produced by L. casei. The highest LA concentration (30 g L^-1 ) was obtained with L. plantarum when the hydrolysate was supplemented with 5 g L^-1 yeast extract and the pH was controlled with CaCO₃ . The concentration of acetic acid differed depending on the concentration of CaCO₃ added, producing its maximum value with 27 g L^-1 CaCO₃ . The proposed kinetic model was able to predict the evolution of substrates and products depending on the variation of the pH in the hydrolysate, according to the amount of CaCO₃ added.

CONCLUSIONS

ESBPP can be revalorised to produce LA. A pure LA stream or a mixture of LA and acetic acid, depending on the pH control method of the fermentation, can be produced. Thus, this control is of great interest depending on the destination of the effluent. © 2020 Society of Chemical Industry.

Synergisms of machine learning and constraint-based modeling of metabolism for analysis and optimization of fermentation parameters

Recent noteworthy advances in the development of high-performing microbial and mammalian strains have enabled the sustainable production of bio-economically valuable substances such as bio-compounds, biofuels, and biopharmaceuticals. However, to obtain an industrially viable mass-production scheme, much time and effort are required. The robust and rational design of fermentation processes requires analysis and optimization of different extracellular conditions and medium components, which have a massive effect on growth and productivity. In this regard, knowledge- and data-driven modeling methods have received much attention. Constraint-based modeling (CBM) is a knowledge-driven mathematical approach that has been widely used in fermentation analysis and optimization due to its capabilities of predicting the cellular phenotype from genotype through high-throughput means. On the other hand, machine learning (ML) is a data-driven statistical method that identifies the data patterns within sophisticated biological systems and processes, where there is inadequate knowledge to represent underlying mechanisms. Furthermore, ML models are becoming a viable complement to constraint-based models in a reciprocal manner when one is used as a pre-step of another. As a result, more predictable model is produced. This review highlights the applications of CBM and ML independently and the combination of these two approaches for analyzing and optimizing fermentation parameters. This article is protected by copyright. All rights reserved.

Modelling Cell Metabolism: A Review on Constraint-Based Steady-State and Kinetic Approaches

Studying cell metabolism serves a plethora of objectives such as the enhancement of bioprocess performance, and advancement in the understanding of cell biology, of drug target discovery, and in metabolic therapy. Remarkable successes in these fields emerged from heuristics approaches, for instance, with the introduction of effective strategies for genetic modifications, drug developments and optimization of bioprocess management. However, heuristics approaches have showed significant shortcomings, such as to describe regulation of metabolic pathways and to extrapolate experimental conditions. In the specific case of bioprocess management, such shortcomings limit their capacity to increase product quality, while maintaining desirable productivity and reproducibility levels. For instance, since heuristics approaches are not capable of prediction of the cellular functions under varying experimental conditions, they may lead to sub-optimal processes. Also, such approaches used for bioprocess control often fail in regulating a process under unexpected variations of external conditions. Therefore, methodologies inspired by the systematic mathematical formulation of cell metabolism have been used to address such drawbacks and achieve robust reproducible results. Mathematical modelling approaches are effective for both the characterization of the cell physiology, and the estimation of metabolic pathways utilization, thus allowing to characterize a cell population metabolic behavior. In this article, we present a review on methodology used and promising mathematical modelling approaches, focusing primarily to investigate metabolic events and regulation. Proceeding from a topological representation of the metabolic networks, we first present the metabolic modelling approaches that investigate cell metabolism at steady state, complying to the constraints imposed by mass conservation law and thermodynamics of reactions reversibility. Constraint-based models (CBMs) are reviewed highlighting the set of assumed optimality functions for reaction pathways. We explore models simulating cell growth dynamics, by expanding flux balance models developed at steady state. Then, discussing a change of metabolic modelling paradigm, we describe dynamic kinetic models that are based on the mathematical representation of the mechanistic description of nonlinear enzyme activities. In such approaches metabolic pathway regulations are considered explicitly as a function of the activity of other components of metabolic networks and possibly far from the metabolic steady state. We have also assessed the significance of metabolic model parameterization in kinetic models, summarizing a standard parameter estimation procedure frequently employed in kinetic metabolic modelling literature. Finally, some optimization practices used for the parameter estimation are reviewed.

Generic Workflow for the Setup of Mechanistic Process Models

Process models, consisting of transferable and applicable knowledge, can be used for various tasks such as process development and optimization, and for predicting and controlling critical process variables. In this regard, mechanistic process models, describing the mechanism of a system with a distinct model structure and characteristic parameters, are very promising.The development of a reliable and applicable model is usually the critical step, before model simulation and application show beneficial effects. In this chapter, a workflow for the generation of mechanistic process models is presented and applied on a typical cell culture process. The workflow includes the definition of critical reactions and the identification of their kinetics. By an iterative approach different reactions and kinetics are tested and model quality is assessed, leading to a final, target-oriented model.

Recent advances in systems and synthetic biology approaches for developing novel cell-factories in non-conventional yeasts

Microbial bioproduction of chemicals, proteins, and primary metabolites from cheap carbon sources is currently an advancing area in industrial research. The model yeast, Saccharomyces cerevisiae, is a well-established biorefinery host that has been used extensively for commercial manufacturing of bioethanol from myriad carbon sources. However, its Crabtree-positive nature often limits the use of this organism for the biosynthesis of commercial molecules that do not belong in the fermentative pathway. To avoid extensive strain engineering of S. cerevisiae for the production of metabolites other than ethanol, non-conventional yeasts can be selected as hosts based on their natural capacity to produce desired commodity chemicals. Non-conventional yeasts like Kluyveromyces marxianus, K. lactis, Yarrowia lipolytica, Pichia pastoris, Scheffersomyces stipitis, Hansenula polymorpha, and Rhodotorula toruloides have been considered as potential industrial eukaryotic hosts owing to their desirable phenotypes such as thermotolerance, assimilation of a wide range of carbon sources, as well as ability to secrete high titers of protein and lipid. However, the advanced metabolic engineering efforts in these organisms are still lacking due to the limited availability of systems and synthetic biology methods like in silico models, well-characterised genetic parts, and optimized genome engineering tools. This review provides an insight into the recent advances and challenges of systems and synthetic biology as well as metabolic engineering endeavours towards the commercial usage of non-conventional yeasts. Particularly, the approaches in emerging non-conventional yeasts for the production of enzymes, therapeutic proteins, lipids, and metabolites for commercial applications are extensively discussed here. Various attempts to address current limitations in designing novel cell factories have been highlighted that include the advances in the fields of genome-scale metabolic model reconstruction, flux balance analysis, 'omics'-data integration into models, genome-editing toolkit development, and rewiring of cellular metabolisms for desired chemical production. Additionally, the understanding of metabolic networks using ¹³C-labelling experiments as well as the utilization of metabolomics in deciphering intracellular fluxes and reactions have also been discussed here. Application of cutting-edge nuclease-based genome editing platforms like CRISPR/Cas9, and its optimization towards efficient strain engineering in non-conventional yeasts have also been described. Additionally, the impact of the advances in promising non-conventional yeasts for efficient commercial molecule synthesis has been meticulously reviewed. In the future, a cohesive approach involving systems and synthetic biology will help in widening the horizon of the use of unexplored non-conventional yeast species towards industrial biotechnology.

Model Parameterization with Quantitative Proteomics: Case Study with Trehalose Metabolism in Saccharomyces cerevisiae

When Saccharomyces cerevisiae undergoes heat stress it stimulates several changes that are necessary for its survival, notably in carbon metabolism. Notable changes include increase in trehalose production and glycolytic flux. The increase in glycolytic flux has been postulated to be due to the regulatory effects in upper glycolysis, but this has not been confirmed. Additionally, trehalose is a useful industrial compound for its protective properties. A model of trehalose metabolism in S. cerevisiae was constructed using Convenient Modeller, a software that uses a combination of convenience kinetics and a genetic algorithm. The model was parameterized with quantitative omics under standard conditions and validated using data collected under heat stress conditions. The completed model was used to show that feedforward activation of pyruvate kinase by fructose 1,6-bisphosphate during heat stress contributes to the increase in metabolic flux. We were also able to demonstrate in silico that overexpression of enzymes involved in production and degradation of trehalose can lead to higher trehalose yield in the cell. By integrating quantitative proteomics with metabolic modelling, we were able to confirm that the flux increase in trehalose metabolic pathways during heat stress is due to regulatory effects and not purely changes in enzyme expression. The overexpression of enzymes involved in trehalose metabolism is a potential approach to be exploited for trehalose production without need for increasing temperature.

Safety and efficacy of a novel home-use device for light-potentiated (LED) skin treatment

Skin structure and function results from a dynamic interplay between dermal and epidermal cell types. Optimizing skin health through an effective and long-lasting skin care regime therefore requires a global approach, encompassing various mechanisms to stimulate this interplay beyond the action scope of a classical topical solution. This study evaluates the impact of a novel home-use device combining a topical serum, light-emitting diodes and massage on the clinical signs of extrinsic skin aging. The innovative principle relies on potentiating the effect of active ingredients contained in the topical serum with visible and near infra-red photons to prevent extracellular matrix degradation and promote its reconstruction. After in vitro and ex vivo investigations, a clinical study assessed the safety and efficacy of a daily treatment with the home-use device for 28 days. A significant increases in skin density and radiance while reducing the wrinkles was obtained with no side effects.

Predicting By-Product Gradients of Baker’s Yeast Production at Industrial Scale: A Practical Simulation Approach

Scaling up bioprocesses is one of the most crucial steps in the commercialization of bioproducts. While it is known that concentration and shear rate gradients occur at larger scales, it is often too risky, if feasible at all, to conduct validation experiments at such scales. Using computational fluid dynamics equipped with mechanistic biochemical engineering knowledge of the process, it is possible to simulate such gradients. In this work, concentration profiles for the by-products of baker’s yeast production are investigated. By applying a mechanistic black-box model, concentration heterogeneities for oxygen, glucose, ethanol, and carbon dioxide are evaluated. The results suggest that, although at low concentrations, ethanol is consumed in more than 90% of the tank volume, which prevents cell starvation, even when glucose is virtually depleted. Moreover, long exposure to high dissolved carbon dioxide levels is predicted. Two biomass concentrations, i.e., 10 and 25 g/L, are considered where, in the former, ethanol production is solely because of overflow metabolism while, in the latter, 10% of the ethanol formation is due to dissolved oxygen limitation. This method facilitates the prediction of the living conditions of the microorganism and its utilization to address the limitations via change of strain or bioreactor design or operation conditions. The outcome can also be of value to design a representative scale-down reactor to facilitate strain studies.

Using optimal control to understand complex metabolic pathways

BACKGROUND

Optimality principles have been used to explain the structure and behavior of living matter at different levels of organization, from basic phenomena at the molecular level, up to complex dynamics in whole populations. Most of these studies have assumed a single-criteria approach. Such optimality principles have been justified from an evolutionary perspective. In the context of the cell, previous studies have shown how dynamics of gene expression in small metabolic models can be explained assuming that cells have developed optimal adaptation strategies. Most of these works have considered rather simplified representations, such as small linear pathways, or reduced networks with a single branching point, and a single objective for the optimality criteria.

RESULTS

Here we consider the extension of this approach to more realistic scenarios, i.e. biochemical pathways of arbitrary size and structure. We first show that exploiting optimality principles for these networks poses great challenges due to the complexity of the associated optimal control problems. Second, in order to surmount such challenges, we present a computational framework which has been designed with scalability and efficiency in mind, including mechanisms to avoid the most common pitfalls. Third, we illustrate its performance with several case studies considering the central carbon metabolism of S. cerevisiae and B. subtilis. In particular, we consider metabolic dynamics during nutrient shift experiments.

CONCLUSIONS

We show how multi-objective optimal control can be used to predict temporal profiles of enzyme activation and metabolite concentrations in complex metabolic pathways. Further, we also show how to consider general cost/benefit trade-offs. In this study we have considered metabolic pathways, but this computational framework can also be applied to analyze the dynamics of other complex pathways, such as signal transduction or gene regulatory networks.

A model-driven approach towards rational microbial bioprocess optimization

Due to sustainability concerns, bio-based production capitalizing on microbes as cell factories is in demand to synthesize valuable products. Nevertheless, the nonhomogenous variations of the extracellular environment in bioprocesses often challenge the biomass growth and the bioproduction yield. To enable a more rational bioprocess optimization, we have established a model-driven approach that systematically integrates experiments with modeling, executed from flask to bioreactor scale, and using ferulic acid to vanillin bioconversion as a case study. The impacts of mass transfer and aeration on the biomass growth and bioproduction performances were examined using minimal small-scale experiments. An integrated model coupling the cell factory kinetics with the three-dimensional computational hydrodynamics of bioreactor was developed to better capture the spatiotemporal distributions of bioproduction. Full-factorial predictions were then performed to identify the desired operating conditions. A bioconversion yield of 94% was achieved, which is one of the highest for recombinant Escherichia coli using ferulic acid as the precursor.