Predicting genetic interactions from Boolean models of biological networks

Patient-specific Boolean models of signalling networks guide personalised treatments

Prostate cancer is the second most occurring cancer in men worldwide. To better understand the mechanisms of tumorigenesis and possible treatment responses, we developed a mathematical model of prostate cancer which considers the major signalling pathways known to be deregulated. We personalised this Boolean model to molecular data to reflect the heterogeneity and specific response to perturbations of cancer patients. A total of 488 prostate samples were used to build patient-specific models and compared to available clinical data. Additionally, eight prostate cell line-specific models were built to validate our approach with dose-response data of several drugs. The effects of single and combined drugs were tested in these models under different growth conditions. We identified 15 actionable points of interventions in one cell line-specific model whose inactivation hinders tumorigenesis. To validate these results, we tested nine small molecule inhibitors of five of those putative targets and found a dose-dependent effect on four of them, notably those targeting HSP90 and PI3K. These results highlight the predictive power of our personalised Boolean models and illustrate how they can be used for precision oncology.

Unifying the mechanism of mitotic exit control in a spatiotemporal logical model

The transition from mitosis into the first gap phase of the cell cycle in budding yeast is controlled by the Mitotic Exit Network (MEN). The network interprets spatiotemporal cues about the progression of mitosis and ensures that release of Cdc14 phosphatase occurs only after completion of key mitotic events. The MEN has been studied intensively; however, a unified understanding of how localisation and protein activity function together as a system is lacking. In this paper, we present a compartmental, logical model of the MEN that is capable of representing spatial aspects of regulation in parallel to control of enzymatic activity. We show that our model is capable of correctly predicting the phenotype of the majority of mutants we tested, including mutants that cause proteins to mislocalise. We use a continuous time implementation of the model to demonstrate that Cdc14 Early Anaphase Release (FEAR) ensures robust timing of anaphase, and we verify our findings in living cells. Furthermore, we show that our model can represent measured cell-cell variation in Spindle Position Checkpoint (SPoC) mutants. This work suggests a general approach to incorporate spatial effects into logical models. We anticipate that the model itself will be an important resource to experimental researchers, providing a rigorous platform to test hypotheses about regulation of mitotic exit.

Conceptual and computational framework for logical modelling of biological networks deregulated in diseases

Mathematical models can serve as a tool to formalize biological knowledge from diverse sources, to investigate biological questions in a formal way, to test experimental hypotheses, to predict the effect of perturbations and to identify underlying mechanisms. We present a pipeline of computational tools that performs a series of analyses to explore a logical model's properties. A logical model of initiation of the metastatic process in cancer is used as a transversal example. We start by analysing the structure of the interaction network constructed from the literature or existing databases. Next, we show how to translate this network into a mathematical object, specifically a logical model, and how robustness analyses can be applied to it. We explore the visualization of the stable states, defined as specific attractors of the model, and match them to cellular fates or biological read-outs. With the different tools we present here, we explain how to assign to each solution of the model a probability and how to identify genetic interactions using mutant phenotype probabilities. Finally, we connect the model to relevant experimental data: we present how some data analyses can direct the construction of the network, and how the solutions of a mathematical model can also be compared with experimental data, with a particular focus on high-throughput data in cancer biology. A step-by-step tutorial is provided as a Supplementary Material and all models, tools and scripts are provided on an accompanying website: https://github.com/sysbio-curie/Logical_modelling_pipeline.

RMut: R package for a Boolean sensitivity analysis against various types of mutations

There have been many in silico studies based on a Boolean network model to investigate network sensitivity against gene or interaction mutations. However, there are no proper tools to examine the network sensitivity against many different types of mutations, including user-defined ones. To address this issue, we developed RMut, which is an R package to analyze the Boolean network-based sensitivity by efficiently employing not only many well-known node-based and edgetic mutations but also novel user-defined mutations. In addition, RMut can specify the mutation area and the duration time for more precise analysis. RMut can be used to analyze large-scale networks because it is implemented in a parallel algorithm using the OpenCL library. In the first case study, we observed that the real biological networks were most sensitive to overexpression/state-flip and edge-addition/-reverse mutations among node-based and edgetic mutations, respectively. In the second case study, we showed that edgetic mutations can predict drug-targets better than node-based mutations. Finally, we examined the network sensitivity to double edge-removal mutations and found an interesting synergistic effect. Taken together, these findings indicate that RMut is a flexible R package to efficiently analyze network sensitivity to various types of mutations. RMut is available at https://github.com/csclab/RMut.

Next generation of network medicine: interdisciplinary signaling approaches

