Dataset Information

Analysis of complex metabolic behavior through pathway decomposition.

ABSTRACT:

Background

Understanding complex systems through decomposition into simple interacting components is a pervasive paradigm throughout modern science and engineering. For cellular metabolism, complexity can be reduced by decomposition into pathways with particular biochemical functions, and the concept of elementary flux modes provides a systematic way for organizing metabolic networks into such pathways. While decomposition using elementary flux modes has proven to be a powerful tool for understanding and manipulating cellular metabolism, its utility, however, is severely limited since the number of modes in a network increases exponentially with its size.

Results

Here, we present a new method for decomposition of metabolic flux distributions into elementary flux modes. Our method can easily operate on large, genome-scale networks since it does not require all relevant modes of the metabolic network to be generated. We illustrate the utility of our method for metabolic engineering of Escherichia coli and for understanding the survival of Mycobacterium tuberculosis (MTB) during infection.

Conclusions

Our method can achieve computational time improvements exceeding 2000-fold and requires only several seconds to generate elementary mode decompositions on genome-scale networks. These improvements arise from not having to generate all relevant elementary modes prior to initiating the decomposition. The decompositions from our method are useful for understanding complex flux distributions and debugging genome-scale models.

SUBMITTER: Ip K

PROVIDER: S-EPMC3135541 | biostudies-literature | 2011 Jun

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Analysis of complex metabolic behavior through pathway decomposition.

Ip Kuhn K Colijn Caroline C Lun Desmond S DS

BMC systems biology 20110603

<h4>Background</h4>Understanding complex systems through decomposition into simple interacting components is a pervasive paradigm throughout modern science and engineering. For cellular metabolism, complexity can be reduced by decomposition into pathways with particular biochemical functions, and the concept of elementary flux modes provides a systematic way for organizing metabolic networks into such pathways. While decomposition using elementary flux modes has proven to be a powerful tool for ...[more]

PMID: 21639889

Similar Datasets

Project description:The quantification of experimentally-induced alterations in biological pathways remains a major challenge in systems biology. One example of this is the quantitative characterization of alterations in defined, established metabolic pathways from complex metabolomic data. At present, the disruption of a given metabolic pathway is inferred from metabolomic data by observing an alteration in the level of one or more individual metabolites present within that pathway. Not only is this approach open to subjectivity, as metabolites participate in multiple pathways, but it also ignores useful information available through the pairwise correlations between metabolites. This extra information may be incorporated using a higher-level approach that looks for alterations between a pair of correlation networks. In this way experimentally-induced alterations in metabolic pathways can be quantitatively defined by characterizing group differences in metabolite clustering. Taking this approach increases the objectivity of interpreting alterations in metabolic pathways from metabolomic data.We present and justify a new technique for comparing pairs of networks--in our case these networks are based on the same set of nodes and there are two distinct types of weighted edges. The algorithm is based on the Generalized Singular Value Decomposition (GSVD), which may be regarded as an extension of Principle Components Analysis to the case of two data sets. We show how the GSVD can be interpreted as a technique for reordering the two networks in order to reveal clusters that are exclusive to only one. Here we apply this algorithm to a new set of metabolomic data from the prefrontal cortex (PFC) of a translational model relevant to schizophrenia, rats treated subchronically with the N-methyl-D-Aspartic acid (NMDA) receptor antagonist phencyclidine (PCP). This provides us with a means to quantify which predefined metabolic pathways (Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolite pathway database) were altered in the PFC of PCP-treated rats. Several significant changes were discovered, notably: 1) neuroactive ligands active at glutamate and GABA receptors are disrupted in the PFC of PCP-treated animals, 2) glutamate dysfunction in these animals was not limited to compromised glutamatergic neurotransmission but also involves the disruption of metabolic pathways linked to glutamate; and 3) a specific series of purine reactions Xanthine ? Hypoxyanthine ? Inosine ? IMP ? adenylosuccinate is also disrupted in the PFC of PCP-treated animals.Network reordering via the GSVD provides a means to discover statistically validated differences in clustering between a pair of networks. In practice this analytical approach, when applied to metabolomic data, allows us to quantify the alterations in metabolic pathways between two experimental groups. With this new computational technique we identified metabolic pathway alterations that are consistent with known results. Furthermore, we discovered disruption in a novel series of purine reactions that may contribute to the PFC dysfunction and cognitive deficits seen in schizophrenia.

