Comparative performance of different scale-down simulators of substrate gradients in Penicillium chrysogenum cultures: the need of a biological systems response analysis

In a 54 m³ large‐scale penicillin fermentor, the cells experience substrate gradient cycles at the timescales of global mixing time about 20–40 s. Here, we used an intermittent feeding regime (IFR) and a two‐compartment reactor (TCR) to mimic these substrate gradients at laboratory‐scale continuous cultures. The IFR was applied to simulate substrate dynamics experienced by the cells at full scale at timescales of tens of seconds to minutes (30 s, 3 min and 6 min), while the TCR was designed to simulate substrate gradients at an applied mean residence time (τc) of 6 min. A biological systems analysis of the response of an industrial high‐yielding P. chrysogenum strain has been performed in these continuous cultures. Compared to an undisturbed continuous feeding regime in a single reactor, the penicillin productivity (q_PenG) was reduced in all scale‐down simulators. The dynamic metabolomics data indicated that in the IFRs, the cells accumulated high levels of the central metabolites during the feast phase to actively cope with external substrate deprivation during the famine phase. In contrast, in the TCR system, the storage pool (e.g. mannitol and arabitol) constituted a large contribution of carbon supply in the non‐feed compartment. Further, transcript analysis revealed that all scale‐down simulators gave different expression levels of the glucose/hexose transporter genes and the penicillin gene clusters. The results showed that q_PenG did not correlate well with exposure to the substrate regimes (excess, limitation and starvation), but there was a clear inverse relation between q_PenG and the intracellular glucose level.

Stoichiometry and kinetics of single and mixed substrate uptake in Aspergillus niger

In its natural environment, the filamentous fungus Aspergillus niger grows on decaying fruits and plant material, thereby enzymatically degrading the lignocellulosic constituents (lignin, cellulose, hemicellulose, and pectin) into a mixture of mono- and oligosaccharides. To investigate the kinetics and stoichiometry of growth of this fungus on lignocellulosic sugars, we carried out batch cultivations on six representative monosaccharides (glucose, xylose, mannose, rhamnose, arabinose, and galacturonic acid) and a mixture of these. Growth on these substrates was characterized in terms of biomass yields, oxygen/biomass ratios, and specific conversion rates. Interestingly, in combination, some of the carbon sources were consumed simultaneously and some sequentially. With a previously developed protocol, a sequential chemostat cultivation experiment was performed on a feed mixture of the six substrates. We found that the uptake of glucose, xylose, and mannose could be described with a Michaelis–Menten-type kinetics; however, these carbon sources seem to be competing for the same transport systems, while the uptake of arabinose, galacturonic acid, and rhamnose appeared to be repressed by the presence of other substrates.

Euler-Lagrange computational fluid dynamics for (bio)reactor scale down: An analysis of organism lifelines

The trajectories, referred to as lifelines, of individual microorganisms in an industrial scale fermentor under substrate limiting conditions were studied using an Euler‐Lagrange computational fluid dynamics approach. The metabolic response to substrate concentration variations along these lifelines provides deep insight in the dynamic environment inside a large‐scale fermentor, from the point of view of the microorganisms themselves. We present a novel methodology to evaluate this metabolic response, based on transitions between metabolic “regimes” that can provide a comprehensive statistical insight in the environmental fluctuations experienced by microorganisms inside an industrial bioreactor. These statistics provide the groundwork for the design of representative scale‐down simulators, mimicking substrate variations experimentally. To focus on the methodology we use an industrial fermentation of Penicillium chrysogenum in a simplified representation, dealing with only glucose gradients, single‐phase hydrodynamics, and assuming no limitation in oxygen supply, but reasonably capturing the relevant timescales. Nevertheless, the methodology provides useful insight in the relation between flow and component fluctuation timescales that are expected to hold in physically more thorough simulations. Microorganisms experience substrate fluctuations at timescales of seconds, in the order of magnitude of the global circulation time. Such rapid fluctuations should be replicated in truly industrially representative scale‐down simulators.

Transport and metabolism of fumaric acid in Saccharomyces cerevisiae in aerobic glucose-limited chemostat culture

Currently, research is being focused on the industrial-scale production of fumaric acid and other relevant organic acids from renewable feedstocks via fermentation, preferably at low pH for better product recovery. However, at low pH a large fraction of the extracellular acid is present in the undissociated form, which is lipophilic and can diffuse into the cell. There have been no studies done on the impact of high extracellular concentrations of fumaric acid under aerobic conditions in S. cerevisiae, which is a relevant issue to study for industrial-scale production. In this work we studied the uptake and metabolism of fumaric acid in S. cerevisiae in glucose-limited chemostat cultures at a cultivation pH of 3.0 (pH < pK). Steady states were achieved with different extracellular levels of fumaric acid, obtained by adding different amounts of fumaric acid to the feed medium. The experiments were carried out with the wild-type S. cerevisiae CEN.PK 113-7D and an engineered S. cerevisiae ADIS 244 expressing a heterologous dicarboxylic acid transporter (DCT-02) from Aspergillus niger, to examine whether it would be capable of exporting fumaric acid. We observed that fumaric acid entered the cells most likely via passive diffusion of the undissociated form. Approximately two-thirds of the fumaric acid in the feed was metabolized together with glucose. From metabolic flux analysis, an increased ATP dissipation was observed only at high intracellular concentrations of fumarate, possibly due to the export of fumarate via an ABC transporter. The implications of our results for the industrial-scale production of fumaric acid are discussed.

Comparative fluxome and metabolome analysis for overproduction of succinate in Escherichia coli

An aerobic succinate-producing Escherichia coli mutant was compared to its wild-type by quantitatively analyzing both the metabolome and fluxome, during glucose-limited steady-state and succinate excess dynamic conditions, in order to identify targets for further strain engineering towards more efficient succinate production. The mutant had four functional mutations under the conditions investigated: increased expression of a succinate exporter (DcuC), deletion of a succinate importer (Dct), deletion of succinate dehydrogenase (SUCDH) and expression of a PEP carboxylase (PPC) with increased capacity due to a point mutation. The steady-state and dynamic patterns of the intracellular metabolite levels and fluxes in response to changes were used to locate the quantitative differences in the physiology/metabolism of the mutant strain. Unexpectedly the mutant had a higher energy efficiency, indicated by a much lower rate of oxygen consumption, under glucose-limited conditions, caused by the deletion of the transcription factors IclR and ArcA. Furthermore the mutant had a much lower uptake capacity for succinate (26-fold) and oxygen (17-fold under succinate excess) compared to the wild-type strain. The mutant strain produced 7.9 mmol.CmolX(-1).h(-1) succinate during chemostat cultivation, showing that the choice of the applied genetic modifications was a successful strategy. Furthermore, the applied genetic modifications resulted in multiple large changes in metabolite levels (FBP, pyruvate, 6PG, NAD(+) /NADH ratio, α-ketogluarate) corresponding to large changes in fluxes. Compared to the wild-type a considerable flux shift occurred from the tricarboxylic acid (TCA) cycle to the oxidative part of the pentose phosphate pathway, including an inversion of the pyruvate kinase flux. The mutant responded very differently to excess of succinate, with a remarkable possible reversal of the TCA cycle. The mutant and the wild-type both showed homeostatic behaviour with respect to the energy charge. In contrast, large changes in redox ratios (NAD(+) /NADH) occurred in the wild-type, while the mutant showed even larger changes. This large redox change can be associated to the reversal of flux directions. The observed large flexibility in the central metabolism following genetic (deletions) and environmental (substrate excess) perturbations of the mutant, indicates that introducing a more efficient succinate exporter could result in an even higher succinate production rate.

