Search

Soil microorganisms mineralize lignin-derived aromatic carbon sources using oxidative catabolic pathways, such as the β-ketoadipate pathway. Although this aromatic pathway is one of the best-studied pathways in biochemistry, the complete pathway, including its regulation by aromatic carbon sources, has not been integrated into the metabolic network. In particular, information about the in vivo operation (e.g., kinetics and flux capacity) of the pathway is lacking. In this contribution, we use kinetic modeling and thermodynamic analysis to evaluate the in vivo operation of this key aromatic multi-step pathway. The resulting ab initio deterministic model of benzoate degradation via the β-ketoadipate (ortho-cleavage) pathway in Pseudomonas putida KT2440 is presented. The kinetic model includes mechanistic rate expressions for the enzymes and transport processes. The design and experimental validation of the model are driven by data generated from short-term perturbation experiments in a benzoate-limited continuous culture. The results of rigorous modeling of the in vivo dynamics provide strong support for flux regulation by the benzoate transporter and the enzymes forming and cleaving catechol. Revisiting the β-ketoadipate pathway might be valuable for applications in different fields, such as biochemistry and metabolic engineering, that use lignin monomers as a carbon source., Highlights • We describe a kinetic model for the β-ketoadipate pathway. • Short term metabolic responses were tracked on metabolite level by rapid sampling. • The model captures steady state and dynamic conditions of the β-ketoadipate pathway. • Thermodynamic analysis revealed regulation points of the pathway. • The results are discussed in the context of metabolic network operation.

The concepts called STED/RESOLFT superresolve features by a light-driven transfer of closely packed molecules between two different states, typically a nonfluorescent "off" state and a fluorescent "on" state at well-defined coordinates on subdiffraction scales. For this, the applied light intensity must be sufficient to guarantee the state difference for molecules spaced at the resolution sought. Relatively high intensities have therefore been applied throughout the imaging to obtain the highest resolutions. At regions where features are far enough apart that molecules could be separated with lower intensity, the excess intensity just adds to photobleaching. Here, we introduce DyMIN (standing for Dynamic Intensity Minimum) scanning, generalizing and expanding on earlier concepts of RESCue and MINFIELD to reduce sample exposure. The principle of DyMIN is that it only uses as much on/off-switching light as needed to image at the desired resolution. Fluorescence can be recorded at those positions where fluorophores are found within a subresolution neighborhood. By tuning the intensity (and thus resolution) during the acquisition of each pixel/voxel, we match the size of this neighborhood to the structures being imaged. DyMIN is shown to lower the dose of STED light on the scanned region up to ∼20-fold under common biological imaging conditions, and >100-fold for sparser 2D and 3D samples. The bleaching reduction can be converted into accordingly brighter images at

The intra- and extracellular concentrations of 16 metabolites were measured in chemostat (D = 0.1 h−1) anaerobic cultures of the yeast Saccharomyces cerevisiae CEN.PK-113-7D growing on minimal medium. Two independent sampling workflows were employed: (i) conventional cold methanol quenching and (ii) a differential approach. Metabolites were quantified in different sample fractions (total, extracellular, quenching supernatant, methanol/water extract and pellet) in order to derive their mass balance. The differential method in combination with absolute metabolite quantification by gas-chromatography with isotope dilution mass spectrometry (GC–IDMS) was used as a benchmark to assess quality of the cold methanol quenching procedure. Quantitative comparison of metabolite concentrations in all fractions collected by different quenching techniques indicates asystematic loss of the total mass of various metabolites in course of the cold methanol quenching. Pellet resulting from the cold methanol quenching besides biomass contains considerable amounts of precipitated inorganic salts from the fermentation media. Quantitative analysis has revealed significant co-precipitation of polar extracellular metabolites together with these salts. This phenomenon is especially significant for metabolites with large extracellular mass-fraction. We report that the co-precipitation is a hitherto neglected phenomenon and concluded that its degree strongly linked to culturing conditions (i.e. media composition) and chemical properties of the particular metabolite. Thus, intracellular metabolite levels measured from samples collected by cold methanol quenching might be uncertain and variably biased due to corruption by described phenomena.

The super-resolution microscopy called RESOLFT relying on fluorophore switching between longlived states, stands out by its coordinate-targeted sequential sample interrogation using low light levels. While RESOLFT has been shown to discern nanostructures in living cells, the reversibly photoswitchable green fluorescent protein (rsEGFP) employed in these experiments was switched rather slowly and recording lasted tens of minutes. We now report on the generation of rsEGFP2 providing faster switching and the use of this protein to demonstrate 25-250 times faster recordings.DOI:http://dx.doi.org/10.7554/eLife.00248.001. peerReviewed

We investigate the cooperative effect of molecular tilt and defocus on fluorophore localization by centroid calculation in far-field superresolution microscopy based on stochastic single molecule switching. If tilt angle and defocus are unknown, the localization contains systematic errors up to about ±125 nm. When imaging rotation-impaired fluorophores of unknown random orientation, the average localization accuracy in three-dimensional samples is typically limited to about ±32 nm, restricting the attainable resolution accordingly.