Project description:BackgroundTerpenoids are of high interest as chemical building blocks and pharmaceuticals. In microbes, terpenoids can be synthesized via the methylerythritol phosphate (MEP) or mevalonate (MVA) pathways. Although the MEP pathway has a higher theoretical yield, metabolic engineering has met with little success because the regulation of the pathway is poorly understood.ResultsWe applied metabolic control analysis to the MEP pathway in Escherichia coli expressing a heterologous isoprene synthase gene (ispS). The expression of ispS led to the accumulation of isopentenyl pyrophosphate (IPP)/dimethylallyl pyrophosphate (DMAPP) and severely impaired bacterial growth, but the coexpression of ispS and isopentenyl diphosphate isomerase (idi) restored normal growth and wild-type IPP/DMAPP levels. Targeted proteomics and metabolomics analysis provided a quantitative description of the pathway, which was perturbed by randomizing the ribosome binding site in the gene encoding 1-deoxyxylulose 5-phosphate synthase (Dxs). Dxs has a flux control coefficient of 0.35 (i.e., a 1% increase in Dxs activity resulted in a 0.35% increase in pathway flux) in the isoprene-producing strain and therefore exerted significant control over the flux though the MEP pathway. At higher dxs expression levels, the intracellular concentration of 2-C-methyl-D-erythritol-2,4-cyclopyrophosphate (MEcPP) increased substantially in contrast to the other MEP pathway intermediates, which were linearly dependent on the abundance of Dxs. This indicates that 4-hydroxy-3-methylbut-2-en-1-yl diphosphate synthase (IspG), which consumes MEcPP, became saturated and therefore limited the flux towards isoprene. The higher intracellular concentrations of MEcPP led to the efflux of this intermediate into the growth medium.DiscussionThese findings show the importance of Dxs, Idi and IspG and metabolite export for metabolic engineering of the MEP pathway and will facilitate further approaches for the microbial production of valuable isoprenoids.

