Model Based Targeting of IL-6-Induced Inflammatory Responses in Cultured Primary Hepatocytes to Improve Application of the JAK Inhibitor Ruxolitinib.

IL-6 is a central mediator of the immediate induction of hepatic acute phase proteins (APP) in the liver during infection and after injury, but increased IL-6 activity has been associated with multiple pathological conditions. In hepatocytes, IL-6 activates JAK1-STAT3 signaling that induces the negative feedback regulator SOCS3 and expression of APPs. While different inhibitors of IL-6-induced JAK1-STAT3-signaling have been developed, understanding their precise impact on signaling dynamics requires a systems biology approach. Here we present a mathematical model of IL-6-induced JAK1-STAT3 signaling that quantitatively links physiological IL-6 concentrations to the dynamics of IL-6-induced signal transduction and expression of target genes in hepatocytes. The mathematical model consists of coupled ordinary differential equations (ODE) and the model parameters were estimated by a maximum likelihood approach, whereas identifiability of the dynamic model parameters was ensured by the Profile Likelihood. Using model simulations coupled with experimental validation we could optimize the long-term impact of the JAK-inhibitor Ruxolitinib, a therapeutic compound that is quickly metabolized. Model-predicted doses and timing of treatments helps to improve the reduction of inflammatory APP gene expression in primary mouse hepatocytes close to levels observed during regenerative conditions. The concept of improved efficacy of the inhibitor through multiple treatments at optimized time intervals was confirmed in primary human hepatocytes. Thus, combining quantitative data generation with mathematical modeling suggests that repetitive treatment with Ruxolitinib is required to effectively target excessive inflammatory responses without exceeding doses recommended by the clinical guidelines.

T160­phosphorylated CDK2 defines threshold for HGF dependent proliferation in primary hepatocytes.

Liver regeneration is a tightly controlled process mainly achieved by proliferation of usually quiescent hepatocytes. The specific molecular mechanisms ensuring cell division only in response to proliferative signals such as hepatocyte growth factor (HGF) are not fully understood. Here, we combined quantitative time-resolved analysis of primary mouse hepatocyte proliferation at the single cell and at the population level with mathematical modeling. We showed that numerous G1/S transition components are activated upon hepatocyte isolation whereas DNA replication only occurs upon additional HGF stimulation. In response to HGF, Cyclin:CDK complex formation was increased, p21 rather than p27 was regulated, and Rb expression was enhanced. Quantification of protein levels at the restriction point showed an excess of CDK2 over CDK4 and limiting amounts of the transcription factor E2F-1. Analysis with our mathematical model revealed that T160 phosphorylation of CDK2 correlated best with growth factor-dependent proliferation, which we validated experimentally on both the population and the single cell level. In conclusion, we identified CDK2 phosphorylation as a gate-keeping mechanism to maintain hepatocyte quiescence in the absence of HGF.

Dynamic mathematical modeling of IL13-induced signaling in Hodgkin and primary mediastinal B-cell lymphoma allows prediction of therapeutic targets.

Primary mediastinal B-cell lymphoma (PMBL) and classical Hodgkin lymphoma (cHL) share a frequent constitutive activation of JAK (Janus kinase)/STAT signaling pathway. Because of complex, nonlinear relations within the pathway, key dynamic properties remained to be identified to predict possible strategies for intervention. We report the development of dynamic pathway models based on quantitative data collected on signaling components of JAK/STAT pathway in two lymphoma-derived cell lines, MedB-1 and L1236, representative of PMBL and cHL, respectively. We show that the amounts of STAT5 and STAT6 are higher whereas those of SHP1 are lower in the two lymphoma cell lines than in normal B cells. Distinctively, L1236 cells harbor more JAK2 and less SHP1 molecules per cell than MedB-1 or control cells. In both lymphoma cell lines, we observe interleukin-13 (IL13)-induced activation of IL4 receptor α, JAK2, and STAT5, but not of STAT6. Genome-wide, 11 early and 16 sustained genes are upregulated by IL13 in both lymphoma cell lines. Specifically, the known STAT-inducible negative regulators CISH and SOCS3 are upregulated within 2 hours in MedB-1 but not in L1236 cells. On the basis of this detailed quantitative information, we established two mathematical models, MedB-1 and L1236 model, able to describe the respective experimental data. Most of the model parameters are identifiable and therefore the models are predictive. Sensitivity analysis of the model identifies six possible therapeutic targets able to reduce gene expression levels in L1236 cells and three in MedB-1. We experimentally confirm reduction in target gene expression in response to inhibition of STAT5 phosphorylation, thereby validating one of the predicted targets.