In biotechnology, the heterologous biosynthesis of isoprenoid compounds in Escherichia coli is a field of great interest and growth. In order to achieve higher isoprenoid yields in heterologous E. coli strains, it is necessary to quantify the pathway intermediates and adjust gene expression. In this study, we developed a precise and sensitive nonradioactive method for the simultaneous quantification of the isoprenoid precursors farnesyl diphosphate (FPP) and geranylgeranyl diphosphate (GGPP) in recombinant and wild-type E. coli cells. The method is based on the dephosphorylation of FPP and GGPP into the respective alcohols and involves their in situ extraction followed by separation and detection using gas chromatography-mass spectrometry. The integration of a geranylgeranyl diphosphate synthase gene into the E. coli chromosome leads to the accumulation of GGPP, generating quantities as high as those achieved with a multicopy expression vector.

The biosynthesis of natural products in a fast growing and easy to manipulate heterologous host system, such as Escherichia coli, is of increasing interest in biotechnology. This procedure allows the investigation of complex natural product biosynthesis and facilitates the engineering of pathways. Here we describe the cloning and the heterologous expression of tocochromanol (vitamin E) biosynthesis genes in E. coli. Tocochromanols are synthesized solely in photosynthetic organisms (cyanobacteria, algae, and higher green plants). For recombinant tocochromanol biosynthesis, the genes encoding hydroxyphenylpyruvate dioxygenase (hpd), geranylgeranylpyrophosphate synthase (crtE), geranylgeranylpyrophosphate reductase (ggh), homogentisate phytyltransferase (hpt), and tocopherol-cyclase (cyc) were cloned in a stepwise fashion and expressed in E. coli. Recombinant E. coli cells were cultivated and analyzed for tocochromanol compounds and their biosynthesis precursors. The expression of only hpd from Pseudomonas putida or crtE from Pantoea ananatis resulted in the accumulation of 336 mg L(-1) homogentisate and 84 microg L(-1) geranylgeranylpyrophosphate in E. coli cultures. Simultaneous expression of hpd, crtE, and hpt from Synechocystis sp. under the control of single tac-promoter resulted in the production of methyl-6-geranylgeranyl-benzoquinol (67.9 microg g(-1)). Additional expression of the tocopherol cyclase gene vte1 from Arabidopsis thaliana resulted in the novel formation of a vitamin E compound-delta-tocotrienol (15 microg g(-1))-in E. coli.

A method is proposed for identification of kinetic parameters when diffusion of substrates is limiting in reactions catalyzed by immobilized enzymes. This method overcomes conventional sequential procedures, which assume immobilization does not affect the conformation of the enzyme and, thus, consider intrinsic and inherent kinetics to be the same. The coupled equations describing intraparticle mass transport are solved simultaneously using numerical methods and are used for direct estimation of kinetic parameters by fitting modeling results to time-course measurements in a stirred tank reactor. While most traditional procedures were based on Michaelis-Menten kinetics, the method presented here is applicable to more complex kinetic mechanisms involving multiple state variables, such as ping-pong bi-bi. The method is applied to the kinetic resolution of (R/S)-1-methoxy-2-propanol with vinyl acetate catalyzed by Candida antarctica lipase B. A mathematical model is developed consisting of irreversible ping-pong bi-bi kinetics, including competitive inhibition of both enantiomers. The kinetic model, which fits to experimental data over a wide range of both substrates (5-95%) and temperatures (5-56 degrees C), is used for simulations to study typical behavior of immobilized enzyme systems.

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.

Glucose transport in Saccharomyces cerevisiae relies on a multi-factorial uptake system. The modulation of its efficiency depends on the differential expression of various sets of hexose transport-related proteins whose glucose affinity differs considerably. The expression of three different glucose transport proteins (HXT1, HXT5 and HXT6/7 with low-, intermediate- and high-affinity, respectively) was monitored as a result of modified extracellular glucose concentrations. Cultivation at glucose-limited (continuous) conditions was instantly replaced by a batch (and thus, non-limited) mode. Further, to mimic concentration gradients in large-scale production bioreactors, multiple and rapid transient glucose pulses were applied to chemostat cultivation. Antibodies against the HXT-proteins were used to monitor the proteins' expression levels prior to and after perturbing the external glucose concentrations. HXT5 and HXT6/7 were either expressed during the starvation-like steady-state phases in the chemostat cultivations, whereas HXT1 could not be detected at all. HXT1, however, is subsequently expressed during the excess of glucose in the batch mode, while the HXT5 and HXT6/7 transporters were at least found to decline. These findings coincide well with the transporters' affinity profiles. As a result of repeated and rapid transient glucose pulses during continuous fermentation, especially HXT6/7 pointed out to alter the protein expression pattern.

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.

