Pangenome reconstruction of Lactobacillaceae metabolism predicts species-specific metabolic traits.

Strains across the Lactobacillaceae family form the basis for a trillion-dollar industry. Our understanding of the genomic basis for their key traits is fragmented, however, including the metabolism that is foundational to their industrial uses. Pangenome analysis of publicly available Lactobacillaceae genomes allowed us to generate genome-scale metabolic network reconstructions for 26 species of industrial importance. Their manual curation led to more than 75,000 gene-protein-reaction associations that were deployed to generate 2,446 genome-scale metabolic models. Cross-referencing genomes and known metabolic traits allowed for manual metabolic network curation and validation of the metabolic models. As a result, we provide the first pangenomic basis for metabolism in the Lactobacillaceae family and a collection of predictive computational metabolic models that enable a variety of practical uses.IMPORTANCELactobacillaceae, a bacterial family foundational to a trillion-dollar industry, is increasingly relevant to biosustainability initiatives. Our study, leveraging approximately 2,400 genome sequences, provides a pangenomic analysis of Lactobacillaceae metabolism, creating over 2,400 curated and validated genome-scale models (GEMs). These GEMs successfully predict (i) unique, species-specific metabolic reactions; (ii) niche-enriched reactions that increase organism fitness; (iii) essential media components, offering insights into the global amino acid essentiality of Lactobacillaceae; and (iv) fermentation capabilities across the family, shedding light on the metabolic basis of Lactobacillaceae-based commercial products. This quantitative understanding of Lactobacillaceae metabolic properties and their genomic basis will have profound implications for the food industry and biosustainability, offering new insights and tools for strain selection and manipulation.

Meta-analysis Driven Strain Design for Mitigating Oxidative Stresses Important in Biomanufacturing.

As the availability of data sets increases, meta-analysis leveraging aggregated and interoperable data types is proving valuable. This study leveraged a meta-analysis workflow to identify mutations that could improve robustness to reactive oxygen species (ROS) stresses using an industrially important melatonin production strain as an example. ROS stresses often occur during cultivation and negatively affect strain performance. Cellular response to ROS is also linked to the SOS response and resistance to pH fluctuations, which is important to strain robustness in large-scale biomanufacturing. This work integrated more than 7000 E. coli adaptive laboratory evolution (ALE) mutations across 59 experiments to statistically associate mutated genes to 2 ROS tolerance ALE conditions from 72 unique conditions. Mutant oxyR, fur, iscR, and ygfZ were significantly associated and hypothesized to contribute fitness in ROS stress. Across these genes, 259 total mutations were inspected in conjunction with transcriptomics from 46 iModulon experiments. Ten mutations were chosen for reintroduction based on mutation clustering and coinciding transcriptional changes as evidence of fitness impact. Strains with mutations reintroduced into oxyR, fur, iscR, and ygfZ exhibited increased tolerance to H2O2 and acid stress and reduced SOS response, all of which are related to ROS. Additionally, new evidence was generated toward understanding the function of ygfZ, an uncharacterized gene. This meta-analysis approach utilized aggregated and interoperable multiomics data sets to identify mutations conferring industrially relevant phenotypes with the least drawbacks, describing an approach for data-driven strain engineering to optimize microbial cell factories.

Using the reconstructed genome-scale human metabolic network to study physiology and pathology.

Metabolism plays a key role in many major human diseases. Generation of high-throughput omics data has ushered in a new era of systems biology. Genome-scale metabolic network reconstructions provide a platform to interpret omics data in a biochemically meaningful manner. The release of the global human metabolic network, Recon 1, in 2007 has enabled new systems biology approaches to study human physiology, pathology and pharmacology. There are currently more than 20 publications that utilize Recon 1, including studies of cancer, diabetes, host-pathogen interactions, heritable metabolic disorders and off-target drug binding effects. In this mini-review, we focus on the reconstruction of the global human metabolic network and four classes of its application. We show that computational simulations for numerous pathologies have yielded clinically relevant results, many corroborated by existing or newly generated experimental data.