Determination of the Cytosolic NADPH/NADP Ratio in Saccharomyces cerevisiae using Shikimate Dehydrogenase as Sensor Reaction

Eukaryotic metabolism is organised in complex networks of enzyme catalysed reactions which are distributed over different organelles. To quantify the compartmentalised reactions, quantitative measurements of relevant physiological variables in different compartments are needed, especially of cofactors. NADP(H) are critical components in cellular redox metabolism. Currently, available metabolite measurement methods allow whole cell measurements. Here a metabolite sensor based on a fast equilibrium reaction is introduced to monitor the cytosolic NADPH/NADP ratio in Saccharomyces cerevisiae: NADP + shikimate ⇄ NADPH + H(+) + dehydroshikimate. The cytosolic NADPH/NADP ratio was determined by measuring the shikimate and dehydroshikimate concentrations (by GC-MS/MS). The cytosolic NADPH/NADP ratio was determined under batch and chemostat (aerobic, glucose-limited, D = 0.1 h(-1)) conditions, to be 22.0 ± 2.6 and 15.6 ± 0.6, respectively. These ratios were much higher than the whole cell NADPH/NADP ratio (1.05 ± 0.08). In response to a glucose pulse, the cytosolic NADPH/NADP ratio first increased very rapidly and restored the steady state ratio after 3 minutes. In contrast to this dynamic observation, the whole cell NADPH/NADP ratio remained nearly constant. The novel cytosol NADPH/NADP measurements provide new insights into the thermodynamic driving forces for NADP(H)-dependent reactions, like amino acid synthesis, product pathways like fatty acid production or the mevalonate pathway.

Determination of the in vivo NAD:NADH ratio in Saccharomyces cerevisiae under anaerobic conditions, using alcohol dehydrogenase as sensor reaction

With the current quantitative metabolomics techniques, only whole-cell concentrations of NAD and NADH can be quantified. These measurements cannot provide information on the in vivo redox state of the cells, which is determined by the ratio of the free forms only. In this work we quantified free NAD:NADH ratios in yeast under anaerobic conditions, using alcohol dehydrogenase (ADH) and the lumped reaction of glyceraldehyde-3-phosphate dehydrogenase and 3-phosphoglycerate kinase as sensor reactions. We showed that, with an alternative accurate acetaldehyde determination method, based on rapid sampling, instantaneous derivatization with 2,4 diaminophenol hydrazine (DNPH) and quantification with HPLC, the ADH-catalysed oxidation of ethanol to acetaldehyde can be applied as a relatively fast and simple sensor reaction to quantify the free NAD:NADH ratio under anaerobic conditions. We evaluated the applicability of ADH as a sensor reaction in the yeast Saccharomyces cerevisiae, grown in anaerobic glucose-limited chemostats under steady-state and dynamic conditions. The results found in this study showed that the cytosolic redox status (NAD:NADH ratio) of yeast is at least one order of magnitude lower, and is thus much more reduced, under anaerobic conditions compared to aerobic glucose-limited steady-state conditions. The more reduced state of the cytosol under anaerobic conditions has major implications for (central) metabolism. Accurate determination of the free NAD:NADH ratio is therefore of importance for the unravelling of in vivo enzyme kinetics and to judge accurately the thermodynamic reversibility of each redox reaction.

A fast sensor for in vivo quantification of cytosolic phosphate in Saccharomyces cerevisiae

Eukaryotic metabolism consists of a complex network of enzymatic reactions and transport processes which are distributed over different subcellular compartments. Currently, available metabolite measurement protocols allow to measure metabolite whole cell amounts which hinder progress to describe the in vivo dynamics in different compartments, which are driven by compartment specific concentrations. Phosphate (Pi) is an essential component for: (1) the metabolic balance of upper and lower glycolytic flux; (2) Together with ATP and ADP determines the phosphorylation energy. Especially, the cytosolic Pi has a critical role in disregulation of glycolysis in tps1 knockout. Here we developed a method that enables us to monitor the cytosolic Pi concentration in S. cerevisiae using an equilibrium sensor reaction: maltose + Pi < = > glucose + glucose-1-phosphate. The required enzyme, maltose phosphorylase from L. sanfranciscensis was overexpressed in S. cerevisiae. With this reaction in place, the cytosolic Pi concentration was obtained from intracellular glucose, G1P and maltose concentrations. The cytosolic Pi concentration was determined in batch and chemostat (D = 0.1 h(-1) ) conditions, which was 17.88 µmol/gDW and 25.02 µmol/gDW, respectively under Pi-excess conditions. Under Pi-limited steady state (D = 0.1 h(-1) ) conditions, the cytosolic Pi concentration dropped to only 17.7% of the cytosolic Pi in Pi-excess condition (4.42 µmol/gDW vs. 25.02 µmol/gDW). In response to a Pi pulse, the cytosolic Pi increased very rapidly, together with the concentration of sugar phosphates. Main sources of the rapid Pi increase are vacuolar Pi (and not the polyPi), as well as Pi uptake from the extracellular space. The temporal increase of cytosolic Pi increases the driving force of GAPDH reaction of the lower glycolytic reactions. The novel cytosol specific Pi concentration measurements provide new insight into the thermodynamic driving force for ATP hydrolysis, GAPDH reaction, and Pi transport over the plasma and vacuolar membranes.

Development of tools for quantitative intracellular metabolomics of Aspergillus niger chemostat cultures

In view of the high citric acid production capacity of Aspergillus niger, it should be well suited as a cell factory for the production of other relevant acids as succinic, fumaric, itaconic and malic. Quantitative metabolomics is an important omics tool in a synthetic biology approach to develop A. niger for the production of these acids. Such studies require well defined and tightly controlled cultivation conditions and proper rapid sampling, sample processing and analysis methods. In this study we present the development of a chemostat for homogeneous steady state cultivation of A. niger, equipped with a new dedicated rapid sampling device. A quenching method for quantitative metabolomics in A. niger based on cold methanol was evaluated using balances and optimized with the aim of avoiding metabolite leakage during sample processing. The optimization was based on measurements of the intermediates of the glycolysis, TCA and PPP pathways and amino acids, using a balance approach. Leakage was found to be absent at -20 °C for a 40 % (v/v) methanol concentration in water. Under these conditions the average metabolite recovery was close to 100 %. When comparing A. niger and Penicillium chrysogenum metabolomes, under the same cultivation conditions, similar metabolite fingerprints were found in both fungi, except for the intracellular citrate level which is higher for A. niger.