Immunoaffinity enrichment of proteotypic peptides, coupled with selected reaction monitoring, enables indirect protein quantification. However the lack of suitable antibodies limits its widespread application. We developed a method in which multi-specific antibodies are used to enrich groups of peptides, thus facilitating multiplexed quantitative protein assays. We tested this strategy in a pharmacokinetic experiment by targeting a group of homologous drug transforming proteins in human hepatocytes. Our results indicate the generic applicability of this method to any biological system.

As part of a research project focusing on the development of a sustainable biocatalytic process for production of chiral secondary alcohols, the pervaporative separation of ethanol from ethanol/ethyl acetate/1-methoxy-2-propanol/1-methoxy-2-propyl acetate-mixtures through a commercial PVA-based membrane was investigated. Separation behavior of this mixture was studied in a range of mol fractions (10–70%), temperatures relevant for biocatalytic conversions (35–55 °C) and downstream pressures (35–200 mbar). Pervaporation of the non-diluted multicomponent mixture was shown to be strongly influenced by interactions between the permeants and the membrane. Investigation of these interactions contributed to the understanding of the mass transport mechanism of this mixture. Overall, high fluxes were obtained, but small differences between the fastest permeating species were found. The fastest permeating species was ethanol, ethyl acetate or 1-methoxy-2-propanol depending on the feed composition.

This paper examines the validity of the linlog approach, which was recently developed in our laboratory, by comparison of two different kinetic models for the metabolic network of Escherichia coli. The first model is a complete mechanistic model; the second is an approximative model in which linlog kinetics are applied. The parameters of the linlog model (elasticities) are derived from the mechanistic model. Three different optimization cases are examined. In all cases, the objective is to calculate the enzyme levels that maximize a certain flux while keeping the total amount of enzyme constant and preventing large changes of metabolite concentrations. For an average variation of metabolite levels of 10% and individual changes of a factor 2, the predicted enzyme levels, metabolite concentrations and fluxes of both models are highly similar. This similarity holds for changes in enzyme level of a factor 4-6 and for changes in fluxes up to a factor 6. In all three cases, the predicted optimal enzyme levels could neither have been found by intuition-based approaches, nor on basis of flux control coefficients. This demonstrates that kinetic models are essential tools in Metabolic Engineering. In this respect, the linlog approach is a valuable extension of MCA, since it allows construction of kinetic models, based on MCA parameters, that can be used for constrained optimization problems and are valid for large changes of metabolite and enzyme levels.

We have investigated the role and the kinetic properties of the Hxt5 glucose transporter of Saccharomyces cerevisiae. The HXT5 gene was not expressed during growth of the yeast cells in rich medium with glucose or raffinose. However, it became strongly induced during nitrogen or carbon starvation. We have constructed yeast strains constitutively expressing only Hxt5, Hxt1 (low affinity) or Hxt7 (high affinity), but no other glucose transporters. Aerobic fed-batch cultures at quasi steady-state conditions, and aerobic and anaerobic chemostat cultures at steady-state conditions of these strains were used for estimation of the kinetic properties of the individual transporters under in vivo conditions, by investigating the dynamic responses of the strains to changes in extracellular glucose concentration. The Km value and the growth properties of the HXT5 single expression strain indicate that Hxt5 is a transporter with intermediate affinity.

Glucose fatty acid monoesters were produced in a stirred-tank reactor in combination with a membrane-pervaporation separation unit for the regeneration of the solvent ethyl methylketone (EMK) on a preparative scale. Yields up to 79% (169 g product) were achieved in the synthesis of 6-O-stearoyl-D-glucose in a mainly solid phase system in the presence of a small amount of EMK maintaining a catalytic liquid phase as well as forming a suitable azeotropic mixture with the reaction water. High residual activity was found using lipase B from Candida antarctica under reduced pressure at 60 °C. EMK could be easily removed from the reaction mixture and is regarded as biocompatible in the preparation of food additives. Des monoesters d'acides gras de glucose sont produits dans un reacteur agite en combinaison avec une unite de separation par pervaporation pour regenerer le solvant ethylmethylcetone (EMC) a une echelle preparative. Des rendements jusqu'a 79% (169 g de produit) sont obtenus pour la synthese du 6-O-stearoyl-D-glucose en presence d'une petite quantite de solvant EMC biocompatible et utilisable dans l'industrie alimentaire.

Escherichia coli fed-batch cultivations at 22 m3 scale were compared to corresponding laboratory scale processes and cultivations using a scale-down reactor furnished with a high-glucose concentration zone to mimic the conditions in a feed zone of the large bioreactor. Formate accumulated in the large reactor, indicating the existence of oxygen limitation zones. It is suggested that the reduced biomass yield at large scale partly is due to repeated production/re-assimilation of acetate from overflow metabolism and mixed acid fermentation products due to local moving zones with oxygen limitation. The conditions that generated mixed-acid fermentation in the scale-down reactor also induced a number of stress responses, monitored by analysis of mRNA of selected stress induced genes. The stress responses were relaxed when the cells returned to the substrate limited and oxygen sufficient compartment of the reactor. Corresponding analysis in the large reactor showed that the concentration of mRNA of four stress induced genes was lowest at the sampling port most distant from the feed zone. It is assumed that repeated induction/relaxation of stress responses in a large bioreactor may contribute to altered physiological properties of the cells grown in large-scale bioreactor. Flow cytometric analysis revealed reduced damage with respect to cytoplasmic membrane potential and integrity in cells grown in the dynamic environments of the large scale reactor and the scale-down reactor.

