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.

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

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.

Unsupervised Identification of Isotope-Labeled Peptides.

In vivo isotopic labeling coupled with high-resolution proteomics is used to investigate primary metabolism in techniques such as stable isotope probing (protein-SIP) and peptide-based metabolic flux analysis (PMFA). Isotopic enrichment of carbon substrates and intracellular metabolism determine the distribution of isotopes within amino acids. The resulting amino acid mass distributions (AMDs) are convoluted into peptide mass distributions (PMDs) during protein synthesis. With no a priori knowledge on metabolic fluxes, the PMDs are therefore unknown. This complicates labeled peptide identification because prior knowledge on PMDs is used in all available peptide identification software. An automated framework for the identification and quantification of PMDs for nonuniformly labeled samples is therefore lacking. To unlock the potential of peptide labeling experiments for high-throughput flux analysis and other complex labeling experiments, an unsupervised peptide identification and quantification method was developed that uses discrete deconvolution of mass distributions of identified peptides to inform on the mass distributions of otherwise unidentifiable peptides. Uniformly (13)C-labeled Escherichia coli protein was used to test the developed feature reconstruction and deconvolution algorithms. The peptide identification was validated by comparing MS(2)-identified peptides to peptides identified from PMDs using unlabeled E. coli protein. Nonuniformly labeled Glycine max protein was used to demonstrate the technology on a representative sample suitable for flux analysis. Overall, automatic peptide identification and quantification were comparable or superior to manual extraction, enabling proteomics-based technology for high-throughput flux analysis studies.

Metabolic flux analysis using ¹³C peptide label measurements.

¹³C metabolic flux analysis (MFA) has become the experimental method of choice to investigate the cellular metabolism of microbes, cell cultures and plant seeds. Conventional steady-state MFA utilizes isotopic labeling measurements of amino acids obtained from protein hydrolysates. To retain spatial information in conventional steady-state MFA, tissues or subcellular fractions must be dissected or biochemically purified. In contrast, peptides retain their identity in complex protein extracts, and may therefore be associated with a specific time of expression, tissue type and subcellular compartment. To enable 'single-sample' spatially and temporally resolved steady-state flux analysis, we investigated the suitability of peptide mass distributions (PMDs) as an alternative to amino acid label measurements. PMDs are the discrete convolution of the mass distributions of the constituent amino acids of a peptide. We investigated the requirements for the unique deconvolution of PMDs into amino acid mass distributions (AAMDs), the influence of peptide sequence length on parameter sensitivity, and how AAMD and flux estimates that are determined through deconvolution compare to estimates from a conventional GC-MS measurement-based approach. Deconvolution of PMDs of the storage protein ß-conglycinin of soybean (Glycine max) resulted in good AAMD and flux estimates if fluxes were directly fitted to PMDs. Unconstrained deconvolution resulted in inferior AAMD and flux estimates. PMD measurements do not include amino acid backbone fragments, which increase the information content in GC-MS-derived analyses. Nonetheless, the resulting flux maps were of comparable quality due to the precision of Orbitrap quantification and the larger number of peptide measurements.

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.

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.

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.

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.

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.

Design of substrate label for steady state flux measurements in plant systems using the metabolic network of Brassica napus embryos.

Steady state metabolic flux analysis using (13)C labeled substrates is of growing importance in plant physiology and metabolic engineering. The quality of the flux estimates in (13)C metabolic flux analysis depend on the: (i) network structure; (ii) flux values; (iii) design of the labeling substrate; and (iv) label measurements performed. Whereas the first two parameters are facts of nature, the latter two are to some extent controlled by the experimenter, yet they have received little attention in most plant studies. Using the metabolic flux map of developing Brassica napus (Rapeseed) embryos, this study explores the value of optimal substrate label designs obtained with different statistical criteria and addresses the applicability of different optimal designs to biological questions. The results demonstrate the value of optimizing the choice of labeled substrates and show that substrate combinations commonly used in bacterial studies can be far from optimal for mapping fluxes in plant systems. The value of performing additional experiments and the inclusion of measurements is also evaluated.

Stringent control of FLP recombinase in Escherichia coli.

Site specific recombinases are invaluable tools in molecular biology, and are emerging as powerful recorders of cellular events in synthetic biology. We have developed a stringently controlled FLP recombinase system in Escherichia coli using an arabinose inducible promoter combined with a weak ribosome binding site.

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.

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.

Optimal design of isotope labeling experiments.

Stable isotope labeling experiments (ILE) constitute a powerful methodology for estimating metabolic fluxes. An optimal label design for such an experiment is necessary to maximize the precision with which fluxes can be determined. But often, precision gained in the determination of one flux comes at the expense of the precision of other fluxes, and an appropriate label design therefore foremost depends on the question the investigator wants to address. One could liken ILE to shadows that metabolism casts on products. Optimal label design is the placement of the lamp; creating clear shadows for some parts of metabolism and obscuring others.An optimal isotope label design is influenced by: (1) the network structure; (2) the true flux values; (3) the available label measurements; and, (4) commercially available substrates. The first two aspects are dictated by nature and constrain any optimal design. The second two aspects are suitable design parameters. To create an optimal label design, an explicit optimization criterion needs to be formulated. This usually is a property of the flux covariance matrix, which can be augmented by weighting label substrate cost. An optimal design is found by using such a criterion as an objective function for an optimizer. This chapter uses a simple elementary metabolite units (EMU) representation of the TCA cycle to illustrate the process of experimental design of isotope labeled substrates.

Nitric oxide methods in seed biology.

