Topology Of Metabolic Networks Research Articles

Teasing out missing reactions in genome-scale metabolic networks through hypergraph learning

GEnome-scale Metabolic models (GEMs) are powerful tools to predict cellular metabolism and physiological states in living organisms. However, due to our imperfect knowledge of metabolic processes, even highly curated GEMs have knowledge gaps (e.g., missing reactions). Existing gap-filling methods typically require phenotypic data as input to tease out missing reactions. We still lack a computational method for rapid and accurate gap-filling of metabolic networks before experimental data is available. Here we present a deep learning-based method — CHEbyshev Spectral HyperlInk pREdictor (CHESHIRE) — to predict missing reactions in GEMs purely from metabolic network topology. We demonstrate that CHESHIRE outperforms other topology-based methods in predicting artificially removed reactions over 926 high- and intermediate-quality GEMs. Furthermore, CHESHIRE is able to improve the phenotypic predictions of 49 draft GEMs for fermentation products and amino acids secretions. Both types of validation suggest that CHESHIRE is a powerful tool for GEM curation to reveal unknown links between reactions and observed metabolic phenotypes.

The topology of genome-scale metabolic reconstructions unravels independent modules and high network flexibility.

The topology of metabolic networks is recognisably modular with modules weakly connected apart from sharing a pool of currency metabolites. Here, we defined modules as sets of reversible reactions isolated from the rest of metabolism by irreversible reactions except for the exchange of currency metabolites. Our approach identifies topologically independent modules under specific conditions associated with different metabolic functions. As case studies, the E.coli iJO1366 and Human Recon 2.2 genome-scale metabolic models were split in 103 and 321 modules respectively, displaying significant correlation patterns in expression data. Finally, we addressed a fundamental question about the metabolic flexibility conferred by reversible reactions: "Of all Directed Topologies (DTs) defined by fixing directions to all reversible reactions, how many are capable of carrying flux through all reactions?". Enumeration of the DTs for iJO1366 model was performed using an efficient depth-first search algorithm, rejecting infeasible DTs based on mass-imbalanced and loopy flux patterns. We found the direction of 79% of reversible reactions must be defined before all directions in the network can be fixed, granting a high degree of flexibility.

Shiny GATOM: omics-based identification of regulated metabolic modules in atom transition networks.

Multiple high-throughput omics techniques provide different angles on systematically quantifying and studying metabolic regulation of cellular processes. However, an unbiased analysis of such data and, in particular, integration of multiple types of data remains a challenge. Previously, for this purpose we developed GAM web-service for integrative metabolic network analysis. Here we describe an updated pipeline GATOM and the corresponding web-service Shiny GATOM, which takes as input transcriptional and/or metabolomic data and finds a metabolic subnetwork most regulated between the two conditions of interest. GATOM features a new metabolic network topology based on atom transition, which significantly improves interpretability of the analysis results. To address computational challenges arising with the new network topology, we introduce a new variant of the maximum weight connected subgraph problem and provide a corresponding exact solver. To make the used networks up-to-date we upgraded the KEGG-based network construction pipeline and developed one based on the Rhea database, which allows analysis of lipidomics data. Finally, we simplified local installation, providing R package mwcsr for solving relevant graph optimization problems and R package gatom, which implements the GATOM pipeline. The web-service is available at https://ctlab.itmo.ru/shiny/gatom and https://artyomovlab.wustl.edu/shiny/gatom.

SCOUR: a stepwise machine learning framework for predicting metabolite-dependent regulatory interactions

