A Modeling Approach to Explain Mutually Exclusive and Co-Occurring Genetic Alterations in Bladder Tumorigenesis

Super multiple primary lung cancers harbor high-frequency BRAF and low-frequency EGFR mutations in the MAPK pathway

The incidence of multiple primary lung cancer (MPLC) is increasing, with some of our surgical patients exhibiting numerous lesions. We defined lung cancer with five or more primary lesions as super MPLCs. Elucidating the genomic characteristics of this special MPLC subtype can help reduce disease burden and understand tumor evolution. In our cohort of synchronous super early-stage MPLCs (PUMCH-ssesMPLC), whole-exome sequencing on 130 resected malignant specimens from 18 patients provided comprehensive super-MPLC genomic landscapes. Mutations are enriched in PI3k-Akt and MAPK pathways. Their BRAF mutation frequency (31.5%) is significantly higher than MPLC with fewer lesions and early-stage single-lesion cancer, while EGFR mutations are significantly fewer (13.8%). As lesion counts increase, BRAF mutations gradually become dominant. Also, invasive lesions more tend to have classic super-MPLC mutation patterns. High-frequency BRAF mutations, especially Class II, and low-frequency EGFR mutations could be a reason for the limited effectiveness of targeted therapy in super-MPLC patients.

PhysiBoSS 2.0: a sustainable integration of stochastic Boolean and agent-based modelling frameworks

In systems biology, mathematical models and simulations play a crucial role in understanding complex biological systems. Different modelling frameworks are employed depending on the nature and scales of the system under study. For instance, signalling and regulatory networks can be simulated using Boolean modelling, whereas multicellular systems can be studied using agent-based modelling. Herein, we present PhysiBoSS 2.0, a hybrid agent-based modelling framework that allows simulating signalling and regulatory networks within individual cell agents. PhysiBoSS 2.0 is a redesign and reimplementation of PhysiBoSS 1.0 and was conceived as an add-on that expands the PhysiCell functionalities by enabling the simulation of intracellular cell signalling using MaBoSS while keeping a decoupled, maintainable and model-agnostic design. PhysiBoSS 2.0 also expands the set of functionalities offered to the users, including custom models and cell specifications, mechanistic submodels of substrate internalisation and detailed control over simulation parameters. Together with PhysiBoSS 2.0, we introduce PCTK, a Python package developed for handling and processing simulation outputs, and generating summary plots and 3D renders. PhysiBoSS 2.0 allows studying the interplay between the microenvironment, the signalling pathways that control cellular processes and population dynamics, suitable for modelling cancer. We show different approaches for integrating Boolean networks into multi-scale simulations using strategies to study the drug effects and synergies in models of cancer cell lines and validate them using experimental data. PhysiBoSS 2.0 is open-source and publicly available on GitHub with several repositories of accompanying interoperable tools.

A computational model of the DNA damage-induced IKK/ NF-κB pathway reveals a critical dependence on irradiation dose and PARP-1

The activation of IKK/NF-κB by genotoxic stress is a crucial process in the DNA damage response. Due to the anti-apoptotic impact of NF-κB, it can affect cell-fate decisions upon DNA damage and therefore interfere with tumor therapy-induced cell death. Here, we developed a dynamical model describing IKK/NF-κB signaling that faithfully reproduces quantitative time course data and enables a detailed analysis of pathway regulation. The approach elucidates a pathway topology with two hubs, where the first integrates signals from two DNA damage sensors and the second forms a coherent feedforward loop. The analyses reveal a critical role of the sensor protein PARP-1 in the pathway regulation. Introducing a method for calculating the impact of changes in individual components on pathway activity in a time-resolved manner, we show how irradiation dose influences pathway activation. Our results give a mechanistic understanding relevant for the interpretation of experimental and clinical studies.

Reducing Boolean networks with backward equivalence

BACKGROUND

Boolean Networks (BNs) are a popular dynamical model in biology where the state of each component is represented by a variable taking binary values that express, for instance, activation/deactivation or high/low concentrations. Unfortunately, these models suffer from the state space explosion, i.e., there are exponentially many states in the number of BN variables, which hampers their analysis.

RESULTS

We present Boolean Backward Equivalence (BBE), a novel reduction technique for BNs which collapses system variables that, if initialized with same value, maintain matching values in all states. A large-scale validation on 86 models from two online model repositories reveals that BBE is effective, since it is able to reduce more than 90% of the models. Furthermore, on such models we also show that BBE brings notable analysis speed-ups, both in terms of state space generation and steady-state analysis. In several cases, BBE allowed the analysis of models that were originally intractable due to the complexity. On two selected case studies, we show how one can tune the reduction power of BBE using model-specific information to preserve all dynamics of interest, and selectively exclude behavior that does not have biological relevance.

CONCLUSIONS

BBE complements existing reduction methods, preserving properties that other reduction methods fail to reproduce, and vice versa. BBE drops all and only the dynamics, including attractors, originating from states where BBE-equivalent variables have been initialized with different activation values The remaining part of the dynamics is preserved exactly, including the length of the preserved attractors, and their reachability from given initial conditions, without adding any spurious behaviours. Given that BBE is a model-to-model reduction technique, it can be combined with further reduction methods for BNs.

Target Control of Asynchronous Boolean Networks

We study the target control of asynchronous Boolean networks, to identify interventions that can drive the dynamics of a given Boolean network from any initial state to the desired target attractor. Based on the application time, the control can be realised with three types of perturbations, including instantaneous, temporary and permanent perturbations. We develop efficient methods to compute the target control for a given target attractor with these three types of perturbations. We compare our methods with the stable motif-based control method on a variety of real-life biological networks to evaluate their performance. We show that our methods scale well for large Boolean networks and they are able to identify a rich set of solutions with a small number of perturbations.

Cyclic attractors of nonexpanding q-ary networks

Discrete dynamical systems in which model components take on categorical values have been successfully applied to biological networks to study their global dynamic behavior. Boolean models in particular have been used extensively. However, multi-state models have also emerged as effective computational tools for the analysis of complex mechanisms underlying biological networks. Models in which variables assume more than two discrete states provide greater resolution, but this scheme introduces discontinuities. In particular, variables can increase or decrease by more than one unit in one time step. This can be corrected, without changing fixed points of the system, by applying an additional rule to each local activation function. On the other hand, if one is interested in cyclic attractors of their system, then this rule can potentially introduce new cyclic attractors that were not observed previously. This article makes some advancements in understanding the state space dynamics of multi-state network models with synchronous, sequential, or block-sequential update schedules and establishes conditions under which no new cyclic attractors are added to networks when the additional rule is applied. Our analytical results have the potential to be incorporated into modeling software and aid researchers in their analyses of biological multi-state networks.

