ViSiBooL-visualization and simulation of Boolean networks with temporal constraints

Boolean Modeling of Biological Network Applied to Protein-Protein Interaction Network of Autism Patients

Cellular molecules interact with one another in a structured manner, defining a regulatory network topology that describes cellular mechanisms. Genetic mutations alter these networks' pathways, generating complex disorders such as autism spectrum disorder (ASD). Boolean models have assisted in understanding biological system dynamics since Kauffman's 1969 discovery, and various analytical tools for regulatory networks have been developed. This study examined the protein-protein interaction network created in our previous publication of four ASD patients using the SPIDDOR R package, a Boolean model-based method. The aim is to examine how patients' genetic variations in INTS6L, USP9X, RSK4, FGF5, FLNA, SUMF1, and IDS affect mTOR and Wnt cell signaling convergence. The Boolean network analysis revealed abnormal activation levels of essential proteins such as β-catenin, MTORC1, RPS6, eIF4E, Cadherin, and SMAD. These proteins affect gene expression, translation, cell adhesion, shape, and migration. Patients 1 and 2 showed consistent patterns of increased β-catenin activity and decreased MTORC1, RPS6, and eIF4E activity. However, patient 2 had an independent decrease in Cadherin and SMAD activity due to the FLNA mutation. Patients 3 and 4 have an abnormal activation of the mTOR pathway, which includes the MTORC1, RPS6, and eIF4E genes. The shared mTOR pathway behavior in these patients is explained by a shared mutation in two closely related proteins (SUMF1 and IDS). Diverse activities in β-catenin, MTORC1, RPS6, eIF4E, Cadherin, and SMAD contributed to the reported phenotype in these individuals. Furthermore, it unveiled the potential therapeutic options that could be suggested to these individuals.

GatekeepR: an R Shiny application for the identification of nodes with high dynamic impact in Boolean networks

MOTIVATION

Boolean networks can serve as straightforward models for understanding processes such as gene regulation, and employing logical rules. These rules can either be derived from existing literature or by data-driven approaches. However, in the context of large networks, the exhaustive search for intervention targets becomes challenging due to the exponential expansion of a Boolean network's state space and the multitude of potential target candidates, along with their various combinations. Instead, we can employ the logical rules and resultant interaction graph as a means to identify targets of specific interest within larger-scale models. This approach not only facilitates the screening process but also serves as a preliminary filtering step, enabling the focused investigation of candidates that hold promise for more profound dynamic analysis. However, applying this method requires a working knowledge of R, thus restricting the range of potential users. We, therefore, aim to provide an application that makes this method accessible to a broader scientific community.

RESULTS

Here, we introduce GatekeepR, a graphical, web-based R Shiny application that enables scientists to screen Boolean network models for possible intervention targets whose perturbation is likely to have a large impact on the system's dynamics. This application does not require a local installation or knowledge of R and provides the suggested targets along with additional network information and visualizations in an intuitive, easy-to-use manner. The Supplementary Material describes the underlying method for identifying these nodes along with an example application in a network modeling pancreatic cancer.

AVAILABILITY AND IMPLEMENTATION

https://www.github.com/sysbio-bioinf/GatekeepR https://abel.informatik.uni-ulm.de/shiny/GatekeepR/.

A systems biology approach to define mechanisms, phenotypes, and drivers in PanNETs with a personalized perspective

Pancreatic neuroendocrine tumors (PanNETs) are a rare tumor entity with largely unpredictable progression and increasing incidence in developed countries. Molecular pathways involved in PanNETs development are still not elucidated, and specific biomarkers are missing. Moreover, the heterogeneity of PanNETs makes their treatment challenging and most approved targeted therapeutic options for PanNETs lack objective responses. Here, we applied a systems biology approach integrating dynamic modeling strategies, foreign classifier tailored approaches, and patient expression profiles to predict PanNETs progression as well as resistance mechanisms to clinically approved treatments such as the mammalian target of rapamycin complex 1 (mTORC1) inhibitors. We set up a model able to represent frequently reported PanNETs drivers in patient cohorts, such as Menin-1 (MEN1), Death domain associated protein (DAXX), Tuberous Sclerosis (TSC), as well as wild-type tumors. Model-based simulations suggested drivers of cancer progression as both first and second hits after MEN1 loss. In addition, we could predict the benefit of mTORC1 inhibitors on differentially mutated cohorts and hypothesize resistance mechanisms. Our approach sheds light on a more personalized prediction and treatment of PanNET mutant phenotypes.

Review and assessment of Boolean approaches for inference of gene regulatory networks

