Aarash Bordbar

Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0

Over the past decade, a growing community of researchers has emerged around the use of constraint-based reconstruction and analysis (COBRA) methods to simulate, analyze and predict a variety of metabolic phenotypes using genome-scale models. The COBRA Toolbox, a MATLAB package for implementing COBRA methods, was presented earlier. Here we present a substantial update of this in silico toolbox. Version 2.0 of the COBRA Toolbox expands the scope of computations by including in silico analysis methods developed since its original release. New functions include (i) network gap filling, (ii) (13)C analysis, (iii) metabolic engineering, (iv) omics-guided analysis and (v) visualization. As with the first version, the COBRA Toolbox reads and writes systems biology markup language-formatted models. In version 2.0, we improved performance, usability and the level of documentation. A suite of test scripts can now be used to learn the core functionality of the toolbox and validate results. This toolbox lowers the barrier of entry to use powerful COBRA methods.

Constraint-based models predict metabolic and associated cellular functions

The prediction of cellular function from a genotype is a fundamental goal in biology. For metabolism, constraint-based modelling methods systematize biochemical, genetic and genomic knowledge into a mathematical framework that enables a mechanistic description of metabolic physiology. The use of constraint-based approaches has evolved over ~30 years, and an increasing number of studies have recently combined models with high-throughput data sets for prospective experimentation. These studies have led to validation of increasingly important and relevant biological predictions. As reviewed here, these recent successes have tangible implications in the fields of microbial evolution, interaction networks, genetic engineering and drug discovery.

Genomic landscapes of Chinese hamster ovary cell lines as revealed by the Cricetulus griseus draft genome

Chinese hamster ovary (CHO) cells, first isolated in 1957, are the preferred production host for many therapeutic proteins. Although genetic heterogeneity among CHO cell lines has been well documented, a systematic, nucleotide-resolution characterization of their genotypic differences has been stymied by the lack of a unifying genomic resource for CHO cells. Here we report a 2.4-Gb draft genome sequence of a female Chinese hamster, Cricetulus griseus, harboring 24,044 genes. We also resequenced and analyzed the genomes of six CHO cell lines from the CHO-K1, DG44 and CHO-S lineages. This analysis identified hamster genes missing in different CHO cell lines, and detected >3.7 million single-nucleotide polymorphisms (SNPs), 551,240 indels and 7,063 copy number variations. Many mutations are located in genes with functions relevant to bioprocessing, such as apoptosis. The details of this genetic diversity highlight the value of the hamster genome as the reference upon which CHO cells can be studied and engineered for protein production.

Large-scale in silico modeling of metabolic interactions between cell types in the human brain

Metabolic interactions between multiple cell types are difficult to model using existing approaches. Here we present a workflow that integrates gene expression data, proteomics data and literature-based manual curation to model human metabolism within and between different types of cells. Transport reactions are used to account for the transfer of metabolites between models of different cell types via the interstitial fluid. We apply the method to create models of brain energy metabolism that recapitulate metabolic interactions between astrocytes and various neuron types relevant to Alzheimers disease. Analysis of the models identifies genes and pathways that may explain observed experimental phenomena, including the differential effects of the disease on cell types and regions of the brain. Constraint-based modeling can thus contribute to the study and analysis of multicellular metabolic processes in the human tissue microenvironment and provide detailed mechanistic insight into high-throughput data analysis.A workflow is presented that integrates gene expression data, proteomic data, and literature-based manual curation to construct multicellular, tissue-specific models of human brain energy metabolism that recapitulate metabolic interactions between astrocytes and various neuron types. Three analyses are applied for gene identification, analysis of omics data, and analysis of physiological states. First, we identify glutamate decarboxylase as a target that may contribute to cell-type and regional specificity in Alzheimer’s disease. Second, the decreased metabolic rate seen in affected brain regions in Alzheimer’s disease is consistent with a suppression of central metabolic gene expression in histopathologically normal neurons. Third, we identify pathways in cholinergic neurons that couple mitochondrial metabolism and cytosolic acetylcholine production, and subsequently find that cholinergic neurotransmission accounts for ∼3% of brain neurotransmission. Constraint-based modeling can thus contribute to the study and analysis of multicellular metabolic processes in human tissues, and provide detailed mechanistic insight into high-throughput data analysis.

Insight into human alveolar macrophage and M. tuberculosis interactions via metabolic reconstructions.

Metabolic coupling of Mycobacterium tuberculosis to its host is foundational to its pathogenesis. Computational genome‐scale metabolic models have shown utility in integrating ‐omic as well as physiologic data for systemic, mechanistic analysis of metabolism. To date, integrative analysis of host–pathogen interactions using in silico mass‐balanced, genome‐scale models has not been performed. We, therefore, constructed a cell‐specific alveolar macrophage model, iAB‐AMØ‐1410, from the global human metabolic reconstruction, Recon 1. The model successfully predicted experimentally verified ATP and nitric oxide production rates in macrophages. This model was then integrated with an M. tuberculosis H37Rv model, iNJ661, to build an integrated host–pathogen genome‐scale reconstruction, iAB‐AMØ‐1410‐Mt‐661. The integrated host–pathogen network enables simulation of the metabolic changes during infection. The resulting reaction activity and gene essentiality targets of the integrated model represent an altered infectious state. High‐throughput data from infected macrophages were mapped onto the host–pathogen network and were able to describe three distinct pathological states. Integrated host–pathogen reconstructions thus form a foundation upon which understanding the biology and pathophysiology of infections can be developed.

