John A. Morgan

Flux balance analysis of primary metabolism in Chlamydomonas reinhardtii

BackgroundPhotosynthetic organisms convert atmospheric carbon dioxide into numerous metabolites along the pathways to make new biomass. Aquatic photosynthetic organisms, which fix almost half of global inorganic carbon, have great potential: as a carbon dioxide fixation method, for the economical production of chemicals, or as a source for lipids and starch which can then be converted to biofuels. To harness this potential through metabolic engineering and to maximize production, a more thorough understanding of photosynthetic metabolism must first be achieved. A model algal species, C. reinhardtii, was chosen and the metabolic network reconstructed. Intracellular fluxes were then calculated using flux balance analysis (FBA).ResultsThe metabolic network of primary metabolism for a green alga, C. reinhardtii, was reconstructed using genomic and biochemical information. The reconstructed network accounts for the intracellular localization of enzymes to three compartments and includes 484 metabolic reactions and 458 intracellular metabolites. Based on BLAST searches, one newly annotated enzyme (fructose-1,6-bisphosphatase) was added to the Chlamydomonas reinhardtii database. FBA was used to predict metabolic fluxes under three growth conditions, autotrophic, heterotrophic and mixotrophic growth. Biomass yields ranged from 28.9 g per mole C for autotrophic growth to 15 g per mole C for heterotrophic growth.ConclusionThe flux balance analysis model of central and intermediary metabolism in C. reinhardtii is the first such model for algae and the first model to include three metabolically active compartments. In addition to providing estimates of intracellular fluxes, metabolic reconstruction and modelling efforts also provide a comprehensive method for annotation of genome databases. As a result of our reconstruction, one new enzyme was annotated in the database and several others were found to be missing; implying new pathways or non-conserved enzymes. The use of FBA to estimate intracellular fluxes also provides flux values that can be used as a starting point for rational engineering of C. reinhardtii. From these initial estimates, it is clear that aerobic heterotrophic growth on acetate has a low yield on carbon, while mixotrophically and autotrophically grown cells are significantly more carbon efficient.

Mapping photoautotrophic metabolism with isotopically nonstationary 13C flux analysis

Understanding in vivo regulation of photoautotrophic metabolism is important for identifying strategies to improve photosynthetic efficiency or re-route carbon fluxes to desirable end products. We have developed an approach to reconstruct comprehensive flux maps of photoautotrophic metabolism by computational analysis of dynamic isotope labeling measurements and have applied it to determine metabolic pathway fluxes in the cyanobacterium Synechocystis sp. PCC6803. Comparison to a theoretically predicted flux map revealed inefficiencies in photosynthesis due to oxidative pentose phosphate pathway and malic enzyme activity, despite negligible photorespiration. This approach has potential to fill important gaps in our understanding of how carbon and energy flows are systemically regulated in cyanobacteria, plants, and algae.

Flux Balance Analysis of Photoautotrophic Metabolism

Photosynthesis is the principal process responsible for fixation of inorganic carbon dioxide into organic molecules with sunlight as the energy source. Potentially, many chemicals could be inexpensively produced by photosynthetic organisms. Mathematical modeling of photoautotrophic metabolism is therefore important to evaluate maximum theoretical product yields and to deeply understand the interactions between biochemical energy, carbon fixation, and assimilation pathways. Flux balance analysis based on linear programming is applied to photoautotrophic metabolism. The stoichiometric network of a model photosynthetic prokaryote, Synechocystis sp. PCC 6803, has been reconstructed from genomic data and biochemical literature and coupled with a model of the photophosphorylation processes. Flux map topologies for the hetero‐, auto‐, and mixotrophic modes of metabolism under conditions of optimal growth were determined and compared. The roles of important metabolic reactions such as the glyoxylate shunt and the transhydrogenase reaction were analyzed. We also theoretically evaluated the effect of gene deletions or additions on biomass yield and metabolic flux distributions.

Metabolic flux analysis of CHO cell metabolism in the late non-growth phase.

