Construction and elementary mode analysis of a metabolic model for Shewanella oneidensis MR-1.

A stoichiometric model describing the central metabolism of Shewanella oneidensis MR-1 wild-type and derivative strains was developed and used in elementary mode analysis (EMA). Shewanella oneidensis MR-1 can anaerobically respire a diverse pool of electron acceptors, and may be applied in several biotechnology settings, including bioremediation of toxic metals, electricity generation in microbial fuel cells, and whole-cell biocatalysis. The metabolic model presented here was adapted and verified by comparing the growth phenotypes of 13 single- and 1 double-knockout strains, while considering respiration via aerobic, anaerobic fumarate, and anaerobic metal reduction (Mtr) pathways, and utilizing acetate, n-acetylglucosamine (NAG), or lactate as carbon sources. The gene ppc, which encodes phosphoenolpyruvate carboxylase (Ppc), was determined to be necessary for aerobic growth on NAG and lactate, while not essential for growth on acetate. This suggests that Ppc is the only active anaplerotic enzyme when cultivated on lactate and NAG. The application of regulatory and substrate limitations to EMA has enabled creation of metabolic models that better reflect biological conditions, and significantly reduce the solution space for each condition, facilitating rapid strain optimization. This wild-type model can be easily adapted to include utilization of different carbon sources or secretion of different metabolic products, and allows the prediction of single- and multiple-knockout strains that are expected to operate under defined conditions with increased efficiency when compared to wild type cells.

Cell cycle-dependent protein secretion by Saccharomyces cerevisiae.

Synchronized Saccharomyces cerevisiae cell populations were used to examine secretion rates of a heterologous protein as a function of cell cycle position. The synchronization procedure had a profound effect on the type and quality of data obtained. When cell synchrony was induced by cell cycle-arresting drugs, a significant physiological perturbation of cells was observed that obscured representative secretion data. In contrast, synchronization with centrifugal elutriation resulted in synchronized first-generation daughter cells with undetectable perturbation of the physiological state. The synchronized cells did not secrete significant amounts of protein until they reached cell division, suggesting that the secretion process in these cells is strongly cell cycle dependent. However, the maximum secretion rate of the synchronized culture (7-14 molecules/cell/second) was significantly lower than that of an asynchronous culture (29-51 molecules/cell/second). This result indicates that young daughter cells isolated in the synchronization process exhibit different protein secretion behavior than older mother cells that are absent in the synchronized cell population but present in the asynchronous culture.

Glucose uptake rates of single E. coli cells grown in glucose-limited chemostat cultures.

To evaluate the extent to which single-cell glucose uptake rates determine the overall specific growth rate of a culture, dilute chemostat cultures of Escherichia coli BL21 were grown in defined medium under glucose limitation. The glucose uptake dynamics of the cell population was examined at the single-cell level using the fluorescent glucose analog, 2-NBDG. Between dilution rates of 0.12 h(-1) and 0.40 h(-1), mean cellular protein content and steady-state, extracellular glucose concentrations increased with increasing dilution rate. However, the distribution of 2-NBDG uptake rates in the population remained constant over the range of dilution rates studied. This indicates that the growth of cells in continuous culture is not limited by the maximum rate of uptake of glucose but by the availability of glucose for transport. The work represents an example of how quantitative flow cytometry can be applied to gain detailed insight into microbial growth physiology.

A sensitive method to detect initiation of growth in Streptococcus gordonii using ribosomal RNA operon-reporter gene fusions.

A system for studying the early growth response of Streptococcus gordonii to environmental stimuli has been developed. A reporter gene, encoding alpha-amylase, has been integrated into an rRNA operon to monitor changes in cellular physiology associated with the initiation of growth. Two such strains with single integrants have been characterized during the transition from lag phase to exponential growth. Synthesis of the reporter is correlated to growth initiation in both strains, and the reporter enzyme is detectable with sufficient sensitivity. Comparison of the expression profiles of the two rrn operons containing the reporter gene suggests that they are differentially expressed over the course of growth.