BackgroundThe topology of metabolic networks is both well-studied and remarkably well-conserved across many species. The regulation of these networks, however, is much more poorly characterized, though it is known to be divergent across organisms—two characteristics that make it difficult to model metabolic networks accurately. While many computational methods have been built to unravel transcriptional regulation, there have been few approaches developed for systems-scale analysis and study of metabolic regulation. Here, we present a stepwise machine learning framework that applies established algorithms to identify regulatory interactions in metabolic systems based on metabolic data: stepwise classification of unknown regulation, or SCOUR.ResultsWe evaluated our framework on both noiseless and noisy data, using several models of varying sizes and topologies to show that our approach is generalizable. We found that, when testing on data under the most realistic conditions (low sampling frequency and high noise), SCOUR could identify reaction fluxes controlled only by the concentration of a single metabolite (its primary substrate) with high accuracy. The positive predictive value (PPV) for identifying reactions controlled by the concentration of two metabolites ranged from 32 to 88% for noiseless data, 9.2 to 49% for either low sampling frequency/low noise or high sampling frequency/high noise data, and 6.6–27% for low sampling frequency/high noise data, with results typically sufficiently high for lab validation to be a practical endeavor. While the PPVs for reactions controlled by three metabolites were lower, they were still in most cases significantly better than random classification.ConclusionsSCOUR uses a novel approach to synthetically generate the training data needed to identify regulators of reaction fluxes in a given metabolic system, enabling metabolomics and fluxomics data to be leveraged for regulatory structure inference. By identifying and triaging the most likely candidate regulatory interactions, SCOUR can drastically reduce the amount of time needed to identify and experimentally validate metabolic regulatory interactions. As high-throughput experimental methods for testing these interactions are further developed, SCOUR will provide critical impact in the development of predictive metabolic models in new organisms and pathways.

IsoSolve: An Integrative Framework to Improve Isotopic Coverage and Consolidate Isotopic Measurements by Mass Spectrometry and/or Nuclear Magnetic Resonance.