Chinese hamster ovary (CHO) cell cultures are commonly used for production of recombinant human therapeutic proteins. Often the goal of such a process is to separate the growth phase of the cells, from the non‐growth phase where ideally the cells are diverting resources to produce the protein of interest. Characterizing the way that the cells use nutrients in terms of metabolic fluxes as a function of culture conditions can provide a deeper understanding of the cell biology offering guidance for process improvements. To evaluate the fluxes, metabolic flux analysis of the CHO cell culture in the non‐growth phase was performed by a combination of steady‐state isotopomer balancing and stoichiometric modeling. Analysis of the glycolytic pathway and pentose phosphate pathway (PPP) indicated that almost all of the consumed glucose is diverted towards PPP with a high NADPH production; with even recycle from PPP to G6P in some cases. Almost all of the pyruvate produced from glycolysis entered the TCA cycle with little or no lactate production. Comparison of the non‐growth phase against previously reported fluxes from growth phase cultures indicated marked differences in the fluxes, in terms of the split between glycolysis and PPP, and also around the pyruvate node. Possible reasons for the high NADPH production are also discussed. Evaluation of the fluxes indicated that the medium strength, carbon dioxide level, and temperature with dissolved oxygen have statistically significant impacts on different nodes of the flux network. Biotechnol. Bioeng. 2011; 108:82–92.

Targeted metabolomic analysis of Escherichia coli by desorption electrospray ionization and extractive electrospray ionization mass spectrometry

Desorption electrospray ionization (DESI) was utilized to monitor the presence of targeted central carbon metabolites within bacterial cell extracts and the quench supernatant of Escherichia coli. The targeted metabolites were identified through tandem mass spectrometry (MS/MS) product ion scans using collision-induced dissociation in the negative ion mode. Picogram detection limits were achieved for a majority of the metabolites during MS/MS analysis of standard metabolite solutions. In a [U-(13)C]glucose pulse experiment, where uniformly labeled glucose was fed to E. coli, the corresponding fragment ions from labeled metabolites in extracts were generally observed. There was evidence of matrix effects including moderate suppression by other metabolites within the spectra of the labeled and unlabeled extracts. To improve the specificity and sensitivity of detection, optimized in situ ambient chemical reactions using DESI and extractive electrospray ionization (EESI) were carried out for targeted compounds. This study provides the first indication of the potential to perform in situ targeted metabolomics of a bacterial sample via ambient ionization mass spectrometry.

Metabolic cartography: experimental quantification of metabolic fluxes from isotopic labelling studies

For the past decade, flux maps have provided researchers with an in-depth perspective on plant metabolism. As a rapidly developing field, significant headway has been made recently in computation, experimentation, and overall understanding of metabolic flux analysis. These advances are particularly applicable to the study of plant metabolism. New dynamic computational methods such as non-stationary metabolic flux analysis are finding their place in the toolbox of metabolic engineering, allowing more organisms to be studied and decreasing the time necessary for experimentation, thereby opening new avenues by which to explore the vast diversity of plant metabolism. Also, improved methods of metabolite detection and measurement have been developed, enabling increasingly greater resolution of flux measurements and the analysis of a greater number of the multitude of plant metabolic pathways. Methods to deconvolute organelle-specific metabolism are employed with increasing effectiveness, elucidating the compartmental specificity inherent in plant metabolism. Advances in metabolite measurements have also enabled new types of experiments, such as the calculation of metabolic fluxes based on (13)CO(2) dynamic labelling data, and will continue to direct plant metabolic engineering. Newly calculated metabolic flux maps reveal surprising and useful information about plant metabolism, guiding future genetic engineering of crops to higher yields. Due to the significant level of complexity in plants, these methods in combination with other systems biology measurements are necessary to guide plant metabolic engineering in the future.

Integrating cybernetic modeling with pathway analysis provides a dynamic, systems-level description of metabolic control

Cybernetic modeling strives to uncover the inbuilt regulatory programs of biological systems and leverage them toward computational prediction of metabolic dynamics. Because of its focus on incorporating the global aims of metabolism, cybernetic modeling provides a systems‐oriented approach for describing regulatory inputs and inferring the impact of regulation within biochemical networks. Combining cybernetic control laws with concepts from metabolic pathway analysis has culminated in a systematic strategy for constructing cybernetic models, which was previously lacking. The newly devised framework relies upon the simultaneous application of local controls that maximize the net flux through each elementary flux mode and global controls that modulate the activities of these modes to optimize the overall nutritional state of the cell. The modeling concepts are illustrated using a simple linear pathway and a larger network representing anaerobic E. coli central metabolism. The E. coli model successfully describes the metabolic shift that occurs upon deleting the pta‐ackA operon that is responsible for fermentative acetate production. The model also furnishes predictions that are consistent with experimental results obtained from additional knockout strains as well as strains expressing heterologous genes. Because of the stabilizing influence of the included control variables, the resulting cybernetic models are more robust and reliable than their predecessors in simulating the network response to imposed genetic and environmental perturbations. Biotechnol. Bioeng. 2008;100: 542–559.

