Data integration strategies for whole-cell modeling.

Data makes the world go round-and high quality data is a prerequisite for precise models, especially for whole-cell models (WCM). Data for WCM must be reusable, contain information about the exact experimental background, and should-in its entirety-cover all relevant processes in the cell. Here, we review basic requirements to data for WCM and strategies how to combine them. As a species-specific resource, we introduce the Yeast Cell Model Data Base (YCMDB) to illustrate requirements and solutions. We discuss recent standards for data as well as for computational models including the modeling process as data to be reported. We outline strategies for constructions of WCM despite their inherent complexity.

A transcriptome-wide analysis deciphers distinct roles of G1 cyclins in temporal organization of the yeast cell cycle.

Oscillating gene expression is crucial for correct timing and progression through cell cycle. In Saccharomyces cerevisiae, G1 cyclins Cln1-3 are essential drivers of the cell cycle and have an important role for temporal fine-tuning. We measured time-resolved transcriptome-wide gene expression for wild type and cyclin single and double knockouts over cell cycle with and without osmotic stress. Clustering of expression profiles, peak time detection of oscillating genes, integration with transcription factor network dynamics, and assignment to cell cycle phases allowed us to quantify the effect of genetic or stress perturbations on the duration of cell cycle phases. Cln1 and Cln2 showed functional differences, especially affecting later phases. Deletion of Cln3 led to a delay of START followed by normal progression through later phases. Our data and network analysis suggest mutual effects of cyclins with the transcriptional regulators SBF and MBF.

Computational Modeling of Lipid Metabolism in Yeast.

Lipid metabolism is essential for all major cell functions and has recently gained increasing attention in research and health studies. However, mathematical modeling by means of classical approaches such as stoichiometric networks and ordinary differential equation systems has not yet provided satisfactory insights, due to the complexity of lipid metabolism characterized by many different species with only slight differences and by promiscuous multifunctional enzymes. Here, we present an object-oriented stochastic model approach as a way to cope with the complex lipid metabolic network. While all lipid species are treated objects in the model, they can be modified by the respective converting reactions based on reaction rules, a hybrid method that integrates benefits of agent-based and classical stochastic simulation. This approach allows to follow the dynamics of all lipid species with different fatty acids, different degrees of saturation and different headgroups over time and to analyze the effect of parameter changes, potential mutations in the catalyzing enzymes or provision of different precursors. Applied to yeast metabolism during one cell cycle period, we could analyze the distribution of all lipids to the various membranes in time-dependent manner. The presented approach allows to efficiently treat the complexity of cellular lipid metabolism and to derive conclusions on the time- and location-dependent distributions of lipid species and their properties such as saturation. It is widely applicable, easily extendable and will provide further insights in healthy and diseased states of cell metabolism.

Global characterization of in vivo enzyme catalytic rates and their correspondence to in vitro kcat measurements.

Turnover numbers, also known as kcat values, are fundamental properties of enzymes. However, kcat data are scarce and measured in vitro, thus may not faithfully represent the in vivo situation. A basic question that awaits elucidation is: how representative are kcat values for the maximal catalytic rates of enzymes in vivo? Here, we harness omics data to calculate kmax(vivo), the observed maximal catalytic rate of an enzyme inside cells. Comparison with kcat values from Escherichia coli, yields a correlation ofr(2)= 0.62 in log scale (p < 10(-10)), with a root mean square difference of 0.54 (3.5-fold in linear scale), indicating that in vivo and in vitro maximal rates generally concur. By accounting for the degree of saturation of enzymes and the backward flux dictated by thermodynamics, we further refine the correspondence between kmax(vivo) and kcat values. The approach we present here characterizes the quantitative relationship between enzymatic catalysis in vitro and in vivo and offers a high-throughput method for extracting enzyme kinetic constants from omics data.

Value or waste: Perfusion imaging following radiofrequency ablation - early experience.

