In silico strategy to rationally engineer metabolite production: A case study for threonine in Escherichia coli

Genetic engineering of metabolic pathways is a standard strategy to increase the production of metabolites of economic interest. However, such flux increases could very likely lead to undesirable changes in metabolite concentrations, producing deleterious perturbations on other cellular processes. These negative effects could be avoided by implementing a balanced increase of enzyme concentrations according to the Universal Method [Kacser and Acerenza (1993) Eur J Biochem 216:361–367]. Exact application of the method usually requires modification of many reactions, which is difficult to achieve in practice. Here, improvement of threonine production via pyruvate kinase deletion in Escherichia coli is used as a case study to demonstrate a partial application of the Universal Method, which includes performing sensitivity analysis. Our analysis predicts that manipulating a few reactions is sufficient to obtain an important increase in threonine production without major perturbations of metabolite concentrations. Biotechnol. Bioeng. 2009;103: 609–620. © 2009 Wiley Periodicals, Inc.     

Observing local and global properties of metabolic pathways: 'load points' and 'choke points' in the metabolic networks

MOTIVATION: The local and global aspects of metabolic network analyses allow us to identify enzymes or reactions that are crucial for the survival of the organism(s), therefore directing us towards the discovery of potential drug targets. RESULTS: We demonstrate a new method ('load points') to rank the enzymes/metabolites in the metabolic network and propose a model to determine and rank the biochemical lethality in metabolic networks (enzymes/metabolites) through 'choke points'. Based on an extended form of the graph theory model of metabolic networks, metabolite structural information was used to calculate the k-shortest paths between metabolites (the presence of more than one competing path between substrate and product). On the basis of these paths and connectivity information, load points were calculated and used to empirically rank the importance of metabolites/enzymes in the metabolic network. The load point analysis emphasizes the role that the biochemical structure of a metabolite, rather than its connectivity (hubs), plays in the conversion pathway. In order to identify potential drug targets (based on the biochemical lethality of metabolic networks), the concept of choke points and load points was used to find enzymes (edges) which uniquely consume or produce a particular metabolite (nodes). A non-pathogenic bacterial strain Bacillus subtilis 168 (lactic acid producing bacteria) and a related pathogenic bacterial strain Bacillus anthracis Sterne (avirulent but toxigenic strain, producing the toxin Anthrax) were selected as model organisms. The choke point strategy was implemented on the pathogen bacterial network of B.anthracis Sterne. Potential drug targets are proposed based on the analysis of the top 10 choke points in the bacterial network. A comparative study between the reported top 10 bacterial choke points and the human metabolic network was performed. Further biological inferences were made on results obtained by performing a homology search against the human genome. AVAILABILITY: The load and choke point modules are introduced in the Pathway Hunter Tool (PHT), the basic version of which is available on http://www.pht.uni-koeln.de.     

Metabolite damage and repair in metabolic engineering design

The necessarily sharp focus of metabolic engineering and metabolic synthetic biology on pathways and their fluxes has tended to divert attention from the damaging enzymatic and chemical side-reactions that pathway metabolites can undergo. Although historically overlooked and underappreciated, such metabolite damage reactions are now known to occur throughout metabolism and to generate (formerly enigmatic) peaks detected in metabolomics datasets. It is also now known that metabolite damage is often countered by dedicated repair enzymes that undo or prevent it. Metabolite damage and repair are highly relevant to engineered pathway design: metabolite damage reactions can reduce flux rates and product yields, and repair enzymes can provide robust, host-independent solutions. Herein, after introducing the core principles of metabolite damage and repair, we use case histories to document how damage and repair processes affect efficient operation of engineered pathways – particularly those that are heterologous, non-natural, or cell-free. We then review how metabolite damage reactions can be predicted, how repair reactions can be prospected, and how metabolite damage and repair can be built into genome-scale metabolic models. Lastly, we propose a versatile ‘plug and play’ set of well-characterized metabolite repair enzymes to solve metabolite damage problems known or likely to occur in metabolic engineering and synthetic biology projects.     

Parametric analysis of metabolic fluxes of α-amylase and protease-producing Bacillus subtilis

Parametric analysis was applied for a metabolic flux model for the fed-batch culture of Bacillus subtilis producing recombinant α-amylase and protease. The metabolic flux model was formulated as a linear programming problem consisting of 49 reactions (decision variables) and 50 metabolites (equality constraints). This study was aimed to determine the response of the metabolic fluxes and objective function value of minimizing the difference between ATP consumption and ATP production (ATP balance). With regard to intracellular metabolite accumulation, the objective function value was least sensitive to variation in succinate and most sensitive to variation in malate. Amongst the variations in the accumulation rates of extracellular metabolites, the objective function value was least sensitive to variation in glutamate and most sensitive to variation in starch hydrolysis and triglyceride synthesis. A 10% variation in metabolite accumulation rates caused a maximum of 13.8% variation (standard error = 3.8%) in the objective function value.     