In the last decade, network approaches have transformed our understanding of biological systems. Network analyses and visualizations have allowed us to identify essential molecules and modules in biological systems, and improved our understanding of how changes in cellular processes can lead to complex diseases, such as cancer, infectious and neurodegenerative diseases. "Network medicine" involves unbiased large-scale network-based analyses of diverse data describing interactions between genes, diseases, phenotypes, drug targets, drug transport, drug side-effects, disease trajectories and more. In terms of drug discovery, network medicine exploits our understanding of the network connectivity and signaling system dynamics to help identify optimal, often novel, drug targets. Contrary to initial expectations, however, network approaches have not yet delivered a revolution in molecular medicine. In this review, we propose that a key reason for the limited impact, so far, of network medicine is a lack of quantitative multi-disciplinary studies involving scientists from different backgrounds. To support this argument, we present existing approaches from structural biology, 'omics' technologies (e.g., genomics, proteomics, lipidomics) and computational modeling that point towards how multi-disciplinary efforts allow for important new insights. We also highlight some breakthrough studies as examples of the potential of these approaches, and suggest ways to make greater use of the power of interdisciplinarity. This review reflects discussions held at an interdisciplinary signaling workshop which facilitated knowledge exchange from experts from several different fields, including in silico modelers, computational biologists, biochemists, geneticists, molecular and cell biologists as well as cancer biologists and pharmacologists.

A new ETV6-NTRK3 cell line model reveals MALAT1 as a novel therapeutic target - a short report

BACKGROUND

Previously, the chromosomal translocation t(12;15)(p13;q25) has been found to recurrently occur in both solid tumors and leukemias. This translocation leads to ETV6-NTRK3 (EN) gene fusions resulting in ectopic expression of the NTRK3 neurotropic tyrosine receptor kinase moiety as well as oligomerization through the donated ETV6-sterile alpha motif domain. As yet, no in vitro cell line model carrying this anomaly is available. Here we genetically characterized the acute promyelocytic leukemia (APL) cell line AP-1060 and, by doing so, revealed the presence of a t(12;15)(p13;q25). Subsequently, we evaluated its suitability as a model for this important clinical entity.

METHODS

Spectral karyotyping, fluorescence in situ hybridization (FISH), and genomic and transcriptomic microarray-based profiling were used to screen for the presence of EN fusions. qRT-PCR was used for quantitative expression analyses. Responses to AZ-23 (NTRK) and wortmannin (PI3K) inhibitors, as well as to arsenic trioxide (ATO), were assessed using colorimetric assays. An AZ-23 microarray screen was used to define the EN targetome, which was parsed bioinformatically. MAPK1 and MALAT1 activation were assayed using Western blotting and RNA-FISH, respectively, whereas an AML patient cohort was used to assess the clinical occurrence of MALAT1 activation.

RESULTS

An EN fusion was detected in AP1060 cells which, accordingly, turned out to be hypersensitive to AZ-23. We also found that AZ-23 can potentiate the effect of ATO and inhibit the phosphorylation of its canonical target MAPK1. The AZ-23 microarray screen highlighted a novel EN target, MALAT1, which also proved sensitive to wortmannin. Finally, we found that MALAT1 was massively up-regulated in a subset of AML patients.

CONCLUSIONS

From our data we conclude that AP-1060 may serve as a first publicly available preclinical model for EN. In addition, we conclude that these EN-positive cells are sensitive to the NTRK inhibitor AZ-23 and that this inhibitor may potentiate the therapeutic efficacy of ATO. Our data also highlight a novel AML EN target, MALAT1, which was so far only conspicuous in solid tumors.

Logic Modeling in Quantitative Systems Pharmacology

Here we present logic modeling as an approach to understand deregulation of signal transduction in disease and to characterize a drug's mode of action. We discuss how to build a logic model from the literature and experimental data and how to analyze the resulting model to obtain insights of relevance for systems pharmacology. Our workflow uses the free tools OmniPath (network reconstruction from the literature), CellNOpt (model fit to experimental data), MaBoSS (model analysis), and Cytoscape (visualization).

Identification of bifurcation transitions in biological regulatory networks using Answer-Set Programming

Background

Numerous cellular differentiation processes can be captured using discrete qualitative models of biological regulatory networks. These models describe the temporal evolution of the state of the network subject to different competing transitions, potentially leading the system to different attractors. This paper focusses on the formal identification of states and transitions that are crucial for preserving or pre-empting the reachability of a given behaviour.

Methods

In the context of non-deterministic automata networks, we propose a static identification of so-called bifurcations, i.e., transitions after which a given goal is no longer reachable. Such transitions are naturally good candidates for controlling the occurrence of the goal, notably by modulating their propensity. Our method combines Answer-Set Programming with static analysis of reachability properties to provide an under-approximation of all the existing bifurcations.

Results

