Basic and applied uses of genome‐scale metabolic network reconstructions of Escherichia coli

The genome‐scale model (GEM) of metabolism in the bacterium Escherichia coli K‐12 has been in development for over a decade and is now in wide use. GEM‐enabled studies of E. coli have been primarily focused on six applications: (1) metabolic engineering, (2) model‐driven discovery, (3) prediction of cellular phenotypes, (4) analysis of biological network properties, (5) studies of evolutionary processes, and (6) models of interspecies interactions. In this review, we provide an overview of these applications along with a critical assessment of their successes and limitations, and a perspective on likely future developments in the field. Taken together, the studies performed over the past decade have established a genome‐scale mechanistic understanding of genotype–phenotype relationships in E. coli metabolism that forms the basis for similar efforts for other microbial species. Future challenges include the expansion of GEMs by integrating additional cellular processes beyond metabolism, the identification of key constraints based on emerging data types, and the development of computational methods able to handle such large‐scale network models with sufficient accuracy.     

Applications of genome‐scale metabolic reconstructions

The availability and utility of genome‐scale metabolic reconstructions have exploded since the first genome‐scale reconstruction was published a decade ago. Reconstructions have now been built for a wide variety of organisms, and have been used toward five major ends: (1) contextualization of high‐throughput data, (2) guidance of metabolic engineering, (3) directing hypothesis‐driven discovery, (4) interrogation of multi‐species relationships, and (5) network property discovery. In this review, we examine the many uses and future directions of genome‐scale metabolic reconstructions, and we highlight trends and opportunities in the field that will make the greatest impact on many fields of biology.     

Improving the phenotype predictions of a yeast genome‐scale metabolic model by incorporating enzymatic constraints

Genome‐scale metabolic models (GEMs) are widely used to calculate metabolic phenotypes. They rely on defining a set of constraints, the most common of which is that the production of metabolites and/or growth are limited by the carbon source uptake rate. However, enzyme abundances and kinetics, which act as limitations on metabolic fluxes, are not taken into account. Here, we present GECKO, a method that enhances a GEM to account for enzymes as part of reactions, thereby ensuring that each metabolic flux does not exceed its maximum capacity, equal to the product of the enzyme's abundance and turnover number. We applied GECKO to a Saccharomyces cerevisiae GEM and demonstrated that the new model could correctly describe phenotypes that the previous model could not, particularly under high enzymatic pressure conditions, such as yeast growing on different carbon sources in excess, coping with stress, or overexpressing a specific pathway. GECKO also allows to directly integrate quantitative proteomics data; by doing so, we significantly reduced flux variability of the model, in over 60% of metabolic reactions. Additionally, the model gives insight into the distribution of enzyme usage between and within metabolic pathways. The developed method and model are expected to increase the use of model‐based design in metabolic engineering.     

Dissecting the energy metabolism in Mycoplasma pneumoniae through genome‐scale metabolic modeling

Mycoplasma pneumoniae, a threatening pathogen with a minimal genome, is a model organism for bacterial systems biology for which substantial experimental information is available. With the goal of understanding the complex interactions underlying its metabolism, we analyzed and characterized the metabolic network of M. pneumoniae in great detail, integrating data from different omics analyses under a range of conditions into a constraint‐based model backbone. Iterating model predictions, hypothesis generation, experimental testing, and model refinement, we accurately curated the network and quantitatively explored the energy metabolism. In contrast to other bacteria, M. pneumoniae uses most of its energy for maintenance tasks instead of growth. We show that in highly linear networks the prediction of flux distributions for different growth times allows analysis of time‐dependent changes, albeit using a static model. By performing an in silico knock‐out study as well as analyzing flux distributions in single and double mutant phenotypes, we demonstrated that the model accurately represents the metabolism of M. pneumoniae. The experimentally validated model provides a solid basis for understanding its metabolic regulatory mechanisms.     

A comprehensive genome‐scale reconstruction of Escherichia coli metabolism—2011

The initial genome‐scale reconstruction of the metabolic network of Escherichia coli K‐12 MG1655 was assembled in 2000. It has been updated and periodically released since then based on new and curated genomic and biochemical knowledge. An update has now been built, named iJO1366, which accounts for 1366 genes, 2251 metabolic reactions, and 1136 unique metabolites. iJO1366 was (1) updated in part using a new experimental screen of 1075 gene knockout strains, illuminating cases where alternative pathways and isozymes are yet to be discovered, (2) compared with its predecessor and to experimental data sets to confirm that it continues to make accurate phenotypic predictions of growth on different substrates and for gene knockout strains, and (3) mapped to the genomes of all available sequenced E. coli strains, including pathogens, leading to the identification of hundreds of unannotated genes in these organisms. Like its predecessors, the iJO1366 reconstruction is expected to be widely deployed for studying the systems biology of E. coli and for metabolic engineering applications.     