Fast "Feast/Famine" Cycles for Studying Microbial Physiology Under Dynamic Conditions: A Case Study with Saccharomyces cerevisiae

Microorganisms are constantly exposed to rapidly changing conditions, under natural as well as industrial production scale environments, especially due to large-scale substrate mixing limitations. In this work, we present an experimental approach based on a dynamic feast/famine regime (400 s) that leads to repetitive cycles with moderate changes in substrate availability in an aerobic glucose cultivation of Saccharomyces cerevisiae. After a few cycles, the feast/famine produced a stable and repetitive pattern with a reproducible metabolic response in time, thus providing a robust platform for studying the microorganism's physiology under dynamic conditions. We found that the biomass yield was slightly reduced (-5%) under the feast/famine regime, while the averaged substrate and oxygen consumption as well as the carbon dioxide production rates were comparable. The dynamic response of the intracellular metabolites showed specific differences in comparison to other dynamic experiments (especially stimulus-response experiments, SRE). Remarkably, the frequently reported ATP paradox observed in single pulse experiments was not present during the repetitive perturbations applied here. We found that intracellular dynamic accumulations led to an uncoupling of the substrate uptake rate (up to 9-fold change at 20 s.) Moreover, the dynamic profiles of the intracellular metabolites obtained with the feast/famine suggest the presence of regulatory mechanisms that resulted in a delayed response. With the feast famine setup many cellular states can be measured at high frequency given the feature of reproducible cycles. The feast/famine regime is thus a versatile platform for systems biology approaches, which can help us to identify and investigate metabolite regulations under realistic conditions (e.g., large-scale bioreactors or natural environments).

Fatal attraction in glycolysis: how Saccharomyces cerevisiae manages sudden transitions to high glucose

In the model eukaryote Saccharomyces cerevisiae, it has long been known that a functional trehalose pathway is indispensable for transitions to high glucose conditions. Upon addition of glucose, cells with a defect in trehalose 6-phosphate synthase (Tps1), the first committed step in the trehalose pathway, display what we have termed an imbalanced glycolytic state; in this state the flux through the upper part of glycolysis outpaces that through the lower part of glycolysis. As a consequence, the intermediate fructose 1,6-bisphosphate (FBP) accumulates at low concentrations of ATP and inorganic phosphate (P_i). Despite significant research efforts, a satisfactory understanding of the regulatory role that trehalose metabolism plays during such transitions has remained infamously unresolved. In a recent study, we demonstrate that the startup of glycolysis exhibits two dynamic fates: a proper, functional, steady state or the imbalanced state described above. Both states are stable, attracting states, and the probability distribution of initial states determines the fate of a yeast cell exposed to glucose. Trehalose metabolism steers the dynamics of glycolysis towards the proper functional state through its ATP hydrolysis activity; a mechanism that ensures that the demand and supply of ATP is balanced with P_i availability under dynamic conditions. [van Heerden et al. Science (2014), DOI: 10.1126/science.1245114.]

Lost in transition: start-up of glycolysis yields subpopulations of nongrowing cells

Cells need to adapt to dynamic environments. Yeast that fail to cope with dynamic changes in the abundance of glucose can undergo growth arrest. We show that this failure is caused by imbalanced reactions in glycolysis, the essential pathway in energy metabolism in most organisms. The imbalance arises largely from the fundamental design of glycolysis, making this state of glycolysis a generic risk. Cells with unbalanced glycolysis coexisted with vital cells. Spontaneous, nongenetic metabolic variability among individual cells determines which state is reached and, consequently, which cells survive. Transient ATP (adenosine 5'-triphosphate) hydrolysis through futile cycling reduces the probability of reaching the imbalanced state. Our results reveal dynamic behavior of glycolysis and indicate that cell fate can be determined by heterogeneity purely at the metabolic level.

Quantitative metabolomics using ID-MS

Quantitative intracellular metabolite measurements are essential for systems biology and modeling of cellular metabolism. The MS-based quantification is error prone because (1) several sampling processing steps have to be performed, (2) the sample contains a complex mixture of partly compounds with the same mass and similar retention time, and (3) especially salts influence the ionization efficiency. Therefore internal standards are required, best for each measured compound. The use of labeled biomass, (13)C extract, is a valuable tool, reducing the standard deviations of intracellular concentration measurements significantly (especially regarding technical reproducibility). Using different platforms, i.e., LC-MS and GC-MS, a large number of different metabolites can be quantified (currently about 110).

Glucose-methanol co-utilization in Pichia pastoris studied by metabolomics and instationary ¹³C flux analysis

Background

Several studies have shown that the utilization of mixed carbon feeds instead of methanol as sole carbon source is beneficial for protein production with the methylotrophic yeast Pichia pastoris. In particular, growth under mixed feed conditions appears to alleviate the metabolic burden related to stress responses triggered by protein overproduction and secretion. Yet, detailed analysis of the metabolome and fluxome under mixed carbon source metabolizing conditions are missing. To obtain a detailed flux distribution of central carbon metabolism, including the pentose phosphate pathway under methanol-glucose conditions, we have applied metabolomics and instationary ¹³C flux analysis in chemostat cultivations.

Results

Instationary ¹³C-based metabolic flux analysis using GC-MS and LC-MS measurements in time allowed for an accurate mapping of metabolic fluxes of glycolysis, pentose phosphate and methanol assimilation pathways. Compared to previous results from NMR-derived stationary state labelling data (proteinogenic amino acids, METAFoR) more fluxes could be determined with higher accuracy. Furthermore, using a thermodynamic metabolic network analysis the metabolite measurements and metabolic flux directions were validated. Notably, the concentration of several metabolites of the upper glycolysis and pentose phosphate pathway increased under glucose-methanol feeding compared to the reference glucose conditions, indicating a shift in the thermodynamic driving forces. Conversely, the extracellular concentrations of all measured metabolites were lower compared with the corresponding exometabolome of glucose-grown P. pastoris cells.

The instationary ¹³C flux analysis resulted in fluxes comparable to previously obtained from NMR datasets of proteinogenic amino acids, but allowed several additional insights. Specifically, i) in vivo metabolic flux estimations were expanded to a larger metabolic network e.g. by including trehalose recycling, which accounted for about 1.5% of the glucose uptake rate; ii) the reversibility of glycolytic/gluconeogenesis, TCA cycle and pentose phosphate pathways reactions was estimated, revealing a significant gluconeogenic flux from the dihydroxyacetone phosphate/glyceraldehydes phosphate pool to glucose-6P. The origin of this finding could be carbon recycling from the methanol assimilatory pathway to the pentose phosphate pool. Additionally, high exchange fluxes of oxaloacetate with aspartate as well as malate indicated amino acid pool buffering and the activity of the malate/Asp shuttle; iii) the ratio of methanol oxidation vs utilization appeared to be lower (54 vs 79% assimilated methanol directly oxidized to CO₂).