Boolean descriptions of gene regulatory networks can provide an insight into interactions between genes. Boolean networks hold predictive power, are easy to understand, and can be used to simulate the observed networks in different scenarios. We review fundamental and state-of-the-art methods for inference of Boolean networks. We introduce a methodology for a straightforward evaluation of Boolean inference approaches based on the generation of evaluation datasets, application of selected inference methods, and evaluation of performance measures to guide the selection of the best method for a given inference problem. We demonstrate this procedure on inference methods REVEAL (REVerse Engineering ALgorithm), Best-Fit Extension, MIBNI (Mutual Information-based Boolean Network Inference), GABNI (Genetic Algorithm-based Boolean Network Inference) and ATEN (AND/OR Tree ENsemble algorithm), which infers Boolean descriptions of gene regulatory networks from discretised time series data. Boolean inference approaches tend to perform better in terms of dynamic accuracy, and slightly worse in terms of structural correctness. We believe that the proposed methodology and provided guidelines will help researchers to develop Boolean inference approaches with a good predictive capability while maintaining structural correctness and biological relevance.

Discrete Logic Modeling of Cell Signaling Pathways

Cell signaling pathways often crosstalk generating complex biological behaviors observed in different cellular contexts. Frequently, laboratory experiments focus on a few putative regulators, alone unable to predict the molecular mechanisms behind the observed phenotypes. Here, systems biology complements these approaches by giving a holistic picture to complex signaling crosstalk. In particular, Boolean network models are a meaningful tool to study large network behaviors and can cope with incomplete kinetic information. By introducing a model describing pathways involved in hematopoietic stem cell maintenance, we present a general approach on how to model cell signaling pathways with Boolean network models.

Unraveling the Molecular Tumor-Promoting Regulation of Cofilin-1 in Pancreatic Cancer

Cofilin-1 (CFL1) overexpression in pancreatic cancer correlates with high invasiveness and shorter survival. Besides a well-documented role in actin remodeling, additional cellular functions of CFL1 remain poorly understood. Here, we unraveled molecular tumor-promoting functions of CFL1 in pancreatic cancer. For this purpose, we first show that a knockdown of CFL1 results in reduced growth and proliferation rates in vitro and in vivo, while apoptosis is not induced. By mechanistic modeling we were able to predict the underlying regulation. Model simulations indicate that an imbalance in actin remodeling induces overexpression and activation of CFL1 by acting on transcription factor 7-like 2 (TCF7L2) and aurora kinase A (AURKA). Moreover, we could predict that CFL1 impacts proliferation and apoptosis via the signal transducer and activator of transcription 3 (STAT3). These initial model-based regulations could be substantiated by studying protein levels in pancreatic cancer cell lines and human datasets. Finally, we identified the surface protein CD44 as a promising therapeutic target for pancreatic cancer patients with high CFL1 expression.

Awakening the HSC: Dynamic Modeling of HSC Maintenance Unravels Regulation of the TP53 Pathway and Quiescence

Hematopoietic stem cells (HSCs) provide all types of blood cells during the entire life of the organism. HSCs are mainly quiescent and can eventually enter the cell cycle to differentiate. HSCs are maintained and tightly regulated in a particular environment. The stem cell niche regulates dormancy and awakening. Deregulations of this interplay can lead to hematopoietic failure and diseases. In this paper, we present a Boolean network model that recapitulates HSC regulation in virtue of external signals coming from the niche. This Boolean network integrates and summarizes the current knowledge of HSC regulation and is based on extensive literature research. Furthermore, dynamic simulations suggest a novel systemic regulation of TP53 in homeostasis. Thereby, our model indicates that TP53 activity is balanced depending on external stimulations, engaging a regulatory mechanism involving ROS regulators and RAS activated transcription factors. Finally, we investigated different mouse models and compared them to in silico knockout simulations. Here, the model could recapitulate in vivo observed behaviors and thus sustains our results.

Concepts in Boolean network modeling: What do they all mean?

Boolean network models are one of the simplest models to study complex dynamic behavior in biological systems. They can be applied to unravel the mechanisms regulating the properties of the system or to identify promising intervention targets. Since its introduction by Stuart Kauffman in 1969 for describing gene regulatory networks, various biologically based networks and tools for their analysis were developed. Here, we summarize and explain the concepts for Boolean network modeling. We also present application examples and guidelines to work with and analyze Boolean network models.

Representing dynamic biological networks with multi-scale probabilistic models

Dynamic models analyzing gene regulation and metabolism face challenges when adapted to modeling signal transduction networks. During signal transduction, molecular reactions and mechanisms occur in different spatial and temporal frames and involve feedbacks. This impedes the straight-forward use of methods based on Boolean networks, Bayesian approaches, and differential equations. We propose a new approach, ProbRules, that combines probabilities and logical rules to represent the dynamics of a system across multiple scales. We demonstrate that ProbRules models can represent various network motifs of biological systems. As an example of a comprehensive model of signal transduction, we provide a Wnt network that shows remarkable robustness under a range of phenotypical and pathological conditions. Its simulation allows the clarification of controversially discussed molecular mechanisms of Wnt signaling by predicting wet-lab measurements. ProbRules provides an avenue in current computational modeling by enabling systems biologists to integrate vast amounts of available data on different scales.

Automatic Screening for Perturbations in Boolean Networks

