23
views
0
recommends
+1 Recommend
0 collections
    0
    shares
      • Record: found
      • Abstract: found
      • Article: found
      Is Open Access

      Flux Design: In silico design of cell factories based on correlation of pathway fluxes to desired properties

      research-article

      Read this article at

      Bookmark
          There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.

          Abstract

          Background

          The identification of genetic target genes is a key step for rational engineering of production strains towards bio-based chemicals, fuels or therapeutics. This is often a difficult task, because superior production performance typically requires a combination of multiple targets, whereby the complex metabolic networks complicate straightforward identification. Recent attempts towards target prediction mainly focus on the prediction of gene deletion targets and therefore can cover only a part of genetic modifications proven valuable in metabolic engineering. Efficient in silico methods for simultaneous genome-scale identification of targets to be amplified or deleted are still lacking.

          Results

          Here we propose the identification of targets via flux correlation to a chosen objective flux as approach towards improved biotechnological production strains with optimally designed fluxes. The approach, we name Flux Design, computes elementary modes and, by search through the modes, identifies targets to be amplified (positive correlation) or down-regulated (negative correlation). Supported by statistical evaluation, a target potential is attributed to the identified reactions in a quantitative manner. Based on systems-wide models of the industrial microorganisms Corynebacterium glutamicum and Aspergillus niger, up to more than 20,000 modes were obtained for each case, differing strongly in production performance and intracellular fluxes. For lysine production in C. glutamicum the identified targets nicely matched with reported successful metabolic engineering strategies. In addition, simulations revealed insights, e.g. into the flexibility of energy metabolism. For enzyme production in A.niger flux correlation analysis suggested a number of targets, including non-obvious ones. Hereby, the relevance of most targets depended on the metabolic state of the cell and also on the carbon source.

          Conclusions

          Objective flux correlation analysis provided a detailed insight into the metabolic networks of industrially relevant prokaryotic and eukaryotic microorganisms. It was shown that capacity, pathway usage, and relevant genetic targets for optimal production partly depend on the network structure and the metabolic state of the cell which should be considered in future metabolic engineering strategies. The presented strategy can be generally used to identify priority sorted amplification and deletion targets for metabolic engineering purposes under various conditions and thus displays a useful strategy to be incorporated into efficient strain and bioprocess optimization.

          Related collections

          Most cited references46

          • Record: found
          • Abstract: found
          • Article: not found

          Analysis of optimality in natural and perturbed metabolic networks.

          An important goal of whole-cell computational modeling is to integrate detailed biochemical information with biological intuition to produce testable predictions. Based on the premise that prokaryotes such as Escherichia coli have maximized their growth performance along evolution, flux balance analysis (FBA) predicts metabolic flux distributions at steady state by using linear programming. Corroborating earlier results, we show that recent intracellular flux data for wild-type E. coli JM101 display excellent agreement with FBA predictions. Although the assumption of optimality for a wild-type bacterium is justifiable, the same argument may not be valid for genetically engineered knockouts or other bacterial strains that were not exposed to long-term evolutionary pressure. We address this point by introducing the method of minimization of metabolic adjustment (MOMA), whereby we test the hypothesis that knockout metabolic fluxes undergo a minimal redistribution with respect to the flux configuration of the wild type. MOMA employs quadratic programming to identify a point in flux space, which is closest to the wild-type point, compatibly with the gene deletion constraint. Comparing MOMA and FBA predictions to experimental flux data for E. coli pyruvate kinase mutant PB25, we find that MOMA displays a significantly higher correlation than FBA. Our method is further supported by experimental data for E. coli knockout growth rates. It can therefore be used for predicting the behavior of perturbed metabolic networks, whose growth performance is in general suboptimal. MOMA and its possible future extensions may be useful in understanding the evolutionary optimization of metabolism.
            Bookmark
            • Record: found
            • Abstract: found
            • Article: not found

            Global analysis of mRNA decay and abundance in Escherichia coli at single-gene resolution using two-color fluorescent DNA microarrays.

            Much of the information available about factors that affect mRNA decay in Escherichia coli, and by inference in other bacteria, has been gleaned from study of less than 25 of the approximately 4,300 predicted E. coli messages. To investigate these factors more broadly, we examined the half-lives and steady-state abundance of known and predicted E. coli mRNAs at single-gene resolution by using two-color fluorescent DNA microarrays. An rRNA-based strategy for normalization of microarray data was developed to permit quantitation of mRNA decay after transcriptional arrest by rifampicin. We found that globally, mRNA half-lives were similar in nutrient-rich media and defined media in which the generation time was approximately tripled. A wide range of stabilities was observed for individual mRNAs of E. coli, although approximately 80% of all mRNAs had half-lives between 3 and 8 min. Genes having biologically related metabolic functions were commonly observed to have similar stabilities. Whereas the half-lives of a limited number of mRNAs correlated positively with their abundance, we found that overall, increased mRNA stability is not predictive of increased abundance. Neither the density of putative sites of cleavage by RNase E, which is believed to initiate mRNA decay in E. coli, nor the free energy of folding of 5' or 3' untranslated region sequences was predictive of mRNA half-life. Our results identify previously unsuspected features of mRNA decay at a global level and also indicate that generalizations about decay derived from the study of individual gene transcripts may have limited applicability.
              Bookmark
              • Record: found
              • Abstract: found
              • Article: not found

              On the frequency of protein glycosylation, as deduced from analysis of the SWISS-PROT database.

              The SWISS-PROT protein sequence data bank contains at present nearly 75,000 entries, almost two thirds of which include the potential N-glycosylation consensus sequence, or sequon, NXS/T (where X can be any amino acid but proline) and thus may be glycoproteins. The number of proteins filed as glycoproteins is however considerably smaller, 7942, of which 749 have been characterized with respect to the total number of their carbohydrate units and sites of attachment of the latter to the protein, as well as the nature of the carbohydrate-peptide linking group. Of these well characterized glycoproteins, about 90% carry either N-linked carbohydrate units alone or both N- and O-linked ones, attached at 1297 N-glycosylation sites (1.9 per glycoprotein molecule) and the rest are O-glycosylated only. Since the total number of sequons in the well characterized glycoproteins is 1968, their rate of occupancy is 2/3. Assuming that the same number of N-linked units and rate of sequon occupancy occur in all sequon containing proteins and that the proportion of solely O-glycosylated proteins (ca. 10%) will also be the same as among the well characterized ones, we conclude that the majority of sequon containing proteins will be found to be glycosylated and that more than half of all proteins are glycoproteins.
                Bookmark

                Author and article information

                Journal
                BMC Syst Biol
                BMC Systems Biology
                BioMed Central
                1752-0509
                2009
                25 December 2009
                : 3
                : 120
                Affiliations
                [1 ]Institute of Biochemical Engineering, Technische Universität Braunschweig, Gaußstr 17, 38106 Braunschweig, Germany
                Article
                1752-0509-3-120
                10.1186/1752-0509-3-120
                2808316
                20035624
                0e4b9b54-c44c-47d1-ae9e-d324c7fa411c
                Copyright ©2009 Melzer et al; licensee BioMed Central Ltd.

                This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

                History
                : 4 June 2009
                : 25 December 2009
                Categories
                Research article

                Quantitative & Systems biology
                Quantitative & Systems biology

                Comments

                Comment on this article