Search

A powerful approach for the optimization of industrial bioprocesses is to perform detailed simulations integrating large-scale computational fluid dynamics (CFD) and cellular reaction dynamics (CRD). However, complex metabolic kinetic models containing a large number of equations pose formidable challenges in CFD-CRD coupling and computation time afterward. This necessitates to formulate a relatively simple but yet representative model structure. Such a kinetic model should be able to reproduce metabolic responses for short-term (mixing time scale of tens of seconds) and long-term (fed-batch cultivation of hours/days) dynamics in industrial bioprocesses. In this paper, we used Penicillium chrysogenum as a model system and developed a metabolically structured kinetic model for growth and production. By lumping the most important intracellular metabolites in 5 pools and 4 intracellular enzyme pools, linked by 10 reactions, we succeeded in maintaining the model structure relatively simple, while providing informative insight into the state of the organism. The performance of this 9-pool model was validated with a periodic glucose feast-famine cycle experiment at the minute time scale. Comparison of this model and a reported black box model for this strain shows the necessity of employing a structured model under feast-famine conditions. This proposed model provides deeper insight into the in vivo kinetics and, most importantly, can be straightforwardly integrated into a computational fluid dynamic framework for simulating complete fermentation performance and cell population dynamics in large scale and small scale fermentors. Biotechnol. Bioeng. 2017;114: 1733-1743. © 2017 Wiley Periodicals, Inc.

A powerful approach for the optimization of industrial bioprocesses is to perform detailed simulations integrating large-scale computational fluid dynamics (CFD) and cellular reaction dynamics (CRD). However, complex metabolic kinetic models containing a large number of equations pose formidable challenges in CFD-CRD coupling and computation time afterward. This necessitates to formulate a relatively simple but yet representative model structure. Such a kinetic model should be able to reproduce metabolic responses for short-term (mixing time scale of tens of seconds) and long-term (fed-batch cultivation of hours/days) dynamics in industrial bioprocesses. In this paper, we used Penicillium chrysogenum as a model system and developed a metabolically structured kinetic model for growth and production. By lumping the most important intracellular metabolites in 5 pools and 4 intracellular enzyme pools, linked by 10 reactions, we succeeded in maintaining the model structure relatively simple, while providing informative insight into the state of the organism. The performance of this 9-pool model was validated with a periodic glucose feast-famine cycle experiment at the minute time scale. Comparison of this model and a reported black box model for this strain shows the necessity of employing a structured model under feast-famine conditions. This proposed model provides deeper insight into the in vivo kinetics and, most importantly, can be straightforwardly integrated into a computational fluid dynamic framework for simulating complete fermentation performance and cell population dynamics in large scale and small scale fermentors. Biotechnol. Bioeng. 2017;114: 1733-1743. © 2017 Wiley Periodicals, Inc.

This contribution addresses the identification of metabolic fluxes and metabolite concentrations in mammalian cells from transient (13)C-labeling experiments. Whilst part I describes experimental set-up and acquisition of required metabolite and (13)C-labeling data, part II focuses on setting up network models and the estimation of intracellular fluxes. Metabolic fluxes were determined in glycolysis, pentose-phosphate pathway (PPP), and citric acid cycle (TCA) in a hepatoma cell line grown in aerobic batch cultures. In glycolytic and PPP metabolite pools isotopic stationarity was observed within 30 min, whereas in the TCA cycle the labeling redistribution did not reach isotopic steady state even within 180 min. In silico labeling dynamics were in accordance with in vivo (13)C-labeling data. Split ratio between glycolysis and PPP was 57%:43%; intracellular glucose concentration was estimated at 101.6 nmol per 10(6) cells. In contrast to isotopic stationary (13)C-flux analysis, transient (13)C-flux analysis can also be applied to industrially relevant mammalian cell fed-batch and batch cultures.

The novel concept of isotopic dynamic 13C metabolic flux analysis (ID-13C MFA) enables integrated analysis of isotopomer data from isotopic transient and/or isotopic stationary phase of a 13C labeling experiment, short-time experiments, and an extended range of applications of 13C MFA. In the presented work, an experimental and computational framework consisting of short-time 13C labeling, an integrated rapid sampling procedure, a LC-MS analytical method, numerical integration of the system of isotopomer differential equations, and estimation of metabolic fluxes was developed and applied to determine intracellular fluxes in glycolysis, pentose phosphate pathway (PPP), and citric acid cycle (TCA) in Escherichia coli grown in aerobic, glucose-limited chemostat culture at a dilution rate of D = 0.10 h−1. Intracellular steady state concentrations were quantified for 12 metabolic intermediates. A total of 90 LC-MS mass isotopomers were quantified at sampling times t = 0, 91, 226, 346, 589 s and at isotopic stationary conditions. Isotopic stationarity was reached within 10 min in glycolytic and PPP metabolites. Consistent flux solutions were obtained by ID-13C MFA using isotopic dynamic and isotopic stationary 13C labeling data and by isotopic stationary 13C MFA (IS-13C MFA) using solely isotopic stationary data. It is demonstrated that integration of dynamic 13C labeling data increases the sensitivity of flux estimation, particularly at the glucose-6-phosphate branch point. The identified split ratio between glycolysis and PPP was 55%:44%. These results were confirmed by IS-13C MFA additionally using labeling data in proteinogenic amino acids (GC-MS) obtained after 5 h from sampled biomass. Biotechnol. Bioeng. 2008;99: 1170–1185. © 2007 Wiley Periodicals, Inc.

