Large-scale evaluation of in silico gene deletions in Saccharomyces cerevisiae.

A large-scale in silico evaluation of gene deletions in Saccharomyces cerevisiae was conducted using a genome-scale reconstructed metabolic model. The effect of 599 single gene deletions on cell viability was simulated in silico and compared to published experimental results. In 526 cases (87.8%), the in silico results were in agreement with experimental observations when growth on synthetic complete medium was simulated. Viable phenotypes were correctly predicted in 89.4% (496 out of 555) and lethal phenotypes were correctly predicted in 68.2% (30 out of 44) of the cases considered. The in silico evaluation was solely based on the topological properties of the metabolic network which is based on well-established reaction stoichiometry. No interaction or regulatory information was accounted for in the in silico model. False predictions were analyzed on a case-by-case basis for four possible inadequacies of the in silico model: (1) incomplete media composition, (2) substitutable biomass components, (3) incomplete biochemical information, and (4) missing regulation. This analysis eliminated a number of false predictions and suggested a number of experimentally testable hypotheses. A genome-scale in silico model can thus be used to systematically reconcile existing data and fill in our knowledge gaps about an organism.

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.

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

BACKGROUND: Genome-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. RESULTS: To 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. CONCLUSION: The 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.

The genomic sequence of the Chinese hamster ovary (CHO)-K1 cell line.

Chinese hamster ovary (CHO)-derived cell lines are the preferred host cells for the production of therapeutic proteins. Here we present a draft genomic sequence of the CHO-K1 ancestral cell line. The assembly comprises 2.45 Gb of genomic sequence, with 24,383 predicted genes. We associate most of the assembled scaffolds with 21 chromosomes isolated by microfluidics to identify chromosomal locations of genes. Furthermore, we investigate genes involved in glycosylation, which affect therapeutic protein quality, and viral susceptibility genes, which are relevant to cell engineering and regulatory concerns. Homologs of most human glycosylation-associated genes are present in the CHO-K1 genome, although 141 of these homologs are not expressed under exponential growth conditions. Many important viral entry genes are also present in the genome but not expressed, which may explain the unusual viral resistance property of CHO cell lines. We discuss how the availability of this genome sequence may facilitate genome-scale science for the optimization of biopharmaceutical protein production.

Genome-scale constraint-based modeling of Geobacter metallireducens.

BACKGROUND: Geobacter metallireducens was the first organism that can be grown in pure culture to completely oxidize organic compounds with Fe(III) oxide serving as electron acceptor. Geobacter species, including G. sulfurreducens and G. metallireducens, are used for bioremediation and electricity generation from waste organic matter and renewable biomass. The constraint-based modeling approach enables the development of genome-scale in silico models that can predict the behavior of complex biological systems and their responses to the environments. Such a modeling approach was applied to provide physiological and ecological insights on the metabolism of G. metallireducens. RESULTS: The genome-scale metabolic model of G. metallireducens was constructed to include 747 genes and 697 reactions. Compared to the G. sulfurreducens model, the G. metallireducens metabolic model contains 118 unique reactions that reflect many of G. metallireducens' specific metabolic capabilities. Detailed examination of the G. metallireducens model suggests that its central metabolism contains several energy-inefficient reactions that are not present in the G. sulfurreducens model. Experimental biomass yield of G. metallireducens growing on pyruvate was lower than the predicted optimal biomass yield. Microarray data of G. metallireducens growing with benzoate and acetate indicated that genes encoding these energy-inefficient reactions were up-regulated by benzoate. These results suggested that the energy-inefficient reactions were likely turned off during G. metallireducens growth with acetate for optimal biomass yield, but were up-regulated during growth with complex electron donors such as benzoate for rapid energy generation. Furthermore, several computational modeling approaches were applied to accelerate G. metallireducens research. For example, growth of G. metallireducens with different electron donors and electron acceptors were studied using the genome-scale metabolic model, which provided a fast and cost-effective way to understand the metabolism of G. metallireducens. CONCLUSION: We have developed a genome-scale metabolic model for G. metallireducens that features both metabolic similarities and differences to the published model for its close relative, G. sulfurreducens. Together these metabolic models provide an important resource for improving strategies on bioremediation and bioenergy generation.