A multi-tissue type genome-scale metabolic network for analysis of whole-body systems physiology

BackgroundGenome-scale metabolic reconstructions provide a biologically meaningful mechanistic basis for the genotype-phenotype relationship. The global human metabolic network, termed Recon 1, has recently been reconstructed allowing the systems analysis of human metabolic physiology and pathology. Utilizing high-throughput data, Recon 1 has recently been tailored to different cells and tissues, including the liver, kidney, brain, and alveolar macrophage. These models have shown utility in the study of systems medicine. However, no integrated analysis between human tissues has been done.ResultsTo describe tissue-specific functions, Recon 1 was tailored to describe metabolism in three human cells: adipocytes, hepatocytes, and myocytes. These cell-specific networks were manually curated and validated based on known cellular metabolic functions. To study intercellular interactions, a novel multi-tissue type modeling approach was developed to integrate the metabolic functions for the three cell types, and subsequently used to simulate known integrated metabolic cycles. In addition, the multi-tissue model was used to study diabetes: a pathology with systemic properties. High-throughput data was integrated with the network to determine differential metabolic activity between obese and type II obese gastric bypass patients in a whole-body context.ConclusionThe multi-tissue type modeling approach presented provides a platform to study integrated metabolic states. As more cell and tissue-specific models are released, it is critical to develop a framework in which to study their interdependencies.

Model-driven multi-omic data analysis elucidates metabolic immunomodulators of macrophage activation.

Macrophages are central players in immune response, manifesting divergent phenotypes to control inflammation and innate immunity through release of cytokines and other signaling factors. Recently, the focus on metabolism has been reemphasized as critical signaling and regulatory pathways of human pathophysiology, ranging from cancer to aging, often converge on metabolic responses. Here, we used genome‐scale modeling and multi‐omics (transcriptomics, proteomics, and metabolomics) analysis to assess metabolic features that are critical for macrophage activation. We constructed a genome‐scale metabolic network for the RAW 264.7 cell line to determine metabolic modulators of activation. Metabolites well‐known to be associated with immunoactivation (glucose and arginine) and immunosuppression (tryptophan and vitamin D3) were among the most critical effectors. Intracellular metabolic mechanisms were assessed, identifying a suppressive role for de‐novo nucleotide synthesis. Finally, underlying metabolic mechanisms of macrophage activation are identified by analyzing multi‐omic data obtained from LPS‐stimulated RAW cells in the context of our flux‐based predictions. Our study demonstrates metabolisms role in regulating activation may be greater than previously anticipated and elucidates underlying connections between activation and metabolic effectors.

Using the reconstructed genome‐scale human metabolic network to study physiology and pathology

Abstract.  Bordbar A, Palsson BO (University of California San Diego, La Jolla, CA, USA). Using the reconstructed genome‐scale human metabolic network to study physiology and pathology (Key Symposium). J Intern Med 2012; 271: 131–141.

Multiscale Modeling of Metabolism and Macromolecular Synthesis in E. coli and Its Application to the Evolution of Codon Usage

Biological systems are inherently hierarchal and multiscale in time and space. A major challenge of systems biology is to describe biological systems as a computational model, which can be used to derive novel hypothesis and drive experiments leading to new knowledge. The constraint-based reconstruction and analysis approach has been successfully applied to metabolism and to the macromolecular synthesis machinery assembly. Here, we present the first integrated stoichiometric multiscale model of metabolism and macromolecular synthesis for Escherichia coli K12 MG1655, which describes the sequence-specific synthesis and function of almost 2000 gene products at molecular detail. We added linear constraints, which couple enzyme synthesis and catalysis reactions. Comparison with experimental data showed improvement of growth phenotype prediction with the multiscale model over E. coli’s metabolic model alone. Many of the genes covered by this integrated model are well conserved across enterobacters and other, less related bacteria. We addressed the question of whether the bias in synonymous codon usage could affect the growth phenotype and environmental niches that an organism can occupy. We created two classes of in silico strains, one with more biased codon usage and one with more equilibrated codon usage than the wildtype. The reduced growth phenotype in biased strains was caused by tRNA supply shortage, indicating that expansion of tRNA gene content or tRNA codon recognition allow E. coli to respond to changes in codon usage bias. Our analysis suggests that in order to maximize growth and to adapt to new environmental niches, codon usage and tRNA content must co-evolve. These results provide further evidence for the mutation-selection-drift balance theory of codon usage bias. This integrated multiscale reconstruction successfully demonstrates that the constraint-based modeling approach is well suited to whole-cell modeling endeavors.

Identified metabolic signature for assessing red blood cell unit quality is associated with endothelial damage markers and clinical outcomes

There has been interest in determining whether older red blood cell (RBC) units have negative clinical effects. Numerous observational studies have shown that older RBC units are an independent factor for patient mortality. However, recently published randomized clinical trials have shown no difference of clinical outcome for patients receiving old or fresh RBCs. An overlooked but essential issue in assessing RBC unit quality and ultimately designing the necessary clinical trials is a metric for what constitutes an old or fresh RBC unit.