Stable-isotope labeling experiments are widely used to investigate the topology and functioning of metabolic networks. Label incorporation into metabolites can be quantified using a broad range of mass spectrometry (MS) and nuclear magnetic resonance (NMR) spectroscopy methods, but in general, no single approach can completely cover isotopic space, even for small metabolites. The number of quantifiable isotopic species could be increased and the coverage of isotopic space improved by integrating measurements obtained by different methods; however, this approach has remained largely unexplored because no framework able to deal with partial, heterogeneous isotopic measurements has yet been developed. Here, we present a generic computational framework based on symbolic calculus that can integrate any isotopic data set by connecting measurements to the chemical structure of the molecules. As a test case, we apply this framework to isotopic analyses of amino acids, which are ubiquitous to life, central to many biological questions, and can be analyzed by a broad range of MS and NMR methods. We demonstrate how this integrative framework helps to (i) clarify and improve the coverage of isotopic space, (ii) evaluate the complementarity and redundancy of different techniques, (iii) consolidate isotopic data sets, (iv) design experiments, and (v) guide future analytical developments. This framework, which can be applied to any labeled element, isotopic tracer, metabolite, and analytical platform, has been implemented in IsoSolve (available at https://github.com/MetaSys-LISBP/IsoSolve and https://pypi.org/project/IsoSolve), an open-source software that can be readily integrated into data analysis pipelines.

Metabolic Network Topology of Alzheimer's Disease and Dementia with Lewy Bodies Generated Using Fluorodeoxyglucose Positron Emission Tomography.

Background:Alzheimer’s disease (AD) and dementia with Lewy bodies (DLB) are often misdiagnosed with each other because of similar symptoms including progressive memory loss. The metabolic network topology that describes inter-regional metabolic connections can be generated using fluorodeoxyglucose positron emission tomography (FDG-PET) data with the graph-theoretical method. We hypothesized that different metabolic connectivity underlies the symptoms of AD patients, DLB patients, and cognitively normal (CN) individuals.Objective:This study aimed to generate metabolic connectivity using FDG-PET data and assess the network topology to differentiate AD patients, DLB patients, and CN individuals.Methods:This study included 45 AD patients, 18 DLB patients, and 142 CN controls. We analyzed FDG-PET data using the graph-theoretical method and generated the network topology in AD patients, DLB patients, and CN individuals. We statistically assessed the topology with global and nodal parameters.Results:The whole metabolic network was preserved in CN; however, diffusely decreased connection was found in AD and partially but more deeply decreased connection was observed in DLB. The metabolic topology revealed that the right posterior cingulate and the left transverse temporal gyrus were significantly different between AD and DLB.Conclusion:The present findings indicate that metabolic connectivity decreased in both AD and DLB, compared with CN. DLB was characterized restricted but deeper stereotyped network disruption compared with AD. The right posterior cingulate and the left transverse temporal gyrus are significant regions in the metabolic connectivity for differentiating AD from DLB.

13C metabolic flux analysis on roles of malate transporter in lipid accumulation of Mucor circinelloides

BackgroundMitochondrial and cytoplasmic malate transporter proteins are responsible for transmembrane transport of malate, thereby linking malate metabolism in various subcellular regions of the cell. These transporters play an important role in fatty acid biosynthesis of oleaginous microorganisms. Our previous studies have found that lipid content of the recombinant Mucor circinelloides (M. circinelloides) strain with mitochondrial malate transporter (mt) gene overexpression was increased by 70%, while that of strain with mt gene knockout was decreased by 27%. However, the mechanism of malate transporter promoting the transport of mitochondrial malate and citrate related to lipid accumulation is not clear. Therefore, 13C-labeled glucose metabolic flux analysis was carried out to identify the metabolic network topology and estimate intracellular fluxes of genetically engineered M. circinelloides strains for the purpose of better understanding the roles of malate transporters in citrate transport systems and lipid accumulation.ResultsThe metabolic flux distribution analysis suggested that tricarboxylic acid (TCA) cycle flux ratio of mt-overexpression strains was decreased compared to that of the control strain, but in contrast, glyoxylic acid (GOX) cycle flux ratio was increased. Accordingly, the mt-knockout strain showed an opposite phenomenon with a higher TCA cycle flux ratio and a lower GOX cycle flux ratio than the control strain. GOX cycle might be more effective than TCA cycle in producing malate and oxaloacetate replenishment. Moreover, a relatively higher flux ratio of the pentose phosphate (PP) pathway was obtained in mt-overexpression strains, but no significant difference in the malic enzyme flux between recombinant strains and the control strain. Our results confirmed that PP pathway might play an important role for supplying NADPH and malic enzyme is not a limiting factor for fatty acid synthesis in oleaginous fungus M. circinelloides strains.ConclusionIntracellular metabolic flux information suggested that mt-overexpression strains had higher flux in PP pathway and GOX cycle, lower flux in TCA cycle, and no difference in malic enzyme cycle. Together, the role of malate transporter was assumed to further participate in transporting cycle of acetyl-CoA and drive PP pathway to supply NADPH required for lipid accumulation in recombinant M. circinelloides strains.

Topological assessment of metabolic networks reveals evolutionary information

Evolutionary information was inferred from the topology of metabolic networks corresponding to 17 plant species belonging to major plant lineages Chlorophytes, Bryophytes, Lycophytes and Angiosperms. The plant metabolic networks were built using the substrate-product network modeling based on the metabolic reactions available on the PlantCyc database (version 9.5), from which their local topological properties such as degree, in-degree, out-degree, clustering coefficient, hub-score, authority-score, local efficiency, betweenness and eigencentrality were measured. The topological measurements corresponding to each metabolite within the networks were considered as a set of metabolic characters to compound a feature vector representing each plant. Our results revealed that some local topological characters are able to discern among plant kinships, since similar phylogenies were found when comparing dendrograms obtained by topological metrics to the one obtained by DNA sequences of chloroplast genes. Furthermore, we also found that even a smaller number of metabolic characters is able to separate among major clades with high bootstrap support (BS > 95), while for some suborders a bigger content has been required.

Methodology for the Validation of Isotopic Analyses by Mass Spectrometry in Stable-Isotope Labeling Experiments.

Stable-isotope labeling experiments (ILEs) are widely used to investigate the topology and operation of metabolic networks. The quality of isotopic data collected in ILEs is of utmost importance to ensure reliable biological interpretations, but current evaluation approaches are limited due to a lack of suitable reference material and relevant evaluation criteria. In this work, we present a complete methodology to evaluate mass spectrometry (MS) methods used for quantitative isotopic studies of metabolic systems. This methodology, based on a biological sample containing metabolites with controlled labeling patterns, exploits different quality metrics specific to isotopic analyses (accuracy and precision of isotopologue masses, abundances, and mass shifts and isotopic working range). We applied this methodology to evaluate a novel LC-MS method for the analysis of amino acids, which was tested on high resolution (Orbitrap operating in full scan mode) and low resolution (triple quadrupole operating in multiple reaction monitoring mode) mass spectrometers. Results show excellent accuracy and precision over a large working range and revealed matrix-specific as well as mode-specific characteristics. The proposed methodology can identify reliable (and unreliable) isotopic data in an easy and straightforward way and efficiently supports the identification of sources of systematic biases as well as of the main factors that influence the overall accuracy and precision of measurements. This approach is generic and can be used to validate isotopic analyses on different matrices, analytical platforms, labeled elements, or classes of metabolites. It is expected to strengthen the reliability of isotopic measurements and thereby the biological value of ILEs.

SynBioEcoli: a comprehensive metabolism network of engineered E. coli in three dimensional visualization

BackgroundA comprehensive metabolism network of engineered E. coli is very important in systems biology and metabolomics studies. Many tools focus on two‐dimensional space to display pathways in metabolic network. However, the usage of three‐dimensional visualization may help to understand better the intricate topology of metabolic and regulatory networks.MethodsWe manually curated large amount of experimental data (including pathways, reactions and metabolites) from literature related with different types of engineered E. coli and then utilized a novel technology of three dimensional visualization to develop a comprehensive metabolic network named SynBioEcoli.ResultsSynBioEoli contains 740 biosynthetic pathways, 3,889 metabolic reactions, 2,255 chemical compoundsmanually curated from about 11,000 metabolism publications related with different types of engineered E. coli. Furthermore, SynBioEcoli integrates with various informatics techniques.ConclusionsSynBioEcoli could be regarded as a comprehensive knowledgebase of engineered E. coli and represents the next generation cellular metabolism network visualization technology. It could be accessed via web browsers (such as Google Chrome) supporting WebGL, at http://www.rxnfinder.org/synbioecoli/.

Cutting the Gordian Knot: Identifiability of anaplerotic reactions in Corynebacterium glutamicum by means of (13) C-metabolic flux analysis.

Corynebacterium glutamicum is the major workhorse for the microbial production of several amino and organic acids. As long as these derive from tricarboxylic acid cycle intermediates, the activity of anaplerotic reactions is pivotal for a high biosynthetic yield. To determine single anaplerotic activities (13) C-Metabolic Flux Analysis ((13) C-MFA) has been extensively used for C. glutamicum, however with different network topologies, inconsistent or poorly determined anaplerotic reaction rates. Therefore, in this study we set out to investigate whether a focused isotopomer model of the anaplerotic node can at all admit a unique solution for all fluxes. By analyzing different scenarios of active anaplerotic reactions, we show in full generality that for C. glutamicum only certain anaplerotic deletion mutants allow to uniquely determine the anaplerotic fluxes from (13) C-isotopomer data. We stress that the result of this analysis for different assumptions on active enzymes is directly transferable to other compartment-free organisms. Our results demonstrate that there exist biologically relevant metabolic network topologies for which the flux distribution cannot be inferred by classical (13) C-MFA.

13C-metabolic flux analysis of lipid accumulation in the oleaginous fungus Mucor circinelloides

The oleaginous fungus Mucor circinelloides is of industrial interest because it can produce high levels of polyunsaturated fatty acid γ-linolenic acid. M. circinelloides CBS 277.49 is able to accumulate less than 15% of cell dry weight as lipids, while M. circinelloides WJ11 can accumulate lipid up to 36%. In order to better understand the mechanisms behind the differential lipid accumulation in these two strains, tracer experiments with 13C-glucose were performed with the growth of M. circinelloides and subsequent gas chromatography–mass spectrometric detection of 13C-patterns in proteinogenic amino acids was carried out to identify the metabolic network topology and estimate intracellular fluxes. Our results showed that the high oleaginous strain WJ11 had higher flux of pentose phosphate pathway and malic enzyme, lower flux in tricarboxylic acid cycle, higher flux in glyoxylate cycle and ATP: citrate lyase, together, it might provide more NADPH and substrate acetyl-CoA for fatty acid synthesis.

NetCooperate: a network-based tool for inferring host-microbe and microbe-microbe cooperation.

BackgroundHost-microbe and microbe-microbe interactions are often governed by the complex exchange of metabolites. Such interactions play a key role in determining the way pathogenic and commensal species impact their host and in the assembly of complex microbial communities. Recently, several studies have demonstrated how such interactions are reflected in the organization of the metabolic networks of the interacting species, and introduced various graph theory-based methods to predict host-microbe and microbe-microbe interactions directly from network topology. Using these methods, such studies have revealed evolutionary and ecological processes that shape species interactions and community assembly, highlighting the potential of this reverse-ecology research paradigm.ResultsNetCooperate is a web-based tool and a software package for determining host-microbe and microbe-microbe cooperative potential. It specifically calculates two previously developed and validated metrics for species interaction: the Biosynthetic Support Score which quantifies the ability of a host species to supply the nutritional requirements of a parasitic or a commensal species, and the Metabolic Complementarity Index which quantifies the complementarity of a pair of microbial organisms’ niches. NetCooperate takes as input a pair of metabolic networks, and returns the pairwise metrics as well as a list of potential syntrophic metabolic compounds.ConclusionsThe Biosynthetic Support Score and Metabolic Complementarity Index provide insight into host-microbe and microbe-microbe metabolic interactions. NetCooperate determines these interaction indices from metabolic network topology, and can be used for small- or large-scale analyses. NetCooperate is provided as both a web-based tool and an open-source Python module; both are freely available online at http://elbo.gs.washington.edu/software_netcooperate.html.

Predicting whole genome protein interaction networks from primary sequence data in model and non-model organisms using ENTS

BackgroundThe large-scale identification of physical protein-protein interactions (PPIs) is an important step toward understanding how biological networks evolve and generate emergent phenotypes. However, experimental identification of PPIs is a laborious and error-prone process, and current methods of PPI prediction tend to be highly conservative or require large amounts of functional data that may not be available for newly-sequenced organisms.ResultsIn this study we demonstrate a random-forest based technique, ENTS, for the computational prediction of protein-protein interactions based only on primary sequence data. Our approach is able to efficiently predict interactions on a whole-genome scale for any eukaryotic organism, using pairwise combinations of conserved domains and predicted subcellular localization of proteins as input features. We present the first predicted interactome for the forest tree Populus trichocarpa in addition to the predicted interactomes for Saccharomyces cerevisiae, Homo sapiens, Mus musculus, and Arabidopsis thaliana. Comparing our approach to other PPI predictors, we find that ENTS performs comparably to or better than a number of existing approaches, including several that utilize a variety of functional information for their predictions. We also find that the predicted interactions are biologically meaningful, as indicated by similarity in functional annotations and enrichment of co-expressed genes in public microarray datasets. Furthermore, we demonstrate some of the biological insights that can be gained from these predicted interaction networks. We show that the predicted interactions yield informative groupings of P. trichocarpa metabolic pathways, literature-supported associations among human disease states, and theory-supported insight into the evolutionary dynamics of duplicated genes in paleopolyploid plants.ConclusionWe conclude that the ENTS classifier will be a valuable tool for the de novo annotation of genome sequences, providing initial clues about regulatory and metabolic network topology, and revealing relationships that are not immediately obvious from traditional homology-based annotations.

Novel approach to engineer strains for simultaneous sugar utilization

Use of lignocellulosic biomass as a second generation feedstock in the biofuels industry is a pressing challenge. Among other difficulties in using lignocellulosic biomass, one major challenge is the optimal utilization of both 6-carbon (glucose) and 5-carbon (xylose) sugars by industrial microorganisms. Most industrial microorganisms preferentially utilize glucose over xylose owing to the regulatory phenomenon of carbon catabolite repression (CCR). Microorganisms that can co-utilize glucose and xylose are of considerable interest to the biofuels industry due to their ability to simplify the fermentation processes. However, elimination of CCR in microorganisms is challenging due to the multiple coordinating mechanisms involved. We report a novel algorithm, SIMUP, which finds metabolic engineering strategies to force co-utilization of two sugars, without targeting the regulatory pathways of CCR. Mutants of Escherichia coli based on SIMUP algorithm showed predicted growth phenotypes and co-utilized glucose and xylose; however, consumed the sugars slower than the wild-type. Some solutions identified by the algorithm were based on stoichiometric imbalance and were not obvious from the metabolic network topology. Furthermore, sequencing studies on the genes involved in CCR showed that the mechanism for co-utilization of the sugars could be different from previously known mechanisms.

MultiMetEval: Comparative and Multi-Objective Analysis of Genome-Scale Metabolic Models

Comparative metabolic modelling is emerging as a novel field, supported by the development of reliable and standardized approaches for constructing genome-scale metabolic models in high throughput. New software solutions are needed to allow efficient comparative analysis of multiple models in the context of multiple cellular objectives. Here, we present the user-friendly software framework Multi-Metabolic Evaluator (MultiMetEval), built upon SurreyFBA, which allows the user to compose collections of metabolic models that together can be subjected to flux balance analysis. Additionally, MultiMetEval implements functionalities for multi-objective analysis by calculating the Pareto front between two cellular objectives. Using a previously generated dataset of 38 actinobacterial genome-scale metabolic models, we show how these approaches can lead to exciting novel insights. Firstly, after incorporating several pathways for the biosynthesis of natural products into each of these models, comparative flux balance analysis predicted that species like Streptomyces that harbour the highest diversity of secondary metabolite biosynthetic gene clusters in their genomes do not necessarily have the metabolic network topology most suitable for compound overproduction. Secondly, multi-objective analysis of biomass production and natural product biosynthesis in these actinobacteria shows that the well-studied occurrence of discrete metabolic switches during the change of cellular objectives is inherent to their metabolic network architecture. Comparative and multi-objective modelling can lead to insights that could not be obtained by normal flux balance analyses. MultiMetEval provides a powerful platform that makes these analyses straightforward for biologists. Sources and binaries of MultiMetEval are freely available from https://github.com/PiotrZakrzewski/MetEval/downloads.

Reconstruction of genome-scale metabolic models for 126 human tissues using mCADRE

BackgroundHuman tissues perform diverse metabolic functions. Mapping out these tissue-specific functions in genome-scale models will advance our understanding of the metabolic basis of various physiological and pathological processes. The global knowledgebase of metabolic functions categorized for the human genome (Human Recon 1) coupled with abundant high-throughput data now makes possible the reconstruction of tissue-specific metabolic models. However, the number of available tissue-specific models remains incomplete compared with the large diversity of human tissues.ResultsWe developed a method called metabolic Context-specificity Assessed by Deterministic Reaction Evaluation (mCADRE). mCADRE is able to infer a tissue-specific network based on gene expression data and metabolic network topology, along with evaluation of functional capabilities during model building. mCADRE produces models with similar or better functionality and achieves dramatic computational speed up over existing methods. Using our method, we reconstructed draft genome-scale metabolic models for 126 human tissue and cell types. Among these, there are models for 26 tumor tissues along with their normal counterparts, and 30 different brain tissues. We performed pathway-level analyses of this large collection of tissue-specific models and identified the eicosanoid metabolic pathway, especially reactions catalyzing the production of leukotrienes from arachidnoic acid, as potential drug targets that selectively affect tumor tissues.ConclusionsThis large collection of 126 genome-scale draft metabolic models provides a useful resource for studying the metabolic basis for a variety of human diseases across many tissues. The functionality of the resulting models and the fast computational speed of the mCADRE algorithm make it a useful tool to build and update tissue-specific metabolic models.

Random phenotypic variation of yeast (Saccharomyces cerevisiae) single-gene knockouts fits a double pareto-lognormal distribution.

BackgroundDistributed robustness is thought to influence the buffering of random phenotypic variation through the scale-free topology of gene regulatory, metabolic, and protein-protein interaction networks. If this hypothesis is true, then the phenotypic response to the perturbation of particular nodes in such a network should be proportional to the number of links those nodes make with neighboring nodes. This suggests a probability distribution approximating an inverse power-law of random phenotypic variation. Zero phenotypic variation, however, is impossible, because random molecular and cellular processes are essential to normal development. Consequently, a more realistic distribution should have a y-intercept close to zero in the lower tail, a mode greater than zero, and a long (fat) upper tail. The double Pareto-lognormal (DPLN) distribution is an ideal candidate distribution. It consists of a mixture of a lognormal body and upper and lower power-law tails.Objective and MethodsIf our assumptions are true, the DPLN distribution should provide a better fit to random phenotypic variation in a large series of single-gene knockout lines than other skewed or symmetrical distributions. We fit a large published data set of single-gene knockout lines in Saccharomyces cerevisiae to seven different probability distributions: DPLN, right Pareto-lognormal (RPLN), left Pareto-lognormal (LPLN), normal, lognormal, exponential, and Pareto. The best model was judged by the Akaike Information Criterion (AIC).ResultsPhenotypic variation among gene knockouts in S. cerevisiae fits a double Pareto-lognormal (DPLN) distribution better than any of the alternative distributions, including the right Pareto-lognormal and lognormal distributions.Conclusions and SignificanceA DPLN distribution is consistent with the hypothesis that developmental stability is mediated, in part, by distributed robustness, the resilience of gene regulatory, metabolic, and protein-protein interaction networks. Alternatively, multiplicative cell growth, and the mixing of lognormal distributions having different variances, may generate a DPLN distribution.

NetSeed: a network-based reverse-ecology tool for calculating the metabolic interface of an organism with its environment

NetSeed is a web tool and Perl module for analyzing the topology of metabolic networks and calculating the set of exogenously acquired compounds. NetSeed is based on the seed detection algorithm, developed and validated in previous studies. The NetSeed web-based tool, open-source Perl module, examples and documentation are freely available online at: http://depts.washington.edu/elbogs/NetSeed.

A compendium of inborn errors of metabolism mapped onto the human metabolic network

Inborn errors of metabolism (IEMs) are hereditary metabolic defects, which are encountered in almost all major metabolic pathways occurring in man. Many IEMs are screened for in neonates through metabolomic analysis of dried blood spot samples. To enable the mapping of these metabolomic data onto the published human metabolic reconstruction, we added missing reactions and pathways involved in acylcarnitine (AC) and fatty acid oxidation (FAO) metabolism. Using literary data, we reconstructed an AC/FAO module consisting of 352 reactions and 139 metabolites. When this module was combined with the human metabolic reconstruction, the synthesis of 39 acylcarnitines and 22 amino acids, which are routinely measured, was captured and 235 distinct IEMs could be mapped. We collected phenotypic and clinical features for each IEM enabling comprehensive classification. We found that carbohydrate, amino acid, and lipid metabolism were most affected by the IEMs, while the brain was the most commonly affected organ. Furthermore, we analyzed the IEMs in the context of metabolic network topology to gain insight into common features between metabolically connected IEMs. While many known examples were identified, we discovered some surprising IEM pairs that shared reactions as well as clinical features but not necessarily causal genes. Moreover, we could also re-confirm that acetyl-CoA acts as a central metabolite. This network based analysis leads to further insight of hot spots in human metabolism with respect to IEMs. The presented comprehensive knowledge base of IEMs will provide a valuable tool in studying metabolic changes involved in inherited metabolic diseases.