BACKGROUND: Radiofrequency ablation (RFA) is an evolving technique in treatment of hepatic malignant tumors. By heating local tissue it leads to coagulative necrotic areas around the ablation probe. Temperature falls with increasing distance to the probe, risking incomplete necrosis at the margins of the RFA-induced lesion. Therefore, immediate non-invasive and precise detection of incomplete ablation is necessary for early enlargement of the ablation if needed. OBJECTIVES: This in vivo pig study compares early experiences of immediate post-interventional computed tomography (CT) perfusion volume analysis to macroscopic and CT image evaluation in healthy pig liver. MATERIAL AND METHODS: RFA was performed in vivo in healthy pig livers. Different CT perfusion algorithms (Maximum slope analysis and Patlak plot) were used to quantify three different perfusion parameters. Data points were acquired from rectangular grids. These grids were semiautomatically overlayed to macroscopic images documented after liver explantation. Each data point was visually assigned to zones defined as "inner" and "outer necrotic zone", "margin" or "vital tissue". RESULTS: Significant differences between necrotic zones and vital tissue are shown for equivalent blood volume (p <  0.0001), arterial flow (p <  0.01) and flow extraction product (p <  0.001). Looking at equivalent blood volume and flow extraction product, there were also significant differences (EquivBV: p <  0.0001, FE: p <  0.001) between margins, necrotic and vital areas. CONCLUSIONS: In a porcine model these early results could show that all of the used CT perfusion parameters allowed discrimination of necrosis from vital tissue after RFA at high levels of significance. In addition, the parameters EquivBV and FE that give an estimate of the tissue blood volume and the permeability, were able to precisely discern different zones also seen macroscopically. From this data CT perfusion analysis could be precise tool for measurement and visualization of ablated liver lesions and for immediate detection of incomplete ablation areas.

Dynamic metabolic models in context: biomass backtracking.

Mathematical modeling has proven to be a powerful tool to understand and predict functional and regulatory properties of metabolic processes. High accuracy dynamic modeling of individual pathways is thereby opposed by simplified but genome scale constraint based approaches. A method that links these two powerful techniques would greatly enhance predictive power but is so far lacking. We present biomass backtracking, a workflow that integrates the cellular context in existing dynamic metabolic models via stoichiometrically exact drain reactions based on a genome scale metabolic model. With comprehensive examples, for different species and environmental contexts, we show the importance and scope of applications and highlight the improvement compared to common boundary formulations in existing metabolic models. Our method allows for the contextualization of dynamic metabolic models based on all available information. We anticipate this to greatly increase their accuracy and predictive power for basic research and also for drug development and industrial applications.

New types of experimental data shape the use of enzyme kinetics for dynamic network modeling.

Since the publication of Leonor Michaelis and Maude Menten's paper on the reaction kinetics of the enzyme invertase in 1913, molecular biology has evolved tremendously. New measurement techniques allow in vivo characterization of the whole genome, proteome or transcriptome of cells, whereas the classical enzyme essay only allows determination of the two Michaelis-Menten parameters V and K(m). Nevertheless, Michaelis-Menten kinetics are still commonly used, not only in the in vitro context of enzyme characterization but also as a rate law for enzymatic reactions in larger biochemical reaction networks. In this review, we give an overview of the historical development of kinetic rate laws originating from Michaelis-Menten kinetics over the past 100 years. Furthermore, we briefly summarize the experimental techniques used for the characterization of enzymes, and discuss web resources that systematically store kinetic parameters and related information. Finally, describe the novel opportunities that arise from using these data in dynamic mathematical modeling. In this framework, traditional in vitro approaches may be combined with modern genome-scale measurements to foster thorough understanding of the underlying complex mechanisms.

Single-cell analysis of ribonucleotide reductase transcriptional and translational response to DNA damage.

The ribonucleotide reductase (RNR) enzyme catalyzes an essential step in the production of deoxyribonucleotide triphosphates (dNTPs) in cells. Bulk biochemical measurements in synchronized Saccharomyces cerevisiae cells suggest that RNR mRNA production is maximal in late G(1) and S phases; however, damaged DNA induces RNR transcription throughout the cell cycle. But such en masse measurements reveal neither cell-to-cell heterogeneity in responses nor direct correlations between transcript and protein expression or localization in single cells which may be central to function. We overcame these limitations by simultaneous detection of single RNR transcripts and also Rnr proteins in the same individual asynchronous S. cerevisiae cells, with and without DNA damage by methyl methanesulfonate (MMS). Surprisingly, RNR subunit mRNA levels were comparably low in both damaged and undamaged G(1) cells and highly induced in damaged S/G(2) cells. Transcript numbers became correlated with both protein levels and localization only upon DNA damage in a cell cycle-dependent manner. Further, we showed that the differential RNR response to DNA damage correlated with variable Mec1 kinase activity in the cell cycle in single cells. The transcription of RNR genes was found to be noisy and non-Poissonian in nature. Our results provide vital insight into cell cycle-dependent RNR regulation under conditions of genotoxic stress.