Igor G. L. Libourel

Metabolic flux analysis in plants: coping with complexity.

Theory and experience in metabolic engineering both show that metabolism operates at the network level. In plants, this complexity is compounded by a high degree of compartmentation and the synthesis of a very wide array of secondary metabolic products. A further challenge to understanding and predicting plant metabolic function is posed by our ignorance about the structure of metabolic networks even in well-studied systems. Metabolic flux analysis (MFA) provides tools to measure and model the functioning of metabolism, and is making significant contributions to coping with their complexity. This review gives an overview of different MFA approaches, the measurements required to implement them and the information they yield. The application of MFA methods to plant systems is then illustrated by several examples from the recent literature. Next, the challenges that plant metabolism poses for MFA are discussed together with ways that these can be addressed. Lastly, new developments in MFA are described that can be expected to improve the range and reliability of plant MFA in the coming years.

Metabolic Flux Analysis in Plants: From Intelligent Design to Rational Engineering

Metabolic flux analysis (MFA) is a rapidly developing field concerned with the quantification and understanding of metabolism at the systems level. The application of MFA has produced detailed maps of flow through metabolic networks of a range of plant systems. These maps represent detailed metabolic phenotypes, contribute significantly to our understanding of metabolism in plants, and have led to the discovery of new metabolic routes. The presentation of thorough statistical evaluation with current flux maps has set a new standard for the quality of quantitative flux studies. In microbial systems, powerful methods have been developed for the reconstruction of metabolic networks from genomic and transcriptomic data, pathway analysis, and predictive modeling. This review brings together the recent developments in quantitative MFA and predictive modeling. The application of predictive tools to high quality flux maps in particular promises to be important in the rational metabolic engineering of plants.

Nitrite reduces cytoplasmic acidosis under Anoxia.

The ameliorating effect of nitrate on the acidification of the cytoplasm during short-term anoxia was investigated in maize (Zea mays) root segments. Seedlings were grown in the presence or absence of nitrate, and changes in the cytoplasmic and vacuolar pH in response to the imposition of anoxia were measured by in vivo 31P nuclear magnetic resonance spectroscopy. Soluble ions and metabolites released to the suspending medium by the anoxic root segments were measured by high-performance liquid chromatography and 1H nuclear magnetic resonance spectroscopy, and volatile metabolites were measured by gas chromatography and gas chromatography-mass spectrometry. The beneficial effect of nitrate on cytoplasmic pH regulation under anoxia occurred despite limited metabolism of nitrate under anoxia, and modest effects on the ions and metabolites, including fermentation end products, released from the anoxic root segments. Interestingly, exposing roots grown and treated in the absence of nitrate to micromolar levels of nitrite during anoxia had a beneficial effect on the cytoplasmic pH that was comparable to the effect observed for roots grown and treated in the presence of nitrate. It is argued that nitrate itself is not directly responsible for improved pH regulation under anoxia, contrary to the usual assumption, and that nitrite rather than nitrate should be the focus for further work on the beneficial effect of nitrate on flooding tolerance.

Genome-wide metabolic network reconstruction of the picoalga Ostreococcus

The green picoalga Ostreococcus is emerging as a simple plant model organism, and two species, O. lucimarinus and O. tauri, have now been sequenced and annotated manually. To evaluate the completeness of the metabolic annotation of both species, metabolic networks of O. lucimarinus and O. tauri were reconstructed from the KEGG database, thermodynamically constrained, elementally balanced, and functionally evaluated. The draft networks contained extensive gaps and, in the case of O. tauri, no biomass components could be produced due to an incomplete Calvin cycle. To find and remove gaps from the networks, an extensive reference biochemical reaction database was assembled using a stepwise approach that minimized the inclusion of microbial reactions. Gaps were then removed from both Ostreococcus networks using two existing gap-filling methodologies. In the first method, a bottom-up approach, a minimal list of reactions was added to each model to enable the production of all metabolites included in our biomass equation. In the second method, a top-down approach, all reactions in the reference database were added to the target networks and subsequently trimmed away based on the sequence alignment scores of identified orthologues. Because current gap-filling methods do not produce unique solutions, a quality metric that includes a weighting for phylogenetic distance and sequence similarity was developed to distinguish between gap-filling results automatically. The draft O. lucimarinus and O. tauri networks required the addition of 56 and 70 reactions, respectively, in order to produce the same biomass precursor metabolites that were produced by our plant reference database.