A variational principle for computing nonequilibrium fluxes and potentials in genome-scale biochemical networks.

We derive a convex optimization problem on a steady-state nonequilibrium network of biochemical reactions, with the property that energy conservation and the second law of thermodynamics both hold at the problem solution. This suggests a new variational principle for biochemical networks that can be implemented in a computationally tractable manner. We derive the Lagrange dual of the optimization problem and use strong duality to demonstrate that a biochemical analogue of Tellegen's theorem holds at optimality. Each optimal flux is dependent on a free parameter that we relate to an elementary kinetic parameter when mass action kinetics is assumed.

Characterization of metabolism in the Fe(III)-reducing organism Geobacter sulfurreducens by constraint-based modeling.

Geobacter sulfurreducens is a well-studied representative of the Geobacteraceae, which play a critical role in organic matter oxidation coupled to Fe(III) reduction, bioremediation of groundwater contaminated with organics or metals, and electricity production from waste organic matter. In order to investigate G. sulfurreducens central metabolism and electron transport, a metabolic model which integrated genome-based predictions with available genetic and physiological data was developed via the constraint-based modeling approach. Evaluation of the rates of proton production and consumption in the extracellular and cytoplasmic compartments revealed that energy conservation with extracellular electron acceptors, such as Fe(III), was limited relative to that associated with intracellular acceptors. This limitation was attributed to lack of cytoplasmic proton consumption during reduction of extracellular electron acceptors. Model-based analysis of the metabolic cost of producing an extracellular electron shuttle to promote electron transfer to insoluble Fe(III) oxides demonstrated why Geobacter species, which do not produce shuttles, have an energetic advantage over shuttle-producing Fe(III) reducers in subsurface environments. In silico analysis also revealed that the metabolic network of G. sulfurreducens could synthesize amino acids more efficiently than that of Escherichia coli due to the presence of a pyruvate-ferredoxin oxidoreductase, which catalyzes synthesis of pyruvate from acetate and carbon dioxide in a single step. In silico phenotypic analysis of deletion mutants demonstrated the capability of the model to explore the flexibility of G. sulfurreducens central metabolism and correctly predict mutant phenotypes. These results demonstrate that iterative modeling coupled with experimentation can accelerate the understanding of the physiology of poorly studied but environmentally relevant organisms and may help optimize their practical applications.

Properties of metabolic networks: structure versus function.

Biological data from high-throughput technologies describing the network components (genes, proteins, metabolites) and their associated interactions have driven the reconstruction and study of structural (topological) properties of large-scale biological networks. In this article, we address the relation of the functional and structural properties by using extensively experimentally validated genome-scale metabolic network models to compute observable functional states of a microorganism and compare the "structure versus function" attributes of metabolic networks. It is observed that, functionally speaking, the essentiality of reactions in a node is not correlated with node connectivity as structural analyses of other biological networks have suggested. These findings are illustrated with the analysis of the genome-scale biochemical networks of three species with distinct modes of metabolism. These results also suggest fundamental differences among different biological networks arising out of their representation and functional constraints.

Quantitative analysis of Escherichia coli metabolic phenotypes within the context of phenotypic phase planes.

In silico models of Escherichia coli metabolism have been developed to predict metabolic behavior and propose experimentally testable hypotheses. However, a thorough assessment of the metabolic phenotype requires well-designed experimentation and reproducible experimental techniques. A method for the quantitative analysis of E. coli metabolism in vivo within the framework of in silico phenotypic phase plane analysis is presented. Using this approach, we have quantitatively studied E. coli metabolism in various environmental conditions and nutritional media. Our experimental methodology, in combination with steady-state metabolic models, can be used to study biological properties and evaluate the metabolic capabilities of microbes.

Dynamic simulation of the human red blood cell metabolic network.

We have developed a Mathematica application package to perform dynamic simulations of the red blood cell (RBC) metabolic network. The package relies on, and integrates, many years of mathematical modeling and biochemical work on red blood cell metabolism. The extensive data regarding the red blood cell metabolic network and the previous kinetic analysis of all the individual components makes the human RBC an ideal 'model' system for mathematical metabolic models. The Mathematica package can be used to understand the dynamics and regulatory characteristics of the red blood cell.