Conclusions

In summary, the application of instationary ¹³C-based metabolic flux analysis to P. pastoris provides an experimental framework with improved capabilities to explore the regulation of the carbon and energy metabolism of this yeast, particularly for the case of methanol and multicarbon source metabolism.

Changes in substrate availability in Escherichia coli lead to rapid metabolite, flux and growth rate responses

The interactions between the intracellular metabolome, fluxome and growth rate of Escherichia coli after sudden glycolytic/gluconeogenic substrate shifts are studied based on pulses of different substrates to an aerobic glucose-limited steady-state (dilution rate=0.1h(-1)). After each added glycolytic (glucose) and gluconeogenic (pyruvate and succinate) substrate pulse, no by-products were secreted and a pseudo steady state in flux and metabolites was achieved in about 30-40s. In the pulse experiments a large oxygen uptake capacity of the cells was observed. The in vivo dynamic responses showed massive reorganization and flexibility (1/100-14-fold change) of extra/intracellular metabolic fluxes, matching with large changes in the concentrations of intracellular metabolites, including reversal of reaction rate for pseudo/near equilibrium reactions. The coupling of metabolome and fluxome could be described by Q-linear kinetics. Remarkably, the three different substrate pulses resulted in a very similar increase in growth rate (0.13-0.3h(-1)). Data analysis showed that there must exist as yet unknown mechanisms which couple the protein synthesis rate to changes in central metabolites.

Similar temperature dependencies of glycolytic enzymes: an evolutionary adaptation to temperature dynamics?

BACKGROUND

Temperature strongly affects microbial growth, and many microorganisms have to deal with temperature fluctuations in their natural environment. To understand regulation strategies that underlie microbial temperature responses and adaptation, we studied glycolytic pathway kinetics in Saccharomyces cerevisiae during temperature changes.

RESULTS

Saccharomyces cerevisiae was grown under different temperature regimes and glucose availability conditions. These included glucose-excess batch cultures at different temperatures and glucose-limited chemostat cultures, subjected to fast linear temperature shifts and circadian sinoidal temperature cycles. An observed temperature-independent relation between intracellular levels of glycolytic metabolites and residual glucose concentration for all experimental conditions revealed that it is the substrate availability rather than temperature that determines intracellular metabolite profiles. This observation corresponded with predictions generated in silico with a kinetic model of yeast glycolysis, when the catalytic capacities of all glycolytic enzymes were set to share the same normalized temperature dependency.

CONCLUSIONS

From an evolutionary perspective, such similar temperature dependencies allow cells to adapt more rapidly to temperature changes, because they result in minimal perturbations of intracellular metabolite levels, thus circumventing the need for extensive modification of enzyme levels.

Substrate cycles in Penicillium chrysogenum quantified by isotopic non-stationary flux analysis

Background

Penicillium chrysogenum, the main production strain for penicillin-G, has a high content of intracellular carbohydrates, especially reduced sugars such as mannitol, arabitol, erythritol, as well as trehalose and glycogen. In previous steady state ¹³C wash-in experiments a delay of labeling enrichments in glycolytic intermediates was observed, which suggests turnover of storage carbohydrates. The turnover of storage pools consumes ATP which is expected to reduce the product yield for energy demanding production pathways like penicillin-G.

Results

In this study, a ¹³C labeling wash-in experiment of 1 hour was performed to systematically quantify the intracellular flux distribution including eight substrate cycles. The experiments were performed using a mixed carbon source of 85% Cmol_Glc/Cmol_Glc+EtOH labeled glucose (mixture of 90% [1-¹³C₁] and 10% [U-¹³C₆]) and 15% ethanol [U-¹³C₂]. It was found, that (1) also several extracellular pools are enriched with ¹³C labeling rapidly (trehalose, mannitol, and others), (2) the intra- to extracellular metabolite concentration ratios were comparable for a large set of metabolites while for some carbohydrates (mannitol, trehalose, and glucose) the measured ratios were much higher.

Conclusions

The fast enrichment of several extracellular carbohydrates and a concentration ratio higher than the ratio expected from cell lysis (2%) indicate active (e.g. ATP consuming) transport cycles over the cellular membrane. The flux estimation indicates, that substrate cycles account for about 52% of the gap in the ATP balance based on metabolic flux analysis.

A metabolic network stoichiometry analysis of microbial growth and product formation

Using available biochemical information, metabolic networks have been constructed to describe the biochemistry of growth of Saccharomyces cerevisiae and Candida utilis on a wide variety of carbon substrates. All networks contained only two fitted parameters, the P/O ratio and a maintenance coefficient. It is shown that with a growth-associated maintenance coefficient, K, of 1.37 mol ATP/ C-mol protein for both yeasts and P/O ratios of 1.20 and 1.53 for S. cerevisiae and C. utilis, respectively, measured biomass yields could be described accurately. A metabolic flux analysis of aerobic growth of S. cerevisiae on glucose/ethanol mixtures predicted five different metabolic flux regimes upon transition from 100% glucose to 100% ethanol. The metabolic network constructed for growth of S. cerevisiae on glucose was applied to perform a theoretical exercise on the overproduction of amino acids. It is shown that theoretical operational product yield values can be substantially lower than calculated maximum product yields. A practical case of lysine production was analyzed with respect to theoretical bottlenecks limiting product formation. Predictions of network-derived irreversibility limits for Y(sp) (mu) functions were compared with literature data. The comparisons show that in real systems such irreversibility constraints may be of relevance. It is concluded that analysis of metabolic network stoichiometry is a useful tool to detect metabolic limits and to guide process intensification studies. (c) 1995 John Wiley & Sons, Inc.

Linear constraint relations in biochemical reaction systems: I. Classification of the calculability and the balanceability of conversion rates

Measurements provide the basis for process monitoring and control as well as for model development and validation. Systematic approaches to increase the accuracy and credibility of the empirical data set are therefore of great value. In (bio)chemical conversions, linear conservation relations such as the balance equations for charge, enthalpy, and/or chemical elements, can be employed to relate conversion rates. In a pactical situation, some of these rates will be measured (in effect, be calculated directly from primary measurements of, e.g., concentrations and flow rates), as others can or cannot be calculated from the measured ones. When certain measured rates can also be calculated from other measured rates, the set of equations, the accuracy and credibility of the measured rates can indeed be improved by, respectively, balancing and gross error diagnosis. The balanced conversion rates are more accurate, and form a consistent set of data, which is more suitable for further application (e.g., to calculate nonmeasured rates) than the raw measurements. Such an approach has drawn attention in previous studies. The current study deals mainly with the problem of mathematically classifying the conversion rates into balanceable and calculable rates, given the subset of measured rates. The significance of this problem is illustrated with some examples. It is shown that a simple matrix equation can be derived that contains the vector of measured conversion rates and the redundancy matrix R. Matrix R plays a predominant role in the classification problem. In supplementary articles, significance of the redundancy matrix R for an improved gross error diagnosis approach will be shown. In addition, efficient equations have been derived to calculate the balanceable and/or calculable rates. The method is completely based on matrix algebra (principally different from the graph-theoretical approach), and it is easily implemented into a computer program. (c) 1994 John Wiley & Sons, Inc.

