Lorenza A. D'Alessandro

Heterogeneous kinetics of AKT signaling in individual cells are accounted for by variable protein concentration

In most solid cancers, cells harboring oncogenic mutations represent only a sub-fraction of the entire population. Within this sub-fraction the expression level of mutated proteins can vary significantly due to cellular variability limiting the efficiency of targeted therapy. To address the causes of the heterogeneity, we performed a systematic analysis of one of the most frequently mutated pathways in cancer cells, the phosphatidylinositol 3 kinase (PI3K) signaling pathway. Among others PI3K signaling is activated by the hepatocyte growth factor (HGF) that regulates proliferation of hepatocytes during liver regeneration but also fosters tumor cell proliferation. HGF-mediated responses of PI3K signaling were monitored both at the single cell and cell population level in primary mouse hepatocytes and in the hepatoma cell line Hepa1_6. Interestingly, we observed that the HGF-mediated AKT responses at the level of individual cells is rather heterogeneous. However, the overall average behavior of the single cells strongly resembled the dynamics of AKT activation determined at the cell population level. To gain insights into the molecular cause for the observed heterogeneous behavior of individual cells, we employed dynamic mathematical modeling in a stochastic framework. Our analysis demonstrated that intrinsic noise was not sufficient to explain the observed kinetic behavior, but rather the importance of extrinsic noise has to be considered. Thus, distinct from gene expression in the examined signaling pathway fluctuations of the reaction rates has only a minor impact whereas variability in the concentration of the various signaling components even in a clonal cell population is a key determinant for the kinetic behavior.

Cellular ERK phospho-form profiles with conserved preference for a switch-like pattern.

ERK is a member of the MAPK pathway with essential functions in cell proliferation, differentiation, and survival. Complete ERK activation by the kinase MEK requires dual phosphorylation at T and Y within the activation motif TEY. We show that exposure of primary mouse hepatocytes to hepatocyte growth factor (HGF) results in phosphorylation at the activation motif, but not of other residues nearby. To determine the relative abundances of unphosphorylated ERK and the three ERK phospho-forms pT, pY, and pTpY, we employed an extended one-source peptide/phosphopeptide standard method in combination with nanoUPLC-MS. This method enabled us to determine the abundances of phospho-forms with a relative variability of ≤5% (SD). We observed a switch-like preference of ERK phospho-form abundances toward the active, doubly phosphorylated and the inactive, unphosphorylated form. Interestingly, ERK phospho-form profiles were similar upon growth factor and cytokine stimulation. A screening of several murine and human cell systems revealed that the balance between TY- and pTpY-ERK is conserved while the abundances of pT- and pY-ERK are more variable within cell types. We show that the phospho-form profiles do not change by blocking MEK activity suggesting that cellular phosphatases determine the ERK phospho-form distribution. This study provides novel quantitative insights into multisite phosphorylation.

Unraveling liver complexity from molecular to organ level: Challenges and perspectives

Biological responses are determined by information processing at multiple and highly interconnected scales. Within a tissue the individual cells respond to extracellular stimuli by regulating intracellular signaling pathways that in turn determine cell fate decisions and influence the behavior of neighboring cells. As a consequence the cellular responses critically impact tissue composition and architecture. Understanding the regulation of these mechanisms at different scales is key to unravel the emergent properties of biological systems. In this perspective, a multidisciplinary approach combining experimental data with mathematical modeling is introduced. We report the approach applied within the Virtual Liver Network to analyze processes that regulate liver functions from single cell responses to the organ level using a number of examples. By facilitating interdisciplinary collaborations, the Virtual Liver Network studies liver regeneration and inflammatory processes as well as liver metabolic functions at multiple scales, and thus provides a suitable example to identify challenges and point out potential future application of multi-scale systems biology.

Context-specific flow through the MEK/ERK module produces cell- and ligand-specific patterns of ERK single and double phosphorylation.

ERK phosphorylated on a single residue of the activation motif may be a sign of dysregulated proliferation. ERK phosphorylation patterns In the ERK pathway, the dual-specificity kinase MEK phosphorylates a threonine and a tyrosine residue in ERK, and this dual-phosphorylated form is the fully active kinase. Iwamoto et al. used mass spectrometry, quantitative Western blotting, and mathematical modeling to explore MEK-dependent phosphorylation dynamics of ERK in skin and liver cells exposed to either a cytokine, IL-6, or a growth factor, HGF. Not surprisingly, the different stimuli produced different dynamics of ERK phosphorylation, and skin and liver cells responded differently to the same ligand. The dynamics of the changes in the abundance of the phosphorylated forms of ERK (pT-ERK, pY-ERK, and pTpY-ERK) and the relative distributions of the single- and double-phosphorylated forms of ERK were different. Mathematical modeling indicated that distinct network structures with or without regulated feedback loops produced the different dynamics of ERK phosphorylation and distributions of phosphorylated ERK. This study provides biochemical insight into how a single pathway can produce distinct responses, such as differentiation or proliferation. The same pathway, such as the mitogen-activated protein kinase (MAPK) pathway, can produce different cellular responses, depending on stimulus or cell type. We examined the phosphorylation dynamics of the MAPK kinase MEK and its targets extracellular signal–regulated kinase 1 and 2 (ERK1/2) in primary hepatocytes and the transformed keratinocyte cell line HaCaT A5 exposed to either hepatocyte growth factor or interleukin-6. By combining quantitative mass spectrometry with dynamic modeling, we elucidated network structures for the reversible threonine and tyrosine phosphorylation of ERK in both cell types. In addition to differences in the phosphorylation and dephosphorylation reactions, the HaCaT network model required two feedback mechanisms, which, as the experimental data suggested, involved the induction of the dual-specificity phosphatase DUSP6 and the scaffold paxillin. We assayed and modeled the accumulation of the double-phosphorylated and active form of ERK1/2, as well as the dynamics of the changes in the monophosphorylated forms of ERK1/2. Modeling the differences in the dynamics of the changes in the distributions of the phosphorylated forms of ERK1/2 suggested that different amounts of MEK activity triggered context-specific responses, with primary hepatocytes favoring the formation of double-phosphorylated ERK1/2 and HaCaT A5 cells that produce both the threonine-phosphorylated and the double-phosphorylated form. These differences in phosphorylation distributions explained the threshold, sensitivity, and saturation of the ERK response. We extended the findings of differential ERK phosphorylation profiles to five additional cultured cell systems and matched liver tumor and normal tissue, which revealed context-specific patterns of the various forms of phosphorylated ERK.

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.