Single Cell Transcriptomics to Understand HSC Heterogeneity and Its Evolution upon Aging

Single-cell transcriptomic technologies enable the uncovering and characterization of cellular heterogeneity and pave the way for studies aiming at understanding the origin and consequences of it. The hematopoietic system is in essence a very well adapted model system to benefit from this technological advance because it is characterized by different cellular states. Each cellular state, and its interconnection, may be defined by a specific location in the global transcriptional landscape sustained by a complex regulatory network. This transcriptomic signature is not fixed and evolved over time to give rise to less efficient hematopoietic stem cells (HSC), leading to a well-documented hematopoietic aging. Here, we review the advance of single-cell transcriptomic approaches for the understanding of HSC heterogeneity to grasp HSC deregulations upon aging. We also discuss the new bioinformatics tools developed for the analysis of the resulting large and complex datasets. Finally, since hematopoiesis is driven by fine-tuned and complex networks that must be interconnected to each other, we highlight how mathematical modeling is beneficial for doing such interconnection between multilayered information and to predict how HSC behave while aging.

SNAIL driven by a feed forward loop motif promotes TGFβinduced epithelial to mesenchymal transition

Epithelial to Mesenchymal Transition (EMT) plays an important role in tissue regeneration, embryonic development, and cancer metastasis. Several signaling pathways are known to regulate EMT, among which the modulation of TGFβ(Transforming Growth Factor-β) induced EMT is crucial in several cancer types. Several mathematical models were built to explore the role of core regulatory circuit of ZEB/miR-200, SNAIL/miR-34 double negative feedback loops in modulating TGFβinduced EMT. Different emergent behavior including tristability, irreversible switching, existence of hybrid EMT states were inferred though these models. Some studies have explored the role of TGFβreceptor activation, SMADs nucleocytoplasmic shuttling and complex formation. Recent experiments have revealed that MDM2 along with SMAD complex regulates SNAIL expression driven EMT. Encouraged by this, in the present study we developed a mathematical model for p53/MDM2 dependent TGFβinduced EMT regulation. Inclusion of p53 brings in an additional mechanistic perspective in exploring the EM transition. The network formulated comprises a C1FFL moderating SNAIL expression involving MDM2 and SMAD complex, which functions as a noise filter and persistent detector. The C1FFL was also observed to operate as a coincidence detector driving the SNAIL dependent downstream signaling into phenotypic switching decision. Systems modelling and analysis of the devised network, displayed interesting dynamic behavior, systems response to various inputs stimulus, providing a better understanding of p53/MDM2 dependent TGF-βinduced Epithelial to Mesenchymal Transition.

SteadyCellPhenotype: a web-based tool for the modeling of biological networks with ternary logic

SUMMARY

We introduce SteadyCellPhenotype, a browser-based interface for the analysis of ternary biological networks. It includes tools for deterministically finding all steady states of a network, as well as the simulation and visualization of trajectories with publication quality graphics. Simulations allow us to approximate the size of the basin for attractors and deterministic simulations of trajectories nearby specified points allow us to explore the behavior of the system in that neighborhood.

AVAILABILITY AND IMPLEMENTATION

https://github.com/knappa/steadycellphenotype MIT License.

UPMaBoSS: A Novel Framework for Dynamic Cell Population Modeling

Mathematical modeling aims at understanding the effects of biological perturbations, suggesting ways to intervene and to reestablish proper cell functioning in diseases such as cancer or in autoimmune disorders. This is a difficult task for obvious reasons: the level of details needed to describe the intra-cellular processes involved, the numerous interactions between cells and cell types, and the complex dynamical properties of such populations where cells die, divide and interact constantly, to cite a few. Another important difficulty comes from the spatial distribution of these cells, their diffusion and motility. All of these aspects cannot be easily resolved in a unique mathematical model or with a unique formalism. To cope with some of these issues, we introduce here a novel framework, UPMaBoSS (for Update Population MaBoSS), dedicated to modeling dynamic populations of interacting cells. We rely on the preexisting tool MaBoSS, which enables probabilistic simulations of cellular networks. A novel software layer is added to account for cell interactions and population dynamics, but without considering the spatial dimension. This modeling approach can be seen as an intermediate step towards more complex spatial descriptions. We illustrate our methodology by means of a case study dealing with TNF-induced cell death. Interestingly, the simulation of cell population dynamics with UPMaBoSS reveals a mechanism of resistance triggered by TNF treatment. Relatively easy to encode, UPMaBoSS simulations require only moderate computational power and execution time. To ease the reproduction of simulations, we provide several Jupyter notebooks that can be accessed within the CoLoMoTo Docker image, which contains all software and models used for this study.

An FVS-Based Approach to Attractor Detection in Asynchronous Random Boolean Networks

Boolean networks (BNs)play a crucial role in modeling and analyzing biological systems. One of the central issues in the analysis of BNs is attractor detection, i.e., identification of all possible attractors. This problem becomes more challenging for large asynchronous random Boolean networks (ARBNs)because of the asynchronous and non-deterministic updating scheme. In this paper, we present and formally prove several relations between feedback vertex sets (FVSs)and dynamics of BNs. From these relations, we propose an FVS-based method for detecting attractors in ARBNs. Our approach relies on the principle of removing arcs in the state transition graph to get a candidate set and the reachability property to filter the candidate set. We formally prove the correctness of our method and show its efficiency by conducting experiments on real biological networks and randomly generated N- K networks. The obtained results are very promising since our method can handle large networks whose sizes are up to 101 without using any network reduction technique.

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.

Navigating Multi-Scale Cancer Systems Biology Towards Model-Driven Clinical Oncology and Its Applications in Personalized Therapeutics

Rapid advancements in high-throughput omics technologies and experimental protocols have led to the generation of vast amounts of scale-specific biomolecular data on cancer that now populates several online databases and resources. Cancer systems biology models built using this data have the potential to provide specific insights into complex multifactorial aberrations underpinning tumor initiation, development, and metastasis. Furthermore, the annotation of these single- and multi-scale models with patient data can additionally assist in designing personalized therapeutic interventions as well as aid in clinical decision-making. Here, we have systematically reviewed the emergence and evolution of (i) repositories with scale-specific and multi-scale biomolecular cancer data, (ii) systems biology models developed using this data, (iii) associated simulation software for the development of personalized cancer therapeutics, and (iv) translational attempts to pipeline multi-scale panomics data for data-driven in silico clinical oncology. The review concludes that the absence of a generic, zero-code, panomics-based multi-scale modeling pipeline and associated software framework, impedes the development and seamless deployment of personalized in silico multi-scale models in clinical settings.