Linear constraint relations in biochemical reaction systems: II. Diagnosis and estimation of gross errors

Conservation equations derived from elemental balances, heat balances, and metabolic stoichiometry, can be used to constrain the values of conversion rates of relevant components. In the present work, their use will be discussed for detection and localization of significant errors of the following types: 1.At least one of the primary measurements has a significant error (gross measurement error).2.The system definition is incorrect: a component a.is not included in the system description.b.has a composition different from that specified.3.The specified variances are too small, resulting in a too-sensitive test.The error diagnosis technique presented here, is based on the following: given the conservation equations, for each set of measured rates, a vector of residuals of these equations can be constructed, of which the direction is related to the error source, as its length is a measure of the error size. The similarity of the directions of such a residual vector and certain compare vectors, each corresponding to a specific error source, is considered in a statistical test. If two compare vectors that result from different error sources have (almost) the same direction, errors of these types cannot be distinguished from each other. For each possible error in the primary measurements of flows and concentrations, the compare vector can be constructed a priori, thus allowing analysis beforehand, which errors can be observed. Therefore, the detectability of certain errors likely to occur can be insured by selecting a proper measurement set. The possibility of performing this analysis before experiments are carried out is an important advantage, providing a profound understanding of the detectability of errors. The characteristics of the method with respect to diagnosis of simultaneous errors and error size estimation are discussed and compared to those of the serial elimination method and the serial compensation strategy, published elsewhere. (c) 1994 John Wiley & Sons, Inc.

Detachment of biomass from suspended nongrowing spherical biofilms in airlift reactors

In three-phase internal loop airlift reactors, the detachment of biomass from suspended biofilm pellets in the presence of bare carrier particles was investigated under nongrowth conditions. The detachment rate was dominated by collisions between bare carrier particles and biofilm pellets. The concentration of bare carrier particles and the carrier roughness strongly influenced the detachment rate. A change in flow regime from bubbling to slug flow considerably increased the detachment rate. Otherwise, the superficial gas velocity did not directly affect the detachment rate. The influence of particle size was not clear. The bottom clearance did not affect the detachment rate within the tested range. Other aspects of reactor geometry might be important. The main detachment processes were abrasion and breakage of biofilm pellets. During the detachment process, two phases could be distinguished. In the first phase the detachment was relatively high, and both breakage and abrasion of biofilm pellets occurred. During the second phase, breakage dominated and the detachment rate was lower. The two-phase behavior is explained by differences in strength between the inner and outer biofilm layers, possibly caused by variations in local growth rates during biofilm formation. Differences in growth history might also explain the various detachment rates observed with different biofilm batches. (c) 1995 John Wiley & Sons, Inc.

The role of glucose in ajmalicine production by catharanthus roseus cell cultures

The role of glucose in ajmalicine production by Catharanthus roseus was investigated in the second stage of a two-stage batch process. Activities of tryptophan decar-boxylate (TDC) and anthranilate synthase (AS), two enzymes In the pathway leading to ajmalicine, were higher after induction with 40 g/L glucose than after induction with 60 or 80 g/L glucose. Experiments with different media containing mixtures of glucose and the nonpermeating osmotic agent xylose, and using an already induced culture as inoculum, revealed that a minimum amount of glucose is required to support ajmalicine production after enzyme induction. This requirement was not an osmotic effect. The relation between the glucose concentration and the specific ajmalicine production rate, q(p), was investigated in seven (fed-)batch cultures with constant glucose concentrations: 23, 29, 35, 53, 57, 75, and 98 g/L. In the cultures with a low glucose concentration (23, 29, and 35 g/L) the q(p) was 2.7-times higher than the cultures with 53 and 57 g/L, and almost six times higher than the cultures with a high glucose concentration (75 and 98 g/L). A glucose perturbation experiment (from 53 to 32 g/L) demonstrated that the ajmalicine production rate was adjusted without much delay. A kinetic equation is proposed for the relationship between the glucose concentration and q(p). Differences in enzyme induction and ajmalicine production at different glucose levels could not be explained by the intracellular concentrations of glucose, fructose, sucrose, or starch. (c) 1995 John Wiley & Sons Inc.

Optimization of cold methanol quenching for quantitative metabolomics of Penicillium chrysogenum

A sampling procedure for quantitative metabolomics in Penicillium chrysogenum based on cold aqueous methanol quenching was re-evaluated and optimized to reduce metabolite leakage during sample treatment. The optimization study included amino acids and intermediates of the glycolysis and the TCA-cycle. Metabolite leakage was found to be minimal for a methanol content of the quenching solution (QS) of 40% (v/v) while keeping the temperature of the quenched sample near -20°C. The average metabolite recovery under these conditions was 95.7% (±1.1%). Several observations support the hypothesis that metabolite leakage from quenched mycelia of P. chrysogenum occurs by diffusion over the cell membrane. First, a prolonged contact time between mycelia and the QS lead to a somewhat higher extent of leakage. Second, when suboptimal quenching liquids were used, increased metabolite leakage was found to be correlated with lower molecular weight and with lower absolute net charge. The finding that lowering the methanol content of the quenching liquid reduces metabolite leakage in P. chrysogenum contrasts with recently published quenching studies for two other eukaryotic micro-organisms. This demonstrates that it is necessary to validate and, if needed, optimize the quenching conditions for each particular micro-organism. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11306-011-0367-3) contains supplementary material, which is available to authorized users.

Quantitative metabolomics analysis of amino acid metabolism in recombinant Pichia pastoris under different oxygen availability conditions

Background

Environmental and intrinsic stress factors can result in the global alteration of yeast physiology, as evidenced by several transcriptional studies. Hypoxia has been shown to have a beneficial effect on the expression of recombinant proteins in Pichia pastoris growing on glucose. Furthermore, transcriptional profiling analyses revealed that oxygen availability was strongly affecting ergosterol biosynthesis, central carbon metabolism and stress responses, in particular the unfolded protein response. To contribute to the better understanding of the effect and interplay of oxygen availability and foreign protein secretion on central metabolism, a first quantitative metabolomic analysis of free amino acids pools in a recombinant P. pastoris strain growing under different oxygen availability conditions has been performed.

Results