Modules of co-regulated metabolites in turmeric (Curcuma longa) rhizome suggest the existence of biosynthetic modules in plant specialized metabolism

Turmeric is an excellent example of a plant that produces largenumbers of metabolites from diverse metabolic pathways or networks.It is hypothesized that these metabolic pathways or networkscontain biosynthetic modules, which lead to the formation ofmetabolite modules—groups of metabolites whose productionis co-regulated and biosynthetically linked. To test whethersuch co-regulated metabolite modules do exist in this plant,metabolic profiling analysis was performed on turmeric rhizomesamples that were collected from 16 different growth and developmenttreatments, which had significant impacts on the levels of 249volatile and non-volatile metabolites that were detected. Importantly,one of the many co-regulated metabolite modules that were indeedreadily detected in this analysis contained the three majorcurcuminoids, whereas many other structurally related diarylheptanoidsbelonged to separate metabolite modules, as did groups of terpenoids.The existence of these co-regulated metabolite modules supportedthe hypothesis that the 3-methoxyl groups on the aromatic ringsof the curcuminoids are formed before the formation of the heptanoidbackbone during the biosynthesis of curcumin and also suggestedthe involvement of multiple polyketide synthases with differentsubstrate selectivities in the formation of the array of diarylheptanoidsdetected in turmeric. Similar conclusions about terpenoid biosynthesiscould also be made. Thus, discovery and analysis of metabolitemodules can be a powerful predictive tool in efforts to understandmetabolism in plants. Key words: Biosynthesis, Curcuma longa, curcumin, metabolite module, metabolomics, rhizome, specialized metabolism Received 16 June 2008; Revised 29 September 2008 Accepted 2 October 2008     

Metabolite coupling in genome-scale metabolic networks

Background  

Biochemically detailed stoichiometric matrices have now been reconstructed for various bacteria, yeast, and for the human cardiac mitochondrion based on genomic and proteomic data. These networks have been manually curated based on legacy data and elementally and charge balanced. Comparative analysis of these well curated networks is now possible. Pairs of metabolites often appear together in several network reactions, linking them topologically. This co-occurrence of pairs of metabolites in metabolic reactions is termed herein "metabolite coupling." These metabolite pairs can be directly computed from the stoichiometric matrix, S. Metabolite coupling is derived from the matrix ŜŜ ^T, whose off-diagonal elements indicate the number of reactions in which any two metabolites participate together, where Ŝ is the binary form of S.     

Conservation of high-flux backbone in alternate optimal and near-optimal flux distributions of metabolic networks

Constraint-based flux balance analysis (FBA) has proven successful in predicting the flux distribution of metabolic networks in diverse environmental conditions. FBA finds one of the alternate optimal solutions that maximizes the biomass production rate. Almaas et al. have shown that the flux distribution follows a power law, and it is possible to associate with most metabolites two reactions which maximally produce and consume a given metabolite, respectively. This observation led to the concept of high-flux backbone (HFB) in metabolic networks. In previous work, the HFB has been computed using a particular optima obtained using FBA. In this paper, we investigate the conservation of HFB of a particular solution for a given medium across different alternate optima and near-optima in metabolic networks of E. coli and S. cerevisiae. Using flux variability analysis (FVA), we propose a method to determine reactions that are guaranteed to be in HFB regardless of alternate solutions. We find that the HFB of a particular optima is largely conserved across alternate optima in E. coli, while it is only moderately conserved in S. cerevisiae. However, the HFB of a particular near-optima shows a large variation across alternate near-optima in both organisms. We show that the conserved set of reactions in HFB across alternate near-optima has a large overlap with essential reactions and reactions which are both uniquely consuming (UC) and uniquely producing (UP). Our findings suggest that the structure of the metabolic network admits a high degree of redundancy and plasticity in near-optimal flow patterns enhancing system robustness for a given environmental condition.     

Evolutionary optimization of metabolic pathways. Theoretical reconstruction of the stoichiometry of ATP and NADH producing systems

