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 Ϸ70 -80% of the conditions considered. S ystems biology is commonly viewed as a four-step procedure (1-3): (i) the enumeration of the biological components that make up the biological process of interest, (ii) the reconstruction of the network of interactions among these components, (iii) the in silico simulation of the network function, and (iv) the comparison of computed network properties with actual phenotypic observations. A wealth of available biological data for Saccharomyces cerevisiae (4-7) led to the establishment of the first genome-scale reconstruction of the metabolic network in a eukaryotic cell (8). The initial completion of the first two steps of the systems biology procedure for yeast metabolism has been achieved. Steps iii and iv have not yet been carried out for a eukaryotic organism.Integrated functions of reconstructed metabolic networks can be determined in silico by using a number of analytical approaches (9-11). The relatively young constraint-based approach differs fundamentally from the more traditional kinetic theory-based approaches in that it is not aimed at finding the solution or behavior of the network under certain conditions but rather at eliminating solutions (behaviors) that the network cannot exhibit (Fig. 1). By using this approach, a network of interactions is successively constrained by defining the stoichiometry of the interacting components, the direction of network reactions, and the maximum allowable throughput. In this way, candidate solutions to the network equations are systematically eliminated by the successive application of the governing constraints, i.e., stoichiometry, thermodynamics, and maximal enzymatic rates (12). One thus can define the range of capabilities of the reconstructed network and then, through the use of optimization procedur...