The influence of the high intracellular concentration of macromolecules on cell physiology is increasingly appreciated, but its impact on system-level cellular functions remains poorly quantified. To assess its potential effect, here we develop a flux balance model of Escherichia coli cell metabolism that takes into account a systemslevel constraint for the concentration of enzymes catalyzing the various metabolic reactions in the crowded cytoplasm. We demonstrate that the model's predictions for the relative maximum growth rate of wild-type and mutant E. coli cells in single substratelimited media, and the sequence and mode of substrate uptake and utilization from a complex medium are in good agreement with subsequent experimental observations. These results suggest that molecular crowding represents a bound on the achievable functional states of a metabolic network, and they indicate that models incorporating this constraint can systematically identify alterations in cellular metabolism activated in response to environmental change.flux balance analysis ͉ metabolic networks ͉ systems biology A n important aim of systems biology is the identification of the organizing principles and fundamental constraints that characterize the function of molecular interaction networks and the limits of an organism's phenotypic diversity (1-3). Flux balancebased modeling approaches, combining the constraints imposed by the metabolic network's structure with, e.g., mass-or energyconservation principles (3-6), are especially successful in providing experimentally testable predictions on an organism's metabolic flux state and growth rate. A relative shortcoming of these approaches, however, is that they do not take into account the physical and spatial constraints resulting from the cell's unique intracellular milieu (7-9). For example, Ϸ20-30% of the Escherichia coli cytoplasm is occupied by macromolecules, many of them enzymes, whose cytoplasmic concentration cannot be further increased without drastically affecting protein folding, protein-protein association rates, biochemical reaction kinetics, and transport dynamics within a cell (9, 10). This suggests that constraint-based modeling approaches, such as flux balance analysis (FBA) (3, 11), could be improved if we take into account that the enzymes catalyzing each reaction compete for the available cytoplasmic space (12, 13), potentially limiting the attainable flux rates.Current flux balance-based modeling approaches also have limited ability to predict substrate uptake from the environment. Extensive experimental data indicates that when grown in complex medium bacterial cells use the available substrates either preferentially or simultaneously depending on the growth condition (see, e.g., refs. 14-17). Efforts to model mixed-substrate growth have assumed specific kinetic expressions for substrate uptake and biomass growth rates, and their predictions are formulated in terms of known model parameters (15,18). Similarly, FBA predictions are based on previous knowledge of the...