6-O-β-D(+)-Glucose stearate serving as a model product of glucose fatty acid monoesters was synthesized using lipase B from Candiada antarctica (Chirazyme® L-2) in a mainly solid-phase system in a stirred-tank membrane reactor. Esterification was performed in the presence of a small amount of ethyl methylketone (EMK), maintaining a catalytic liquid phase as well as forming an azeotrope with the reaction water. The azeotrope was evaporated and broken by membrane vapor permeation, then the dried EMK was returned to the reaction medium. The process was optimized by response surface methodology based on five major reaction parameters [time (T), substrate ratio (acyl donor to glucose, S r, temperature (R t), amount of solvent [solvent to substrates, [w/w], S a), and enzyme load (E l)] varied at three levels, resulting in higher yields. Thus, under optimized conditions [T r=58 h; S r=2.7; E l=8.9% (w/w); R t=78°C; S a=1.9], up to 93% yields of glucose stearate were achieved.

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.

In order to produce sophorolipids from whey, thereby lowering the lactose content and biological oxygen demand, a two-step batch cultivation process was developed including medium sterilization by filtration. In the first step, whey was sterilized by a combination of crossflow and sterile filtration. Because the sophorolipid-producing yeast Candida bombicola ATCC 22214 was not able to use lactose as a carbon source directly, the oleaginous yeast Cryptococcus curvatus ATCC 20509 was grown on deproteinized whey concentrates (DWC). With 1:1 diluted DWC-20, lactose was consumed as the carbon source and biomass (24 g/l dry weight content) as well as single-cell oil (SCO, 10 g/l) were produced. The cultivation broth was disrupted with a glass bead mill and it served as medium for growth (29 g cell dry mass/l) and sophorolipid production (12 g/l) of the yeast C. bombicola.

An ion-pair reversed-phase chromatography method was adapted for the simultaneous separation and quantification of adenine and pyridine nucleotide concentrations in cell extracts of Saccharomyces cerevisiae. Microbial extracts including metabolites, macromolecular constituents and inorganic compounds were loaded onto a ODS pre-column in the presence of triethylamine phosphate (TEA-Pi) resulting in a selective binding of the nucleotides and removing of interfering compounds. After washing the enrichment cartridge with 30 mM TEA-Pi buffer, adenine and pyridine nucleotides were eluted with a gradient of Mg(II), the competing hetaeron. This combined cleaning and concentration step leads to remarkable improvement of the detection limit for all nucleotides of interest and column lifetimes. The clean up and separation procedures were performed automatically with a micro robotic-system and a control software package written in PASCAL. The paper reports about the application of the proposed method to separation of adenine and pyridine nucleotides in cells extracts of S. cerevisiae grown anaerobically in a continuous culture (D = 0.1 h-1). Rapidity of analysis, high sensitivity as well as reproducibility of the system and the accurate evaluation of the adenine and pyridine nucleotide concentrations make this method particularly useful for routine analysis.

The explosive 2,4,6-trinitrotoluene (TNT), found as a major contaminant at armament plants from the two world wars, is reduced by a variety of microorganisms when electron donors such as glucose are added. This study shows that the cometabolic reduction of TNT to 2,4,6-triaminotoluene by an undefined anaerobic consortium increased considerably with increasing TNT concentrations and decreased with decreasing concentrations and feeding rates of glucose. The interactions of TNT and its reduction products with montmorillonitic clay and humic acids were investigated in abiotic adsorption experiments and during the microbial reduction of TNT. The results indicate that reduction products of TNT particularly hydroxylaminodinitrotoluenes and 2,4,6-triaminotoluene bind irreversibly to soil components, which would prevent or prolong mineralization of the contaminants. Irreversible binding also hinders a further spread of the contaminants through soil or leaching into the groundwater.

High concentrations of 422 g sophorolipids l−1 were produced using a two-stage cultivation process: first deproteinized whey concentrate (DWC) containing 110 g lactose l−1 was used for cultivation of the yeast Cryptococcus curvatus ATCC 20509, resulting in 34 g dry weight l−1, 20 g single-cell oil l−1 and reducing the chemical oxygen demand (COD) from 159 g l−1 to 35 g oxygen l−1. Afterwards cells were disrupted by passing the cell suspension directly through a high pressure laboratory homogeniser. After autoclavation, the resulting crude cell extract containing the single-cell oil served as substrate for growth of Candida bombicola ATCC 22214 and for sophorolipid production in a second stage. When the single-cell oil was consumed, repeated feeding of 400 g rapeseed oil l−1 was started increasing the yield of sophorolipids to 422 g l−1. A simple technique for product isolation, sedimentation, could be used to harvest the crude sophorolipids. © Rapid Science Ltd. 1998

Experimental data for the transient behavior of Pseudomonas putida grown on phenol in continuous culture after a shift in substrate inlet concentration are presented. Representing these effects by a mathematical model requires the introduction of a model structure including a time delay. The hypothesis of ribosomes as the growth governing compound leads to a closer inspection of the assembly of ribosomes. The `maturation'of complex structures like ribosomes corresponds to a discrete time delay in a mathematical model. The hypothesis of ribosomal assembly being viewed as a polymerization reaction leads to a system of coupled differential equations. This mathematical structure allows the modeling of distributed discrete time delays. Compatibility between dynamic and steady state kinetics must be considered; which leads to restrictions with respect to model structure. The resulting model, which simplifies the biochemical reaction networks in ribosomal synthesis by mathematical formulation of a linear polymerization, is able to predict process dynamics

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.

A continuously stirred tank reactor equipped with a microfiltration membrane unit for effective retention of solid enzyme particles in organic solvent systems was constructed. A microfiltration membrane with pore size 0.1 μm based on a polypropylene matrice was used for enzyme retention and showed no loss of performance or denaturation by the organic solvent. The stereoselective transesterification reaction of (±)-trans-sobrerol with vinyl acetate as both acyl donor and solvent was chosen as a model system. Lipase PS from Pseudomonas sp. was used in its native form without further modification steps and remained active and highly stereoselective at 50 °C when applied to the reaction system for about 600 h of continuous operation. (--)-(1S,5R)-sobrerolacetate and (+)-(1R,5S)-sobrerol could be obtained in the reactor outlet stream at approximately 48% conversion with enantiomeric excesses eeP ≈ 92% and eeS≈ 80%, respectively. Furthermore, theoretical kinetic studies showed that limitations in enantiomeric purity due to the characteristics of an ideal mixed tank may be sufficiently overcome by designing a simple two-stage process.

We present a three-dimensional, multiscale model of vascular tumour growth, which couples nutrient/growth factor transport, blood flow, angiogenesis, vascular remodelling, movement of and interactions between normal and tumour cells, and nutrient-dependent cell cycle dynamics within each cell. We present computational simulations which show how a vascular network may evolve and interact with tumour and healthy cells. We also demonstrate how our model may be combined with experimental data, to predict the spatio-temporal evolution of a vascular tumour.

A method for experimental observation of the diffusivity of yeast cells in a suspension undergoing shear is presented. Local filtrate flux and cake thickness were measured in a radial flow chamber, which has been developed for this purpose. From the simultaneous numerical solution of the Navier-Stokes equation as well as the material balance equation, which govern momentum and material transport in the flat-channel geometry, simulated local flux and cake thickness were obtained for the same operating conditions as in the experiments. The shear-induced diffusion coefficient was then estimated after comparing the measurements with simulated values and minimizing the differences between them.

Conjugation between two Pseudomonas putida species has been studied. Special emphasis has been given to the design of tools for better defined and reproducible experimental conditions. This is an essential prerequisite for any kinetic analysis of the phenomenon. The experimental approach suggested in the paper is a three-stage continuous fermentation, where donor and recipient strains can be cultivated at predefined steady-state growth rates, and conjugation takes place at steady-state conditions in the third stage. For the first time a remarkable influence of agitative speed on conjugation has been experimentally observed.

In the field of metabolomics, GC-MS has rather established itself as a tool for semi-quantitative strategies like metabolic fingerprinting or metabolic profiling. Absolute quantification of intra- or extracellular metabolites is nowadays mostly accomplished by application of diverse LC-MS techniques. Only few groups have so far adopted GC-MS technology for this exceptionally challenging task. Besides numerous and deeply investigated problems related to sample generation, the pronounced matrix effects in biological samples have led to the almost mandatory application of isotope dilution mass spectrometry (IDMS) for the accurate determination of absolute metabolite concentrations. Nevertheless, access to stable isotope labeled internal standards (ILIS), which are in many cases commercially unavailable, is quite laborious and very expensive. Here we present an improved and simplified gas chromatography-isotope dilution mass spectrometry (GC-IDMS) protocol for the absolute determination of intra- and extracellular metabolite levels. Commercially available (13)C-labeled algal cells were used as a convenient source for the preparation of internal standards. Advantages as well as limitations of the described method are discussed.

Transition from glucose limitation to glucose excess in Saccharomyces cerevisiae leads to a fast change from complete to partial oxidative metabolism. To determine rapid changes of in vivo concentrations of yeast metabolites, occurring in the range of seconds after a glucose injection, a fast sampling technique was developed. A harvesting device with negligible dead space and highly efficient inactivation and extraction steps was applied to analyze the behavior of the physiological concentrations of adenine nucleotides as the dynamic response to fast changing glucose concentrations. Adenine nucleotides were selected as appropriate references due to their high turnover rates (in the range of seconds). During the first 5 s after a glucose pulse to a continuous chemostat culture, a significant decrease (30%) in ATP concentration and a coincident rise of ADP and AMP values were found. A further advantage of this technique is the high sampling frequency (5 s) and the long-term sterility ensuring recurrent harvesting.

Background In this paper we apply a novel agent-based simulation method in order to model intracellular reactions in detail. The simulations are performed within a virtual cytoskeleton enriched with further crowding elements, which allows the analysis of molecular crowding effects on intracellular diffusion and reaction rates. The cytoskeleton network leads to a reduction in the mobility of molecules. Molecules can also unspecifically bind to membranes or the cytoskeleton affecting (i) the fraction of unbound molecules in the cytosol and (ii) furthermore reducing the mobility. Binding of molecules to intracellular structures or scaffolds can in turn lead to a microcompartmentalization of the cell. Especially the formation of enzyme complexes promoting metabolic channeling, e.g. in glycolysis, depends on the co-localization of the proteins. Results While the co-localization of enzymes leads to faster reaction rates, the reduced mobility decreases the collision rate of reactants, hence reducing the reaction rate, as expected. This effect is most prominent in diffusion limited reactions. Furthermore, anomalous diffusion can occur due to molecular crowding in the cell. In the context of diffusion controlled reactions, anomalous diffusion leads to fractal reaction kinetics. The simulation framework is used to quantify and separate the effects originating from molecular crowding or the reduced mobility of the reactants. We were able to define three factors which describe the effective reaction rate, namely f diff for the diffusion effect, f volume for the crowding, and f access for the reduced accessibility of the molecules. Conclusions Molecule distributions, reaction rate constants and structural parameters can be adjusted separately in the simulation allowing a comprehensive study of individual effects in the context of a realistic cell environment. As such, the present simulation can help to bridge the gap between in vivo and in vitro kinetics., BMC Systems Biology, 5, ISSN:1752-0509

Agent-based models are rigorous tools for simulating the interactions of individual entities, such as organisms or molecules within cells and assessing their effects on the dynamic behavior of the system as a whole. In context with biopro- cess and biosystems engineering there are several interesting and important appli- cations. This contribution aims at introducing this strategy with the aid of two examples characterized by striking distinctions in the scale of the individual entities and the mode of their interactions. In the first example a structured-segregated model is applied to travel along the lifelines of single cells in the environment of a three-dimensional turbulent field of a stirred bioreactor. The modeling approach is based on an Euler-Lagrange formulation of the system. The strategy permits one to account for the heterogeneity present in real reactors in both the fluid and cellular phases, respectively. The individual response of the cells to local variations in the extracellular concentrations is pictured by a dynamically structured model of the key reactions of the central metabolism. The approach permits analysis of the lifelines of individual cells in space and time. The second application of the individual modeling approach deals with dynamic modeling of signal transduction pathways in individual cells. Usually signal trans- duction networks are portrayed as being wired together in a spatially defined manner. Living circuitry, however, is placed in highly malleable internal architec- ture. Creating a homogenous bag of molecules, a well-mixed system, the dynamic behavior of which is modeled with a set of ordinary differential equations is normally not valid. The dynamics of the MAP kinase and a steroid hormone pathway serve as examples to illustrate how single molecule tracking can be linked with the stochasticity of biochemical reactions, where diffusion and reaction occur in a probabilistic manner. The problem of hindered diffusion caused by macromo- lecular crowding is also taken into account.

The majority of dynamic gene regulatory network (GRN) models are comprised of only a few genes and do not take multiple transcription regulation into account. The models are conceived in this way in order to minimize the number of kinetic parameters. In this paper, we propose a new approach for predicting kinetic parameters from DNA-binding site sequences by correlating the protein–DNA-binding affinities with nucleotide sequence conservation. We present the dynamic modeling of the cra modulon transcription in Escherichia coli during glucose-limited fed-batch cultivation. The concentration of the Cra regulator protein inhibitor, fructose 1,6-bis(phosphate), decreases sharply, eventually leading to the repression of transcription. Total RNA concentration data indicate a strong regulation of transcription through the availability of RNA polymerase. A critical assessment of the results of the model simulations supports this finding. This new approach for the prediction of transcription dynamics may improve the metabolic engineering of gene regulation processes.

Metabolomics emerges to become an important profiling technique in bio(techno)logical systems. In addition to intracellular metabolite concentrations at (quasi) stationary conditions, stimulus-response experiments provide information on the dynamic behavior of metabolic pathways. These data are relevant for bioprocess analysis on the level of metabolism and for application of metabolic engineering principles aiming at a metabolic redesign of producer cells. However, even for the well-studied bacteria Escherichia coli only limited growth-rate dependent intracellular metabolite information is currently available, thereby impeding comprehensive metabolome analysis. Here, we present intracellular metabolite concentration data of representative glycolytic intermediates in E. coli cultivated in glucose-limited chemostats, (i) at systematic variation of growth-rate (D = 0.1, 0.2, 0.3, and 0.4 h(-1)) and (ii) at both steady-state and after a glucose pulse applying a recently introduced integrated sampling procedure and LC-MS analytical method. Whereas intracellular steady-state concentrations of upper part glycolytic intermediates FBP and DHAP+GAP increased 2.3-fold, respectively 2.8-fold, when specific growth-rate is raised from micro = 0.1 h(-1) to micro = 0.4 h(-1), the opposite trend was observed for 2PG+3PG and PEP pools with a decrease by a factor of 2.1, respectively 1.9. In glucose pulse experiments FBP and DHAP+GAP showed a 3.3 (1.8)-fold, respectively 2.8 (2.0)-fold, increase relative to the steady-state level at micro = 0.1 (0.4) h(-1). Also, the dynamics changed with growth-rate for these two metabolite pools. In contrast, 2PG+3PG and PEP were characterized by decreased concentrations in response to a glucose pulse and the relative changes related to steady-state values were significantly smaller compared with FBP and DHAP+GAP. The observed growth-rate dependency of our data clearly indicates the necessity for metabolome studies covering a broader range of physiological growth conditions.

The present work is the first to deal with the determination of cholesterol synthesis rates in primary rat hepatocytes using transient (13)C-flux analysis. The effects of statins on cholesterol biosynthesis and central carbon fluxes were quantified at a therapeutic concentration of 50 nM atorvastatin using carbon-labeled glutamine. The flux through the cholesterol pathway decreased from 0.27 to 0.08 mmol/l(cv)h in response to the administration of the hypolipidemic drug. Isotopic steady state was reached within 4h in the central carbon metabolism but not in the cholesterol pathway, regardless of whether atorvastatin was administered or not. Marked channeling was observed for the symmetrical tricarboxylic acid cycle intermediates, succinate and fumarate. Non-stationary (13)C-based flux identification delivers both intracellular fluxes and intermediate levels, which was for the first time utilized for investigating systems-level effects of the administered drug by quantifying the flux control of the 3-hydroxy-3-methylglutaryl-coenzyme A reductase.

Coupling complex regulatory and metabolic networks for the purpose of dynamic modeling requires knowledge of the quantitative kinetics of the participating reactions as well as the variation of parameters in the context of the physiological state of the system. This chapter aims at demonstrating the integration of the different networks for E. coli exposed to an increasing carbon limitation of a fed-batch process with constant feeding of the carbon and energy source glucose. Starting from a global observation of the response of the bacteria in terms of flux distribution and gene expression in the central metabolism, emphasis is given to the dynamic modeling of regulation phenomena in the catabolism. The cra regulon which is linked to the dynamic response of the metabolite fructose 1,6- bis(phosphate) serves as an example to introduce a new concept, in which the binding constants are estimated from DNA-binding site sequences of the regulatory proteins. By comparison of the nucleotide frequencies within the DNA-binding sites for the individual target genes of the regulon, it is possible to perform a reasonable estimation of the kinetic parameters. Results of these estimations are compared with experimentally observed transcript concentrations measured with the aid of quantitative PCR. In addition it is shown how these outputs of the regulatory networks can be linked to the maximal rates of the enzymes for the metabolic system of interest. The discussion of this issue is embedded within a critical assessment of different conceptual frameworks for modeling the metabolic network, which covers the spectrum of dynamic modeling at different levels of complexity, such as genome scale, modular approaches and reduced models.

The intracellular alarmone guanosine 3',5'-bis(diphosphate) (ppGpp) has been thoroughly investigated over the past 40 years and has become one of the best-known effectors in bacterial physiology. ppGpp is also of great importance for biotechnological applications. Systems biology research, involving experimental and mathematical approaches, has contributed a great deal to uncovering the alarmone's complex regulatory effects. HPLC analysis and UV detection are used to quantify intracellular ppGpp. The samples analyzed also contain other phosphorylated guanine nucleotides and, therefore, are spiked with a standard ppGpp solution. A rapidly growing number of laboratories are turning to synthesizing the nucleotide in vitro involving time-consuming protocols and yielding only low amounts of ppGpp. The current article provides a protocol for the preparation of large quantities of a ribosome extract that contains high ppGpp synthesis activity. The demonstrated upscaling from shaking flask to bioreactor cultivation involves the continuous and refrigerated harvest of the biomass. (13)C NMR analysis enabled the structural characterization of the synthesis product and was complemented by mass spectrometry and methods that are commonly used to identify ppGpp.

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.

Sophorolipids were produced by single-step batch cultivation of Candida bombicola ATCC 22214 on deproteinized whey concentrate and repeated feed of rapeseed oil. A mild sterilization method for whey was developed. High yields of 280 g dry sophorolipids l−1 were obtained from deproteinized whey concentrate containing 100 g lactose l−1 and 300 g rapeseed oil l−1. Surprisingly, the whey lactose was not consumed by the organism. Growth only on the oil was assumed and a high lipase activity of 24 U per g cell dry weight resulted.

Metabolome analysis, the analysis of large sets of intracellular metabolites, has become an important systems analysis method in biotechnological and pharmaceutical research. In metabolic engineering, the integration of metabolome data with fluxome and proteome data into large-scale mathematical models promises to foster rational strategies for strain and cell line improvement. However, the development of reproducible sampling procedures for quantitative analysis of intracellular metabolite concentrations represents a major challenge, accomplishing (i) fast transfer of sample, (ii) efficient quenching of metabolism, (iii) quantitative metabolite extraction, and (iv) optimum sample conditioning for subsequent quantitative analysis. In addressing these requirements, we propose an integrated sampling procedure. Simultaneous quenching and quantitative extraction of intracellular metabolites were realized by short-time exposure of cells to temperatures < or =95 degrees C, where intracellular metabolites are released quantitatively. Based on these findings, we combined principles of heat transfer with knowledge on physiology, for example, turnover rates of energy metabolites, to develop an optimized sampling procedure based on a coiled single tube heat exchanger. As a result, this sampling procedure enables reliable and reproducible measurements through (i) the integration of three unit operations into a one unit operation, (ii) the avoidance of any alteration of the sample due to chemical reagents in quenching and extraction, and (iii) automation. A sampling frequency of 5 s(-)(1) and an overall individual sample processing time faster than 30 s allow observing responses of intracellular metabolite concentrations to extracellular stimuli on a subsecond time scale. Recovery and reliability of the unit operations were analyzed. Impact of sample conditioning on subsequent IC-MS analysis of metabolites was examined as well. The integrated sampling procedure was validated through consistent results from steady-state metabolite analysis of Escherichia coli cultivated in a chemostat at D = 0.1 h(-)(1).

Deproteinized sweet and sour cheese whey concentrates were investigated for their suitability as substrates for the production of single-cell protein with Kluyveromyces marxianus CBS 6556 up to a 100-l scale. An important factor for gaining high cell concentrations was the use of the Crabtree-negative strain K. marxianus CBS 6556. Supplements such as trace elements, ammonium and calcium were required for the complete conversion of sweet whey concentrates into biomass, whereas sour whey concentrates had to be supplemented with ammonium, trace elements and vitamins. After improvement, biomass dry concentrations of up to 50 g l−1 could be reached with Yx/s values of 0.52 for sweet whey and of up to 65 g l−1 with Yx/s values of 0.48 for sour whey concentrates. The chemical oxygen demand of the whey concentrates were reduced by 80%. The cells were used for the analysis of amino acid and ash composition, showing a distinct increase of eight out of ten essential amino acids compared to sweet and sour whey protein and exceeding the World Health Organisation guidelines for valine, leucine, isoleucine, threonine, phenylalanine and tyrosine.

The Candida antarctica lipase B catalyzed kinetic resolution of (R/S)-1-methoxy-2-propyl-acetate was studied as a model system for the biocatalytic production of chiral secondary alcohols. For this purpose, a kinetic model is proposed involving both enantiomers of this reaction using model discrimination and parameter identification. Starting from a ping-pong bi-bi mechanism, a simplified model with sensitive parameters was derived for the R- and S-enantiomer, respectively. It was validated at pH 7.0, using time-course measurements at varying temperatures (30-60 degrees C) and initial substrate conditions (0.05-1.5 M). This model was then used for mechanistic interpretation of the kinetic resolution on a biochemical level. The effect of temperature on kinetic parameters and enantiomeric ratio was investigated and compared to findings from the field of molecular modeling to obtain a better understanding of the reaction system for process design. Values of 21.2 and 9.7 kJmol-1 were determined for the enthalpic (DeltaR-S DeltaH ++ degrees) and the entropic (-T x DeltaR-S DeltaS ++ degrees) contribution of the difference in transition state energy of both enantiomers at 30 degrees C. High enantiomeric ratio's (E of 47-110) especially at lower temperatures, in addition to enzyme activity at a wide pH range, indicate this biotransformation is a promising example for the industrial production of chiral secondary alcohols.

The presumably high potential of a holistic design approach for complex biochemical reaction networks is exemplified here for the network of tryptophan biosynthesis from glucose, a system whose components have been investigated thoroughly before. A dynamic model that combines the behavior of the trp operon gene expression with the metabolic network of central carbon metabolism and tryptophan biosynthesis is investigated. This model is analyzed in terms of metabolic fluxes, metabolic control, and nonlinear optimization. We compare two models for a wild-type strain and another model for a tryptophan producer. An integrated optimization of the whole network leads to a significant increase in tryptophan production rate for all systems under study. This enhancement is well above the increase that can be achieved by an optimization of subsystems. A constant ratio of control coefficients on tryptophan synthesis rate has been identified for the models regarding or disregarding trp operon expression. Although we found some examples where flux control coefficients even contradict the trends of enzyme activity changes in an optimized profile, flux control can be used as an indication for enzymes that have to be taken into account in optimization.

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.