The values obtained indicate significant variations in the intracellular amino acid pools due to different oxygen availability conditions, showing an overall increase of their size under oxygen limitation. Notably, even while foreign protein productivities were relatively low (about 40–80 μg Fab/g_DCW·h), recombinant protein production was found to have a limited but significant impact on the intracellular amino acid pools, which were generally decreased in the producing strain compared with the reference strain. However, observed changes in individual amino acids pools were not correlated with their corresponding relative abundance in the recombinant protein sequence, but to the overall cell protein amino acid compositional variations.

Conclusions

Overall, the results obtained, combined with previous transcriptomic and proteomic analyses provide a systematic metabolic fingerprint of the oxygen availability impact on recombinant protein production in P. pastoris.

Development of quantitative metabolomics for Pichia pastoris

Accurate, reliable and reproducible measurement of intracellular metabolite levels has become important for metabolic studies of microbial cell factories. A first critical step for metabolomic studies is the establishment of an adequate quenching and washing protocol, which ensures effective arrest of all metabolic activity and removal of extracellular metabolites, without causing leakage of metabolites from the cells. Five different procedures based on cold methanol quenching and cell separation by filtration were tested for metabolomics of Pichia pastoris regarding methanol content and temperature of the quenching solution as key parameters. Quantitative evaluation of these protocols was carried out through mass balance analysis, based on metabolite measurements in all sample fractions, those are whole broth, quenched and washed cells, culture filtrate and quenching and washing solution. Finally, the optimal method was used to study the time profiles of free amino acid and central carbon metabolism intermediates in glucose-limited chemostat cultures. Acceptable recoveries (>90%) were obtained for all quenching procedures tested. However, quenching at -27°C in 60% v/v methanol performed slightly better in terms of leakage minimization. We could demonstrate that five residence times under glucose limitation are enough to reach stable intracellular metabolite pools. Moreover, when comparing P. pastoris and S. cerevisiae metabolomes, under the same cultivation conditions, similar metabolite fingerprints were found in both yeasts, except for the lower glycolysis, where the levels of these metabolites in P. pastoris suggested an enzymatic capacity limitation in that part of the metabolism. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11306-011-0308-1) contains supplementary material, which is available to authorized users.

Identification of informative metabolic responses using a minibioreactor: a small step change in the glucose supply rate creates a large metabolic response in Saccharomyces cerevisiae

In this study, a previously developed mini-bioreactor, the Biocurve, was used to identify an informative stimulus-response experiment. The identified stimulus-response experiment was a modest 50% shift-up in glucose uptake rate (qGLC) that unexpectedly resulted in a disproportionate transient metabolic response. The 50% shift-up in qGLC in the Biocurve resulted in a near tripling of the online measured oxygen uptake (qO2) and carbon dioxide production (qCO2) rates, suggesting a considerable mobilization of glycogen and trehalose. The 50% shift-up in qGLC was subsequently studied in detail in a conventional bioreactor (4 l working volume), which confirmed the results obtained with the Biocurve. Especially relevant is the observation that the 50% increase in glucose uptake rate led to a three-fold increase in glycolytic flux, due to mobilization of storage materials. This explains the unexpected ethanol and acetate secretion after the shift-up, in spite of the fact that after the shift-up the qGLC was far less than the critical value. Moreover, these results show that the correct in vivo fluxes in glucose pulse experiments cannot be obtained from the uptake and secretion rates only. Instead, the storage fluxes must also be accurately quantified. Finally, we speculate on the possible role that the transient increase in dissolved CO2 immediately after the 50% shift-up in qGLC could have played a part in triggering glycogen and trehalose mobilization.

Use of sequential-batch fermentations to characterize the impact of mild hypothermic temperatures on the anaerobic stoichiometry and kinetics of Saccharomyces cerevisiae

This work presents a characterization of the stoichiometry and kinetics of anaerobic batch growth of Saccharomyces cerevisiae at cultivation temperatures between 12 and 30°C. To minimize the influence of the inoculum condition and ensure full adaptation to the cultivation temperature, the experiments were carried out in sequencing batch reactors. It was observed that the growth rate obtained in the first batch performed after each temperature shift was 10-30% different compared with the subsequent batches at the same temperature, which were much more reproducible. This indicates that the sequencing batch approach provides accurate and reproducible growth rate data. Data reconciliation was applied to the measured time patterns of substrate, biomass, carbon dioxide and byproducts with the constraint that the elemental conservation relations were satisfied, allowing to obtain consistent best estimates of all uptake and secretion rates. Subsequently, it was attempted to obtain an appropriate model description of the temperature dependency of these rates. It was found that the Ratkowsky model provided a better description of the temperature dependency of growth, uptake and secretion rates than the Arrhenius law. Most interesting was to find that most of the biomass-specific rates have the same temperature dependency, leading to a near temperature independent batch stoichiometry.

pH-dependent uptake of fumaric acid in Saccharomyces cerevisiae under anaerobic conditions

Microbial production of C(4) dicarboxylic acids from renewable resources has gained renewed interest. The yeast Saccharomyces cerevisiae is known as a robust microorganism and is able to grow at low pH, which makes it a suitable candidate for biological production of organic acids. However, a successful metabolic engineering approach for overproduction of organic acids requires an incorporation of a proper exporter to increase the productivity. Moreover, low-pH fermentations, which are desirable for facilitating the downstream processing, may cause back diffusion of the undissociated acid into the cells with simultaneous active export, thereby creating an ATP-dissipating futile cycle. In this work, we have studied the uptake of fumaric acid in S. cerevisiae in carbon-limited chemostat cultures under anaerobic conditions. The effect of the presence of fumaric acid at different pH values (3 to 5) has been investigated in order to obtain more knowledge about possible uptake mechanisms. The experimental results showed that at a cultivation pH of 5.0 and an external fumaric acid concentration of approximately 0.8 mmol · liter(-1), the fumaric acid uptake rate was unexpectedly high and could not be explained by diffusion of the undissociated form across the plasma membrane alone. This could indicate the presence of protein-mediated import. At decreasing pH levels, the fumaric acid uptake rate was found to increase asymptotically to a maximum level. Although this observation is in accordance with protein-mediated import, the presence of a metabolic bottleneck for fumaric acid conversion under anaerobic conditions could not be excluded.

Fast sampling of the cellular metabolome

Obtaining meaningful snapshots of the metabolome of microorganisms requires rapid sampling and immediate quenching of all metabolic activity, to prevent any changes in metabolite levels after sampling. Furthermore, a suitable extraction method is required ensuring complete extraction of metabolites from the cells and inactivation of enzymatic activity, with minimal degradation of labile compounds. Finally a sensitive, high-throughput analysis platform is needed to quantify a large number of metabolites in a small amount of sample. An issue which has often been overlooked in microbial metabolomics is the fact that many intracellular metabolites are also present in significant amounts outside the cells, and may interfere with the endometabolome measurements. Attempts to remove the extracellular metabolites with dedicated quenching methods often induce release of intracellular metabolites into the quenching solution. For eukaryotic microorganisms, leakage can be minimized by adaptation of the quenching method. For prokaryotic cells this had not yet been accomplished, so the application of a differential method whereby metabolites are measured in the culture supernatant as well as in total broth samples, to calculate the intracellular levels by subtraction, seems to be the most suitable approach. Here we present an overview of different sampling, quenching, and extraction methods developed for microbial metabolomics, described in the literature. Detailed protocols are provided for rapid sampling, quenching, and extraction for measurement of metabolites in total broth samples, washed cell samples and supernatant, to be applied for quantitative metabolomics of both eukaryotic and prokaryotic microorganisms.