Rethinking how volatiles are released from plant cells

For plant volatile organic compounds (VOCs) to be emitted, they must cross membrane(s), the aqueous cell wall, and sometimes the cuticle, before moving into the gas phase. It is presumed that VOC movement through each barrier occurs via passive diffusion. However, VOCs, which are primarily nonpolar compounds, will preferentially partition into membranes, making diffusion into aqueous compartments slow. Using Ficks first law, we calculated that to achieve observed VOC emission rates by diffusion alone would necessitate toxic VOC levels in membranes. Here, we propose that biological mechanisms, such as those involved in trafficking other hydrophobic compounds, must contribute to VOC emission. Such parallel biological pathways would lower barrier resistances and, thus, steady-state emission rates could be maintained with significantly reduced intramembrane VOC concentrations.

Systematic development of hybrid cybernetic models: Application to recombinant yeast co-consuming glucose and xylose

The hybrid cybernetic modeling approach of Kim et al. (Kim et al. [2008] Biotechnol. Prog., in press) views the substrate uptake flux in microorganisms as being distributed in a regulated way among different elementary modes (EMs) of a metabolic network, which intracellular fluxes related to the uptake rates by the pseudo‐steady‐state approximation on intracellular metabolites. While the conceptual development has been demonstrated by Kim et al. (Kim et al. [2008] Biotechnol. Prog., in press) using a rather simple example (i.e., Escherichia coli metabolizing a single substrate), its extension to a larger scale network involving multiple substrates results in serious overparameterization (which implies an excessive number of parameters relative to the measurements available to determine them). Through the case study of recombinant Saccharomyces yeast co‐consuming glucose and xylose, we present a systematic way of formulating a minimal order hybrid cybernetic model (HCM) for a general metabolic network. The overparameterization problem mostly arising from a large number of EMs is avoided using a model reduction technique developed by Song and Ramkrishna (Song and Ramkrishna [2009a] Biotechnol. Bioeng. 102(2):554–568) where an original set of EMs is condensed to a much smaller subset. Detailed discussions follow on the issue of determining the minimal set of active modes needed for the description of the simultaneous consumption of multiple substrates. The developed HCM is compared with other metabolic models: macroscopic bioreaction models (Provost et al. [2006] Bioprocess Biosyt. Eng. 29(5–6):349–366), and dynamic flux balance analysis. It is shown that the HCM outperforms the other two as validated using various sets of fermentation data. The difference among the models is more dramatic in a situation such as the sequential utilization of glucose and xylose, which is observed under realistic fermentation conditions. Biotechnol. Bioeng. 2009;103: 984–1002.

Computation of metabolic fluxes and efficiencies for biological carbon dioxide fixation

With rising energy prices and concern over the environmental impact of fossil fuel consumption, the push to develop biomass derived fuels has increased significantly. Although most global carbon fixation occurs via the Calvin Benson Bassham cycle, there are currently five other known pathways for carbon fixation; the goal of this study was to determine the thermodynamic efficiencies of all six carbon fixation pathways for the production of biomass using flux balance analysis. The three chemotrophic pathways, the reductive acetyl-CoA pathway, the 3-hydroxypropionate/4-hydroxybutyrate cycle and the dicarboxylate/4-hydroxybutyrate cycle, were found to be more efficient than photoautotrophic carbon fixation pathways. However, as hydrogen is not freely available, the energetic cost of hydrogen production from sunlight was calculated and included in the overall energy demand, which results in a 5 fold increase in the energy demand of chemoautotrophic carbon fixation. Therefore, when the cost of hydrogen production is included, photoautotrophic pathways are more efficient. However, the energetic cost for the production of 12 metabolic precursors was found to vary widely across the different carbon fixation pathways; therefore, different pathways may be more efficient at producing products from a single precursor than others. The results of this study have significant impact on the selection or design of autotrophic organisms for biofuel or biochemical production. Overall biomass production from solar energy is most efficient in organisms using the reductive TCA cycle, however, products derived from one metabolic precursor may be more efficiently produced using other carbon fixation pathways.