Patient-Specific Modeling of Diffuse Large B-Cell Lymphoma

Personalized medicine aims to tailor treatment to patients based on their individual genetic or molecular background. Especially in diseases with a large molecular heterogeneity, such as diffuse large B-cell lymphoma (DLBCL), personalized medicine has the potential to improve outcome and/or to reduce resistance towards treatment. However, integration of patient-specific information into a computational model is challenging and has not been achieved for DLBCL. Here, we developed a computational model describing signaling pathways and expression of critical germinal center markers. The model integrates the regulatory mechanism of the signaling and gene expression network and covers more than 50 components, many carrying genetic lesions common in DLBCL. Using clinical and genomic data of 164 primary DLBCL patients, we implemented mutations, structural variants and copy number alterations as perturbations in the model using the CoLoMoTo notebook. Leveraging patient-specific genotypes and simulation of the expression of marker genes in specific germinal center conditions allows us to predict the consequence of the modeled pathways for each patient. Finally, besides modeling how genetic perturbations alter physiological signaling, we also predicted for each patient model the effect of rational inhibitors, such as Ibrutinib, that are currently discussed as possible DLBCL treatments, showing patient-dependent variations in effectiveness and synergies.

Logical modelling of in vitro differentiation of human monocytes into dendritic cells unravels novel transcriptional regulatory interactions

Dendritic cells (DCs) are the major specialized antigen-presenting cells, thereby connecting innate and adaptive immunity. Because of their role in establishing adaptive immunity, they constitute promising targets for immunotherapy. Monocytes can differentiate into DCs in vitro in the presence of colony-stimulating factor 2 (CSF2) and interleukin 4 (IL4), activating four signalling pathways (MAPK, JAK/STAT, NFKB and PI3K). However, the downstream transcriptional programme responsible for DC differentiation from monocytes (moDCs) remains unknown. By analysing the scientific literature on moDC differentiation, we established a preliminary logical model that helped us identify missing information regarding the activation of genes responsible for this differentiation, including missing targets for key transcription factors (TFs). Using ChIP-seq and RNA-seq data from the Blueprint consortium, we defined active and inactive promoters, together with differentially expressed genes in monocytes, moDCs and macrophages, which correspond to an alternative cell fate. We then used this functional genomic information to predict novel targets for previously identified TFs. By integrating this information, we refined our model and recapitulated the main established facts regarding moDC differentiation. Prospectively, the resulting model should be useful to develop novel immunotherapies targeting moDCs.

In Silico Logical Modelling to Uncover Cooperative Interactions in Cancer

The multistep development of cancer involves the cooperation between multiple molecular lesions, as well as complex interactions between cancer cells and the surrounding tumour microenvironment. The search for these synergistic interactions using experimental models made tremendous contributions to our understanding of oncogenesis. Yet, these approaches remain labour-intensive and challenging. To tackle such a hurdle, an integrative, multidisciplinary effort is required. In this article, we highlight the use of logical computational models, combined with experimental validations, as an effective approach to identify cooperative mechanisms and therapeutic strategies in the context of cancer biology. In silico models overcome limitations of reductionist approaches by capturing tumour complexity and by generating powerful testable hypotheses. We review representative examples of logical models reported in the literature and their validation. We then provide further analyses of our logical model of Epithelium to Mesenchymal Transition (EMT), searching for additional cooperative interactions involving inputs from the tumour microenvironment and gain of function mutations in NOTCH.

Personalized logical models to investigate cancer response to BRAF treatments in melanomas and colorectal cancers

The study of response to cancer treatments has benefited greatly from the contribution of different omics data but their interpretation is sometimes difficult. Some mathematical models based on prior biological knowledge of signaling pathways facilitate this interpretation but often require fitting of their parameters using perturbation data. We propose a more qualitative mechanistic approach, based on logical formalism and on the sole mapping and interpretation of omics data, and able to recover differences in sensitivity to gene inhibition without model training. This approach is showcased by the study of BRAF inhibition in patients with melanomas and colorectal cancers who experience significant differences in sensitivity despite similar omics profiles. We first gather information from literature and build a logical model summarizing the regulatory network of the mitogen-activated protein kinase (MAPK) pathway surrounding BRAF, with factors involved in the BRAF inhibition resistance mechanisms. The relevance of this model is verified by automatically assessing that it qualitatively reproduces response or resistance behaviors identified in the literature. Data from over 100 melanoma and colorectal cancer cell lines are then used to validate the model's ability to explain differences in sensitivity. This generic model is transformed into personalized cell line-specific logical models by integrating the omics information of the cell lines as constraints of the model. The use of mutations alone allows personalized models to correlate significantly with experimental sensitivities to BRAF inhibition, both from drug and CRISPR targeting, and even better with the joint use of mutations and RNA, supporting multi-omics mechanistic models. A comparison of these untrained models with learning approaches highlights similarities in interpretation and complementarity depending on the size of the datasets. This parsimonious pipeline, which can easily be extended to other biological questions, makes it possible to explore the mechanistic causes of the response to treatment, on an individualized basis.

Modelling of Immune Checkpoint Network Explains Synergistic Effects of Combined Immune Checkpoint Inhibitor Therapy and the Impact of Cytokines in Patient Response

Simple Summary

The future of cancer immunotherapy relies on a combination of individually targeted therapies. However, a lot of experiments are needed to define the most effective combinations of drugs. A computational and modelling approach could help reduce the number of experiments and suggest optimal treatments to test. This article presents a logical model of T cell activation influenced by immune checkpoints, and explores the effect of these checkpoints, suggests mechanisms that would explain why some treatments might be better suited than others. The model includes not only programmed cell death protein 1 (PD1) and cytotoxic T-lymphocyte-associated protein 4 (CTL4) downstream pathways but also those of other immune checkpoints such as T cell immunoglobulin and ITIM (immunoreceptor tyrosine-based inhibition motif) domain (TIGIT), lymphocyte activation gene 3 (LAG3), T cell immunoglobulin and mucin domain-containing protein 3 (TIM3), cluster of differentiation 226 (CD226), inducible T-cell costimulator (ICOS), and tumour necrosis factor receptors (TNFRs).

Abstract