The structural design of ATP and NADH producing systems, such as glycolysis and the citric acid cycle (TCA), is analysed using optimization principles. It is assumed that these pathways combined with oxidative phosphorylation have reached, during their evolution, a high efficiency with respect to ATP production rates. On the basis of kinetic and thermodynamic principles, conclusions are derived concerning the optimal stoichiometry of such pathways. Extending previous investigations, both the concentrations of adenine nucleotides as well as nicotinamide adenine dinucleotides are considered variable quantities. This implies the consideration of the interaction of an ATP and NADH producing system, an ATP consuming system, a system coupling NADH consumption with ATP production and a system consuming NADH decoupled from ATP production. It is examined in what respect real metabolic pathways can be considered optimal by studying a large number of alternative pathways. The kinetics of the individual reactions are described by linear or bilinear functions of reactant concentrations. In this manner, the steady-state ATP production rate can be calculated for any possible ATP and NADH producing pathway. It is shown that most of the possible pathways result in a very low ATP production rate and that the very efficient pathways share common structural properties. Optimization with respect to the ATP production rate is performed by an evolutionary algorithm. The following results of our analysis are in close correspondence to the real design of glycolysis and the TCA cycle. (1) In all efficient pathways the ATP consuming reactions are located near the beginning. (2) In all efficient pathways NADH producing reactions as well as ATP producing reactions are located near the end. (3) The number of NADH molecules produced by the consumption of one energy-rich molecule (glucose) amounts to four in all efficient pathways. A distance measure and a measure for the internal ordering of reactions are introduced to study differences and similarities in the stoichiometries of metabolic pathways.     

The origin of correlations in metabolomics data

A phenomenon observed earlier in the development of metabolomics as a systems biology methodology, consists of a small but significant number of metabolites whose levels are highly correlated between biological replicates. Contrary to initial interpretations, these correlations are not necessarily only between neighboring metabolites in the metabolic network. Most metabolites that participate in common reactions are not correlated in this way, while some non-neighboring metabolites are highly correlated. Here we investigate the origin of such correlations using metabolic control analysis and computer simulation of biochemical networks. A series of cases is identified which lead to high correlation between metabolite pairs in replicate measurement. These are (1) chemical equilibrium, (2) mass conservation, (3) asymmetric control distribution, and (4) unusually high variance in the expression of a single gene. The importance of identifying metabolite correlations within a physiological state and changes of correlation between different states is discussed in the context of systems biology.     

Interpreting Metabolomic Profiles using Unbiased Pathway Models

Human disease is heterogeneous, with similar disease phenotypes resulting from distinct combinations of genetic and environmental factors. Small-molecule profiling can address disease heterogeneity by evaluating the underlying biologic state of individuals through non-invasive interrogation of plasma metabolite levels. We analyzed metabolite profiles from an oral glucose tolerance test (OGTT) in 50 individuals, 25 with normal (NGT) and 25 with impaired glucose tolerance (IGT). Our focus was to elucidate underlying biologic processes. Although we initially found little overlap between changed metabolites and preconceived definitions of metabolic pathways, the use of unbiased network approaches identified significant concerted changes. Specifically, we derived a metabolic network with edges drawn between reactant and product nodes in individual reactions and between all substrates of individual enzymes and transporters. We searched for “active modules”—regions of the metabolic network enriched for changes in metabolite levels. Active modules identified relationships among changed metabolites and highlighted the importance of specific solute carriers in metabolite profiles. Furthermore, hierarchical clustering and principal component analysis demonstrated that changed metabolites in OGTT naturally grouped according to the activities of the System A and L amino acid transporters, the osmolyte carrier SLC6A12, and the mitochondrial aspartate-glutamate transporter SLC25A13. Comparison between NGT and IGT groups supported blunted glucose- and/or insulin-stimulated activities in the IGT group. Using unbiased pathway models, we offer evidence supporting the important role of solute carriers in the physiologic response to glucose challenge and conclude that carrier activities are reflected in individual metabolite profiles of perturbation experiments. Given the involvement of transporters in human disease, metabolite profiling may contribute to improved disease classification via the interrogation of specific transporter activities.     

Functional essentiality from topology features in metabolic networks: a case study in yeast

The relation between the position of mutations in Saccharomyces cerevisiae metabolic network and their lethality is the subject of this work. We represent the topology of the network by a directed graph: nodes are metabolites and arcs represent the reactions; a mutation corresponds to the removal of all the arcs referring to the deleted enzyme. Using publicly available knock-out data, we show that lethality corresponds to the lack of alternative paths in the perturbed network linking the nodes affected by the enzyme deletion. Such feature is at the basis of the recently recognized importance of 'marginal' arcs of metabolic networks.     

A direct cell quenching method for cell-culture based metabolomics

