Anja Sieber

Strong negative feedback from Erk to Raf confers robustness to MAPK signalling

Protein levels within signal transduction pathways vary strongly from cell to cell. Here, we analysed how signalling pathways can still process information quantitatively despite strong heterogeneity in protein levels. We systematically perturbed the protein levels of Erk, the terminal kinase in the MAPK signalling pathway in a panel of human cell lines. We found that the steady‐state phosphorylation of Erk is very robust against perturbations of Erk protein level. Although a multitude of mechanisms exist that may provide robustness against fluctuating protein levels, we found that one single feedback from Erk to Raf‐1 accounts for the observed robustness. Surprisingly, robustness is provided through a fast post‐translational mechanism although variation of Erk levels occurs on a timescale of days.

Network quantification of EGFR signaling unveils potential for targeted combination therapy

The epidermal growth factor receptor (EGFR) signaling network is activated in most solid tumors, and small‐molecule drugs targeting this network are increasingly available. However, often only specific combinations of inhibitors are effective. Therefore, the prediction of potent combinatorial treatments is a major challenge in targeted cancer therapy. In this study, we demonstrate how a model‐based evaluation of signaling data can assist in finding the most suitable treatment combination. We generated a perturbation data set by monitoring the response of RAS/PI3K signaling to combined stimulations and inhibitions in a panel of colorectal cancer cell lines, which we analyzed using mathematical models. We detected that a negative feedback involving EGFR mediates strong cross talk from ERK to AKT. Consequently, when inhibiting MAPK, AKT activity is increased in an EGFR‐dependent manner. Using the model, we predict that in contrast to single inhibition, combined inactivation of MEK and EGFR could inactivate both endpoints of RAS, ERK and AKT. We further could demonstrate that this combination blocked cell growth in BRAF‐ as well as KRAS‐mutated tumor cells, which we confirmed using a xenograft model.

The Two Active X Chromosomes in Female ESCs Block Exit from the Pluripotent State by Modulating the ESC Signaling Network

During early development of female mouse embryos, both X chromosomes are transiently active. X gene dosage is then equalized between the sexes through the process of X chromosome inactivation (XCI). Whether the double dose of X-linked genes in females compared with males leads to sex-specific developmental differences has remained unclear. Using embryonic stem cells with distinct sex chromosome compositions as a model system, we show that two X chromosomes stabilize the naive pluripotent state by inhibiting MAPK and Gsk3 signaling and stimulating the Akt pathway. Since MAPK signaling is required to exit the pluripotent state, differentiation is paused in female cells as long as both X chromosomes are active. By preventing XCI or triggering it precociously, we demonstrate that this differentiation block is released once XX cells have undergone X inactivation. We propose that double X dosage interferes with differentiation, thus ensuring a tight coupling between X chromosome dosage compensation and development.

Drug resistance mechanisms in colorectal cancer dissected with cell type-specific dynamic logic models

Genomic features are used as biomarkers of sensitivity to kinase inhibitors used widely to treat human cancer, but effective patient stratification based on these principles remains limited in impact. Insofar as kinase inhibitors interfere with signaling dynamics, and, in turn, signaling dynamics affects inhibitor responses, we investigated associations in this study between cell-specific dynamic signaling pathways and drug sensitivity. Specifically, we measured 14 phosphoproteins under 43 different perturbed conditions (combinations of 5 stimuli and 7 inhibitors) in 14 colorectal cancer cell lines, building cell line-specific dynamic logic models of underlying signaling networks. Model parameters representing pathway dynamics were used as features to predict sensitivity to a panel of 27 drugs. Specific parameters of signaling dynamics correlated strongly with drug sensitivity for 14 of the drugs, 9 of which had no genomic biomarker. Following one of these associations, we validated a drug combination predicted to overcome resistance to MEK inhibitors by coblockade of GSK3, which was not found based on associations with genomic data. These results suggest that to better understand the cancer resistance and move toward personalized medicine, it is essential to consider signaling network dynamics that cannot be inferred from static genotypes. Cancer Res; 77(12); 3364-75. ©2017 AACR.

Perturbation-response genes reveal signaling footprints in cancer gene expression

Aberrant cell signaling can cause cancer and other diseases and is a focal point of drug research. A common approach is to infer signaling activity of pathways from gene expression. However, mapping gene expression to pathway components disregards the effect of post-translational modifications, and downstream signatures represent very specific experimental conditions. Here we present PROGENy, a method that overcomes both limitations by leveraging a large compendium of publicly available perturbation experiments to yield a common core of Pathway RespOnsive GENes. Unlike pathway mapping methods, PROGENy can (i) recover the effect of known driver mutations, (ii) provide or improve strong markers for drug indications, and (iii) distinguish between oncogenic and tumor suppressor pathways for patient survival. Collectively, these results show that PROGENy accurately infers pathway activity from gene expression in a wide range of conditions.Deregulation of signalling is responsible for many cancer phenotypes. Leveraging available perturbation data, here the authors assess large-scale pathway activity patterns based on consensus downstream readout genes, enabling accurate prediction of the effects of mutations and small molecules.

Analysis of impedance-based cellular growth assays