Peroxisomes as sites for synthesis of polyhydroxyalkanoates in transgenic plants.

Bacterial genes responsible for poly(3-hydroxybutyrate) (PHB) biosynthesis were targeted to plant peroxisomes by adding a carboxy-terminal targeting sequence. The enzymes evidently were transported into peroxisomes, retained their catalytic activity, and reacted with peroxisomally available precursors because PHB synthesis in transgenic plant cells was localized to peroxisomes. Up to 2 mg/g fresh weight PHB was produced in suspension cultures of Black Mexican Sweet maize cells after biolistic transformation with three peroxisomally targeted bacterial genes. An equilibrium effect is proposed to explain the unexpected existence of (R)-3-hydroxybutyryl-CoA in plant peroxisomes.

Production of two phase polyhydroxyalkanoic acid granules in Ralstonia eutropha.

To synthesize layered granules consisting of selected phases of polyhydroxybutyrate (PHB) homopolymer and PH(B-co-V) copolymer, Ralstonia eutropha was grown on fructose and limited quantities (1 g/l) of valeric acid. Exhaustion of the valerate resulted in a carbon source shift and a shift in the composition of polyhydroxyalkanoate (PHA) being synthesized within the cell. The synthesis rates were 0.030 g PH(B-co-V)/l per h and 0.033 g PHB/l per h, giving a copolymer composition of 48% HV. The valerate was exhausted at approximately 12 h at a rate of 0.0894 g/l per h after which only PHB was produced through the remaining 12 h at 0.033 g PHB/l per h from the remaining fructose, which was utilized at a constant rate of 0.0861 g/l per h throughout all 24 h of the experiment. Differential scanning calorimetry (DSC) of isolated granules showed two glass transitions, confirming the presence of two distinct polymer phases within the layered granules. Transmission electron microscopic images stained with RuO4 revealed a heavily stained copolymer core within a lighter stained PHB shell, confirming the expected morphology of granule composition. Thus, biosynthesis can be exploited for the control of domain sizes in layered granules, potentially providing metabolic control over the physical properties of the resultant polymer.

A flow injection flow cytometry system for on-line monitoring of bioreactors.

For direct and on-line study of the physiological states of cell cultures, a robust flow injection system has been designed and interfaced with flow cytometry (FI-FCM). The core of the flow injection system includes a microchamber designed for sample processing. The design of this microchamber allows not only an accurate on-line dilution but also on-line cell fixation, staining, and washing. The flow injection part of the system was tested by monitoring the optical density of a growing E.coli culture on-line using a spectrophotometer. The entire growth curve, from lag phase to stationary phase, was obtained with frequent sampling. The performance of the entire FI-FCM system is demonstrated in three applications. The first is the monitoring of green fluorescent protein fluorophore formation kinetics in E.coli by visualizing the fluorescence evolution after protein synthesis is inhibited. The data revealed a subpopulation of cells that do not become fluorescent. In addition, the data show that single-cell fluorescence is distributed over a wide range and that the fluorescent population contains cells that are capable of reaching significantly higher expression levels than that indicated by the population average. The second application is the detailed flow cytometric evaluation of the batch growth dynamics of E.coli expressing Gfp. The collected single-cell data visualize the batch growth phases and it is shown that a state of balanced growth is never reached by the culture. The third application is the determination of distribution of DNA content of a S. cerevisiae population by automatically staining cells using a DNA-specific stain. Reproducibility of the on-line staining reaction shows that the system is not restricted to measuring the native properties of cells; rather, a wider range of cellular components could be monitored after appropriate sample processing. The system is thus particularly useful because it operates automatically without direct operator supervision for extended time periods.

Dynamics of glucose uptake by single Escherichia coli cells.