Metabolic modeling of microbial strains in silico.

The large volume of genome-scale data that is being produced and made available in databases on the World Wide Web is demanding the development of integrated mathematical models of cellular processes. The analysis of reconstructed metabolic networks as systems leads to the development of an in silico or computer representation of collections of cellular metabolic constituents, their interactions and their integrated function as a whole. The use of quantitative analysis methods to generate testable hypotheses and drive experimentation at a whole-genome level signals the advent of a systemic modeling approach to cellular and molecular biology.

In silico predictions of Escherichia coli metabolic capabilities are consistent with experimental data.

A significant goal in the post-genome era is to relate the annotated genome sequence to the physiological functions of a cell. Working from the annotated genome sequence, as well as biochemical and physiological information, it is possible to reconstruct complete metabolic networks. Furthermore, computational methods have been developed to interpret and predict the optimal performance of a metabolic network under a range of growth conditions. We have tested the hypothesis that Escherichia coli uses its metabolism to grow at a maximal rate using the E. coli MG1655 metabolic reconstruction. Based on this hypothesis, we formulated experiments that describe the quantitative relationship between a primary carbon source (acetate or succinate) uptake rate, oxygen uptake rate, and maximal cellular growth rate. We found that the experimental data were consistent with the stated hypothesis, namely that the E. coli metabolic network is optimized to maximize growth under the experimental conditions considered. This study thus demonstrates how the combination of in silico and experimental biology can be used to obtain a quantitative genotype-phenotype relationship for metabolism in bacterial cells.

Flux-balance analysis of mitochondrial energy metabolism: consequences of systemic stoichiometric constraints.

Mitochondrial metabolism is a critical component in the functioning and maintenance of cellular organs. The stoichiometry of biochemical reaction networks imposes constraints on mitochondrial function. A modeling framework, flux-balance analysis (FBA), was used to characterize the optimal flux distributions for maximal ATP production in the mitochondrion. The model predicted the expected ATP yields for glucose, lactate, and palmitate. Genetic defects that affect mitochondrial functions have been implicated in several human diseases. FBA can characterize the metabolic behavior due to genetic deletions at the metabolic level, and the effect of mutations in the tricarboxylic acid (TCA) cycle on mitochondrial ATP production was simulated. The mitochondrial ATP production is severely affected by TCA-cycle mutations. In addition, the model predicts the secretion of TCA-cycle intermediates, which is observed in clinical studies of mitochondriopathies such as those associated with fumarase deficiency. The model provides a systemic perspective to characterize the effect of stoichiometric constraints and specific metabolic fluxes on mitochondrial function.

Robustness analysis of the Escherichia coli metabolic network.

Genomic, biochemical, and strain-specific data can be assembled to define an in silico representation of the metabolic network for a select group of single cellular organisms. Flux-balance analysis and phenotypic phase planes derived therefrom have been developed and applied to analyze the metabolic capabilities and characteristics of Escherichia coli K-12. These analyses have shown the existence of seven essential reactions in the central metabolic pathways (glycolysis, pentose phosphate pathway, tricarboxylic acid cycle) for the growth in glucose minimal media. The corresponding seven gene products can be grouped into three categories: (1) pentose phosphate pathway genes, (2) three-carbon glycolytic genes, and (3) tricarboxylic acid cycle genes. Here we develop a procedure that calculates the sensitivity of optimal cellular growth to altered flux levels of these essential gene products. The results indicate that the E. coli metabolic network is robust with respect to the flux levels of these enzymes. The metabolic flux in the transketolase and the tricarboxylic acid cycle reactions can be reduced to 15% and 19%, respectively, of the optimal value without significantly influencing the optimal growth flux. The metabolic network also exhibited robustness with respect to the ribose-5-phosphate isomerase, and the ribose-5-phosephate isomerase flux was reduced to 28% of the optimal value without significantly effecting the optimal growth flux. The metabolic network exhibited limited robustness to the three-carbon glycolytic fluxes both increased and decreased. The development presented another dimension to the use of FBA to study the capabilities of metabolic networks.