Theoretical analysis of time-to-peak responses in biological reaction networks.

Processing of information by signaling networks is characterized by properties of the induced kinetics of the activated pathway components. The maximal extent of pathway activation (maximum amplitude) and the time-to-peak-response (position) are key determinants of biological responses that have been linked to specific outcomes. We investigate how the maximum amplitude of pathway activation and its position depend on the input and wiring of a signaling network. For this purpose, we consider a simple reaction A→B that is regulated by a transient input and extended this to include back-reaction and additional partners. In particular, we show that a unique maximum of B(t) exists. Moreover, we prove that the position of the maximum is independent of the applied input but regulated by degradation reactions of B. Indeed, the time-to-peak-response decreases with increasing degradation rate, which we prove for small models and show in simulations for more complex ones. The identified dependencies provide insights into design principles that facilitate the realization dynamical characteristics like constant position of maximal pathway activation and thereby guide the characterization of unknown kinetics within larger protein networks.

Knowledge-based matrix factorization temporally resolves the cellular responses to IL-6 stimulation.

BACKGROUND: External stimulations of cells by hormones, cytokines or growth factors activate signal transduction pathways that subsequently induce a re-arrangement of cellular gene expression. The analysis of such changes is complicated, as they consist of multi-layered temporal responses. While classical analyses based on clustering or gene set enrichment only partly reveal this information, matrix factorization techniques are well suited for a detailed temporal analysis. In signal processing, factorization techniques incorporating data properties like spatial and temporal correlation structure have shown to be robust and computationally efficient. However, such correlation-based methods have so far not be applied in bioinformatics, because large scale biological data rarely imply a natural order that allows the definition of a delayed correlation function. RESULTS: We therefore develop the concept of graph-decorrelation. We encode prior knowledge like transcriptional regulation, protein interactions or metabolic pathways in a weighted directed graph. By linking features along this underlying graph, we introduce a partial ordering of the features (e.g. genes) and are thus able to define a graph-delayed correlation function. Using this framework as constraint to the matrix factorization task allows us to set up the fast and robust graph-decorrelation algorithm (GraDe). To analyze alterations in the gene response in IL-6 stimulated primary mouse hepatocytes, we performed a time-course microarray experiment and applied GraDe. In contrast to standard techniques, the extracted time-resolved gene expression profiles showed that IL-6 activates genes involved in cell cycle progression and cell division. Genes linked to metabolic and apoptotic processes are down-regulated indicating that IL-6 mediated priming renders hepatocytes more responsive towards cell proliferation and reduces expenditures for the energy metabolism. CONCLUSIONS: GraDe provides a novel framework for the decomposition of large-scale 'omics' data. We were able to show that including prior knowledge into the separation task leads to a much more structured and detailed separation of the time-dependent responses upon IL-6 stimulation compared to standard methods. A Matlab implementation of the GraDe algorithm is freely available at http://cmb.helmholtz-muenchen.de/grade.

Comparative proteomic phenotyping of cell lines and primary cells to assess preservation of cell type-specific functions.

Biological experiments are most often performed with immortalized cell lines because they are readily available and can be expanded without limitation. However, cell lines may differ from the in vivo situation in important aspects. Here we introduce a straightforward methodology to compare cell lines to their cognate primary cells and to derive a comparative functional phenotype. We used SILAC (stable isotope labeling by amino acids in cell culture) for quantitative, mass spectrometry-based comparison of the hepatoma cell line Hepa1-6 with primary hepatocytes. The resulting quantitative proteome of 4,063 proteins had an asymmetric distribution, with many proteins down-regulated in the cell line. Bioinformatic analysis of the quantitative proteomics phenotypes revealed that Hepa1-6 cells were deficient in mitochondria, reflecting re-arrangement of metabolic pathways, drastically up-regulate cell cycle-associated functions and largely shut down drug metabolizing enzymes characteristic for the liver. This quantitative knowledge of changes provides an important basis to adapt cell lines to more closely resemble physiological conditions.

