Metabolomics and modelling approaches for systems metabolic engineering

Designing metabolic engineering strategies with genome-scale metabolic flux modeling Jiun Y Yen,1,2 Imen Tanniche,1 Amanda K Fisher,1–3 Glenda E Gillaspy,2 David R Bevan,2,3 Ryan S Senger,1 1Department of Biological Systems Engineering, 2Department of Biochemistry, 3Genomics, Bioinformatics, and Computational Biology Interdisciplinary Program, Virginia Tech, Blacksburg, VA, USA Abstract: New in silico tools that make use of genome-scale metabolic flux modeling are improving the design of metabolic engineering strategies. This review highlights the latest developments in this area, explains the interface between these in silico tools and the experimental implementation tools of metabolic engineers, and provides a way forward so that in silico predictions can better mimic reality and more experimental methods can be considered in simulation studies. The several methodologies for solving genome-scale models (eg, flux balance analysis [FBA], parsimonious FBA, flux variability analysis, and minimization of metabolic adjustment) all have unique advantages and applications. There are two basic approaches to designing metabolic engineering strategies in silico, and both have demonstrated success in the literature. The first involves: 1) making a genetic manipulation in a model; 2) testing for improved performance through simulation; and 3) iterating the process. The second approach has been used in more recently designed in silico tools and involves: 1) comparing metabolic flux profiles of a wild-type and ideally engineered state and 2) designing engineering strategies based on the differences in these flux profiles. Improvements in genome-scale modeling are anticipated in areas such as the inclusion of all relevant cellular machinery, the ability to understand and anticipate the results of combinatorial enrichment experiments, and constructing dynamic and flexible biomass equations that can respond to environmental and genetic manipulations. Keywords: genome-scale modeling, genome-scale modeling, flux balance analysis, flux variability analysis, minimization of metabolic adjustment, metabolic bottleneck, pathway optimization

Spatial Metabolomics and Lipidomics in Kidney Disease.

A dynamic model-based preparation of uniformly-13C-labeled internal standards facilitates quantitative metabolomics analysis of Penicillium chrysogenum

Constraint-based models of metabolism have been used in a variety of studies on drug discovery, metabolic engineering, evolution, and multi-species interactions. These genome-scale models can be generated for any sequenced organism since their main parameters (i.e., reaction stoichiometry) are highly conserved. Their relatively low parameter requirement makes these models easy to develop; however, these models often result in a solution space with multiple possible flux distributions, making it difficult to determine the precise flux state in the cell. Recent research efforts in this modeling field have investigated how additional experimental data, including gene expression, protein expression, metabolite concentrations, and kinetic parameters, can be used to reduce the solution space. This mini-review provides a summary of the data-driven computational approaches that are available for reducing the solution space and thereby improve predictions of intracellular fluxes by constraint-based models.

Metabolic engineering can be defined as purposeful modification of metabolic pathways and other cellular network to achieve desired cellular phenotype and performance. Rational metabolic engineering developed a couple of decades ago changed the way strains have been developed, which had traditionally been performed by random mutagenesis and selection. Now, we are observing another paradigm shift towards systems-level metabolic engineering, powered by the methods and tools developed in the discipline of systems biology. Systems metabolic engineering allows whole-cell-wide metabolic engineering based on the findings of systems biological studies including omics and computational analyses. Not only the metabolic network, but also gene regulatory and signaling networks can be engineered to develop an optimal strain. Also, it is important to consider fermentation and downstream processes during the upstream strain development. In this chapter, the general strategies of systems metabolic engineering are reviewed with relevant examples recently reported.

The advent of model-enabled workflows in systems biology allows for the integration of experimental data types with genome-scale models to discover new features of biology. This work demonstrates such a workflow, aimed at establishing a metabolomics platform applied to study the differences in metabolomes between anaerobic and aerobic growth of Escherichia coli. Constraint-based modeling was utilized to deduce a target list of compounds for downstream method development. An analytical and experimental methodology was developed and tailored to the compound chemistry and growth conditions of interest. This included the construction of a rapid sampling apparatus for use with anaerobic cultures. The resulting genome-scale data sets for anaerobic and aerobic growth were validated by comparison to previous small-scale studies comparing growth of E. coli under the same conditions. The metabolomics data were then integrated with the E. coli genome-scale metabolic model (GEM) via a sensitivity analysis that utilized reaction thermodynamics to reconcile simulated growth rates and reaction directionalities. This analysis highlighted several optimal network usage inconsistencies, including the incorrect use of the beta-oxidation pathway for synthesis of fatty acids. This analysis also identified enzyme promiscuity for the pykA gene, that is critical for anaerobic growth, and which has not been previously incorporated into metabolic models of E coli.