The fluorescent glucose analog, 2-(N-(7-nitrobenz-2-oxa-1,3-diazol-4-yl)amino)-2-deoxyglucose (2-NBDG), was used to measure rates of glucose uptake by single Escherichia coli cells. When cell populations were exposed to the glucose analog, 2-NBDG was actively transported and accumulated in single cells to a steady-state level that depended upon its extracellular concentration, the glucose transport capacity of the cells, and the intracellular degradation rate. The dependence upon substrate concentration could be described according to Michaelis-Menten kinetics with apparent saturation constant KM = 1.75 microM, and maximum 2-NBDG uptake rate= 197 molecules/cell-second. Specificity of glucose transporters to the analog was confirmed by inhibition of uptake of 2-NBDG by D-glucose, 3-o-methyl glucose, and D-glucosamine, and lack of inhibition by L-glucose. Inhibition of 2-NBDG uptake by D-glucose was competitive in nature. The assay for 2-NBDG uptake is extremely sensitive such that the presence of even trace amounts of D-glucose in the culture medium (approximately 0.2 microM) is detectable. The rates of single-cell analog uptake were found to increase proportionally with cell size as measured by microscopy or single-cell light scattering intensity. The assay was used to identify and isolate mutant cells with altered glucose uptake characteristics. A mathematical model was developed to provide a theoretical basis for estimating single-cell glucose uptake rates from single-cell 2-NBDG uptake rates. The assay provides a novel means of estimating the instantaneous rates of nutrient depletion in the growth environment during a batch cultivation.

Comparison of mutant forms of the green fluorescent protein as expression markers in Chinese hamster ovary (CHO) and Saccharomyces cerevisiae cells.

Several green fluorescent protein (Gfp) mutants with increased cellular fluorescence compared to the wildtype protein have recently been generated. We have expressed and compared wildtype Gfp and mutants S65T, F100S/ M154T/V164A, F64L/S65T, and S65A/V68L/S72A under identical growth conditions in CHO and Saccharomyces cerevisiae cells. The results suggest that the last two Gfp mutants are the best candidates as reporter proteins, and they provide a high signal-to-noise ratio in both systems. Single gene copy expression of these mutant forms is easily detectable over background autofluorescence. All Gfps are highly stable within cells, with an estimated 1/2-life between 7 h (wildtype) and 70 h (F100S/M154T/V164A) in S. cerevisiae cells. Although this limits their use in examining rapid cellular events without further modification, Gfp is expected to be a useful marker for monitoring the physiological state of cells in bioreactors using on-line probes.

Selective synchronization of Tetrahymena pyriformis cell populations and cell growth kinetics during the cell cycle.

A selection synchronization technique based on ingestion of tantalum particles has been applied to obtain synchronized cultures of the filter feeding ciliate Tetrahymena pyriformis. Cell concentrations and cell volume distributions of synchronized cell populations have been monitored for more than four average generation times. A simple curve-fitting method has been used to decompose the cell volume distributions of a synchronous population into two populations representing cells before and after cell division. In this way, the time course of the growth of an initially synchronous culture is decomposed into the growth of successive generations. The data indicate that at any given time only two generations of cells are present in significant numbers. The measured cell volume distributions show that T. pyriformis has a complicated growth pattern during the cell cycle. Newborn T. pyriformis cells do not grow significantly at the beginning of the cell cycle. After the nongrowing stage, cells start growing in a possibly exponential rate before cells enter into a second nongrowing stage. The second nongrowing stage lasts until cell division. The presented data demonstrate that growing cell populations can be viewed as the composite of cells belonging to different generations. This concept has important implications for solving corpuscular models of cell growth.

Metabolic modeling of polyhydroxybutyrate biosynthesis.