After the success of the new generation of immune therapies, immune checkpoint receptors have become one important center of attention of molecular oncologists. The initial success and hopes of anti-programmed cell death protein 1 (anti-PD1) and anti-cytotoxic T-lymphocyte-associated protein 4 (anti-CTLA4) therapies have shown some limitations since a majority of patients have continued to show resistance. Other immune checkpoints have raised some interest and are under investigation, such as T cell immunoglobulin and ITIM (immunoreceptor tyrosine-based inhibition motif) domain (TIGIT), inducible T-cell costimulator (ICOS), and T cell immunoglobulin and mucin domain-containing protein 3 (TIM3), which appear as promising targets for immunotherapy. To explore their role and study possible synergetic effects of these different checkpoints, we have built a model of T cell receptor (TCR) regulation including not only PD1 and CTLA4, but also other well studied checkpoints (TIGIT, TIM3, lymphocyte activation gene 3 (LAG3), cluster of differentiation 226 (CD226), ICOS, and tumour necrosis factor receptors (TNFRs)) and simulated different aspects of T cell biology. Our model shows good correspondence with observations from available experimental studies of anti-PD1 and anti-CTLA4 therapies and suggest efficient combinations of immune checkpoint inhibitors (ICI). Among the possible candidates, TIGIT appears to be the most promising drug target in our model. The model predicts that signal transducer and activator of transcription 1 (STAT1)/STAT4-dependent pathways, activated by cytokines such as interleukin 12 (IL12) and interferon gamma (IFNG), could improve the effect of ICI therapy via upregulation of Tbet, suggesting that the effect of the cytokines related to STAT3/STAT1 activity is dependent on the balance between STAT1 and STAT3 downstream signalling.

Simulation of multiple microenvironments shows a pivot role of RPTPs on the control of Epithelial-to-Mesenchymal Transition

Epithelial-to-Mesenchymal Transition (EMT) is a natural and reversible process involved in embryogenesis, wound healing and thought to participate in the process of metastasis. Multiple signals from the microenvironment have been reported to drive EMT. However, the tight control of this process on physiological scenarios and how it is disrupted during cancer progression is not fully understood. Here, we analysed a regulatory network of EMT accounting for 10 key microenvironment signals focusing on the impact of two cell contact signals on the reversibility of EMT and the stability of resulting phenotypes. The analysis showed that the microenvironment is not enough for stabilizing Hybrid and Amoeboid-like phenotypes, requiring intracellular de-regulations as reported during cancer progression. Our simulations demonstrated that RPTP activation by cell contacts have the potential to inhibit the process of EMT and trigger its reversibility under tissue growth and chronic inflammation scenarios. Simulations also showed that hypoxia inhibits the capacity of RPTPs to control EMT. Our analysis further provided a theoretical explanation for the observed correlation between hypoxia and metastasis under chronic inflammation, and predicted that de-regulations in FAT4 signalling may promote Hybrid stabilization. Taken together, we propose a natural control mechanism of EMT that supports the idea that EMT is tightly regulated by the microenvironment.

The impact of ARID1A mutation on molecular characteristics in colorectal cancer

BACKGROUND

ARID1A is a key subunit of the SWItch/Sucrose Non-Fermentable (SWI/SNF) complex which regulates dynamic repositioning of nucleosomes to repair DNA damage. Only small pilot studies have evaluated the role of ARID1A mutation in colorectal cancer (CRC). The aim of the present study was to explore the potential impact of ARID1A mutation on clinicopathological and molecular characteristics in CRC.

METHODS

We used integrated data sets of 7978 CRC cases (one data set from a clinical laboratory improvement amendments [CLIA]-certified laboratory and three independent published data sets). The associations of ARID1A mutation with molecular characteristics including immune profile (the status of microsatellite instability [MSI], tumour mutational burden [TMB], programmed death ligand 1 [PD-L1] and estimated infiltrating immune cells), clinicopathological features and related pathways were analysed using next-generation sequencing, RNA sequencing and immunohistochemistry.

RESULTS

ARID1A mutant samples had more genomically unstable tumour features (MSI-high and TMB-high) and exhibited more characteristics of a T-cell-inflamed microenvironment (PD-L1 expression and high estimated infiltrating cytotoxic T lymphocytes [CTLs]) than ARID1A wild-type samples in the discovery and validation cohorts. Even ARID1A mutant samples without MSI-high status were TMB-high, had high levels of PD-L1 expression and high estimated infiltrating CTLs. ARID1A mutations were more common with right-sided primary and earlier stage tumours. ARID1A mutant tumours mainly had co-occurring gene mutations related to chromatin modifying, DNA repair, WNT signalling and epidermal growth factor receptor inhibitor resistance pathways, and ARID1A mutations strongly regulated DNA repair pathways. Key genes for chemotherapy/radiotherapy sensitivity were suppressed in ARID1A mutant samples.

CONCLUSIONS

Our findings may provide novel insights to develop individualised approaches for treatment of CRC based on ARID1A mutation status.

Analysis, identification and visualization of subgroups in genomics

Motivation

Cancer is a complex and heterogeneous disease involving multiple somatic mutations that accumulate during its progression. In the past years, the wide availability of genomic data from patients’ samples opened new perspectives in the analysis of gene mutations and alterations. Hence, visualizing and further identifying genes mutated in massive sets of patients are nowadays a critical task that sheds light on more personalized intervention approaches.

Results

Here, we extensively review existing tools for visualization and analysis of alteration data. We compare different approaches to study mutual exclusivity and sample coverage in large-scale omics data. We complement our review with the standalone software AVAtar (‘analysis and visualization of alteration data’) that integrates diverse aspects known from different tools into a comprehensive platform. AVAtar supplements customizable alteration plots by a multi-objective evolutionary algorithm for subset identification and provides an innovative and user-friendly interface for the evaluation of concurrent solutions. A use case from personalized medicine demonstrates its unique features showing an application on vaccination target selection.

Availability

AVAtar is available at: https://github.com/sysbio-bioinf/avatar

Contact

hans.kestler@uni-ulm.de, phone: +49 (0) 731 500 24 500, fax: +49 (0) 731 500 24 502

Reconciling qualitative, abstract, and scalable modeling of biological networks

Predicting biological systems' behaviors requires taking into account many molecular and genetic elements for which limited information is available past a global knowledge of their pairwise interactions. Logical modeling, notably with Boolean Networks (BNs), is a well-established approach that enables reasoning on the qualitative dynamics of networks. Several dynamical interpretations of BNs have been proposed. The synchronous and (fully) asynchronous ones are the most prominent, where the value of either all or only one component can change at each step. Here we prove that, besides being costly to analyze, these usual interpretations can preclude the prediction of certain behaviors observed in quantitative systems. We introduce an execution paradigm, the Most Permissive Boolean Networks (MPBNs), which offers the formal guarantee not to miss any behavior achievable by a quantitative model following the same logic. Moreover, MPBNs significantly reduce the complexity of dynamical analysis, enabling to model genome-scale networks.