The effects of alternate optimal solutions in constraint-based genome-scale metabolic models

Constraint-based modeling is a well established modelling methodology used to analyze and study biological networks on both a medium and genome scale. Due to their large size, genome scale models are typically analysed using constraint-based optimization techniques. One widely used method is Flux Balance Analysis (FBA) which, for example, requires a modelling description to include: the definition of a stoichiometric matrix, an objective function and bounds on the values that fluxes can obtain at steady state. The Flux Balance Constraints (FBC) Package extends SBML Level 3 and provides a standardized format for the encoding, exchange and annotation of constraint-based models. It includes support for modelling concepts such as objective functions, flux bounds and model component annotation that facilitates reaction balancing. The FBC package establishes a base level for the unambiguous exchange of genome-scale, constraint-based models, that can be built upon by the community to meet future needs (e. g. by extending it to cover dynamic FBC models).

Summary Constraint-based modeling is a well established modelling methodology used to analyze and study biological networks on both a medium and genome scale. Due to their large size, genome scale models are typically analysed using constraint-based optimization techniques. One widely used method is Flux Balance Analysis (FBA) which, for example, requires a modelling description to include: the definition of a stoichiometric matrix, an objective function and bounds on the values that fluxes can obtain at steady state.The Flux Balance Constraints (FBC) Package extends SBML Level 3 and provides a standardized format for the encoding, exchange and annotation of constraint-based models. It includes support for modelling concepts such as objective functions, flux bounds and model component annotation that facilitates reaction balancing. The FBC package establishes a base level for the unambiguous exchange of genome-scale, constraint-based models, that can be built upon by the community to meet future needs (e. g. by extending it to cover dynamic FBC models).

Increasing demand for petroleum has stimulated industry to develop sustainable production of chemicals and biofuels using microbial cell factories. Fatty acids of chain lengths from C6 to C16 are propitious intermediates for the catalytic synthesis of industrial chemicals and diesel-like biofuels. The abundance of genetic information available for Escherichia coli and specifically, fatty acid metabolism in E. coli, supports this bacterium as a promising host for engineering a biocatalyst for the microbial production of fatty acids. Recent successes rooted in different features of systems metabolic engineering in the strain design of high-yielding medium chain fatty acid producing E. coli strains provide an emerging case study of design methods for effective strain design. Classical metabolic engineering and synthetic biology approaches enabled different and distinct design paths towards a high-yielding strain. Here we highlight a rational strain design process in systems biology, an integrated computational and experimental approach for carboxylic acid production, as an alternative method. Additional challenges inherent in achieving an optimal strain for commercialization of medium chain-length fatty acids will likely require a collection of strategies from systems metabolic engineering. Not only will the continued advancement in systems metabolic engineering result in these highly productive strains more quickly, this knowledge will extend more rapidly the carboxylic acid platform to the microbial production of carboxylic acids with alternate chain-lengths and functionalities.

Constraint-based metabolic models are currently the most comprehensive system-wide models of cellular metabolism. Several challenges arise when building an in silico constraint-based model of an organism that need to be addressed before flux balance analysis (FBA) can be applied for simulations. An algorithm called FBA-Gap is presented here that aids the construction of a working model based on plausible modifications to a given list of reactions that are known to occur in the organism. When applied to a working model, the algorithm gives a hypothesis concerning a minimal medium for sustaining the cell in culture. The utility of the algorithm is demonstrated in creating a new model organism and is applied to four existing working models for generating hypotheses about culture media. In modifying a partial metabolic reconstruction so that biomass may be produced using FBA, the proposed method is more efficient than a previously proposed method in that fewer new reactions are added to complete the model. The proposed method is also more accurate than other approaches in that only biologically plausible reactions and exchange reactions are used.