A mathematical model describing intracellular polyhydroxybutyrate (PHB) synthesis in Alcaligenes eutrophus has been constructed. The model allows investigation of issues such as the existence of rate-limiting enzymatic steps, possible regulatory mechanisms in PHB synthesis, and the effects different types of rate expressions have on model behavior. Simulations with the model indicate that activities of all PHB pathway enzymes influence overall PHB flux and that no single enzymatic step can easily be identified as rate limiting. Simulations also support regulatory roles for both thiolase and reductase, mediated through AcCoA/CoASH and NADPH/NADP+ ratios, respectively. To make the model more realistic, complex rate expressions for enzyme-catalyzed reactions were used which reflect both the reversibility of the reactions and the reaction mechanisms. Use of the complex kinetic expressions dramatically changed the behavior of the system compared to a simple model containing only Michaelis-Menten kinetic expressions; the more complicated model displayed different responses to changes in enzyme activities as well as inhibition of flux by the reaction products CoASH and NADP+. These effects can be attributed to reversible rate expressions, which allow prediction of reaction rates under conditions both near and far from equilibrium.

Quantitating secretion rates of individual cells: design of secretion assays.

To observe events occurring in the microenvironment surrounding individual cells, a mathematical framework has been developed describing the behavior of a compound following its secretion by a single cell. This description is based on the diffusional and binding processes taking place in the vicinity of the cell surface. It allows prediction of the rate of capture and accumulation of a secreted compound around a single cell. This concept provides the basis for the design of two experimental assays for measuring single-cell secretion rates: (1) Cells are immobilized in hydrogel microbeads which contain capture sites for the secreted compound; and (2) artificial receptors are bound directly to the cell surface which are capable of binding molecules secreted by individual cells. This general methodology is developed in the specific case of the model organism Saccharomyces cerevisiae secreting a heterologous protein, but can be applied to any cell/secreted protein combination. Binding studies have shown that approximately 2 x 10(5) of these artificial receptors can be attached to the surface of a single yeast cell. At this surface density of a putative artificial receptor, it is predicted that single-cell secretion rates of 47 molecules/cell/sec of a 150 kDa protein can be detected. Simulations indicate that a microbead loaded with 5 x 10(6) capture antibodies will result in detection of secretion of this protein at rates as low as 4 molecules/cell/sec.

Growth kinetics, nutrient uptake, and expression of the Alcaligenes eutrophus poly(beta-hydroxybutyrate) synthesis pathway in transgenic maize cell suspension cultures.

Transgenic suspension cultures of Black Mexican Sweet maize (Zea mays L.) expressing the Alcaligenes eutrophus genes encoding enzymes of the pathway for biosynthesis of the biodegradable polymer poly(beta-hydroxybutyrate) (PHB) were established as a tool for investigating metabolic regulation of the PHB pathway in plant cells. Cultures were grown in a 2 L modified mammalian cell bioreactor and in shake flasks. Biomass doubling times for transgenic bioreactor cultures (3.42 +/- 0.76 days) were significantly higher than those for untransformed cultures (2.01 +/- 0.33 days). Transgenic expression of the bacterial enzymes beta-ketothiolase (0.140 units/mg protein) and acetoacetyl-CoA reductase (0.636 units/mg protein) was detected by enzyme assays and immunoblots. However, over the first 2 years of cultivation, reductase activity decreased to 0.120 units/mg proteins. Furthermore, the PHB synthase gene, although initially present, was not detectable after 1.5 years of cultivation in suspension culture. These facts suggest that transgenic expression of PHB pathway genes in plant cells may not be stable. A hydroxybutyrate derivative was detected via gas chromatography even after 4 years of cultivation. Although the method used to prepare samples for gas chromatography cannot directly distinguish among PHB polymer, hydroxybutyryl-CoA (HB-CoA), and hydroxybutyric acid, solvent washing experiments indicated that most or all of the signal was non-polymeric, presumably H-CoA. The synthesis of HB-CoA appeared to be linked to substrate growth limitation, with HB-CoA accumulation increasing dramatically and cell growth ceasing upon depletion of ammonium. This suggests that the PHB synthesis pathway in plants is subject to regulatory mechanisms similar to those in prokaryotic cells.