The glucocorticoid receptor controls hepatic dyslipidemia through Hes1.

Aberrant accumulation of lipids in the liver ("fatty liver" or hepatic steatosis) represents a hallmark of the metabolic syndrome and is tightly associated with obesity, type II diabetes, starvation, or glucocorticoid (GC) therapy. While fatty liver has been connected with numerous abnormalities of liver function, the molecular mechanisms of fatty liver development remain largely enigmatic. Here we show that liver-specific disruption of glucocorticoid receptor (GR) action improves the steatotic phenotype in fatty liver mouse models and leads to the induction of transcriptional repressor hairy enhancer of split 1 (Hes1) gene expression. The GR directly interferes with Hes1 promoter activity, triggering the recruitment of histone deacetylase (HDAC) activities to the Hes1 gene. Genetic restoration of hepatic Hes1 levels in steatotic animals normalizes hepatic triglyceride (TG) levels. As glucocorticoid action is increased during starvation, myotonic dystrophy, and Cushing's syndrome, the inhibition of Hes1 through the GR might explain the fatty liver phenotype in these subjects.

Standardizing experimental protocols.

Systems biology aims at understanding the behavior of biological networks by mathematical modeling based on experimental data. However, frequently experimental data is derived from poorly defined cellular systems, the procedures of data generation are insufficiently documented and data processing is arbitrary. For the advancement of systems biology, standardization at multiple levels is essential. Several systems biology consortia have started by focusing on standardization of cellular systems and experimental procedures. Minimum information standards for the description of data sets and common languages for the description of biological pathways as well as for mathematical modeling are being developed. Standardization is required to facilitate data exchange between different research groups and finally the assembly of large integrated models providing novel biological insights.

Dynamic pathway modeling: feasibility analysis and optimal experimental design.

A major challenge in systems biology is to evaluate the feasibility of a biological research project prior to its realization. Since experiments are animals-, cost- and time-consuming, approaches allowing researchers to discriminate alternative hypotheses with a minimal set of experiments are highly desirable. Given a null hypothesis and alternative model, as well as laboratory constraints like observable players, sample size, noise level, and stimulation options, we suggest a method to obtain a list of required experiments in order to significantly reject the null hypothesis model M0 if a specified alternative model MA is realized. For this purpose, we estimate the power to detect a violation of M0 by means of Monte Carlo simulations. Iteratively, the power is maximized over all feasible stimulations of the system using multi-experiment fitting, leading to an optimal combination of experimental settings to discriminate the null hypothesis and alternative model. We prove the importance of simultaneous modeling of combined experiments with quantitative, highly sampled in vivo measurements from the Jak/STAT5 signaling pathway in fibroblasts, stimulated with erythropoietin (Epo). Afterwards we apply the presented iterative experimental design approach to the Jak/STAT3 pathway of primary hepatocytes stimulated with IL-6. Our approach offers the possibility of deciding which scientific questions can be answered based on existing laboratory constraints. To be able to concentrate on feasible questions on account of inexpensive computational simulations yields not only enormous cost and time saving, but also helps to specify realizable, systematic research projects in advance.

Computational processing and error reduction strategies for standardized quantitative data in biological networks.

High-quality quantitative data generated under standardized conditions is critical for understanding dynamic cellular processes. We report strategies for error reduction, and algorithms for automated data processing and for establishing the widely used techniques of immunoprecipitation and immunoblotting as highly precise methods for the quantification of protein levels and modifications. To determine the stoichiometry of cellular components and to ensure comparability of experiments, relative signals are converted to absolute values. A major source for errors in blotting techniques are inhomogeneities of the gel and the transfer procedure leading to correlated errors. These correlations are prevented by randomized gel loading, which significantly reduces standard deviations. Further error reduction is achieved by using housekeeping proteins as normalizers or by adding purified proteins in immunoprecipitations as calibrators in combination with criteria-based normalization. Additionally, we developed a computational tool for automated normalization, validation and integration of data derived from multiple immunoblots. In this way, large sets of quantitative data for dynamic pathway modeling can be generated, enabling the identification of systems properties and the prediction of targets for efficient intervention.