Scale-down of penicillin production in Penicillium chrysogenum

In large-scale production reactors the combination of high broth viscosity and large broth volume leads to insufficient liquid-phase mixing, resulting in gradients in, for example, the concentrations of substrate and oxygen. This often leads to differences in productivity of the full-scale process compared with laboratory scale. In this scale-down study of penicillin production, the influence of substrate gradients on process performance and cell physiology was investigated by imposing an intermittent feeding regime on a laboratory-scale culture of a high yielding strain of Penicillium chrysogenum. It was found that penicillin production was reduced by a factor of two in the intermittently fed cultures relative to constant feed cultivations fed with the same amount of glucose per hour, while the biomass yield was the same. Measurement of the levels of the intermediates of the penicillin biosynthesis pathway, along with the enzyme levels, suggested that the reduction of the flux through the penicillin pathway is mainly the result of a lower influx into the pathway, possibly due to inhibitory levels of adenosine monophosphate and pyrophosphate and lower activating levels of adenosine triphosphate during the zero-substrate phase of each cycle of intermittent feeding.

Integrated multilaboratory systems biology reveals differences in protein metabolism between two reference yeast strains

The field of systems biology is often held back by difficulties in obtaining comprehensive, high-quality, quantitative data sets. In this paper, we undertook an interlaboratory effort to generate such a data set for a very large number of cellular components in the yeast Saccharomyces cerevisiae, a widely used model organism that is also used in the production of fuels, chemicals, food ingredients and pharmaceuticals. With the current focus on biofuels and sustainability, there is much interest in harnessing this species as a general cell factory. In this study, we characterized two yeast strains, under two standard growth conditions. We ensured the high quality of the experimental data by evaluating a wide range of sampling and analytical techniques. Here we show significant differences in the maximum specific growth rate and biomass yield between the two strains. On the basis of the integrated analysis of the high-throughput data, we hypothesize that differences in phenotype are due to differences in protein metabolism.

Effect of iclR and arcA knockouts on biomass formation and metabolic fluxes in Escherichia coli K12 and its implications on understanding the metabolism of Escherichia coli BL21 (DE3)

Background

Gene expression is regulated through a complex interplay of different transcription factors (TFs) which can enhance or inhibit gene transcription. ArcA is a global regulator that regulates genes involved in different metabolic pathways, while IclR as a local regulator, controls the transcription of the glyoxylate pathway genes of the aceBAK operon. This study investigates the physiological and metabolic consequences of arcA and iclR deletions on E. coli K12 MG1655 under glucose abundant and limiting conditions and compares the results with the metabolic characteristics of E. coli BL21 (DE3).

Results

The deletion of arcA and iclR results in an increase in the biomass yield both under glucose abundant and limiting conditions, approaching the maximum theoretical yield of 0.65 c-mole/c-mole glucose under glucose abundant conditions. This can be explained by the lower flux through several CO₂ producing pathways in the E. coli K12 ΔarcAΔiclR double knockout strain. Due to iclR gene deletion, the glyoxylate pathway is activated resulting in a redirection of 30% of the isocitrate molecules directly to succinate and malate without CO₂ production. Furthermore, a higher flux at the entrance of the TCA was noticed due to arcA gene deletion, resulting in a reduced production of acetate and less carbon loss. Under glucose limiting conditions the flux through the glyoxylate pathway is further increased in the ΔiclR knockout strain, but this effect was not observed in the double knockout strain. Also a striking correlation between the glyoxylate flux data and the isocitrate lyase activity was observed for almost all strains and under both growth conditions, illustrating the transcriptional control of this pathway. Finally, similar central metabolic fluxes were observed in E. coli K12 ΔarcA ΔiclR compared to the industrially relevant E. coli BL21 (DE3), especially with respect to the pentose pathway, the glyoxylate pathway, and the TCA fluxes. In addition, a comparison of the genome sequences of the two strains showed that BL21 possesses two mutations in the promoter region of iclR and rare codons are present in arcA implying a lower tRNA acceptance. Both phenomena presumably result in a reduced ArcA and IclR synthesis in BL21, which contributes to the similar physiology as observed in E. coli K12 ΔarcAΔiclR.

Conclusions

The deletion of arcA results in a decrease of repression on transcription of TCA cycle genes under glucose abundant conditions, without significantly affecting the glyoxylate pathway activity. IclR clearly represses transcription of glyoxylate pathway genes under glucose abundance, a condition in which Crp activation is absent. Under glucose limitation, Crp is responsible for the high glyoxylate flux, but IclR still represses transcription. Finally, in E. coli BL21 (DE3), ArcA and IclR are poorly expressed, explaining the similar fluxes observed compared to the ΔarcAΔiclR strain.

Escherichia coli responds with a rapid and large change in growth rate upon a shift from glucose-limited to glucose-excess conditions

Glucose pulse experiments at seconds time scale resolution were performed in aerobic glucose-limited Escherichia coli chemostat cultures. The dynamic responses of oxygen-uptake and growth rate at seconds time scale were determined using a new method based on the dynamic liquid-phase mass balance for oxygen and the pseudo-steady-state ATP balance. Significant fold changes in metabolites (10-1/10) and fluxes (4-1/4) were observed during the short (200 s) period of glucose excess. During glucose excess there was no secretion of by-products and the increased glucose uptake rate led within 40s to a 3.7 fold increase in growth rate. Also within 40-60s a new pseudo-steady-state was reached for both metabolite levels and fluxes. Flux changes of reactions were strongly correlated to the concentrations of involved compounds. Surprisingly the 3.7 fold increase in growth rate and hence protein synthesis rate was not matched by a significant increase in amino acid concentrations. This poses interesting questions for the kinetic factors, which drive protein synthesis by ribosomes.

Cytosolic NADPH balancing in Penicillium chrysogenum cultivated on mixtures of glucose and ethanol

The in vivo flux through the oxidative branch of the pentose phosphate pathway (oxPPP) in Penicillium chrysogenum was determined during growth in glucose/ethanol carbon-limited chemostat cultures, at the same growth rate. Non-stationary (13)C flux analysis was used to measure the oxPPP flux. A nearly constant oxPPP flux was found for all glucose/ethanol ratios studied. This indicates that the cytosolic NADPH supply is independent of the amount of assimilated ethanol. The cofactor assignment in the model of van Gulik et al. (Biotechnol Bioeng 68(6):602-618, 2000) was supported using the published genome annotation of P. chrysogenum. Metabolic flux analysis showed that NADPH requirements in the cytosol remain nearly the same in these experiments due to constant biomass growth. Based on the cytosolic NADPH balance, it is known that the cytosolic aldehyde dehydrogenase in P. chrysogenum is NAD(+) dependent. Metabolic modeling shows that changing the NAD(+)-aldehyde dehydrogenase to NADP(+)-aldehyde dehydrogenase can increase the penicillin yield on substrate.