Effect of amino acid supplementation on titer and glycosylation distribution in hybridoma cell cultures—Systems biology‐based interpretation using genome‐scale metabolic flux balance model and multivariate data analysis

Genome‐scale flux balance analysis (FBA) is a powerful systems biology tool to characterize intracellular reaction fluxes during cell cultures. FBA estimates intracellular reaction rates by optimizing an objective function, subject to the constraints of a metabolic model and media uptake/excretion rates. A dynamic extension to FBA, dynamic flux balance analysis (DFBA), can calculate intracellular reaction fluxes as they change during cell cultures. In a previous study by Read et al. (2013), a series of informed amino acid supplementation experiments were performed on twelve parallel murine hybridoma cell cultures, and this data was leveraged for further analysis (Read et al., Biotechnol Prog. 2013;29:745–753). In order to understand the effects of media changes on the model murine hybridoma cell line, a systems biology approach is applied in the current study. Dynamic flux balance analysis was performed using a genome‐scale mouse metabolic model, and multivariate data analysis was used for interpretation. The calculated reaction fluxes were examined using partial least squares and partial least squares discriminant analysis. The results indicate media supplementation increases product yield because it raises nutrient levels extending the growth phase, and the increased cell density allows for greater culture performance. At the same time, the directed supplementation does not change the overall metabolism of the cells. This supports the conclusion that product quality, as measured by glycoform assays, remains unchanged because the metabolism remains in a similar state. Additionally, the DFBA shows that metabolic state varies more at the beginning of the culture but less by the middle of the growth phase, possibly due to stress on the cells during inoculation. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:1163–1173, 2016     

Identification of anticancer drugs for hepatocellular carcinoma through personalized genome‐scale metabolic modeling

Genome‐scale metabolic models (GEMs) have proven useful as scaffolds for the integration of omics data for understanding the genotype–phenotype relationship in a mechanistic manner. Here, we evaluated the presence/absence of proteins encoded by 15,841 genes in 27 hepatocellular carcinoma (HCC) patients using immunohistochemistry. We used this information to reconstruct personalized GEMs for six HCC patients based on the proteomics data, HMR 2.0, and a task‐driven model reconstruction algorithm (tINIT). The personalized GEMs were employed to identify anticancer drugs using the concept of antimetabolites; i.e., drugs that are structural analogs to metabolites. The toxicity of each antimetabolite was predicted by assessing the in silico functionality of 83 healthy cell type‐specific GEMs, which were also reconstructed with the tINIT algorithm. We predicted 101 antimetabolites that could be effective in preventing tumor growth in all HCC patients, and 46 antimetabolites which were specific to individual patients. Twenty‐two of the 101 predicted antimetabolites have already been used in different cancer treatment strategies, while the remaining antimetabolites represent new potential drugs. Finally, one of the identified targets was validated experimentally, and it was confirmed to attenuate growth of the HepG2 cell line.     

Metabolic reconstruction identifies strain‐specific regulation of virulence in Toxoplasma gondii

Increasingly, metabolic potential is proving to be a critical determinant governing a pathogen's virulence as well as its capacity to expand its host range. To understand the potential contribution of metabolism to strain‐specific infectivity differences, we present a constraint‐based metabolic model of the opportunistic parasite, Toxoplasma gondii. Dominated by three clonal strains (Type I, II, and III demonstrating distinct virulence profiles), T. gondii exhibits a remarkably broad host range. Integrating functional genomic data, our model (which we term as iCS382) reveals that observed strain‐specific differences in growth rates are driven by altered capacities for energy production. We further predict strain‐specific differences in drug susceptibilities and validate one of these predictions in a drug‐based assay, with a Type I strain demonstrating resistance to inhibitors that are effective against a Type II strain. We propose that these observed differences reflect an evolutionary strategy that allows the parasite to extend its host range, as well as result in a subsequent partitioning into discrete strains that display altered virulence profiles across different hosts, different organs, and even cell types.     

Large‐scale molecular genetic analysis in plant‐pathogenic fungi: a decade of genome‐wide functional analysis

Plant‐pathogenic fungi cause diseases to all major crop plants world‐wide and threaten global food security. Underpinning fungal diseases are virulence genes facilitating plant host colonization that often marks pathogenesis and crop failures, as well as an increase in staple food prices. Fungal molecular genetics is therefore the cornerstone to the sustainable prevention of disease outbreaks. Pathogenicity studies using mutant collections provide immense function‐based information regarding virulence genes of economically relevant fungi. These collections are rich in potential targets for existing and new biological control agents. They contribute to host resistance breeding against fungal pathogens and are instrumental in searching for novel resistance genes through the identification of fungal effectors. Therefore, functional analyses of mutant collections propel gene discovery and characterization, and may be incorporated into disease management strategies. In the light of these attributes, mutant collections enhance the development of practical solutions to confront modern agricultural constraints. Here, a critical review of mutant collections constructed by various laboratories during the past decade is provided. We used Magnaporthe oryzae and Fusarium graminearum studies to show how mutant screens contribute to bridge existing knowledge gaps in pathogenicity and fungal–host interactions.     