The ubiquitous signaling molecule nitric oxide (NO) plays an important role in seed biology. Experiments with this biologically important gas require special provisions because NO in aerobic environments is readily converted into other oxides of nitrogen. In this chapter, we describe methods for the application of NO as a gas, and through the use of NO-donor compounds. We included information on the removal or reduction of NO with NO scavengers. Methods for detecting NO using NO-reactive fluorescent probes, and an apparatus incorporating an oxidizer column are also described.

The Arabidopsis aleurone layer responds to nitric oxide, gibberellin, and abscisic acid and is sufficient and necessary for seed dormancy.

Seed dormancy is a common phase of the plant life cycle, and several parts of the seed can contribute to dormancy. Whole seeds, seeds lacking the testa, embryos, and isolated aleurone layers of Arabidopsis (Arabidopsis thaliana) were used in experiments designed to identify components of the Arabidopsis seed that contribute to seed dormancy and to learn more about how dormancy and germination are regulated in this species. The aleurone layer was found to be the primary determinant of seed dormancy. Embryos from dormant seeds, however, had a lesser growth potential than those from nondormant seeds. Arabidopsis aleurone cells were examined by light and electron microscopy, and cell ultrastructure was similar to that of cereal aleurone cells. Arabidopsis aleurone cells responded to nitric oxide (NO), gibberellin (GA), and abscisic acid, with NO being upstream of GA in a signaling pathway that leads to vacuolation of protein storage vacuoles and abscisic acid inhibiting vacuolation. Molecular changes that occurred in embryos and aleurone layers prior to germination were measured, and these data show that both the aleurone layer and the embryo expressed the NO-associated gene AtNOS1, but only the embryo expressed genes for the GA biosynthetic enzyme GA3 oxidase.

Nitric oxide reduces seed dormancy in Arabidopsis.

Dormancy is a property of many mature seeds, and experimentation over the past century has identified numerous chemical treatments that will reduce seed dormancy. Nitrogen-containing compounds including nitrate, nitrite, and cyanide break seed dormancy in a range of species. Experiments are described here that were carried out to further our understanding of the mechanism whereby these and other compounds, such as the nitric oxide (NO) donor sodium nitroprusside (SNP), bring about a reduction in seed dormancy of Arabidopsis thaliana. A simple method was devised for applying the products of SNP photolysis through the gas phase. Using this approach it was shown that SNP, as well as potassium ferricyanide (Fe(III)CN) and potassium ferrocyanide (Fe(II)CN), reduced dormancy of Arabidopsis seeds by generating cyanide (CN). The effects of potassium cyanide (KCN) on dormant seeds were tested and it was confirmed that cyanide vapours were sufficient to break Arabidopsis seed dormancy. Nitrate and nitrite also reduced Arabidopsis seed dormancy and resulted in substantial rates of germination. The effects of CN, nitrite, and nitrate on dormancy were prevented by the NO scavenger c-PTIO. It was confirmed that NO plays a role in reducing seed dormancy by using purified NO gas, and a model to explain how nitrogen-containing compounds may break dormancy in Arabidopsis is presented.

Sodium nitroprusside, cyanide, nitrite, and nitrate break Arabidopsis seed dormancy in a nitric oxide-dependent manner.

The seeds of many plant species are dormant at maturity and dormancy loss is a prerequisite for germination. Numerous environmental and chemical treatments are known to lessen or remove seed dormancy, but the biochemical changes that occur during this change of state are poorly understood. Several lines of research have implicated nitric oxide (NO) as a participant in this process. Here, we show that dormant seeds of Arabidopsis thaliana (L.) Heynh. will germinate following treatment with the NO donor sodium nitroprusside (SNP), cyanide (CN), nitrite or nitrate. In all cases, the NO scavenger c-PTIO effectively promotes the maintenance of seed dormancy. c-PTIO does not, however, inhibit germination of fully after-ripened seeds, and c-PTIO does not interact directly with nitrite, nitrate or CN. We also show that volatile CN effectively breaks dormancy of Arabidopsis seeds, and that CN is the volatile compound in SNP that promotes dormancy loss. Our data support the hypothesis that NO is a signaling molecule that plays an important role in the loss of seed dormancy.

Nitric oxide gas stimulates germination of dormant Arabidopsis seeds: use of a flow-through apparatus for delivery of nitric oxide.

Nitric oxide (NO) is a gaseous free radical that reacts with O(2) in air and aqueous solution. NO donors have been widely used to circumvent the difficulties inherent in working with a reactive gas, but NO donors do not deliver NO at a constant rate for prolonged periods of time. Furthermore, some of the most commonly used NO donors produce additional, bioactive decomposition products. We designed and built an apparatus that allowed for the precise mixing of gaseous NO with air and the delivery of gas through sample vials at fixed rates. This experimental setup has the added advantage that continuous flow of gas over the sample reduces the buildup of volatile breakdown products. To show that this experimental setup was suitable for studies on the dormancy and germination of Arabidopsis thaliana seeds, we introduced vapors from water or sodium nitroprusside (SNP) into the gas stream. Seeds remained dormant when treated with water vapor, but gases generated by SNP increased germination to 90%. When pure NO was mixed with air and passed over dormant seeds, approximately approximately 30% of the seeds germinated. Because nitrite accumulates in aqueous solutions exposed to NO gas, we measured the accumulation of nitrite under our experimental conditions and found that it did not exceed 100 microM. Nitrite or nitrate at concentrations of up to 500 microM did not increase germination of C24 ecotype Arabidopsis seeds to more than 10%. These data support the hypothesis that NO participates in the loss of Arabidopsis seed dormancy, and they show that for some dormant seeds, exposure to exogenous NO is sufficient to trigger germination.