We illustrate our discrete bifurcation analysis on several models of biological systems, for which we identify transitions which impact the reachability of given long-term behaviour. In particular, we apply our implementation on a regulatory network among hundreds of biological species, supporting the scalability of our approach.

Conclusions

Our method allows a formal and scalable identification of transitions which are responsible for the lost of capability to reach a given state. It can be applied to any asynchronous automata networks, which encompass Boolean and multi-valued models. An implementation is provided as part of the Pint software, available at http://loicpauleve.name/pint.

Electronic supplementary material

The online version of this article (doi:10.1186/s13015-017-0110-3) contains supplementary material, which is available to authorized users.

The second European interdisciplinary Ewing sarcoma research summit--A joint effort to deconstructing the multiple layers of a complex disease

Despite multimodal treatment, long term outcome for patients with Ewing sarcoma is still poor. The second “European interdisciplinary Ewing sarcoma research summit” assembled a large group of scientific experts in the field to discuss their latest unpublished findings on the way to the identification of novel therapeutic targets and strategies. Ewing sarcoma is characterized by a quiet genome with presence of an EWSR1-ETS gene rearrangement as the only and defining genetic aberration. RNA-sequencing of recently described Ewing-like sarcomas with variant translocations identified them as biologically distinct diseases. Various presentations adressed mechanisms of EWS-ETS fusion protein activities with a focus on EWS-FLI1. Data were presented shedding light on the molecular underpinnings of genetic permissiveness to this disease uncovering interaction of EWS-FLI1 with recently discovered susceptibility loci. Epigenetic context as a consequence of the interaction between the oncoprotein, cell type, developmental stage, and tissue microenvironment emerged as dominant theme in the discussion of the molecular pathogenesis and inter- and intra-tumor heterogeneity of Ewing sarcoma, and the difficulty to generate animal models faithfully recapitulating the human disease. The problem of preclinical development of biologically targeted therapeutics was discussed and promising perspectives were offered from the study of novel in vitro models. Finally, it was concluded that in order to facilitate rapid pre-clinical and clinical development of novel therapies in Ewing sarcoma, the community needs a platform to maintain knowledge of unpublished results, systems and models used in drug testing and to continue the open dialogue initiated at the first two Ewing sarcoma summits.

Prediction of Genetic Interactions Using Machine Learning and Network Properties

A genetic interaction (GI) is a type of interaction where the effect of one gene is modified by the effect of one or several other genes. These interactions are important for delineating functional relationships among genes and their corresponding proteins, as well as elucidating complex biological processes and diseases. An important type of GI - synthetic sickness or synthetic lethality - involves two or more genes, where the loss of either gene alone has little impact on cell viability, but the combined loss of all genes leads to a severe decrease in fitness (sickness) or cell death (lethality). The identification of GIs is an important problem for it can help delineate pathways, protein complexes, and regulatory dependencies. Synthetic lethal interactions have important clinical and biological significance, such as providing therapeutically exploitable weaknesses in tumors. While near systematic high-content screening for GIs is possible in single cell organisms such as yeast, the systematic discovery of GIs is extremely difficult in mammalian cells. Therefore, there is a great need for computational approaches to reliably predict GIs, including synthetic lethal interactions, in these organisms. Here, we review the state-of-the-art approaches, strategies, and rigorous evaluation methods for learning and predicting GIs, both under general (healthy/standard laboratory) conditions and under specific contexts, such as diseases.

DeDaL: Cytoscape 3 app for producing and morphing data-driven and structure-driven network layouts

Background

Visualization and analysis of molecular profiling data together with biological networks are able to provide new mechanistic insights into biological functions. Currently, it is possible to visualize high-throughput data on top of pre-defined network layouts, but they are not always adapted to a given data analysis task. A network layout based simultaneously on the network structure and the associated multidimensional data might be advantageous for data visualization and analysis in some cases.

Results

We developed a Cytoscape app, which allows constructing biological network layouts based on the data from molecular profiles imported as values of node attributes. DeDaL is a Cytoscape 3 app, which uses linear and non-linear algorithms of dimension reduction to produce data-driven network layouts based on multidimensional data (typically gene expression). DeDaL implements several data pre-processing and layout post-processing steps such as continuous morphing between two arbitrary network layouts and aligning one network layout with respect to another one by rotating and mirroring. The combination of all these functionalities facilitates the creation of insightful network layouts representing both structural network features and correlation patterns in multivariate data. We demonstrate the added value of applying DeDaL in several practical applications, including an example of a large protein-protein interaction network.

Conclusions

DeDaL is a convenient tool for applying data dimensionality reduction methods and for designing insightful data displays based on data-driven layouts of biological networks, built within Cytoscape environment. DeDaL is freely available for downloading at http://bioinfo-out.curie.fr/projects/dedal/.