Hepatocellular carcinoma: a systems biology perspective.

Hepatocellular carcinomas (HCCs) have different etiology and heterogenic genomic alterations lead to high complexity. The molecular features of HCC have largely been studied by gene expression and proteome profiling focusing on the correlations between the expression of specific markers and clinical data. Integration of the increasing amounts of data in databases has facilitated the link of genomic and proteomic profiles of HCC to disease state and clinical outcome. Despite the current knowledge, specific molecular markers remain to be identified and new strategies are required to establish novel-targeted therapies. In the last years, mathematical models reconstructing gene and protein networks based on experimental data of HCC have been developed providing powerful tools to predict candidate interactions and potential targets for therapy. Furthermore, the combination of dynamic and logical mathematical models with quantitative data allows detailed mechanistic insights into system properties. To address effects at the organ level, mathematical models reconstructing the three-dimensional organization of liver lobules were developed. In the future, integration of different modeling approaches capturing the effects at the cellular up to the organ level is required to address the complex properties of HCC and to enable the discovery of new targets for HCC prevention or treatment.

A Systems Biology Study on NFκB Signaling in Primary Mouse Hepatocytes.

The cytokine tumor necrosis factor-alpha (TNFα) is one of the key factors during the priming phase of liver regeneration as well as in hepatocarcinogenesis. TNFα activates the nuclear factor κ-light-chain-enhancer of activated B cells (NFκB) signaling pathway and contributes to the conversion of quiescent hepatocytes to activated hepatocytes that are able to proliferate in response to growth factor stimulation. Different mathematical models have been previously established for TNFα/NFκB signaling in the context of tumor cells. Combining these mathematical models with time-resolved measurements of expression and phosphorylation of TNFα/NFκB pathway constituents in primary mouse hepatocytes revealed that an additional phosphorylation step of the NFκB isoform p65 has to be considered in the mathematical model in order to sufficiently describe the dynamics of pathway activation in the primary cells. Also, we addressed the role of basal protein turnover by experimentally measuring the degradation rate of pivotal players in the absence of TNFα and including this information in the model. To elucidate the impact of variations in the protein degradation rates on TNFα/NFκB signaling on the overall dynamic behavior we used global sensitivity analysis that accounts for parameter uncertainties and showed that degradation and translation of p65 had a major impact on the amplitude and the integral of p65 phosphorylation. Finally, our mathematical model of TNFα/NFκB signaling was able to predict the time-course of the complex formation of p65 and of the inhibitor of NFκB (IκB) in primary mouse hepatocytes, which was experimentally verified. Hence, we here present a mathematical model for TNFα/NFκB signaling in primary mouse hepatocytes that provides an important basis to quantitatively disentangle the complex interplay of multiple factors in liver regeneration and tumorigenesis.

Excemplify: a flexible template based solution, parsing and managing data in spreadsheets for experimentalists.

In systems biology, quantitative experimental data is the basis of building mathematical models. In most of the cases, they are stored in Excel files and hosted locally. To have a public database for collecting, retrieving and citing experimental raw data as well as experimental conditions is important for both experimentalists and modelers. However, the great effort needed in the data handling procedure and in the data submission procedure becomes the crucial limitation for experimentalists to contribute to a database, thereby impeding the database to deliver its benefit. Moreover, manual copy and paste operations which are commonly used in those procedures increase the chance of making mistakes. Excemplify, a web-based application, proposes a flexible and adaptable template-based solution to solve these problems. Comparing to the normal template based uploading approach, which is supported by some public databases, rather than predefining a format that is potentiall impractical, Excemplify allows users to create their own experiment-specific content templates in different experiment stages and to build corresponding knowledge bases for parsing. Utilizing the embedded knowledge of used templates, Excemplify is able to parse experimental data from the initial setup stage and generate following stages spreadsheets automatically. The proposed solution standardizes the flows of data traveling according to the standard procedures of applying the experiment, cuts down the amount of manual effort and reduces the chance of mistakes caused by manual data handling. In addition, it maintains the context of meta-data from the initial preparation manuscript and improves the data consistency. It interoperates and complements RightField and SEEK as well.

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.

P117 HOLISTIC APPROACH TO UNRAVEL FUNCTIONS AND REGULATION OF HGF IN LIVER REGENERATION

P116 INTERACTION OF APOPTOTIC AND CYTOPROTECTIVE BILE ACIDS WITH BIOMEMBRANES T. Sousa, S. Pinto, S.D. Lucas, R. Moreira, A. Coutinho, R.E. Castro, C.M. Rodrigues, M. Prieto, F. Fernandes. Centro de Quimica-Fisica Molecular and Institute of Nanoscience and Nanotechnology, Instituto Superior Tecnico, Instituto de Investigacao do Medicamento (iMed.ULisboa), Faculty of Pharmacy, University of Lisbon, Lisbon, Portugal E-mail: ruieduardocastro@ff.ul.pt