Quantitative analysis of transient gene expression in mammalian cells using the green fluorescent protein.

The green fluorescent protein (Gfp) has been used as a reporter, along with flow cytometric analysis, to follow the dynamics of gene expression in transiently transfected mammalian cells. Gene transfer conditions for lipofection were optimized. The highest fraction of transfectants were obtained when lipid-DNA complexes were formed with 6 microliters lipid and 1 microgram DNA for chinese hamster ovary (CHO) cells and with 9 microliters lipid and 2 micrograms DNA for NIH/3T3 cells. Chinese hamster ovary cells were monitored for Gfp expression and growth for 6 days following transfection. An initial decrease in viability for 36 h was observed after which cell growth followed exponential kinetics with increasing viability. Intracellular accumulation of recombinant protein peaked at 24 h post-transfection and then decreased with first order kinetics at a rate comparable to the specific growth rate. It appears that dilution by growth accounts for the decrease of Gfp in the biomass. Immunofluorescent staining of Gfp and subsequent flow cytometric analysis of transfected cells revealed a linear correlation between the green fluorescence and immunofluorescence. This indicates that green fluorescence is a quantitative measure of intracellular Gfp in single cells in spite of the dynamics of post-translational modifications involved in the conversion of expressed protein into its fluorescent form. A structured model has been formulated to describe the observed kinetics of gene expression and fluorophore formation. The model accurately predicts experimental trends and suggests that the fraction of non-fluorescent Gfp is significant only during the initial period of gene expression.

Production of heteropolymeric polyhydroxyalkanoate in Escherichia coli from a single carbon source.

Poly[beta-hydroxybutyrate-co-beta-hydroxyvalerate] co-polymer, PHBV, is a polyhydroxyalkanoate (PHA) that has greater utility as a biodegradable thermoplastic polyester than poly-beta-hydroxybutyrate, PHB. In order to produce PHBV, a system of pathways is required to produce both hydroxybutyrate (HB) and hydroxyvalerate (HV) monomers from the sources of carbon. A working model for conversion of glucose to PHBV via acetyl- and propionyl-coenzyme A was constructed by expressing the PHA biosynthesis genes from Alcaligenes eutrophus in Escherichia coli strain K-12 under novel growth conditions. When 1 mM valine was added to 1% glucose medium, growth ceased and up to 2.5% of the incorporated monomers were HV; up to 4% were HV when 1 mM threonine was added as well. Threonine dehydratase (TD) converts threonine to alpha-ketobutyrate; TD is required for HV to be incorporated into PHA unless its transaminated reaction product, alpha-aminobutyrate, is added to the medium. Intracellular alpha-ketobutyrate accumulates when valine is added to the medium because valine, which cannot be metabolized to HV by E. coli strain K-12, stimulates TD and inhibits acetolactate synthase. In turn, alpha-ketobutyrate is converted to propionyl-CoA by the E. coli pyruvate dehydrogenase complex. This constitutes a defined system of pathways for synthesis of a heteropolymeric PHA from a single carbon source, which in the future could be transferred to other organisms including plants.

Surface IgG content of murine hybridomas: Direct evidence for variation of antibody secretion rates during the cell cycle.

Previous experiments have shown that population average surface lgG content is correlated with the specific antibody production rates of batch hybridoma cultures. Therefore, surface associated lgG content of single hybridoma cells might indicate antibody secretion rates of individual cells. Moreover, the surface lgG content should reflect the pattern of secretion rates during the cell cycle. To probe for lgG secretion rates during the cellcycle, a double staining procedure has been developed allowing simultaneousflow cytometric analysis of surface lgG content and DNA content of murine hybridoma cells. Crosslinking of the surface associated immunofluorescence with the cell by paraformaldehyde fixation permits subsequent DNA staining without loss of immunofluorescence. The optimized protocol has been used to determine the pattern of the surface lgG fluorescence as a function of the cell cycle position. It is highest during the G2+M cell cycle phase and the experimental data are in excellent agreement with the previously predicted secretion pattern during the cell cycle. (c) 1995 John Wiley & Sons Inc.