Intracellular metabolite determination in the presence of extracellular abundance: Application to the penicillin biosynthesis pathway in Penicillium chrysogenum

Important steps in metabolic pathways are formed by the transport of substrates and products over the cell membrane. The study of in vivo transport kinetics requires accurate quantification of intra- and extracellular levels of the transported compounds. Especially in case of extracellular abundance, the proper determination of intracellular metabolite levels poses challenges. Efficient removal of extracellular substrates and products is therefore important not to overestimate the intracellular amounts. In this study we evaluated two different rapid sampling methods, one combined with cold filtration and the other with centrifugation, for their applicability to determine intracellular amounts of metabolites which are present in high concentrations in the extracellular medium. The filtration-based method combines fast sampling and immediate quenching of cellular metabolism in cold methanol, with rapid and effective removal of all compounds present outside the cells by means of direct filtration and subsequent filtration-based washing. In the centrifugation-based method, removal of the extracellular metabolites from the cells was achieved by means of multiple centrifugation and resuspension steps with the cold quenching solution. The cold filtration method was found to be highly superior to the centrifugation method to determine intracellular amounts of metabolites related to penicillin-G biosynthesis and allowed the quantification of compounds of which the extracellular amounts were 3-4 orders of magnitude higher than the intracellular amounts. Using this method for the first time allowed to measure the intracellular levels of the side chain precursor phenylacetic acid (PAA) and the product penicillin-G of the penicillin biosynthesis pathway, compounds of which the transport mechanism in Penicillium chrysogenum is still far from being sufficiently understood.

Genome-derived minimal metabolic models for Escherichia coli MG1655 with estimated in vivo respiratory ATP stoichiometry

Metabolic network models describing growth of Escherichia coli on glucose, glycerol and acetate were derived from a genome scale model of E. coli. One of the uncertainties in the metabolic networks is the exact stoichiometry of energy generating and consuming processes. Accurate estimation of biomass and product yields requires correct information on the ATP stoichiometry. The unknown ATP stoichiometry parameters of the constructed E. coli network were estimated from experimental data of eight different aerobic chemostat experiments carried out with E. coli MG1655, grown at different dilution rates (0.025, 0.05, 0.1, and 0.3 h(-1)) and on different carbon substrates (glucose, glycerol, and acetate). Proper estimation of the ATP stoichiometry requires proper information on the biomass composition of the organism as well as accurate assessment of net conversion rates under well-defined conditions. For this purpose a growth rate dependent biomass composition was derived, based on measurements and literature data. After incorporation of the growth rate dependent biomass composition in a metabolic network model, an effective P/O ratio of 1.49 +/- 0.26 mol of ATP/mol of O, K(X) (growth dependent maintenance) of 0.46 +/- 0.27 mol of ATP/C-mol of biomass and m(ATP) (growth independent maintenance) of 0.075 +/- 0.015 mol of ATP/C-mol of biomass/h were estimated using a newly developed Comprehensive Data Reconciliation (CDR) method, assuming that the three energetic parameters were independent of the growth rate and the used substrate. The resulting metabolic network model only requires the specific rate of growth, micro, as an input in order to accurately predict all other fluxes and yields.

Metabolome, transcriptome and metabolic flux analysis of arabinose fermentation by engineered Saccharomyces cerevisiae

One of the challenges in strain improvement by evolutionary engineering is to subsequently determine the molecular basis of the improved properties that were enriched from the natural genetic variation during the selective conditions. This study focuses on Saccharomyces cerevisiae IMS0002 which, after metabolic and evolutionary engineering, ferments the pentose sugar arabinose. Glucose- and arabinose-limited anaerobic chemostat cultures of IMS0002 and its non-evolved ancestor were subjected to transcriptome analysis, intracellular metabolite measurements and metabolic flux analysis. Increased expression of the GAL-regulon and deletion of GAL2 in IMS0002 confirmed that the galactose transporter is essential for growth on arabinose. Elevated intracellular concentrations of pentose-phosphate-pathway intermediates and upregulation of TKL2 and YGR043c (encoding transketolase and transaldolase isoenzymes) suggested an involvement of these genes in flux-controlling reactions in arabinose fermentation. Indeed, deletion of these genes in IMS0002 caused a 21% reduction of the maximum specific growth rate on arabinose.

Aerobic batch cultivation in micro bioreactor with integrated electrochemical sensor array

Aerobic batch cultivations of Candida utilis were carried out in two micro bioreactors with a working volume of 100 muL operated in parallel. The dimensions of the micro bioreactors were similar as the wells in a 96-well microtiter plate, to preserve compatibility with the current high-throughput cultivation systems. Each micro bioreactor was equipped with an electrochemical sensor array for the online measurement of temperature, pH, dissolved oxygen, and viable biomass concentration. Furthermore, the CO(2) production rate was obtained from the online measurement of cumulative CO(2) production during the cultivation. The online data obtained by the sensor array and the CO(2) production measurements appeared to be very reproducible for all batch cultivations performed and were highly comparable to measurement results obtained during a similar aerobic batch cultivation carried out in a conventional 4L bench-scale bioreactor. Although the sensor chip certainly needs further improvement on some points, this work clearly shows the applicability of electrochemical sensor arrays for the monitoring of parallel micro-scale fermentations, e.g. using the 96-well microtiterplate format.

Catching prompt metabolite dynamics in Escherichia coli with the BioScope at oxygen rich conditions

The design and application of a BioScope, a mini plug-flow reactor for carrying out pulse response experiments, specifically designed for Escherichia coli is presented. Main differences with the previous design are an increased volume-specific membrane surface for oxygen transfer and significantly decreased sampling intervals. The characteristics of the new device (pressure drop, residence time distribution, plug-flow behavior and O2 mass transfer) were determined and evaluated. Subsequently, 2.8 mM glucose perturbation experiments on glucose-limited aerobic E. coli chemostat cultures were carried out directly in the chemostat as well as in the BioScope (for two time frames: 8 and 40 s). It was ensured that fully aerobic conditions were maintained during the perturbation experiments. To avoid metabolite leakage during quenching, metabolite quantification (glycolytic and TCA-cycle intermediates and nucleotides) was carried out with a differential method, whereby the amounts measured in the filtrate were subtracted from the amounts measured in total broth. The dynamic metabolite profiles obtained from the BioScope perturbations were very comparable with the profiles obtained from the chemostat perturbation. This agreement demonstrates that the BioScope is a promising device for studying in vivo kinetics in E. coli that shows much faster response (< 10 s) in comparison with eukaryotes.