Application of metabolic engineering principles to the rational design of microbial production processes crucially depends on the ability to describe quantitatively the systemic behavior of the central carbon metabolism to redirect carbon fluxes to the product-forming pathways. Despite the importance for several production processes, development of an essential dynamic model for central carbon metabolism of Escherichia coli has been severely hampered by the current lack of kinetic information on the dynamics of the metabolic reactions. Here we present the design and experimental validation of such a dynamic model, which, for the first time, links the sugar transport system (i.e., phosphotransferase system [PTS]) with the reactions of glycolysis and the pentose-phosphate pathway. Experimental observations of intracellular concentrations of metabolites and cometabolites at transient conditions are used to validate the structure of the model and to estimate the kinetic parameters. Further analysis of the detailed characteristics of the system offers the possibility of studying important questions regarding the stability and control of metabolic fluxes. © 2002 Wiley Periodicals, Inc. Biotechnol Bioeng 79: 53–73; 2002.

A high-cell-density fed-batch fermentation for the production of heterologous proteins in Escherichia coli was developed using the positively regulated Escherichia coli rhaBAD promoter. The expression system was improved by reducing of the amount of expensive L-rhamnose necessary for induction of the rhamnose promoter and by increasing the vector stability. Consumption of the inducer L-rhamnose was inhibited by inactivation of L-rhamnulose kinase encoding gene rhaB of Escherichia coli W3110, responsible for the first irreversible step in rhamnose catabolism. Plasmid instability caused by multimerization of the expression vector in the recombination-proficient W3110 was prevented by insertion of the multimer resolution site cer from the ColE1 plasmid into the vector. Fermentation experiments with the optimized system resulted in the production of 100 g x L(-1) cell dry weight and 3.8 g x L(-1) of recombinant L-N-carbamoylase, an enzyme, which is needed for the production of enantiomeric pure amino acids in a two-step reaction from hydantoins.

In the present study, the glucose transport into the yeast Saccharomyces cerevisiae has been investigated. The approach suggested is based on a rapid sampling technique for studying the dynamic response of the yeast to rapid changes in extracellular glucose concentrations. For this purpose a concentrated glucose solution has been injected into a continuous culture at steady state growth conditions resulting in a shift of the extracellular glucose level. Samples have been taken every 5 s for determination of extracellular glucose and intracellular glucose-6-phosphate concentrations. Attempts to fit the experimental observations with simulations from existing models failed. The mechanism then proposed is based on a facilitated diffusion of glucose superimposed by an inhibition of glucose-6-phosphate. The use of the so-called in vivo approach suggested in this article appears to be proper, because the investigations can be performed at defined physiological states of the microbial cultures. Furthermore, the experimental observations are not being corrupted by the preparation of the samples for the transport studies as it happens during radioactive measurements. (c) 1996 John Wiley & Sons, Inc.

A mathematical model of glycolysis in Saccharomyces cerevisiae is presented. The model is based on rate equations for the individual reactions and aims to predict changes in the levels of intra- and extracellular metabolites after a glucose pulse, as described in part I of this study. Kinetic analysis focuses on a time scale of seconds, thereby neglecting biosynthesis of new enzymes. The model structure and experimental observations are related to the aerobic growth of the yeast. The model is based on material balance equations of the key metabolites in the extracellular environment, the cytoplasm and the mitochondria, and includes mechanistically based, experimentally matched rate equations for the individual enzymes. The model includes removal of metabolites from glycolysis and TCC for biosynthesis, and also compartmentation and translocation of adenine nucleotides. The model was verified by in vivo diagnosis of intracellular enzymes, which includes the decomposition of the network of reactions to reduce the number of parameters to be estimated simultaneously. Additionally, sensitivity analysis guarantees that only those parameters are estimated that contribute to systems trajectory with reasonable sensitivity. The model predictions and experimental observations agree reasonably well for most of the metabolites, except for pyruvate and adenine nucleotides. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 55: 592-608, 1997.