Genetic alterations in the 3q26.31-32 locus confer an aggressive prostate cancer phenotype

Large-scale genetic aberrations that underpin prostate cancer development and progression, such as copy-number alterations (CNAs), have been described but the consequences of specific changes in many identified loci is limited. Germline SNPs in the 3q26.31 locus are associated with aggressive prostate cancer, and is the location of NAALADL2, a gene overexpressed in aggressive disease. The closest gene to NAALADL2 is TBL1XR1, which is implicated in tumour development and progression. Using publicly-available cancer genomic data we report that NAALADL2 and TBL1XR1 gains/amplifications are more prevalent in aggressive sub-types of prostate cancer when compared to primary cohorts. In primary disease, gains/amplifications occurred in 15.99% (95% CI: 13.02–18.95) and 14.96% (95% CI: 12.08–17.84%) for NAALADL2 and TBL1XR1 respectively, increasing in frequency in higher Gleason grade and stage tumours. Gains/amplifications result in transcriptional changes and the development of a pro-proliferative and aggressive phenotype. These results support a pivotal role for copy-number gains in this genetic region.

Benjamin Simpson et al. use publicly available cancer genomic data to investigate copy number changes at the 3q26.31–32 locus, which has been associated with aggressive prostate cancer based on single-nucleotide polymorphisms. They find that gains of NAALADL2 and TBL1XR1 in this locus are associated with more aggressive subtypes of prostate cancer and the transcription of pro-proliferative signalling processes.

Integrated multi-omics data analyses for exploring the co-occurring and mutually exclusive gene alteration events in colorectal cancer

Co-occurring and mutually exclusive gene alteration events are helpful for understanding carcinogenesis but systematic screening for such events is quite limited. We conducted pairwise screening tests to identify "hit pairs" in colorectal cancer (CRC) by utilizing the cross-omics data from The Cancer Genome Atlas (TCGA). Numerous hit pairs involving somatic mutations, copy number variations, and DNA methylation were found to occur nonrandomly in CRC, such as KRAS and HOXB6, SMAD4 and PMEPA1. Based on these hit pairs, we identified 32 synthetic lethal pairs and 7,527 co-occurring pairs relating to drug response. Our further biological experiments showed that the co-occurrence of mutant FCGBP and NUDT12 silencing (or mutant TMC3 and RPS6KA6 silencing) with small interfering RNA reduced cell viability. Moreover, novel hit pairs could influence prognosis. The patients who carried concurrent mutations of IRF5 and NEFH, SYNE1 and TTN, or MUC16 and NEFH had worse survival outcomes. Particularly, the presence of mutant SYNE1 and TTN pair not only affects prognosis, but also is related to CRC patients' response to drug treatment. Our "hit pair" genes may provide insights into colorectal carcinogenesis and help open new avenues for CRC therapy.

Basins of Attraction, Commitment Sets, and Phenotypes of Boolean Networks

The attractors of Boolean networks and their basins have been shown to be highly relevant for model validation and predictive modeling, e.g., in systems biology. Yet, there are currently very few tools available that are able to compute and visualize not only attractors but also their basins. In the realm of asynchronous, non-deterministic modeling not only is the repertoire of software even more limited, but also the formal notions for basins of attraction are often lacking. In this setting, the difficulty both for theory and computation arises from the fact that states may be elements of several distinct basins. In this paper, we address this topic by partitioning the state space into sets that are committed to the same attractors. These commitment sets can easily be generalized to sets that are equivalent w.r.t. the long-term behaviors of pre-selected nodes which leads us to the notions of markers and phenotypes which we illustrate in a case study on bladder tumorigenesis. For every concept, we propose equivalent CTL model checking queries and an extension of the state of the art model checking software NuSMV is made available that is capable of computing the respective sets. All notions are fully integrated as three new modules in our Python package PyBoolNet, including functions for visualizing the basins, commitment sets, and phenotypes as quotient graphs and pie charts.

Genetic hallmarks of recurrent/metastatic adenoid cystic carcinoma

BACKGROUNDAdenoid cystic carcinoma (ACC) is a rare malignancy arising in salivary glands and other sites, characterized by high rates of relapse and distant spread. Recurrent/metastatic (R/M) ACCs are generally incurable, due to a lack of active systemic therapies. To improve outcomes, deeper understanding of genetic alterations and vulnerabilities in R/M tumors is needed.METHODSAn integrated genomic analysis of 1,045 ACCs (177 primary, 868 R/M) was performed to identify alterations associated with advanced and metastatic tumors. Intratumoral genetic heterogeneity, germline mutations, and therapeutic actionability were assessed.RESULTSCompared with primary tumors, R/M tumors were enriched for alterations in key Notch (NOTCH1, 26.3% vs. 8.5%; NOTCH2, 4.6% vs. 2.3%; NOTCH3, 5.7% vs. 2.3%; NOTCH4, 3.6% vs. 0.6%) and chromatin-remodeling (KDM6A, 15.2% vs. 3.4%; KMT2C/MLL3, 14.3% vs. 4.0%; ARID1B, 14.1% vs. 4.0%) genes. TERT promoter mutations (13.1% of R/M cases) were mutually exclusive with both NOTCH1 mutations (q = 3.3 × 10-4) and MYB/MYBL1 fusions (q = 5.6 × 10-3), suggesting discrete, alternative mechanisms of tumorigenesis. This network of alterations defined 4 distinct ACC subgroups: MYB+NOTCH1+, MYB+/other, MYBWTNOTCH1+, and MYBWTTERT+. Despite low mutational load, we identified numerous samples with marked intratumoral genetic heterogeneity, including branching evolution across multiregion sequencing.CONCLUSIONThese observations collectively redefine the molecular underpinnings of ACC progression and identify further targets for precision therapies.FUNDINGAdenoid Cystic Carcinoma Research Foundation, Pershing Square Sohn Cancer Research grant, the PaineWebber Chair, Stand Up 2 Cancer, NIH R01 CA205426, the STARR Cancer Consortium, NCI R35 CA232097, the Frederick Adler Chair, Cycle for Survival, the Jayme Flowers Fund, The Sebastian Nativo Fund, NIH K08 DE024774 and R01 DE027738, and MSKCC through NIH/NCI Cancer Center Support Grant (P30 CA008748).

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.