Patients with the renal phosphate–wasting disease X-linked hypophosphatemia (XLH) and Hyp mice, the murine homolog of XLH, are characterized by loss-of-function mutations in phosphate-regulating endopeptidase homolog X-linked (PHEX), leading to excessive secretion of the bone-derived phosphotropic hormone FGF23. The mineralization defect in patients with XLH and Hyp mice is caused by a combination of hypophosphatemia and local accumulation of mineralization-inhibiting molecules in bone. However, the mechanism by which PHEX deficiency regulates bone cell metabolism remains elusive. Here, we used spatial metabolomics by employing matrix-assisted laser desorption/ionization (MALDI) Fourier-transform ion cyclotron resonance mass spectrometry imaging (MSI) of undecalcified bone cryosections to characterize in situ metabolic changes in bones of Hyp mice in a holistic, unbiased manner. We found complex changes in Hyp bone metabolism, including perturbations in pentose phosphate, purine, pyrimidine, and phospholipid metabolism. Importantly, our study identified an upregulation of several biochemical pathways involved in intra- and extracellular production of the mineralization inhibitor pyrophosphate in the bone matrix of Hyp mice. Our data emphasize the utility of MSI–based spatial metabolomics in bone research and provide holistic in situ insights as to how Phex deficiency–induced changes in biochemical pathways in bone cells are linked to impaired bone mineralization.

BackgroundSolventogenic clostridia offer a sustainable alternative to petroleum-based production of butanol--an important chemical feedstock and potential fuel additive or replacement. C. beijerinckii is an attractive microorganism for strain design to improve butanol production because it (i) naturally produces the highest recorded butanol concentrations as a byproduct of fermentation; and (ii) can co-ferment pentose and hexose sugars (the primary products from lignocellulosic hydrolysis). Interrogating C. beijerinckii metabolism from a systems viewpoint using constraint-based modeling allows for simulation of the global effect of genetic modifications.ResultsWe present the first genome-scale metabolic model (iCM925) for C. beijerinckii, containing 925 genes, 938 reactions, and 881 metabolites. To build the model we employed a semi-automated procedure that integrated genome annotation information from KEGG, BioCyc, and The SEED, and utilized computational algorithms with manual curation to improve model completeness. Interestingly, we found only a 34% overlap in reactions collected from the three databases--highlighting the importance of evaluating the predictive accuracy of the resulting genome-scale model. To validate iCM925, we conducted fermentation experiments using the NCIMB 8052 strain, and evaluated the ability of the model to simulate measured substrate uptake and product production rates. Experimentally observed fermentation profiles were found to lie within the solution space of the model; however, under an optimal growth objective, additional constraints were needed to reproduce the observed profiles--suggesting the existence of selective pressures other than optimal growth. Notably, a significantly enriched fraction of actively utilized reactions in simulations--constrained to reflect experimental rates--originated from the set of reactions that overlapped between all three databases (P = 3.52 × 10-9, Fisher's exact test). Inhibition of the hydrogenase reaction was found to have a strong effect on butanol formation--as experimentally observed.ConclusionsMicrobial production of butanol by C. beijerinckii offers a promising, sustainable, method for generation of this important chemical and potential biofuel. iCM925 is a predictive model that can accurately reproduce physiological behavior and provide insight into the underlying mechanisms of microbial butanol production. As such, the model will be instrumental in efforts to better understand, and metabolically engineer, this microorganism for improved butanol production.

The branched chain amino acids (BCAAs), L-valine, L-leucine, and L-isoleucine, have recently been attracting much attention as their potential to be applied in various fields, including animal feed additive, cosmetics, and pharmaceuticals, increased. Strategies for developing microbial strains efficiently producing BCAAs are now in transition toward systems metabolic engineering from random mutagenesis. The metabolism and regulatory circuits of BCAA biosynthesis need to be thoroughly understood for designing system-wide metabolic engineering strategies. Here we review the current knowledge on BCAAs including their biosynthetic pathways, regulations, and export and transport systems. Recent advances in the development of BCAA production strains are also reviewed with a particular focus on L-valine production strain. At the end, the general strategies for developing BCAA overproducers by systems metabolic engineering are suggested.

Chapter 10 - Systems and synthetic metabolic engineering: Challenges and prospects