A crucial step in metabolomic analysis of cellular extracts is the cell quenching process. The conventional method first uses trypsin to detach cells from their growth surface. This inevitably changes the profile of cellular metabolites since the detachment of cells from the extracellular matrix alters their physiology. This conventional method also includes time consuming wash/centrifuge steps after trypsinization, but prior to quenching cell activity. During this time, a considerable portion of intracellular metabolites are lost, rendering the conventional method less than ideal for application to metabolomics. We report here a novel sample preparation method for metabolomics applications using adherent mammalian cells, which eliminates the time consumption and physiological stress of the trypsinization and wash/centrifuge steps. This new method was evaluated in the study of metabolic changes caused by 17α-ethynylestradiol (EE2) in estrogen receptor (ER)-positive MCF-7 and ER-negative MDA-MB-231 human breast cancer cell lines using NMR spectroscopy. The results demonstrated that our direct cell quenching method is rapid, effective, and exhibits greater metabolite retention, providing an increase of approximately a factor of 50 compared to the conventional method.     

Complex Marine Natural Products as Potential Epigenetic and Production Regulators of Antibiotics from a Marine Pseudomonas aeruginosa

Marine microbes are capable of producing secondary metabolites for defense and competition. Factors exerting an impact on secondary metabolite production of microbial communities included bioactive natural products and co-culturing. These external influences may have practical applications such as increased yields or the generation of new metabolites from otherwise silent genes in addition to reducing or limiting the production of undesirable metabolites. In this paper, we discuss the metabolic profiles of a marine Pseudomonas aeruginosa in the presence of a number of potential chemical epigenetic regulators, adjusting carbon sources and co-culturing with other microbes to induce a competitive response. As a result of these stressors certain groups of antibiotics or antimalarial agents were increased most notably when treating P. aeruginosa with sceptrin and co-culturing with another Pseudomonas sp. An interesting cross-talking event between these two Pseudomonas species when cultured together and exposed to sceptrin was observed.     

Analysis of metabolic network based on conservation of molecular structure

Recent work has revealed much about chemical reactions inside hundreds of organisms as well as universal characteristics of metabolic networks, which shed light on the evolution of the networks. However, characteristics of individual metabolites have been neglected. For example, some carbohydrates have structures that are decomposed into small molecules by metabolic reactions, but coenzymes such as ATP are mostly preserved. Such differences in metabolite characteristics are important for understanding the universal characteristics of metabolic networks. To quantify the structure conservation of metabolites, we defined the "structure conservation index" (SCI) for each metabolite as the fraction of metabolite atoms restored to their original positions through metabolic reactions. As expected, coenzymes and coenzyme-like metabolites that have reaction loops in the network show a higher SCI. Using the index, we found that the sum of metabolic fluxes is negatively correlated with the structure preservation of metabolite. Also, we found that each reaction path around high SCI metabolites changes independently, while changes in reaction paths involving low SCI metabolites coincide through evolution processes. These correlations may provide a clue to universal properties of metabolic networks.     

Optimizing metabolic pathways by screening for feasible synthetic reactions

Background

Reconstruction of genome-scale metabolic networks has resulted in models capable of reproducing experimentally observed biomass yield/growth rates and predicting the effect of alterations in metabolism for biotechnological applications. The existing studies rely on modifying the metabolic network of an investigated organism by removing or inserting reactions taken either from evolutionary similar organisms or from databases of biochemical reactions (e.g., KEGG). A potential disadvantage of these knowledge-driven approaches is that the result is biased towards known reactions, as such approaches do not account for the possibility of including novel enzymes, together with the reactions they catalyze.

Results

Here, we explore the alternative of increasing biomass yield in three model organisms, namely Bacillus subtilis, Escherichia coli, and Hordeum vulgare, by applying small, chemically feasible network modifications. We use the predicted and experimentally confirmed growth rates of the wild-type networks as reference values and determine the effect of inserting mass-balanced, thermodynamically feasible reactions on predictions of growth rate by using flux balance analysis.

Conclusions

While many replacements of existing reactions naturally lead to a decrease or complete loss of biomass production ability, in all three investigated organisms we find feasible modifications which facilitate a significant increase in this biological function. We focus on modifications with feasible chemical properties and a significant increase in biomass yield. The results demonstrate that small modifications are sufficient to substantially alter biomass yield in the three organisms. The method can be used to predict the effect of targeted modifications on the yield of any set of metabolites (e.g., ethanol), thus providing a computational framework for synthetic metabolic engineering.     

Prediction of Rate‐Limiting Reactions for Growth‐Associated Production Using a Constraint‐Based Approach