Control of Intracellular Molecular Networks Using Algebraic Methods

Many problems in biology and medicine have a control component. Often, the goal might be to modify intracellular networks, such as gene regulatory networks or signaling networks, in order for cells to achieve a certain phenotype, what happens in cancer. If the network is represented by a mathematical model for which mathematical control approaches are available, such as systems of ordinary differential equations, then this problem might be solved systematically. Such approaches are available for some other model types, such as Boolean networks, where structure-based approaches have been developed, as well as stable motif techniques. However, increasingly many published discrete models are mixed-state or multistate, that is, some or all variables have more than two states, and thus the development of control strategies for multistate networks is needed. This paper presents a control approach broadly applicable to general multistate models based on encoding them as polynomial dynamical systems over a finite algebraic state set, and using computational algebra for finding appropriate intervention strategies. To demonstrate the feasibility and applicability of this method, we apply it to a recently developed multistate intracellular model of E2F-mediated bladder cancerous growth and to a model linking intracellular iron metabolism and oncogenic pathways. The control strategies identified for these published models are novel in some cases and represent new hypotheses, or are supported by the literature in others as potential drug targets. Our Macaulay2 scripts to find control strategies are publicly available through GitHub at https://github.com/luissv7/multistatepdscontrol.

Algorithms for the Sequential Reprogramming of Boolean Networks

Cellular reprogramming, a technique that opens huge opportunities in modern and regenerative medicine, heavily relies on identifying key genes to perturb. Most of the existing computational methods for controlling which attractor (steady state) the cell will reach focus on finding mutations to apply to the initial state. However, it has been shown, and is proved in this article, that waiting between perturbations so that the update dynamics of the system prepares the ground, allows for new reprogramming strategies. To identify such sequential perturbations, we consider a qualitative model of regulatory networks, and rely on Binary Decision Diagrams to model their dynamics and the putative perturbations. Our method establishes a set identification of sequential perturbations, whether permanent (mutations) or only temporary, to achieve the existential or inevitable reachability of an arbitrary state of the system. We apply an implementation for temporary perturbations on models from the literature, illustrating that we are able to derive sequential perturbations to achieve trans-differentiation.

Connecting signaling and metabolic pathways in EGF receptor-mediated oncogenesis of glioblastoma

As malignant transformation requires synchronization of growth-driving signaling (S) and metabolic (M) pathways, defining cancer-specific S-M interconnected networks (SMINs) could lead to better understanding of oncogenic processes. In a systems-biology approach, we developed a mathematical model for SMINs in mutated EGF receptor (EGFRvIII) compared to wild-type EGF receptor (EGFRwt) expressing glioblastoma multiforme (GBM). Starting with experimentally validated human protein-protein interactome data for S-M pathways, and incorporating proteomic data for EGFRvIII and EGFRwt GBM cells and patient transcriptomic data, we designed a dynamic model for EGFR-driven GBM-specific information flow. Key nodes and paths identified by in silico perturbation were validated experimentally when inhibition of signaling pathway proteins altered expression of metabolic proteins as predicted by the model. This demonstrated capacity of the model to identify unknown connections between signaling and metabolic pathways, explain the robustness of oncogenic SMINs, predict drug escape, and assist identification of drug targets and the development of combination therapies.

Complex and highly dynamic interconnected networks allow cancer to take different routes and circumvent chemotherapy. Therefore, understanding these context-specific networks and their dynamics of molecular interactions driven by different oncogenic signaling and metabolic pathways is very much needed to predict drug targets and the effect of therapeutics. We incorporated high-throughput transcriptome and proteome data into mathematical models to deduce properties of cancer cells through systems biology approach. Here we report the development, testing and validation of an integrated systems biology model of information flow between signaling and metabolic pathways to understand the regulation of the interconnection between them in cancer. Our model efficiently identified unique connections and key nodes important in signaling-metabolic information flow. We predicted some potential novel targets before performing actual drug tests. We have successfully applied this model to identify the interconnections altered in the constitutive signaling of the mutated EGFR by comparing EGF-dependent and wild-type EGFR signaling in glioblastoma multiforme.

Identifying mutual exclusivity across cancer genomes: computational approaches to discover genetic interaction and reveal tumor vulnerability

Systematic sequencing of cancer genomes has revealed prevalent heterogeneity, with patients harboring various combinatorial patterns of genetic alteration. In particular, a phenomenon that a group of genes exhibits mutually exclusive patterns has been widespread across cancers, covering a broad spectrum of crucial cancer pathways. Recently, there is considerable evidence showing that, mutual exclusivity reflects alternative functions in tumor initiation and progression, or suggests adverse effects of their concurrence. Given its importance, numerous computational approaches have been proposed to study mutual exclusivity using genomic profiles alone, or by integrating networks and phenotypes. Some of them have been routinely used to explore genetic associations, which lead to a deeper understanding of carcinogenic mechanisms and reveals unexpected tumor vulnerabilities. Here, we present an overview of mutual exclusivity from the perspective of cancer genome. We describe the common hypothesis underlying mutual exclusivity, summarize the strategies for the identification of significant mutually exclusive patterns, compare the performance of representative algorithms from simulated data sets and discuss their common confounders.

Controlling large Boolean networks with single-step perturbations

Motivation

The control of Boolean networks has traditionally focussed on strategies where the perturbations are applied to the nodes of the network for an extended period of time. In this work, we study if and how a Boolean network can be controlled by perturbing a minimal set of nodes for a single-step and letting the system evolve afterwards according to its original dynamics. More precisely, given a Boolean network (BN), we compute a minimal subset Cmin of the nodes such that BN can be driven from any initial state in an attractor to another ‘desired’ attractor by perturbing some or all of the nodes of Cmin for a single-step. Such kind of control is attractive for biological systems because they are less time consuming than the traditional strategies for control while also being financially more viable. However, due to the phenomenon of state-space explosion, computing such a minimal subset is computationally inefficient and an approach that deals with the entire network in one-go, does not scale well for large networks.

Results

We develop a ‘divide-and-conquer’ approach by decomposing the network into smaller partitions, computing the minimal control on the projection of the attractors to these partitions and then composing the results to obtain Cmin for the whole network. We implement our method and test it on various real-life biological networks to demonstrate its applicability and efficiency.

Supplementary information

Supplementary data are available at Bioinformatics online.

Personalization of Logical Models With Multi-Omics Data Allows Clinical Stratification of Patients