Dynamics of maize endosperm development and DNA endoreduplication.

Endosperm development in Zea mays is characterized by a period of intense mitotic activity followed by a period in which mitosis is essentially eliminated and the cell cycle becomes one of alternating S and G phases, leading to endoreduplication of the nuclear DNA. The endosperm represents a significant contribution to the grain yield of maize; thus, methods that facilitate the study of cellular kinetics may be useful in discerning cellular and molecular components of grain yield. Two mathematical models have been developed to describe the kinetics of endosperm growth. The first describes the kinetics of mitosis during endosperm development; the second describes the kinetics of DNA endoreduplication during endosperm development. The mitotic model is a modification of standard growth curves. The endoreduplication model is composed of six differential equations that represent the progression of nuclei from one DNA content to another during the endoreduplication process. Total nuclei number per endosperm and the number of 3C, 6C, 12C, 24C, 48C, and 96C nuclei per endosperm (C is the haploid DNA content per nucleus) for inbred W64A from 8 to 18 days after pollination were determined by flow cytometry. The results indicate that the change in number of nuclei expressed as a function of the number of days after pollination is the same from one yearly crop to another. These data were used in the model to determine the endosperm growth rate, the maximum nuclei number per endosperm, and transition rates from one C value to the next higher C value. The kinetics of endosperm development are reasonably well represented by the models. Thus, the models provide a means to quantify the complex pattern of endosperm development.

Tracking of individual cell cohorts in asynchronous Saccharomyces cerevisiae populations.

A novel flow cytometric procedure has been developed with the aim to obtain the growth properties of individual Saccharomyces cerevisiae cells in asynchronous culture. The method is based on labeling of the cell surface with FITC-conjugated concanavalin A and detection of the single-cell fluorescence with flow cytometry after cell exposure to growth conditions. Because the formation of new cell wall material in budded cells is restricted to the bud tip, exposure of the stained cells to growth conditions results in three cell types: (i) stained cells, (ii) partially stained cells, and (iii) unstained cells. Analysis of the staining pattern over time permits the determination of the specific growth rate of the cell population, the length of the budded cell cycle phase, and the growth pattern during the cell cycle of newly formed, partially stained daughter cells. The procedure has been tested with yeast cell populations growing at different rates. The data suggest an exponential increase in the size of individual cells during the cell cycle, as reflected by the forward angle light scattering (FALS) signals. It has been found that the apparent single-cell specific cell size growth rates, determined by FALS intensity, are significantly lower than the specific growth rates of the overall population. This could indicate that the tracking of a cohort of cells is significantly perturbed by a distribution of staining levels of daughter cells at cell division and that FALS may not be a good indicator of the cell size.

Multistaged corpuscular models of microbial growth: Monte Carlo simulations.

A new framework is developed by extending the existing population balance framework for modeling the growth of microbial populations. The new class of multistaged corpuscular models allows further structuring of the microbial life cycle into separate phases or stages and thus facilitates the incorporation of cell cycle phenomena to population models. These multistaged models consist of systems of population balance equations coupled by appropriate boundary conditions. The specific form of the equations depend on the assumed forms for the transition rate functions, the growth rate functions, and the partitioning function, which determines how the biological material is distributed at division. A growth model for ciliated protozoa is formulated to demonstrate the proposed framework. To obtain a solution to the system of the partial integro differential equations that results from such formulation, we adopted a Monte Carlo simulation technique which is very stable, versatile, and insensitive to the complexity of the model. The theory and implementation of the Monte Carlo simulation algorithm is analyzed and results from the simulation of the ciliate growth model are presented. The proposed approach seems to be promising for integrating single-cell mechanisms into population models.