A common approach to address biological questions in systems biology is to simulate regulatory mechanisms using dynamic models. Among others, Boolean networks can be used to model the dynamics of regulatory processes in biology. Boolean network models allow simulating the qualitative behavior of the modeled processes. A central objective in the simulation of Boolean networks is the computation of their long-term behavior—so-called attractors. These attractors are of special interest as they can often be linked to biologically relevant behaviors. Changing internal and external conditions can influence the long-term behavior of the Boolean network model. Perturbation of a Boolean network by stripping a component of the system or simulating a surplus of another element can lead to different attractors. Apparently, the number of possible perturbations and combinations of perturbations increases exponentially with the size of the network. Manually screening a set of possible components for combinations that have a desired effect on the long-term behavior can be very time consuming if not impossible. We developed a method to automatically screen for perturbations that lead to a user-specified change in the network's functioning. This method is implemented in the visual simulation framework ViSiBool utilizing satisfiability (SAT) solvers for fast exhaustive attractor search.

A Boolean network of the crosstalk between IGF and Wnt signaling in aging satellite cells

Aging is a complex biological process, which determines the life span of an organism. Insulin-like growth factor (IGF) and Wnt signaling pathways govern the process of aging. Both pathways share common downstream targets that allow competitive crosstalk between these branches. Of note, a shift from IGF to Wnt signaling has been observed during aging of satellite cells. Biological regulatory networks necessary to recreate aging have not yet been discovered. Here, we established a mathematical in silico model that robustly recapitulates the crosstalk between IGF and Wnt signaling. Strikingly, it predicts critical nodes following a shift from IGF to Wnt signaling. These findings indicate that this shift might cause age-related diseases.

ATLANTIS - Attractor Landscape Analysis Toolbox for Cell Fate Discovery and Reprogramming

Boolean modelling of biological networks is a well-established technique for abstracting dynamical biomolecular regulation in cells. Specifically, decoding linkages between salient regulatory network states and corresponding cell fate outcomes can help uncover pathological foundations of diseases such as cancer. Attractor landscape analysis is one such methodology which converts complex network behavior into a landscape of network states wherein each state is represented by propensity of its occurrence. Towards undertaking attractor landscape analysis of Boolean networks, we propose an Attractor Landscape Analysis Toolbox (ATLANTIS) for cell fate discovery, from biomolecular networks, and reprogramming upon network perturbation. ATLANTIS can be employed to perform both deterministic and probabilistic analyses. It has been validated by successfully reconstructing attractor landscapes from several published case studies followed by reprogramming of cell fates upon therapeutic treatment of network. Additionally, the biomolecular network of HCT-116 colorectal cancer cell line has been screened for therapeutic evaluation of drug-targets. Our results show agreement between therapeutic efficacies reported by ATLANTIS and the published literature. These case studies sufficiently highlight the in silico cell fate prediction and therapeutic screening potential of the toolbox. Lastly, ATLANTIS can also help guide single or combinatorial therapy responses towards reprogramming biomolecular networks to recover cell fates.

Stability of Signaling Pathways during Aging-A Boolean Network Approach

Biological pathways are thought to be robust against a variety of internal and external perturbations. Fail-safe mechanisms allow for compensation of perturbations to maintain the characteristic function of a pathway. Pathways can undergo changes during aging, which may lead to changes in their stability. Less stable or less robust pathways may be consequential to or increase the susceptibility of the development of diseases. Among others, NF- κ B signaling is a crucial pathway in the process of aging. The NF- κ B system is involved in the immune response and dealing with various internal and external stresses. Boolean networks as models of biological pathways allow for simulation of signaling behavior. They can help to identify which proposed mechanisms are biologically representative and which ones function but do not mirror physical processes-for instance, changes of signaling pathways during the aging process. Boolean networks can be inferred from time-series of gene expression data. This allows us to get insights into the changes of behavior of pathways such as NF- κ B signaling in aged organisms in comparison to young ones.

Simulation and visualization of multiple KEGG pathways using BioNSi

Motivation: Many biologists are discouraged from using network simulation tools because these require manual, often tedious network construction. This situation calls for building new tools or extending existing ones with the ability to import biological pathways previously deposited in databases and analyze them, in order to produce novel biological insights at the pathway level.

Results: We have extended a network simulation tool (BioNSi), which now allows merging of multiple pathways from the KEGG pathway database into a single, coherent network, and visualizing its properties. Furthermore, the enhanced tool enables loading experimental expression data into the network and simulating its dynamics under various biological conditions or perturbations. As a proof of concept, we tested two sets of published experimental data, one related to inflammatory bowel disease condition and the other to breast cancer treatment. We predict some of the major observations obtained following these laboratory experiments, and provide new insights that may shed additional light on these results.

Tool requirements: Cytoscape 3.x, JAVA 8

Availability: The tool is freely available at http://bionsi.wix.com/bionsi, where a complete user guide and a step-by-step manual can also be found.