Logical models of cancer pathways are typically built by mining the literature for relevant experimental observations. They are usually generic as they apply for large cohorts of individuals. As a consequence, they generally do not capture the heterogeneity of patient tumors and their therapeutic responses. We present here a novel framework, referred to as PROFILE, to tailor logical models to a particular biological sample such as a patient tumor. This methodology permits to compare the model simulations to individual clinical data, i.e., survival time. Our approach focuses on integrating mutation data, copy number alterations (CNA), and expression data (transcriptomics or proteomics) to logical models. These data need first to be either binarized or set between 0 and 1, and can then be incorporated in the logical model by modifying the activity of the node, the initial conditions or the state transition rates. The use of MaBoSS, a tool based on Monte-Carlo kinetic algorithm to perform stochastic simulations on logical models results in model state probabilities, and allows for a semi-quantitative study of the model phenotypes and perturbations. As a proof of concept, we use a published generic model of cancer signaling pathways and molecular data from METABRIC breast cancer patients. For this example, we test several combinations of data incorporation and discuss that, with these data, the most comprehensive patient-specific cancer models are obtained by modifying the nodes' activity of the model with mutations, in combination or not with CNA data, and altering the transition rates with RNA expression. We conclude that these model simulations show good correlation with clinical data such as patients' Nottingham prognostic index (NPI) subgrouping and survival time. We observe that two highly relevant cancer phenotypes derived from personalized models, Proliferation and Apoptosis, are biologically consistent prognostic factors: patients with both high proliferation and low apoptosis have the worst survival rate, and conversely. Our approach aims to combine the mechanistic insights of logical modeling with multi-omics data integration to provide patient-relevant models. This work leads to the use of logical modeling for precision medicine and will eventually facilitate the choice of patient-specific drug treatments by physicians.

Molecular characterization of low grade and high grade bladder cancer

Background

Bladder cancer (BC) is the 9^th most common cancer diagnosis worldwide. Low grade (LG) represents 70% of all BCs, characterized by recurrence and rare ability (10–15%) to progress to high grade (HG) and invade. The remaining 30% is high grade (HG), fast invasive BC, which is resistant to therapy. Identifying biomarkers for predicting those tumors able to progress is a key goal for patient outcome improvement. This study focuses on the most promising prognostic markers.

Materials and methods

TP53 and FGFR3 mutational status, Survivin, CK19, CK20, E-cadherin and CD44 gene expression analysis were performed on 66 BCs.

Results

Survivin was found associated to tumor grade (p<0.05). Moreover, Survivin correlated with CD44 in TP53 wild type (p = 0.0242) and FGFR3 wild type (p = 0.0036) tumors. In particular the Survivin-CD44 correlation was associated to HG FGFR3 wild type BCs (p = 0.0045). Unsupervised hierarchical clustering based on gene expression data identified four distinct molecular groups reflecting the patient histology (p = 0.038).

Conclusion

We suggest Survivin, both as a biomarker associated to G3 BCs but negatively related to TP53 mutational status, and as a potential novel therapeutic target.

A pan-cancer atlas of cancer hallmark-associated candidate driver lncRNAs

Substantial cancer genome sequencing efforts have discovered many important driver genes contributing to tumorigenesis. However, very little is known about the genetic alterations of long non‐coding RNAs (lncRNAs) in cancer. Thus, there is a need for systematic surveys of driver lncRNAs. Through integrative analysis of 5918 tumors across 11 cancer types, we revealed that lncRNAs have undergone dramatic genomic alterations, many of which are mutually exclusive with well‐known cancer genes. Using the hypothesis of functional redundancy of mutual exclusivity, we developed a computational framework to identify driver lncRNAs associated with different cancer hallmarks. Applying it to pan‐cancer data, we identified 378 candidate driver lncRNAs whose genomic features highly resemble the known cancer driver genes (e.g. high conservation and early replication). We further validated the candidate driver lncRNAs involved in ‘Tissue Invasion and Metastasis’ in lung adenocarcinoma and breast cancer, and also highlighted their potential roles in improving clinical outcomes. In summary, we have generated a comprehensive landscape of cancer candidate driver lncRNAs that could act as a starting point for future functional explorations, as well as the identification of biomarkers and lncRNA‐based target therapy.

Estimating Attractor Reachability in Asynchronous Logical Models

Logical models are well-suited to capture salient dynamical properties of regulatory networks. For networks controlling cell fate decisions, cell fates are associated with model attractors (stable states or cyclic attractors) whose identification and reachability properties are particularly relevant. While synchronous updates assume unlikely instantaneous or identical rates associated with component changes, the consideration of asynchronous updates is more realistic but, for large models, may hinder the analysis of the resulting non-deterministic concurrent dynamics. This complexity hampers the study of asymptotical behaviors, and most existing approaches suffer from efficiency bottlenecks, being generally unable to handle cyclical attractors and quantify attractor reachability. Here, we propose two algorithms providing probability estimates of attractor reachability in asynchronous dynamics. The first algorithm, named Firefront, exhaustively explores the state space from an initial state, and provides quasi-exact evaluations of the reachability probabilities of model attractors. The algorithm progresses in breadth, propagating the probabilities of each encountered state to its successors. Second, Avatar is an adapted Monte Carlo approach, better suited for models with large and intertwined transient and terminal cycles. Avatar iteratively explores the state space by randomly selecting trajectories and by using these random walks to estimate the likelihood of reaching an attractor. Unlike Monte Carlo simulations, Avatar is equipped to avoid getting trapped in transient cycles and to identify cyclic attractors. Firefront and Avatar are validated and compared to related methods, using as test cases logical models of synthetic and biological networks. Both algorithms are implemented as new functionalities of GINsim 3.0, a well-established software tool for logical modeling, providing executable GUI, Java API, and scripting facilities.

Logical Modeling and Analysis of Cellular Regulatory Networks With GINsim 3.0

The logical formalism is well adapted to model large cellular networks, in particular when detailed kinetic data are scarce. This tutorial focuses on this well-established qualitative framework. Relying on GINsim (release 3.0), a software implementing this formalism, we guide the reader step by step toward the definition, the analysis and the simulation of a four-node model of the mammalian p53-Mdm2 network.

The Discovery of Mutated Driver Pathways in Cancer: Models and Algorithms

The pathogenesis of cancer in human is still poorly understood. With the rapid development of high-throughput sequencing technologies, huge volumes of cancer genomics data have been generated. Deciphering that data poses great opportunities and challenges to computational biologists. One of such key challenges is to distinguish driver mutations, genes as well as pathways from passenger ones. Mutual exclusivity of gene mutations (each patient has no more than one mutation in the gene set) has been observed in various cancer types and thus has been used as an important property of a driver gene set or pathway. In this article, we aim to review the recent development of computational models and algorithms for discovering driver pathways or modules in cancer with the focus on mutual exclusivity-based ones.