Metabolic flux balance analysis and the in silico analysis of Escherichia coli K-12 gene deletions.

BACKGROUND: Genome sequencing and bioinformatics are producing detailed lists of the molecular components contained in many prokaryotic organisms. From this 'parts catalogue' of a microbial cell, in silico representations of integrated metabolic functions can be constructed and analyzed using flux balance analysis (FBA). FBA is particularly well-suited to study metabolic networks based on genomic, biochemical, and strain specific information. RESULTS: Herein, we have utilized FBA to interpret and analyze the metabolic capabilities of Escherichia coli. We have computationally mapped the metabolic capabilities of E. coli using FBA and examined the optimal utilization of the E. coli metabolic pathways as a function of environmental variables. We have used an in silico analysis to identify seven gene products of central metabolism (glycolysis, pentose phosphate pathway, TCA cycle, electron transport system) essential for aerobic growth of E. coli on glucose minimal media, and 15 gene products essential for anaerobic growth on glucose minimal media. The in silico tpi-, zwf, and pta- mutant strains were examined in more detail by mapping the capabilities of these in silico isogenic strains. CONCLUSIONS: We found that computational models of E. coli metabolism based on physicochemical constraints can be used to interpret mutant behavior. These in silica results lead to a further understanding of the complex genotype-phenotype relation.

The Escherichia coli MG1655 in silico metabolic genotype: its definition, characteristics, and capabilities.

The Escherichia coli MG1655 genome has been completely sequenced. The annotated sequence, biochemical information, and other information were used to reconstruct the E. coli metabolic map. The stoichiometric coefficients for each metabolic enzyme in the E. coli metabolic map were assembled to construct a genome-specific stoichiometric matrix. The E. coli stoichiometric matrix was used to define the system's characteristics and the capabilities of E. coli metabolism. The effects of gene deletions in the central metabolic pathways on the ability of the in silico metabolic network to support growth were assessed, and the in silico predictions were compared with experimental observations. It was shown that based on stoichiometric and capacity constraints the in silico analysis was able to qualitatively predict the growth potential of mutant strains in 86% of the cases examined. Herein, it is demonstrated that the synthesis of in silico metabolic genotypes based on genomic, biochemical, and strain-specific information is possible, and that systems analysis methods are available to analyze and interpret the metabolic phenotype.

Extension of osmolality-induced podia is observed from fluorescently labeled hematopoietic cell lines in hyperosmotic medium.

BACKGROUND: Since the description of long podia extended by hematopoietic cells and cell lines, the reliable elicitation of podia extensions is needed to study these podia systemically. In this study, hyperosmotic stress was considered as an elicitor. METHODS: Using two fluorescent membrane dyes PKH2 and PKH26, and an automated fluorescence microscopy system, morphological changes of seven human cell lines (six hematopoietic, one fibrosarcoma) at different osmolalities were monitored. Presence of surface molecules on the hyperosmolality-induced podia (osmopodia) was examined. RESULTS: In hyperosmotic medium, cells shrank rapidly, followed by osmopodia extension. Cells exhibited variable number (up to five) and length (up to longer than 100 microm) of osmopodia in about 1 h. Dead cells did not extend podia. Frequency, length, and number of podia were variable among cell lines studied. CD44 and CD45 were not present on the osmopodia, although they were present on the cell surface, showing that osmopodia characteristics differ from the podia observed previously in isotonic media. The osmopodia extension process was shown to be reversible upon repeated osmolality changes. CONCLUSIONS: Osmopodia extended by human hematopoietic cell lines display a newly observed cellular morphology and provide a tool for investigation of dynamic cellular response to environmental changes.

Intercellular adhesion can be visualized using fluorescently labeled fibrosarcoma HT1080 cells cocultured with hematopoietic cell lines or CD34(+) enriched human mobilized peripheral blood cells.

