Cells react to extracellular perturbations with complex and intertwined responses. Systematic identification of the regulatory mechanisms that control these responses is still a challenge and requires tailored analyses integrating different types of molecular data. Here we acquired time-resolved metabolomics measurements in yeast under salt and pheromone stimulation and developed a machine learning approach to explore regulatory associations between metabolism and signal transduction. Existing phosphoproteomics measurements under the same conditions and kinase-substrate regulatory interactions were used to in silico estimate the enzymatic activity of signalling kinases. Our approach identified informative associations between kinases and metabolic enzymes capable of predicting metabolic changes. We extended our analysis to two studies containing transcriptomics, phosphoproteomics and metabolomics measurements across a comprehensive panel of kinases/phosphatases knockouts and time-resolved perturbations to the nitrogen metabolism. Changes in activity of transcription factors, kinases and phosphatases were estimated in silico and these were capable of building predictive models to infer the metabolic adaptations of previously unseen conditions across different dynamic experiments. Time-resolved experiments were significantly more informative than genetic perturbations to infer metabolic adaptation. This difference may be due to the indirect nature of the associations and of general cellular states that can hinder the identification of causal relationships. This work provides a novel genome-scale integrative analysis to propose putative transcriptional and post-translational regulatory mechanisms of metabolic processes.

MOTIVATION: Simple forms of mutualism between microorganisms are widespread in nature. Nevertheless, the role played by the environmental nutrient composition in mediating cross-feeding in microbial ecosystems is still poorly understood. RESULTS: Here, we use mixed-integer bilevel linear programming to investigate the cost of sharing metabolic resources in microbial communities. The algorithm infers an optimal combination of nutrients that can selectively sustain synergistic growth for a pair of species and guarantees minimum cost of cross-fed metabolites. To test model-based predictions, we selected a pair of Escherichia coli single gene knockouts auxotrophic, respectively, for arginine and leucine: ΔargB and ΔleuB and we experimentally verified that model-predicted medium composition significantly favors mutualism. Moreover, mass spectrometry profiling of exchanged metabolites confirmed the predicted cross-fed metabolites, supporting our constraint based modeling approach as a promising tool for engineering microbial consortia. AVAILABILITY AND IMPLEMENTATION: The software is freely available as a matlab script in the Supplementary materials. CONTACT: [email protected] SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Understanding the structure and function of complex gene regulatory networks using classical genetic assays is an error-prone procedure that frequently generates ambiguous outcomes. Even some of the best-characterized gene networks contain interactions whose validity is not conclusively proven. Founded on dynamic experimental data, mechanistic mathematical models are able to offer detailed insights that would otherwise require prohibitively large numbers of genetic experiments. Here we attempt mechanistic modeling of the transcriptional network formed by the four GATA-factor proteins, a well-studied system of central importance for nitrogen-source regulation of transcription in the yeast Saccharomyces cerevisiae. To resolve ambiguities in the network organization, we encoded a set of five interactions hypothesized in the literature into a set of 32 mathematical models, and employed Bayesian model selection to identify the most plausible set of interactions based on dynamic gene expression data. The top-ranking model was validated on newly generated GFP reporter dynamic data and was subsequently used to gain a better understanding of how yeast cells organize their transcriptional response to dynamic changes of nitrogen sources. Our work constitutes a necessary and important step towards obtaining a holistic view of the yeast nitrogen regulation mechanisms; on the computational side, it provides a demonstration of how powerful Monte Carlo techniques can be creatively combined and used to address the great challenges of large-scale dynamical system inference.

MOTIVATION: A common problem in understanding a biochemical system is to infer its correct structure or topology. This topology consists of all relevant state variables-usually molecules and their interactions. Here we present a method called topological augmentation to infer this structure in a statistically rigorous and systematic way from prior knowledge and experimental data. RESULTS: Topological augmentation starts from a simple model that is unable to explain the experimental data and augments its topology by adding new terms that capture the experimental behavior. This process is guided by representing the uncertainty in the model topology through stochastic differential equations whose trajectories contain information about missing model parts. We first apply this semiautomatic procedure to a pharmacokinetic model. This example illustrates that a global sampling of the parameter space is critical for inferring a correct model structure. We also use our method to improve our understanding of glutamine transport in yeast. This analysis shows that transport dynamics is determined by glutamine permeases with two different kinds of kinetics. Topological augmentation can not only be applied to biochemical systems, but also to any system that can be described by ordinary differential equations. AVAILABILITY AND IMPLEMENTATION: Matlab code and examples are available at: http://www.csb.ethz.ch/tools/index

BACKGROUND: Current approaches to parameter estimation are often inappropriate or inconvenient for the modelling of complex biological systems. For systems described by nonlinear equations, the conventional approach is to first numerically integrate the model, and then, in a second a posteriori step, check for consistency with experimental constraints. Hence, only single parameter sets can be considered at a time. Consequently, it is impossible to conclude that the "best" solution was identified or that no good solution exists, because parameter spaces typically cannot be explored in a reasonable amount of time. RESULTS: We introduce a novel approach based on semidefinite programming to directly identify consistent steady state concentrations for systems consisting of mass action kinetics, i.e., polynomial equations and inequality constraints. The duality properties of semidefinite programming allow to rigorously certify infeasibility for whole regions of parameter space, thus enabling the simultaneous multi-dimensional analysis of entire parameter sets. CONCLUSION: Our algorithm reduces the computational effort of parameter estimation by several orders of magnitude, as illustrated through conceptual sample problems. Of particular relevance for systems biology, the approach can discriminate between structurally different candidate models by proving inconsistency with the available data.

BACKGROUND: Quantitative knowledge of intracellular fluxes is important for a comprehensive characterization of metabolic networks and their functional operation. In contrast to direct assessment of metabolite concentrations, in vivo metabolite fluxes must be inferred indirectly from measurable quantities in 13C experiments. The required experience, the complicated network models, large and heterogeneous data sets, and the time-consuming set-up of highly controlled experimental conditions largely restricted metabolic flux analysis to few expert groups. A conceptual simplification of flux analysis is the analytical determination of metabolic flux ratios exclusively from MS data, which can then be used in a second step to estimate absolute in vivo fluxes. RESULTS: Here we describe the user-friendly software package FiatFlux that supports flux analysis for non-expert users. In the first module, ratios of converging fluxes are automatically calculated from GC-MS-detected 13C-pattern in protein-bound amino acids. Predefined fragmentation patterns are automatically identified and appropriate statistical data treatment is based on the comparison of redundant information in the MS spectra. In the second module, absolute intracellular fluxes may be calculated by a 13C-constrained flux balancing procedure that combines experimentally determined fluxes in and out of the cell and the above flux ratios. The software is preconfigured to derive flux ratios and absolute in vivo fluxes from [1-13C] and [U-13C]glucose experiments and GC-MS analysis of amino acids for a variety of microorganisms. CONCLUSION: FiatFlux is an intuitive tool for quantitative investigations of intracellular metabolism by users that are not familiar with numerical methods or isotopic tracer experiments. The aim of this open source software is to enable non-specialists to adapt the software to their specific scientific interests, including other 13C-substrates, labeling mixtures, and organisms.