Therapeutic target discovery using Boolean network attractors: improvements of kali

In a previous article, an algorithm for identifying therapeutic targets in Boolean networks modelling pathological mechanisms was introduced. In the present article, the improvements made on this algorithm, named kali, are described. These improvements are (i) the possibility to work on asynchronous Boolean networks, (ii) a finer assessment of therapeutic targets and (iii) the possibility to use multivalued logic. kali assumes that the attractors of a dynamical system, such as a Boolean network, are associated with the phenotypes of the modelled biological system. Given a logic-based model of pathological mechanisms, kali searches for therapeutic targets able to reduce the reachability of the attractors associated with pathological phenotypes, thus reducing their likeliness. kali is illustrated on an example network and used on a biological case study. The case study is a published logic-based model of bladder tumorigenesis from which kali returns consistent results. However, like any computational tool, kali can predict but cannot replace human expertise: it is a supporting tool for coping with the complexity of biological systems in the field of drug discovery.

Inferring hidden causal relations between pathway members using reduced Google matrix of directed biological networks

Signaling pathways represent parts of the global biological molecular network which connects them into a seamless whole through complex direct and indirect (hidden) crosstalk whose structure can change during development or in pathological conditions. We suggest a novel methodology, called Googlomics, for the structural analysis of directed biological networks using spectral analysis of their Google matrices, using parallels with quantum scattering theory, developed for nuclear and mesoscopic physics and quantum chaos. We introduce analytical "reduced Google matrix" method for the analysis of biological network structure. The method allows inferring hidden causal relations between the members of a signaling pathway or a functionally related group of genes. We investigate how the structure of hidden causal relations can be reprogrammed as a result of changes in the transcriptional network layer during cancerogenesis. The suggested Googlomics approach rigorously characterizes complex systemic changes in the wiring of large causal biological networks in a computationally efficient way.

CM-viewer: Visualizing interaction network of co-mutated and mutually exclusively mutated cancer genes

Cancer genes usually play a crucial role in regulating cell growth. Normal cells transform into malignant tumors by the acquisition of accumulated genetic mutations that enable them to evade normal growth control. It is therefore important to understand the relationships between mutations during cancer development and progression. Although cancer genes with co-occurring and mutually exclusive mutations have already been studied on different scales, there is no timely updated interaction network available for co-mutated and mutually exclusively mutated cancer genes. Therefore, we firstly downloaded 567 cancer genes from COSMIC (catalogue of somatic mutations in cancer) cancer gene census. Secondly, somatic mutations of 71 cancer genomics projects were downloaded from the ICGC (International Cancer Genome Consortium) data portal. Thirdly, mutated cancer genes and affected donors were extracted from the ICGC data to form a mutation matrix where rows are genes, columns are donors, 1 denotes occurrence, and 0 denotes absence of mutation. Afterwards, co-mutated and mutually exclusively mutated cancer gene pairs were identified using DISCOVER (discrete independence statistic controlling for observations with varying event rates). Finally, CM-viewer was developed to visualize the interaction network of cancer genes with co-occurring and mutually exclusive mutations. It is an online visualization tool as well as a biological database. It promises to understand how gene mutations contribute to tumorigenesis and to identify key biomarkers and drug targets for cancer. CM-viewer is freely available at http://www.zhounan.org/comutgene.

Logical model specification aided by model-checking techniques: application to the mammalian cell cycle regulation

MOTIVATION

Understanding the temporal behaviour of biological regulatory networks requires the integration of molecular information into a formal model. However, the analysis of model dynamics faces a combinatorial explosion as the number of regulatory components and interactions increases.

RESULTS

We use model-checking techniques to verify sophisticated dynamical properties resulting from the model regulatory structure in the absence of kinetic assumption. We demonstrate the power of this approach by analysing a logical model of the molecular network controlling mammalian cell cycle. This approach enables a systematic analysis of model properties, the delineation of model limitations, and the assessment of various refinements and extensions based on recent experimental observations. The resulting logical model accounts for the main irreversible transitions between cell cycle phases, the sequential activation of cyclins, and the inhibitory role of Skp2, and further emphasizes the multifunctional role for the cell cycle inhibitor Rb.

AVAILABILITY AND IMPLEMENTATION

The original and revised mammalian cell cycle models are available in the model repository associated with the public modelling software GINsim (http://ginsim.org/node/189).

CONTACT

thieffry@ens.fr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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).

Discovery of cancer common and specific driver gene sets

Cancer is known as a disease mainly caused by gene alterations. Discovery of mutated driver pathways or gene sets is becoming an important step to understand molecular mechanisms of carcinogenesis. However, systematically investigating commonalities and specificities of driver gene sets among multiple cancer types is still a great challenge, but this investigation will undoubtedly benefit deciphering cancers and will be helpful for personalized therapy and precision medicine in cancer treatment. In this study, we propose two optimization models to de novo discover common driver gene sets among multiple cancer types (ComMDP) and specific driver gene sets of one certain or multiple cancer types to other cancers (SpeMDP), respectively. We first apply ComMDP and SpeMDP to simulated data to validate their efficiency. Then, we further apply these methods to 12 cancer types from The Cancer Genome Atlas (TCGA) and obtain several biologically meaningful driver pathways. As examples, we construct a common cancer pathway model for BRCA and OV, infer a complex driver pathway model for BRCA carcinogenesis based on common driver gene sets of BRCA with eight cancer types, and investigate specific driver pathways of the liquid cancer lymphoblastic acute myeloid leukemia (LAML) versus other solid cancer types. In these processes more candidate cancer genes are also found.

Targeted therapies in the treatment of urothelial cancers

Progress has been slow in systemic management of locally advanced and metastatic bladder cancer over the past 20 years. However, the recent approval of immunotherapy with atezolizumab and nivolumab for second-line salvage therapy may usher in an era of more rapid improvement. Systemic treatment is suboptimal and is an area of substantial unmet medical need. The recent findings from The Cancer Genome Atlas project revealed promising pathways that may be amenable to targeted therapies. Promising results with treatment using vascular endothelial growth factor inhibitors such as ramucirumab, sunitinib or bevacizumab, and human epidermal growth factor receptor 2 targeted therapies, epidermal growth factor receptor inhibitors, and fibroblast growth factor receptor inhibitors, are undergoing clinical trials and are discussed later.