BACKGROUND: Intercellular contacts between adjacent cells migrating over each other are important in many cellular processes. However, it has been difficult to visualize and identify dynamic intercellular adhesions between migrating cells in situ. METHODS: Two fluorescent membrane dyes, PKH2 and PKH26 for staining HT1080 and hematopoietic cells and cell lines, and an automated fluorescence microscopy system were used to monitor intercellular adhesion. RESULTS: Cellular extensions connecting two or more adjacent cells were visualized, showing the intercellular adhesion between migrating cells for minutes and up to hours. After cells adhered to each other, followed by cell migration in different directions, cellular extensions were dragged from the pivotal contact points in different focal planes. CD34(+)-enriched mobilized peripheral blood cells and six hematopoietic cell lines showed intercellular connections in cocultures with HT1080. However, the frequency of intercellular connections was variable in different cocultures. A cell density of about 3.1 x 10(4) cells/cm(2) for both cell lines in cocultures provided an adequate number of cells in each field of view, showing up to four intercellular connections per 100 total cells plated. DISCUSSION: The tools derived from this study will open new areas of investigation for understanding the mechanism of the intercellular adhesion process.

Theory for the systemic definition of metabolic pathways and their use in interpreting metabolic function from a pathway-oriented perspective.

Cellular metabolism is most often described and interpreted in terms of the biochemical reactions that make up the metabolic network. Genomics is providing near complete information regarding the genes/gene products participating in cellular metabolism for a growing number of organisms. As the true functional units of metabolic systems are its pathways, the time has arrived to define metabolic pathways in the context of whole-cell metabolism for the analysis of the structural design and capabilities of the metabolic network. In this study, we present the theoretical foundations for the identification of the unique set of systemically independent biochemical pathways, termed extreme pathways, based on system stochiometry and limited thermodynamics. These pathways represent the edges of the steady-state flux cone derived from convex analysis, and they can be used to represent any flux distribution achievable by the metabolic network. An algorithm is presented to determine the set of extreme pathways for a system of any complexity and a classification scheme is introduced for the characterization of these pathways. The property of systemic independence is discussed along with its implications for issues related to metabolic regulation and the evolution of cellular metabolic networks. The underlying pathway structure that is determined from the set of extreme pathways now provides us with the ability to analyse, interpret, and perhaps predict metabolic function from a pathway-based perspective in addition to the traditional reaction-based perspective. The algorithm and classification scheme developed can be used to describe the pathway structure in annotated genomes to explore the capabilities of an organism.

Assessment of the metabolic capabilities of Haemophilus influenzae Rd through a genome-scale pathway analysis.

The annotated full DNA sequence is becoming available for a growing number of organisms. This information along with additional biochemical and strain-specific data can be used to define metabolic genotypes and reconstruct cellular metabolic networks. The first free-living organism for which the entire genomic sequence was established was Haemophilus influenzae. Its metabolic network is reconstructed herein and contains 461 reactions operating on 367 intracellular and 84 extracellular metabolites. With the metabolic reaction network established, it becomes necessary to determine its underlying pathway structure as defined by the set of extreme pathways. The H. influenzae metabolic network was subdivided into six subsystems and the extreme pathways determined for each subsystem based on stoichiometric, thermodynamic, and systems-specific constraints. Positive linear combinations of these pathways can be taken to determine the extreme pathways for the complete system. Since these pathways span the capabilities of the full system, they could be used to address a number of important physiological questions. First, they were used to reconcile and curate the sequence annotation by identifying reactions whose function was not supported in any of the extreme pathways. Second, they were used to predict gene products that should be co-regulated and perhaps co-expressed. Third, they were used to determine the composition of the minimal substrate requirements needed to support the production of 51 required metabolic products such as amino acids, nucleotides, phospholipids, etc. Fourth, sets of critical gene deletions from core metabolism were determined in the presence of the minimal substrate conditions and in more complete conditions reflecting the environmental niche of H. influenzae in the human host. In the former case, 11 genes were determined to be critical while six remained critical under the latter conditions. This study represents an important milestone in theoretical biology, namely the establishment of the first extreme pathway structure of a whole genome.