Photobioreactor design for isotopic non‐stationary 13C‐metabolic flux analysis (INST 13C‐MFA) under photoautotrophic conditions

Adaptive metabolic behavior of photoautotrophic microorganisms toward genetic and environmental perturbations can be interpreted in a quantitative depiction of carbon flow through a biochemical reaction network using isotopic non‐stationary ¹³C‐metabolic flux analysis (INST ¹³C‐MFA). To evaluate ¹³C‐metabolic flux maps for Chlamydomonas reinhardtii, an original experimental framework was designed allowing rapid, reliable collection of high‐quality isotopomer data against time. It involved (i) a short‐time ¹³C labeling injection device based on mixing control in a torus‐shaped photobioreactor with plug‐flow hydrodynamics allowing a sudden step‐change in the ¹³C proportion in the substrate feed and (ii) a rapid sampling procedure using an automatic fast filtration method coupled to a manual rapid liquid nitrogen quenching step. ¹³C‐substrate labeling enrichment was controlled through the total dissolved inorganic carbon concentration in the pulsed solution. First results were obtained from steady‐state continuous culture measurements allowing the characterization of the kinetics of label incorporation into light‐limited growing cells cultivated in a photobioreactor operating at the maximal biomass productivity for an incident photon flux density of 200 µmol m^?2 s^?1. ¹³C label incorporation was measured for 21 intracellular metabolites using IC‐MS/MS in 58 samples collected across a labeling experiment duration of 7 min. The fastest labeling rate was observed for 2/3‐phosphoglycerate with an apparent isotopic stationary state reached after 300 s. The labeling rate was consistent with the optimized mixing time of about 4.9 s inside the reactor and the shortest reliable sampling period assessed at 5 s. Biotechnol. Bioeng. 2012; 109: 3030–3040. © 2012 Wiley Periodicals, Inc.     

代谢网络定量分析研究进展

综述了代谢工程中代谢控制分析、代谢通量分析、生化系统理论、途径分析、控制论模型等定量分析方法的基本理论,以实例说明了这些方法的应用,并对代谢分析方法的发展进行了展望。     

Combining mechanistic and data‐driven approaches to gain process knowledge on the control of the metabolic shift to lactate uptake in a fed‐batch CHO process

A growing body of knowledge is available on the cellular regulation of overflow metabolism in mammalian hosts of recombinant protein production. However, to develop strategies to control the regulation of overflow metabolism in cell culture processes, the effect of process parameters on metabolism has to be well understood. In this study, we investigated the effect of pH and temperature shift timing on lactate metabolism in a fed‐batch Chinese hamster ovary (CHO) process by using a Design of Experiments (DoE) approach. The metabolic switch to lactate consumption was controlled in a broad range by the proper timing of pH and temperature shifts. To extract process knowledge from the large experimental dataset, we proposed a novel methodological concept and demonstrated its usefulness with the analysis of lactate metabolism. Time‐resolved metabolic flux analysis and PLS‐R VIP were combined to assess the correlation of lactate metabolism and the activity of the major intracellular pathways. Whereas the switch to lactate uptake was mainly triggered by the decrease in the glycolytic flux, lactate uptake was correlated to TCA activity in the last days of the cultivation. These metabolic interactions were visualized on simple mechanistic plots to facilitate the interpretation of the results. Taken together, the combination of knowledge‐based mechanistic modeling and data‐driven multivariate analysis delivered valuable insights into the metabolic control of lactate production and has proven to be a powerful tool for the analysis of large metabolic datasets. © 2015 American Institute of Chemical Engineers Biotechnol. Prog., 31:1657–1668, 2015     

Knowledge‐guided laboratory evolution of protein thermolability

In rare but nevertheless important cases it is of practical interest to decrease the thermostability of an enzyme, that is, to increase thermolability in a controlled manner. In the present model study, this unconventional goal has been reached by applying directed evolution to the lipase from Pseudomonas aeruginosa (PAL). By utilizing the B‐factor iterative test (B‐FIT), previously developed to increase the thermostability of enzymes, it was possible to reduce the value from 71.6°C in the case of wild type (WT‐PAL) to 35.6°C (best mutant) without affecting the catalytic profile in terms of substrate acceptance or enantioselectivity at room temperature. Accordingly, saturation mutagenesis was performed at sites in PAL, which on the basis of its X‐ray structure, have the lowest B‐factors indicative of high rigidity. Focused mutations were introduced which can be expected to decrease rigidity, the ensuing increased flexibility leading to higher thermolability without changing the actual catalytic profile. Biotechnol. Bioeng. 2009;102: 1712–1717. © 2008 Wiley Periodicals, Inc.