The goal of this work was to obtain rapid sampling technique to measure transient metabolites in vivo. First, a pulse of glucose was added to a culture of the yeast Saccharomyces cerevisiae growing aerobically under glucose limitation. Next, samples were removed at 2 to 5 s intervals and quenched using methods that depend on the metabolite measured. Extracellular glucose, excreted products, as well as glycolytic intermediates (G6P, F6P, FBP, GAP, 3-PG, PEP, Pyr) and cometabolites (ATP, ADP, AMP, NAD(+), NADH) were measured using enzymatic or HPLC methods. Significant differences between the adenine nucleotide concentrations in the cytoplasm and mitochondria indicated the importance of compartmentation for the regulation of the glycolysis. Changes in the intra- and extracellular levels of metabolites confirmed that glycolysis is regulated on a time scale of seconds. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 55: 305-316, 1997.

An experimental set-up for acquiring metabolite and transient 13C-labeling data in mammalian cells is presented. An efficient sampling procedure was established for hepatic cells cultured in six-well plates as a monolayer attached to collagen, which allowed simultaneous quenching of metabolism and extraction of the intracellular intermediates of interest. Extracellular concentrations of glucose, amino acids, lactate, pyruvate, and urea were determined by GC–MS procedures and were used for estimation of metabolic uptake and excretion rates. Sensitive LC–MS and GC–MS methods were used to quantify the intracellular intermediates of tricarboxylic acid cycle, glycolysis, and pentose phosphate pathway and for the determination of isotopomer fractions of the respective metabolites. Mass isotopomer fractions were determined in a transient 13C-labeling experiment using 13C-labeled glucose as substrate. The absolute amounts of intracellular metabolites were obtained from a non-labeled experiment carried out in exactly the same way as the 13C-labeling experiment, except that the media contained naturally labeled glucose only. Estimation of intracellular metabolic fluxes from the presented data is addressed in part II of this contribution. Biotechnol. Bioeng. 2008;100: 344–354. © 2007 Wiley Periodicals, Inc.

Knowledge of concentrations of intracellular metabolites is important for quantitative analysis of metabolic networks. As far as the very fast response of intracellular metabolites in the millisecond range is concerned, the frequently used pulse technique shows an inherent limitation. The time span between the disturbance and the first sample is constrained by the time necessary to obtain a homogeneous distribution of the pertubation within the bioreactor. For determination of rapid changes, a novel sampling technique based on the stopped-flow method has been developed. A continuous stream of biosuspension leaving the bioreactor is being mixed with a glucose solution in a turbulent mixing chamber. Through computer-aided activation of sequentially positioned three-way valves, different residence times and thus reaction times can be verified. The application of this new sampling method is illustrated with examples including measurements of adenine nucleotides and glucose-6-phosphate in Saccharomyces cerevisiae as well as measurements related to the PTS system in Escherichia coli.

This study provides a mathematical model of T7 RNA polymerase (T7 RNAP) kinetics under in vitro conditions targeted at application of this model to simulation of dynamic transcription performance. A functional dependence of transcript synthesis rate is derived based on: (a) essential reactant concentrations, including T7 RNAP and its promoter, substrate nucleotides, and the inhibitory byproduct inorganic pyrophosphate; (b) a distinction among vector characteristics such as recognition sequences regulating transcription initiation and termination, respectively; and (c) specific properties of the nucleotide sequence including both transcript length and nucleotide composition. Inactivation kinetics showed a half-life of T7 RNAP activity of 50 min under the conditions applied in vitro using the isolated enzyme. Model parameters and their precision are estimated using dynamic simulation and nonlinear regression analysis. The particular novelty of this model is its capability to incorporate linear genomic sequence information for simulation of nonlinear in vitro transcription kinetics.

The nonaqueous lipase-catalyzed kinetic resolution of racemic 2-methyl-1-pentanol in a continuously operated fixed bed reactor was theoretically and experimentally studied. A 2-dimensional mathematical model was developed to predict the performance of the heterogeneous enantioselective biotransformation and to optimize the productivity and effective stereoselectivity of the process by taking pore diffusion, solid(SINGLEBOND)liquid mass transfer, convection, and axial dispersion into account. Experimental investigations were conducted with lipase from Pseudomonas sp. immobilized by ionic binding in the pores of the anion exchange resin Duolite A 568. As determined in prior initial rate experiments, the specific activities of the immobilized lipase were maximal at a water activity of approximately a(w) = 0.67 and revealed a significant dependence on the amount of enzyme bound to the carrier material, with enantiomeric ratios slightly increasing with increasing water activities. Continuous resolution processes were carried out at a wide range of enzyme loadings. By controlled immobilization according to the theoretically evaluated optimal enzyme loading, the continuous racemic resolution could be optimized to obtain (R)-2-methyl-1-pentanol at high productivity and enantiopurity (ee > 95%). The steady-state characteristics of the system could be generally predicted by the model, despite the necessity to reevaluate the kinetic properties of the immobilized lipase to account for the complex non-aqueous microkinetics in a heterogeneous environment. Further model extension introducing competitive inhibition of free water in solutions as well as water diffusion and adsorption to the biocatalyst were useful in providing a more accurate description of the experimental results. (c) 1996 John Wiley & Sons, Inc.