MOTIVATION Impedance-based technologies are advancing methods for measuring proliferation of adherent cell cultures non-invasively and in real time. The analysis of the resulting data has so far been hampered by inappropriate computational methods and the lack of systematic data to evaluate the characteristics of the assay. RESULTS We used a commercially available system for impedance-based growth measurement (xCELLigence) and compared the reported cell index with data from microscopy. We found that the measured signal correlates linearly with the cell number throughout the time of an experiment with sufficient accuracy in subconfluent cell cultures. The resulting growth curves for various colon cancer cells could be well described with the empirical Richards growth model, which allows for extracting quantitative parameters (such as characteristic cycle times). We found that frequently used readouts like the cell index at a specific time or the area under the growth curve cannot be used to faithfully characterize growth inhibition. We propose to calculate the average growth rate of selected time intervals to accurately estimate time-dependent IC50 values of drugs from growth curves. CONTACT nils.bluethgen@charite.de SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

An immediate–late gene expression module decodes ERK signal duration

The RAF‐MEK‐ERK signalling pathway controls fundamental, often opposing cellular processes such as proliferation and apoptosis. Signal duration has been identified to play a decisive role in these cell fate decisions. However, it remains unclear how the different early and late responding gene expression modules can discriminate short and long signals. We obtained both protein phosphorylation and gene expression time course data from HEK293 cells carrying an inducible construct of the proto‐oncogene RAF. By mathematical modelling, we identified a new gene expression module of immediate–late genes (ILGs) distinct in gene expression dynamics and function. We find that mRNA longevity enables these ILGs to respond late and thus translate ERK signal duration into response amplitude. Despite their late response, their GC‐rich promoter structure suggested and metabolic labelling with 4SU confirmed that transcription of ILGs is induced immediately. A comparative analysis shows that the principle of duration decoding is conserved in PC12 cells and MCF7 cells, two paradigm cell systems for ERK signal duration. Altogether, our findings suggest that ILGs function as a gene expression module to decode ERK signal duration.

Modelling signalling networks from perturbation data

Motivation Intracellular signalling is realized by complex signalling networks, which are almost impossible to understand without network models, especially if feedbacks are involved. Modular Response Analysis (MRA) is a convenient modelling method to study signalling networks in various contexts. Results We developed the software package STASNet (STeady-STate Analysis of Signalling Networks) that provides an augmented and extended version of MRA suited to model signalling networks from incomplete perturbation schemes and multi-perturbation data. Using data from the Dialogue on Reverse Engineering Assessment and Methods challenge, we show that predictions from STASNet models are among the top-performing methods. We applied the method to study the effect of SHP2, a protein that has been implicated in resistance to targeted therapy in colon cancer, using a novel dataset from the colon cancer cell line Widr and a SHP2-depleted derivative. We find that SHP2 is required for mitogen-activated protein kinase signalling, whereas AKT signalling only partially depends on SHP2. Availability and implementation An R-package is available at https://github.com/molsysbio/STASNet. Supplementary information Supplementary data are available at Bioinformatics online.

Comparative Network Reconstruction using mixed integer programming

Motivation Signal‐transduction networks are often aberrated in cancer cells, and new anti‐cancer drugs that specifically target oncogenes involved in signaling show great clinical promise. However, the effectiveness of such targeted treatments is often hampered by innate or acquired resistance due to feedbacks, crosstalks or network adaptations in response to drug treatment. A quantitative understanding of these signaling networks and how they differ between cells with different oncogenic mutations or between sensitive and resistant cells can help in addressing this problem. Results Here, we present Comparative Network Reconstruction (CNR), a computational method to reconstruct signaling networks based on possibly incomplete perturbation data, and to identify which edges differ quantitatively between two or more signaling networks. Prior knowledge about network topology is not required but can straightforwardly be incorporated. We extensively tested our approach using simulated data and applied it to perturbation data from a BRAF mutant, PTPN11 KO cell line that developed resistance to BRAF inhibition. Comparing the reconstructed networks of sensitive and resistant cells suggests that the resistance mechanism involves re‐establishing wild‐type MAPK signaling, possibly through an alternative RAF‐isoform. Availability and implementation CNR is available as a python module at https://github.com/NKI‐CCB/cnr. Additionally, code to reproduce all figures is available at https://github.com/NKI‐CCB/CNR‐analyses. Supplementary information Supplementary data are available at Bioinformatics online.

Dissecting cancer resistance to therapies with cell-type-specific dynamic logic models

Therapies targeting specific molecular processes, in particular kinases, are major strategies to treat cancer. Genomic features are commonly used as biomarkers for drug sensitivity, but our ability to stratify patients based on these features is still limited. As response to kinase inhibitors is a dynamic process affecting largely signal transduction, we investigated the association between cell-specific dynamic signaling pathways and drug sensitivity. We measured 14 phosphoproteins under 43 different perturbed conditions (combination of 5 stimuli and 7 inhibitors) for 14 colorectal cancer cell-lines, and built cell-line-specific dynamic logic models of the underlying signaling network. Model parameters, representing pathway dynamics, were used as features to predict sensitivity to a panel of 27 drugs. This analysis revealed associations between cell-specific signaling pathways and drug sensitivity for 14 of the drugs, 9 of which have no genomic biomarker. Following one of these associations, we validated a drug combination predicted to overcome resistance to MEK inhibitors by co-blockade of GSK3. These results underscore the value of perturbation-based studies to find biomarkers and combination therapies complementing those based on a static genomic characterization.