Identification of a rate‐limiting step in pathways is a key challenge in metabolic engineering. Although the prediction of rate‐limiting steps using a kinetic model is a powerful approach, there are several technical hurdles for developing a kinetic model. In this study, an in silico screening algorithm of key enzyme for metabolic engineering is developed to identify the possible rate‐limiting reactions for the growth‐coupled target production using a stoichiometric model without any experimental data and kinetic parameters. In this method, for each reaction, an upper‐bound flux constraint is imposed and the target production is predicted by linear programming. When the constraint decreases the target production at the optimal growth state, the reaction is thought to be a possible rate‐limiting step. For validation, this method is applied to the production of succinate or 1,4‐butanediol (1,4‐BDO) in Escherichia coli, in which the experimental engineering for eliminating rate‐limiting steps has been previously reported. In succinate production from glycerol, nine reactions including phosphoenolpyruvate carboxylase are predicted as the rate‐limiting steps. In 1,4‐BDO production from glucose, eight reactions including pyruvate dehydrogenase are predicted as the rate‐limiting steps. These predictions include experimentally identified rate‐limiting steps, which would contribute to metabolic engineering as a practical tool for screening candidates of rate‐limiting reactions.     

Systematic phenome analysis of Escherichia coli multiple‐knockout mutants reveals hidden reactions in central carbon metabolism

Central carbon metabolism is a basic and exhaustively analyzed pathway. However, the intrinsic robustness of the pathway might still conceal uncharacterized reactions. To test this hypothesis, we constructed systematic multiple‐knockout mutants involved in central carbon catabolism in Escherichia coli and tested their growth under 12 different nutrient conditions. Differences between in silico predictions and experimental growth indicated that unreported reactions existed within this extensively analyzed metabolic network. These putative reactions were then confirmed by metabolome analysis and in vitro enzymatic assays. Novel reactions regarding the breakdown of sedoheptulose‐7‐phosphate to erythrose‐4‐phosphate and dihydroxyacetone phosphate were observed in transaldolase‐deficient mutants, without any noticeable changes in gene expression. These reactions, triggered by an accumulation of sedoheptulose‐7‐phosphate, were catalyzed by the universally conserved glycolytic enzymes ATP‐dependent phosphofructokinase and aldolase. The emergence of an alternative pathway not requiring any changes in gene expression, but rather relying on the accumulation of an intermediate metabolite may be a novel mechanism mediating the robustness of these metabolic networks.     

Reconstruction of the central carbon metabolism of Aspergillus niger.

The topology of central carbon metabolism of Aspergillus niger was identified and the metabolic network reconstructed, by integrating genomic, biochemical and physiological information available for this microorganism and other related fungi. The reconstructed network may serve as a valuable database for annotation of genes identified in future genome sequencing projects on aspergilli. Based on the metabolic reconstruction, a stoichiometric model was set up that includes 284 metabolites and 335 reactions, of which 268 represent biochemical conversions and 67 represent transport processes between the different intracellular compartments and between the cell and the extracellular medium. The stoichiometry of the metabolic reactions was used in combination with biosynthetic requirements for growth and pseudo-steady state mass balances over intracellular metabolites for the quantification of metabolic fluxes using metabolite balancing. This framework was employed to perform an in silico characterisation of the phenotypic behaviour of A. niger grown on different carbon sources. The effects on growth of single reaction deletions were assessed and essential biochemical reactions were identified for different carbon sources. Furthermore, application of the stoichiometric model for assessing the metabolic capabilities of A. niger to produce metabolites was evaluated by using succinate production as a case study.     

K-FIT: An accelerated kinetic parameterization algorithm using steady-state fluxomic data

Kinetic models predict the metabolic flows by directly linking metabolite concentrations and enzyme levels to reaction fluxes. Robust parameterization of organism-level kinetic models that faithfully reproduce the effect of different genetic or environmental perturbations remains an open challenge due to the intractability of existing algorithms. This paper introduces Kinetics-based Fluxomics Integration Tool (K-FIT), a robust kinetic parameterization workflow that leverages a novel decomposition approach to identify steady-state fluxes in response to genetic perturbations followed by a gradient-based update of kinetic parameters until predictions simultaneously agree with the fluxomic data in all perturbed metabolic networks. The applicability of K-FIT to large-scale models is demonstrated by parameterizing an expanded kinetic model for E. coli (307 reactions and 258 metabolites) using fluxomic data from six mutants. The achieved thousand-fold speed-up afforded by K-FIT over meta-heuristic approaches is transformational enabling follow-up robustness of inference analyses and optimal design of experiments to inform metabolic engineering strategies.