Towards multidimensional genome annotation.

Our information about the gene content of organisms continues to grow as more genomes are sequenced and gene products are characterized. Sequence-based annotation efforts have led to a list of cellular components, which can be thought of as a one-dimensional annotation. With growing information about component interactions, facilitated by the advancement of various high-throughput technologies, systemic, or two-dimensional, annotations can be generated. Knowledge about the physical arrangement of chromosomes will lead to a three-dimensional spatial annotation of the genome and a fourth dimension of annotation will arise from the study of changes in genome sequences that occur during adaptive evolution. Here we discuss all four levels of genome annotation, with specific emphasis on two-dimensional annotation methods.

k-Cone analysis: determining all candidate values for kinetic parameters on a network scale.

The absence of comprehensive measured kinetic values and the observed inconsistency in the available in vitro kinetic data has hindered the formulation of network-scale kinetic models of biochemical reaction networks. To meet this challenge we present an approach to construct a convex space, termed the k-cone, which contains all the allowable numerical values of the kinetic constants in large-scale biochemical networks. The definition of the k-cone relies on the incorporation of in vivo concentration data and a simplified approach to represent enzyme kinetics within an established constraint-based modeling approach. The k-cone approach was implemented to define the allowable combination of numerical values for a full kinetic model of human red blood cell metabolism and to study its correlated kinetic parameters. The k-cone approach can be used to determine consistency between in vitro measured kinetic values and in vivo concentration and flux measurements when used in a network-scale kinetic model. k-Cone analysis was successful in determining whether in vitro measured kinetic values used in the reconstruction of a kinetic-based model of Saccharomyces cerevisiae central metabolism could reproduce in vivo measurements. Further, the k-cone can be used to determine which numerical values of in vitro measured parameters are required to be changed in a kinetic model if in vivo measured values are not reproduced. k-Cone analysis could identify what minimum number of in vitro determined kinetic parameters needed to be adjusted in the S. cerevisiae model to be consistent with the in vivo data. Applying the k-cone analysis a priori to kinetic model development may reduce the time and effort involved in model building and parameter adjustment. With the recent developments in high-throughput profiling of metabolite concentrations at a whole-cell scale and advances in metabolomics technologies, the k-cone approach presented here may hold the promise for kinetic characterization of metabolic networks as well as other biological functions at a whole-cell level.

The convex basis of the left null space of the stoichiometric matrix leads to the definition of metabolically meaningful pools.

The stoichiometric matrix, S, represents a mapping of reaction rate vectors into a space of concentration time derivatives. The left null space of the stoichiometric matrix contains the dynamic invariants: a combination of concentration variables, referred to as metabolic pools, whose total concentration does not change over time. By analogy to the traditional reaction map formed by S, a compound map can be derived from -S(T). The analogy to flux analysis of the (right) null space of S enables us to classify the metabolic pools into three categories: Type A that contains chemical elements and their combinations in the form of certain moieties, Type B that contains such moieties in addition to cofactors carrying such moieties that are internal to the network, and Type C that contains only the cofactors. A convex formulation of the basis for the left null space allows us to directly classify the metabolic pools into these three categories. Type B metabolic pools include conservation pools that form conjugates of moiety-occupied and moiety-vacant concentration states of metabolites and cofactors. Type B metabolic pools thus describe the various states of moiety exchange between the primary substrates and the cofactors that capture properties like energy and redox potential. The convex basis gives clear insight into this exchange for glycolytic pathway in human red blood cell, including the identification of high and low energy pools that form conjugates. Examples suggest that pool maps may be more appropriate for signaling pathways than flux maps. The analysis of the left null space of the stoichiometric matrix allows us to define the achievable states of the cell and their physiological relevance.