Analysis of isotopic labeling in peptide fragments by tandem mass spectrometry

Phenotype in multicellular organisms is the consequence of dynamic metabolic events that occur in a spatially dependent fashion. This spatial and temporal complexity presents challenges for investigating metabolism; creating a need for improved methods that effectively probe biochemical events such as amino acid biosynthesis. Isotopic labeling can provide a temporal-spatial recording of metabolic events through, for example, the description of enriched amino acids in the protein pool. Proteins are therefore an important readout of metabolism and can be assessed with modern mass spectrometers. We compared the measurement of isotopic labeling in MS2 spectra obtained from tandem mass spectrometry under either higher energy collision dissociation (HCD) or collision induced dissociation (CID) at varied energy levels. Developing soybean embryos cultured with or without 13C-labeled substrates, and Escherichia coli MG1655 enriched by feeding 7% uniformly labeled glucose served as a source of biological material for protein evaluation. CID with low energies resulted in a disproportionate amount of heavier isotopologues remaining in the precursor isotopic distribution. HCD resulted in fewer quantifiable products; however deviation from predicted distributions were small relative to the CID-based comparisons. Fragment ions have the potential to provide information on the labeling of amino acids in peptides, but our results indicate that without further development the use of this readout in quantitative methods such as metabolic flux analysis is limited.

Quantification of Peptide m/z Distributions from 13C-Labeled Cultures with High-Resolution Mass Spectrometry

Isotopic labeling studies of primary metabolism frequently utilize GC/MS to quantify (13)C in protein-hydrolyzed amino acids. During processing some amino acids are degraded, which reduces the size of the measurement set. The advent of high-resolution mass spectrometers provides a tool to assess molecular masses of peptides with great precision and accuracy and computationally infer information about labeling in amino acids. Amino acids that are isotopically labeled during metabolism result in labeled peptides that contain spatial and temporal information that is associated with the biosynthetic origin of the protein. The quantification of isotopic labeling in peptides can therefore provide an assessment of amino acid metabolism that is specific to subcellular, cellular, or temporal conditions. A high-resolution orbital trap was used to quantify isotope labeling in peptides that were obtained from unlabeled and isotopically labeled soybean embryos and Escherichia coli cultures. Standard deviations were determined by estimating the multinomial variance associated with each element of the m/z distribution. Using the estimated variance, quantification of the m/z distribution across multiple scans was achieved by a nonlinear fitting approach. Observed m/z distributions of uniformly labeled E. coli peptides indicated no significant differences between observed and simulated m/z distributions. Alternatively, amino acid m/z distributions obtained from GC/MS were convolved to simulate peptide m/z distributions but resulted in distinct profiles due to the production of protein prior to isotopic labeling. The results indicate that peptide mass isotopologue measurements faithfully represent mass distributions, are suitable for quantification of isotope-labeling-based studies, and provide additional information over existing methods.

Genome-Scale Metabolic Network Validation of Shewanella oneidensis Using Transposon Insertion Frequency Analysis

Transposon mutagenesis, in combination with parallel sequencing, is becoming a powerful tool for en-masse mutant analysis. A probability generating function was used to explain observed miniHimar transposon insertion patterns, and gene essentiality calls were made by transposon insertion frequency analysis (TIFA). TIFA incorporated the observed genome and sequence motif bias of the miniHimar transposon. The gene essentiality calls were compared to: 1) previous genome-wide direct gene-essentiality assignments; and, 2) flux balance analysis (FBA) predictions from an existing genome-scale metabolic model of Shewanella oneidensis MR-1. A three-way comparison between FBA, TIFA, and the direct essentiality calls was made to validate the TIFA approach. The refinement in the interpretation of observed transposon insertions demonstrated that genes without insertions are not necessarily essential, and that genes that contain insertions are not always nonessential. The TIFA calls were in reasonable agreement with direct essentiality calls for S. oneidensis, but agreed more closely with E. coli essentiality calls for orthologs. The TIFA gene essentiality calls were in good agreement with the MR-1 FBA essentiality predictions, and the agreement between TIFA and FBA predictions was substantially better than between the FBA and the direct gene essentiality predictions.

