Parameter estimation of qualitative biological regulatory networks on high performance computing hardware

Abstract

BackgroundBiological Regulatory Networks (BRNs) are responsible for developmental and maintenance related functions in organisms. These functions are implemented by the dynamics of BRNs and are sensitive to regulations enforced by specific activators and inhibitors. The logical modeling formalism by René Thomas incorporates this sensitivity with a set of logical parameters modulated by available regulators, varying with time. With the increase in complexity of BRNs in terms of number of entities and their interactions, the task of parameters estimation becomes computationally expensive with existing sequential SMBioNET tool. We extend the existing sequential implementation of SMBioNET by using a data decomposition approach using a Java messaging library called MPJ Express. The approach divides the parameters space into different regions and each region is then explored in parallel on High Performance Computing (HPC) hardware.ResultsThe performance of the parallel approach is evaluated on BRNs of different sizes, and experimental results on multicore and cluster computers showed almost linear speed-up. This parallel code can be executed on a wide range of concurrent hardware including laptops equipped with multicore processors, and specialized distributed memory computer systems. To demonstrate the application of parallel implementation, we selected a case study of Hexosamine Biosynthetic Pathway (HBP) in cancer progression to identify potential therapeutic targets against cancer. A set of logical parameters were computed for HBP model that directs the biological system to a state of recovery. Furthermore, the parameters also suggest a potential therapeutic intervention that restores homeostasis. Additionally, the performance of parallel application was also evaluated on a network (comprising of 23 entities) of Fibroblast Growth Factor Signalling in Drosophila melanogaster.ConclusionsQualitative modeling framework is widely used for investigating dynamics of biological regulatory networks. However, computation of model parameters in qualitative modeling is computationally intensive. In this work, we presented results of our Java based parallel implementation that provides almost linear speed-up on both multicore and cluster platforms. The parallel implementation is available at https://psmbionet.github.io.Electronic supplementary materialThe online version of this article (10.1186/s12918-018-0670-y) contains supplementary material, which is available to authorized users.

Decoding the Role of Epigenetics in Breast Cancer Using Formal Modeling and Machine-Learning Methods

“…This study provides an insight into the integrated DNMT1-RUNX3 signaling by taking into account various cancer-related significant upstream and downstream regulators (such as p21, c-myc, p53, and MDM2) and presents the DNMT1-RUNX3 signaling cascades as one consolidated network. To the best of authors' knowledge, this study is one of its kind to illustrate and model the epigenetic-mediated silencing of TSG RUNX3 through the René Thomas framework modeling that has largely been a common practice of systems biology to investigate the dynamics of biological networks (Saeed et al, 2018). Moreover, the model checking technique was adopted to build the interaction graph rendered as a state graph (Figure 6B) by utilizing the literature-driven information of the DNMT1-RUNX3 signaling (Figure 2).…”

“…and the relationship (activation and inhibition) among the chosen entities. The unknown parameters were inferred by encoding wet-laboratory biological observations as propositional calculus or more precisely computation tree logic (CTL) and verified through the model checking technique as previously reported by Saeed et al (2018). Briefly, model checking is an automatic technique based on the exhaustive exploration of the entire state space of a biological system, which therefore allows the analysis and cross-verification of a large number of possible outcomes of a network (Bernot et al, 2004).…”

Deciphering the expression dynamics of ANGPTL8 associated regulatory network in insulin resistance using formal modelling approaches

“…This method for estimation of logical parameters of different BRNs has been used previously in [ 45 , 46 , 51 , 52 , 53 ]. Moreover, we refer our readers to [ 42 , 56 , 57 , 58 ] for a comprehensive literature review of the methodology used for logical parameters estimation.…”

“…2, pp. [47][48][49][50][51][52][53][54][55][56][57][58] This is an open access article published by the IET under the Creative Commons Attribution-NoDerivs License (http://creativecommons.org/licenses/by-nd/3.0/) ANGPTL8 during insulin resistance. In their study, overexpression of ANGPTL8 in mice liver (its normal site of expression) exhibited two responses, i.e.…”

Inferring cancer common and specific gene networks via multi-layer joint graphical model

“…Over the past decades, various computational approaches have been proposed for inferring gene networks from microarrays or bulk RNA-seq data, such as boolean models [4] , information theory-based models [5] , differential equation-based models [6] and Gaussian graphical models [7] . Among these methods, Gaussian graphical models (GGMs) are popular due to their ability in predicting the direct interactions between genes.…”