Systemic metabolic reactions are obtained by singular value decomposition of genome-scale stoichiometric matrices.

Genome-scale metabolic networks can be reconstructed. The systemic biochemical properties of these networks can now be studied. Here, genome-scale reconstructed metabolic networks were analysed using singular value decomposition (SVD). All the individual biochemical conversions contained in a reconstructed metabolic network are described by a stoichiometric matrix (S). SVD of S led to the definition of the underlying modes that characterize the overall biochemical conversions that take place in a network and rank-ordered their importance. The modes were shown to correspond to systemic biochemical reactions and they could be used to identify the groups and clusters of individual biochemical reactions that drive them. Comparative analysis of the Escherichia coli, Haemophilus influenzae, and Helicobacter pylori genome-scale metabolic networks showed that the four dominant modes in all three networks correspond to: (1) the conversion of ATP to ADP, (2) redox metabolism of NADP, (3) proton-motive force, and (4) inorganic phosphate metabolism. The sets of individual metabolic reactions deriving these systemic conversions, however, differed among the three organisms. Thus, we can now define systemic metabolic reactions, or eigen-reactions, for the study of systems biology of metabolism and have a basis for comparing the overall properties of genome-specific metabolic networks.

Saccharomyces cerevisiae phenotypes can be predicted by using constraint-based analysis of a genome-scale reconstructed metabolic network.

Full genome sequences of prokaryotic organisms have led to reconstruction of genome-scale metabolic networks and in silico computation of their integrated functions. The first genome-scale metabolic reconstruction for a eukaryotic cell, Saccharomyces cerevisiae, consisting of 1,175 metabolic reactions and 733 metabolites, has appeared. A constraint-based in silico analysis procedure was used to compute properties of the S. cerevisiae metabolic network. The computed number of ATP molecules produced per pair of electrons donated to the electron transport system (ETS) and energy-maintenance requirements were quantitatively in agreement with experimental results. Computed whole-cell functions of growth and metabolic by-product secretion in aerobic and anaerobic culture were consistent with experimental data, and thus mRNA expression profiles during metabolic shifts were computed. The computed consequences of gene knockouts on growth phenotypes were consistent with experimental observations. Thus, constraint-based analysis of a genome-scale metabolic network for the eukaryotic S. cerevisiae allows for computation of its integrated functions, producing in silico results that were consistent with observed phenotypic functions for approximately 70-80% of the conditions considered.

Identifying constraints that govern cell behavior: a key to converting conceptual to computational models in biology?

Cells must abide by a number of constraints. The environmental constrains of cellular behavior and physicochemical limitations affect cellular processes. To regulate and adapt their functions, cells impose constraints on themselves. Enumerating, understanding, and applying these constraints leads to a constraints-based modeling formalism that has been helpful in converting conceptual models to computational models in biology. The continued success of the constraints-based approach depends upon identification and incorporation of new constraints to more accurately define cellular capabilities. This review considers constraints in terms of environmental, physicochemical, and self-imposed regulatory and evolutionary constraints with the purpose of refining current constraints-based models of cell phenotype.

Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network.

The metabolic network in the yeast Saccharomyces cerevisiae was reconstructed using currently available genomic, biochemical, and physiological information. The metabolic reactions were compartmentalized between the cytosol and the mitochondria, and transport steps between the compartments and the environment were included. A total of 708 structural open reading frames (ORFs) were accounted for in the reconstructed network, corresponding to 1035 metabolic reactions. Further, 140 reactions were included on the basis of biochemical evidence resulting in a genome-scale reconstructed metabolic network containing 1175 metabolic reactions and 584 metabolites. The number of gene functions included in the reconstructed network corresponds to approximately 16% of all characterized ORFs in S. cerevisiae. Using the reconstructed network, the metabolic capabilities of S. cerevisiae were calculated and compared with Escherichia coli. The reconstructed metabolic network is the first comprehensive network for a eukaryotic organism, and it may be used as the basis for in silico analysis of phenotypic functions.