Sequence-based Network Completion Reveals the Integrality of Missing Reactions in Metabolic Networks.

Background: Genome-scale draft metabolic networks are incomplete, even for well studied organisms. Results: Reactions selected by minimizing flux through unlikely reactions resulted in networks of superior quality. Conclusion: Genome-scale models have many network completion solutions but require the addition of unsupported reactions to be functional. Significance: Metabolic networks guide synthetic biology efforts, and the quality of networks determines their predictive power. Genome-scale metabolic models are central in connecting genotypes to metabolic phenotypes. However, even for well studied organisms, such as Escherichia coli, draft networks do not contain a complete biochemical network. Missing reactions are referred to as gaps. These gaps need to be filled to enable functional analysis, and gap-filling choices influence model predictions. To investigate whether functional networks existed where all gap-filling reactions were supported by sequence similarity to annotated enzymes, four draft networks were supplemented with all reactions from the Model SEED database for which minimal sequence similarity was found in their genomes. Quadratic programming revealed that the number of reactions that could partake in a gap-filling solution was vast: 3,270 in the case of E. coli, where 72% of the metabolites in the draft network could connect a gap-filling solution. Nonetheless, no network could be completed without the inclusion of orphaned enzymes, suggesting that parts of the biochemistry integral to biomass precursor formation are uncharacterized. However, many gap-filling reactions were well determined, and the resulting networks showed improved prediction of gene essentiality compared with networks generated through canonical gap filling. In addition, gene essentiality predictions that were sensitive to poorly determined gap-filling reactions were of poor quality, suggesting that damage to the network structure resulting from the inclusion of erroneous gap-filling reactions may be predictable.

Open source software to control Bioflo bioreactors.

Bioreactors are designed to support highly controlled environments for growth of tissues, cell cultures or microbial cultures. A variety of bioreactors are commercially available, often including sophisticated software to enhance the functionality of the bioreactor. However, experiments that the bioreactor hardware can support, but that were not envisioned during the software design cannot be performed without developing custom software. In addition, support for third party or custom designed auxiliary hardware is often sparse or absent. This work presents flexible open source freeware for the control of bioreactors of the Bioflo product family. The functionality of the software includes setpoint control, data logging, and protocol execution. Auxiliary hardware can be easily integrated and controlled through an integrated plugin interface without altering existing software. Simple experimental protocols can be entered as a CSV scripting file, and a Python-based protocol execution model is included for more demanding conditional experimental control. The software was designed to be a more flexible and free open source alternative to the commercially available solution. The source code and various auxiliary hardware plugins are publicly available for download from https://github.com/LibourelLab/BiofloSoftware. In addition to the source code, the software was compiled and packaged as a self-installing file for 32 and 64 bit windows operating systems. The compiled software will be able to control a Bioflo system, and will not require the installation of LabVIEW.

Thermodynamic Constraints Improve Metabolic Networks

In pursuit of establishing a realistic metabolic phenotypic space, the reversibility of reactions is thermodynamically constrained in modern metabolic networks. The reversibility constraints follow from heuristic thermodynamic poise approximations that take anticipated cellular metabolite concentration ranges into account. Because constraints reduce the feasible space, draft metabolic network reconstructions may need more extensive reconciliation, and a larger number of genes may become essential. Notwithstanding ubiquitous application, the effect of reversibility constraints on the predictive capabilities of metabolic networks has not been investigated in detail. Instead, work has focused on the implementation and validation of the thermodynamic poise calculation itself. With the advance of fast linear programming-based network reconciliation, the effects of reversibility constraints on network reconciliation and gene essentiality predictions have become feasible and are the subject of this study. Networks with thermodynamically informed reversibility constraints outperformed gene essentiality predictions compared to networks that were constrained with randomly shuffled constraints. Unconstrained networks predicted gene essentiality as accurately as thermodynamically constrained networks, but predicted substantially fewer essential genes. Networks that were reconciled with sequence similarity data and strongly enforced reversibility constraints outperformed all other networks. We conclude that metabolic network analysis confirmed the validity of the thermodynamic constraints, and that thermodynamic poise information is actionable during network reconciliation.