The in vivo distribution of metabolic fluxes in Escherichia coli can be predicted from optimality principles At least two different sets of optimality principles govern the operation of the metabolic network under different environmental conditionsMetabolism during unlimited growth on glucose in batch culture is best described by the nonlinear maximization of ATP yield per unit of flux

The roles of duplicate genes and their contribution to the phenomenon of enzyme dispensability are a central issue in molecular and genome evolution. A comprehensive classification of the mechanisms that may have led to their preservation, however, is currently lacking. In a systems biology approach, we classify here back-up, regulatory, and gene dosage functions for the 105 duplicate gene families of Saccharomyces cerevisiae metabolism. The key tool was the reconciled genome-scale metabolic model iLL672, which was based on the older iFF708. Computational predictions of all metabolic gene knockouts were validated with the experimentally determined phenotypes of the entire singleton yeast library of 4658 mutants under five environmental conditions. iLL672 correctly identified 96%-98% and 73%-80% of the viable and lethal singleton phenotypes, respectively. Functional roles for each duplicate family were identified by integrating the iLL672-predicted in silico duplicate knockout phenotypes, genome-scale carbon-flux distributions, singleton mutant phenotypes, and network topology analysis. The results provide no evidence for a particular dominant function that maintains duplicate genes in the genome. In particular, the back-up function is not favored by evolutionary selection because duplicates do not occur more frequently in essential reactions than singleton genes. Instead of a prevailing role, multigene-encoded enzymes cover different functions. Thus, at least for metabolism, persistence of the paralog fraction in the genome can be better explained with an array of different, often overlapping functional roles.

The aim of this tutorial is to introduce the fundamental concepts of physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) modeling with a special focus on their practical implementation in a typical PBPK model building workflow. To illustrate basic steps in PBPK model building, a PBPK model for ciprofloxacin will be constructed and coupled to a pharmacodynamic model to simulate the antibacterial activity of ciprofloxacin treatment.

Systems biology iteratively combines experimentation with mathematical modeling. However, limited mechanistic knowledge, conflicting hypotheses and scarce experimental data severely hamper the development of predictive mechanistic models in many areas of biology. Even under such high uncertainty, we show here that ensemble modeling, when combined with targeted experimental analysis, can unravel key operating principles in complex cellular pathways. For proof of concept, we develop a library of mechanistically alternative dynamic models for the highly conserved target-of-rapamycin (TOR) pathway of Saccharomyces cerevisiae. In contrast to the prevailing view of a de novo assembly of type 2A phosphatases (PP2As), our integrated computational and experimental analysis proposes a specificity factor, based on Tap42p-Tip41p, for PP2As as the key signaling mechanism that is quantitatively consistent with all available experimental data. Beyond revising our picture of TOR signaling, we expect ensemble modeling to help elucidate other insufficiently characterized cellular circuits.

Large-scale 13 C -flux analysis in yeast

Genome-scale 13^C-flux analysis in Saccharomyces cerevisiae revealed that the apparent dispensability of knockout mutants with metabolic function can be explained by gene inactivity under a particular condition, by network redundancy through duplicated genes or by alternative pathways.