Analysis of metabolic capabilities using singular value decomposition of extreme pathway matrices.

It is now possible to construct genome-scale metabolic networks for particular microorganisms. Extreme pathway analysis is a useful method for analyzing the phenotypic capabilities of these networks. Many extreme pathways are needed to fully describe the functional capabilities of genome-scale metabolic networks, and therefore, a need exists to develop methods to study these large sets of extreme pathways. Singular value decomposition (SVD) of matrices of extreme pathways was used to develop a conceptual framework for the interpretation of large sets of extreme pathways and the steady-state flux solution space they define. The key results of this study were: 1), convex steady-state solution cones describing the potential functions of biochemical networks can be studied using the modes generated by SVD; 2), Helicobacter pylori has a more rigid metabolic network (i.e., a lower dimensional solution space and a more dominant first singular value) than Haemophilus influenzae for the production of amino acids; and 3), SVD allows for direct comparison of different solution cones resulting from the production of different amino acids. SVD was used to identify key network branch points that may identify key control points for regulation. Therefore, SVD of matrices of extreme pathways has proved to be a useful method for analyzing the steady-state solution space of genome-scale metabolic networks.

Monte Carlo sampling can be used to determine the size and shape of the steady-state flux space.

Constraint-based modeling results in a convex polytope that defines a solution space containing all possible steady-state flux distributions. The properties of this polytope have been studied extensively using linear programming to find the optimal flux distribution under various optimality conditions and convex analysis to define its extreme pathways (edges) and elementary modes. The work presented herein further studies the steady-state flux space by defining its hyper-volume. In low dimensions (i.e. for small sample networks), exact volume calculation algorithms were used. However, due to the #P-hard nature of the vertex enumeration and volume calculation problem in high dimensions, random Monte Carlo sampling was used to characterize the relative size of the solution space of the human red blood cell metabolic network. Distributions of the steady-state flux levels for each reaction in the metabolic network were generated to show the range of flux values for each reaction in the polytope. These results give insight into the shape of the high-dimensional solution space. The value of measuring uptake and secretion rates in shrinking the steady-state flux solution space is illustrated through singular value decomposition of the randomly sampled points. The V(max) of various reactions in the network are varied to determine the sensitivity of the solution space to the maximum capacity constraints. The methods developed in this study are suitable for testing the implication of additional constraints on a metabolic network system and can be used to explore the effects of single nucleotide polymorphisms (SNPs) on network capabilities.

Genome-scale metabolic model of Helicobacter pylori 26695.

A genome-scale metabolic model of Helicobacter pylori 26695 was constructed from genome sequence annotation, biochemical, and physiological data. This represents an in silico model largely derived from genomic information for an organism for which there is substantially less biochemical information available relative to previously modeled organisms such as Escherichia coli. The reconstructed metabolic network contains 388 enzymatic and transport reactions and accounts for 291 open reading frames. Within the paradigm of constraint-based modeling, extreme-pathway analysis and flux balance analysis were used to explore the metabolic capabilities of the in silico model. General network properties were analyzed and compared to similar results previously generated for Haemophilus influenzae. A minimal medium required by the model to generate required biomass constituents was calculated, indicating the requirement of eight amino acids, six of which correspond to essential human amino acids. In addition a list of potential substrates capable of fulfilling the bulk carbon requirements of H. pylori were identified. A deletion study was performed wherein reactions and associated genes in central metabolism were deleted and their effects were simulated under a variety of substrate availability conditions, yielding a number of reactions that are deemed essential. Deletion results were compared to recently published in vitro essentiality determinations for 17 genes. The in silico model accurately predicted 10 of 17 deletion cases, with partial support for additional cases. Collectively, the results presented herein suggest an effective strategy of combining in silico modeling with experimental technologies to enhance biological discovery for less characterized organisms and their genomes.