Project description:The quinic acid ulitization (qut) pathway in Aspergillus nidulans is a dispensable carbon utilization pathway that catabolizes quinate to protocatechuate via dehydroquinate and dehydroshikimate(DHS). At the usual in vitro growth pH of 6.5, quinate enters the mycelium by means of a specific permease and is converted into PCA by the sequential action of the enzymes quinate dehydrogenase, 3-dehydroquinase and DHS dehydratase. The extent of control on metabolic flux exerted by the permease and the three pathway enzymes was investigated by applying the techniques of Metabolic Control Analysis. The flux control coefficients for each of the three quinate pathway enzymes were determined empirically, and the flux control coefficient of the quinate permease was inferred by use of the summation theorem. There measurements implied that, under the standard growth conditions used, the values for the flux control coefficients of the components of the quinate pathway were: quinate permease, 0.43; quinate dehydrogenase, 0.36; dehydroquinase, 0.18; DHS dehydratase, <0,03. Attempts to partially decouple quinate permease from the control over flux by measuring flux at pH 3.5 (when a significant percentage of the soluble quinate is protonated and able to enter the mycelium without the aid of a permease) led to an increase of approx. 50% in the flux control coefficient for dehydroquinase. Taken together with the fact that A. nidulans has a very efficient pH homeostasis mechanism, these experiments are consistent with the view that quinate permease exerts a high degree of control over pathway flux under the standard laboratory growth conditions at pH 6.5. The enzymes quinate dehydrogenase and 3-dehydroquinase have previously been overproduced in Escherichia coli, and protocols for their purification published. The remaining qut pathway enzyme DHS dehydratase was overproduced in E. coli and a purification protocol established. The purified DHS dehydratase was shown to have a K(m) of 530 microM for its substrate DHS and a requirement for bivalent metal cations that could be fulfilled by Mg(2+), Mn(2+) or Zn(2+). All three qut pathway enzymes were purified in bulk and their elasticity coefficients with respect to the three quinate pathway intermediates were derived over a range of concentrations in a core tricine/NaOH buffer, augmented with necessary cofactors and bivalent cations as appropriate. Using these empirically determined relative values, in conjunction with the connectivity theorem, the relative ratios of the flux control coefficients for the various quinate pathway enzymes, and how this control shifts between them, was determined over a range of possible metabolic concentrations. These calculations, although clearly subject to caveates about the relationswhip between kinetic measurements in vitro and the situation in vivo, were able to successfully predict the hiearchy of control observed under the standard laboratory growth conditions. The calculations imply that the hierarchy of control exerted by the quinate pathway enzymes is stable and relatively insensitive to changing metabolite concentrations in the ranges most likely to correspond to those found in vivo. The effects of substituting the type I 3-dehydroquinases from Salmonella typhi and the A. nidulans AROM protein (a pentadomain protein catalysing the conversion of 3-deoxy-D-arabinoheptulosonic acid 7-phosphate into 5-enolpyruvylshikimate 3 phosphate), and the Mycobacterium tuberculosis type II 3-dehydroquinase, in the quinate pathway were investigated and found to have an effect. In the case of S. typhi and A. nidulans, overproduction of heterologous dehydroquinase led to a diminuation of pathway flux caused by a lowering of in vivo quinate dehydrogenase levels increased above those of the wild type. We speculate that these changes in qu

Project description:Despite depression being one of the most prevalent and debilitating disorders worldwide, it has been difficult to understand its pathophysiology and to develop more effective treatments. Maladaptive transcriptional regulation within limbic neural circuits, including reward processing regions such as the nucleus accumbens (NAc), in response to chronic stress is thought to be a major contributor to the development of the syndrome. Epigenetic events?in particular, histone writers and erasers?that alter chromatin structure to regulate programs of gene expression have increasingly been associated with depression-related behavioral abnormalities in animal models and in depressed humans examined postmortem. However, very little is known about the ATP-dependent chromatin remodelers that control nucleosome positioning and the packing state of chromatin. Here we show that the ACF complex, part of the ISWI family of chromatin remodelers, is persistently and selectively upregulated in the NAc of mice that are susceptible to chronic social stress, as well as in the NAc of depressed human. We further establish that ACF induction is both necessary and sufficient for susceptibility to stress-induced depressive-like behaviors. Using ChIP-seq, we demonstrate that altered ACF binding after chronic stress is strongly correlated with altered nucleosome positioning, in particular, around the transcriptional start sites of affected genes. These alterations in ACF binding and nucleosome repositioning are associated with repressed expression of a subset of genes in animals that are susceptible to chronic stress. Together, these findings establish that active ATP-dependent chromatin remodeling by the ACF complex is a key regulator in the repression of genes that mediate susceptibility to social stress, and provide novel candidate targets for improved therapeutics of depression and other stress-related disorders. c57bl/6 mice underwent chronic social defeat stress (CSDS), and social interaction test was used to separate animals into control, susceptible and resilient groups. Nucleus accumbens (NAc) tissue was collected 48 hours after the last defeat session, and then Acf1, SNF2H ChIP-seq or H3 MNase-seq were performed based on the control, susceptible, and resilient groups. Three sequencing replicates were performed on each group.

Dataset Information

Analysis of complex metabolic behavior through pathway decomposition.

Background

Results

Conclusions

Publications

Analysis of complex metabolic behavior through pathway decomposition.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets