1
|
Singh N, Khan FM, Bala L, Vera J, Wolkenhauer O, Pützer B, Logotheti S, Gupta SK. Logic-based modeling and drug repurposing for the prediction of novel therapeutic targets and combination regimens against E2F1-driven melanoma progression. BMC Chem 2023; 17:161. [PMID: 37993971 PMCID: PMC10666365 DOI: 10.1186/s13065-023-01082-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2023] [Accepted: 11/08/2023] [Indexed: 11/24/2023] Open
Abstract
Melanoma presents increasing prevalence and poor outcomes. Progression to aggressive stages is characterized by overexpression of the transcription factor E2F1 and activation of downstream prometastatic gene regulatory networks (GRNs). Appropriate therapeutic manipulation of the E2F1-governed GRNs holds the potential to prevent metastasis however, these networks entail complex feedback and feedforward regulatory motifs among various regulatory layers, which make it difficult to identify druggable components. To this end, computational approaches such as mathematical modeling and virtual screening are important tools to unveil the dynamics of these signaling networks and identify critical components that could be further explored as therapeutic targets. Herein, we integrated a well-established E2F1-mediated epithelial-mesenchymal transition (EMT) map with transcriptomics data from E2F1-expressing melanoma cells to reconstruct a core regulatory network underlying aggressive melanoma. Using logic-based in silico perturbation experiments of a core regulatory network, we identified that simultaneous perturbation of Protein kinase B (AKT1) and oncoprotein murine double minute 2 (MDM2) drastically reduces EMT in melanoma. Using the structures of the two protein signatures, virtual screening strategies were performed with the FDA-approved drug library. Furthermore, by combining drug repurposing and computer-aided drug design techniques, followed by molecular dynamics simulation analysis, we identified two potent drugs (Tadalafil and Finasteride) that can efficiently inhibit AKT1 and MDM2 proteins. We propose that these two drugs could be considered for the development of therapeutic strategies for the management of aggressive melanoma.
Collapse
Affiliation(s)
- Nivedita Singh
- Department of Biochemistry, BBDCODS, BBD University, Lucknow, Uttar Pradesh, India
- MRC Laboratory for Molecular Cell Biology, University College London, London, UK
| | - Faiz M Khan
- Department of Systems Biology and Bioinformatics, University of Rostock, Rostock, Germany
| | - Lakshmi Bala
- Department of Biochemistry, BBDCODS, BBD University, Lucknow, Uttar Pradesh, India
| | - Julio Vera
- Department of Dermatology, Universitätsklinikum Erlangen and Friedrich-Alexander Universität Erlangen-Nürnberg (FAU), Erlangen, Germany
- Comprehensive Cancer Center Erlangen-European Metropolitan Area of Nuremberg (CCC ER-EMN), Erlangen, Germany
- Deutsches Zentrum Immuntherapie (DZI), Erlangen, Germany
| | - Olaf Wolkenhauer
- Department of Systems Biology and Bioinformatics, University of Rostock, Rostock, Germany
- Leibniz Institute for Food Systems Biology, Technical University of Munich, Munich, Germany
- Chhattisgarh Swami Vivekanand Technical University, Bhilai, Chhattisgarh, India
- Stellenbosch Institute of Advanced Study, Wallenberg Research Centre, Stellenbosch University, Stellenbosch, South Africa
| | - Brigitte Pützer
- Institute of Experimental Gene Therapy and Cancer Research, Rostock University Medical Center, Rostock, Germany
| | - Stella Logotheti
- Institute of Experimental Gene Therapy and Cancer Research, Rostock University Medical Center, Rostock, Germany
- DNA Damage Laboratory, Physics Department, School of Applied Mathematical and Physical Sciences, National Technical University of Athens (NTUA), Zografou, Athens, Greece
| | - Shailendra K Gupta
- Department of Systems Biology and Bioinformatics, University of Rostock, Rostock, Germany.
- Chhattisgarh Swami Vivekanand Technical University, Bhilai, Chhattisgarh, India.
| |
Collapse
|
2
|
Parmer T, Radicchi F. Dynamical methods for target control of biological networks. ROYAL SOCIETY OPEN SCIENCE 2023; 10:230542. [PMID: 37920567 PMCID: PMC10618059 DOI: 10.1098/rsos.230542] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/01/2023] [Accepted: 09/29/2023] [Indexed: 11/04/2023]
Abstract
Estimating the influence that individual nodes have on one another in a Boolean network is essential to predict and control the system's dynamical behaviour, for example, detecting key therapeutic targets to control pathways in models of biological signalling and regulation. Exact estimation is generally not possible due to the fact that the number of configurations that must be considered grows exponentially with the system size. However, approximate, scalable methods exist in the literature. These methods can be divided into two main classes: (i) graph-theoretic methods that rely on representations of Boolean dynamics into static graphs and (ii) mean-field approaches that describe average trajectories of the system but neglect dynamical correlations. Here, we compare systematically the performance of these state-of-the-art methods on a large collection of real-world gene regulatory networks. We find comparable performance across methods. All methods underestimate the ground truth, with mean-field approaches having a better recall but a worse precision than graph-theoretic methods. Computationally speaking, graph-theoretic methods are faster than mean-field ones in sparse networks, but are slower in dense networks. The preference of which method to use, therefore, depends on a network's connectivity and the relative importance of recall versus precision for the specific application at hand.
Collapse
Affiliation(s)
- Thomas Parmer
- Center for Complex Networks and Systems Research, Luddy School of Informatics, Computing, and Engineering, Indiana University, Bloomington, IN 47408, USA
| | - Filippo Radicchi
- Center for Complex Networks and Systems Research, Luddy School of Informatics, Computing, and Engineering, Indiana University, Bloomington, IN 47408, USA
| |
Collapse
|
3
|
Marku M, Pancaldi V. From time-series transcriptomics to gene regulatory networks: A review on inference methods. PLoS Comput Biol 2023; 19:e1011254. [PMID: 37561790 PMCID: PMC10414591 DOI: 10.1371/journal.pcbi.1011254] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/12/2023] Open
Abstract
Inference of gene regulatory networks has been an active area of research for around 20 years, leading to the development of sophisticated inference algorithms based on a variety of assumptions and approaches. With the ever increasing demand for more accurate and powerful models, the inference problem remains of broad scientific interest. The abstract representation of biological systems through gene regulatory networks represents a powerful method to study such systems, encoding different amounts and types of information. In this review, we summarize the different types of inference algorithms specifically based on time-series transcriptomics, giving an overview of the main applications of gene regulatory networks in computational biology. This review is intended to give an updated reference of regulatory networks inference tools to biologists and researchers new to the topic and guide them in selecting the appropriate inference method that best fits their questions, aims, and experimental data.
Collapse
Affiliation(s)
- Malvina Marku
- CRCT, Université de Toulouse, Inserm, CNRS, Université Toulouse III-Paul Sabatier, Centre de Recherches en Cancérologie de Toulouse, Toulouse, France
| | - Vera Pancaldi
- CRCT, Université de Toulouse, Inserm, CNRS, Université Toulouse III-Paul Sabatier, Centre de Recherches en Cancérologie de Toulouse, Toulouse, France
- Barcelona Supercomputing Center, Barcelona, Spain
| |
Collapse
|
4
|
Argyris GA, Lluch Lafuente A, Tribastone M, Tschaikowski M, Vandin A. Reducing Boolean networks with backward equivalence. BMC Bioinformatics 2023; 24:212. [PMID: 37221494 DOI: 10.1186/s12859-023-05326-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2022] [Accepted: 05/05/2023] [Indexed: 05/25/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Georgios A Argyris
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Lyngby, Denmark
| | - Alberto Lluch Lafuente
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Lyngby, Denmark
| | | | - Max Tschaikowski
- Department of Computer Science, University of Aalborg, Aalborg, Denmark
| | - Andrea Vandin
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Lyngby, Denmark.
- Department of Excellence EMbeDS and Institute of Economics, Sant'Anna School for Advanced Studies, Pisa, Italy.
| |
Collapse
|
5
|
Schnitzer B, Österberg L, Skopa I, Cvijovic M. Multi-scale model suggests the trade-off between protein and ATP demand as a driver of metabolic changes during yeast replicative ageing. PLoS Comput Biol 2022; 18:e1010261. [PMID: 35797415 PMCID: PMC9295998 DOI: 10.1371/journal.pcbi.1010261] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2022] [Revised: 07/19/2022] [Accepted: 05/31/2022] [Indexed: 11/18/2022] Open
Abstract
The accumulation of protein damage is one of the major drivers of replicative ageing, describing a cell’s reduced ability to reproduce over time even under optimal conditions. Reactive oxygen and nitrogen species are precursors of protein damage and therefore tightly linked to ageing. At the same time, they are an inevitable by-product of the cell’s metabolism. Cells are able to sense high levels of reactive oxygen and nitrogen species and can subsequently adapt their metabolism through gene regulation to slow down damage accumulation. However, the older or damaged a cell is the less flexibility it has to allocate enzymes across the metabolic network, forcing further adaptions in the metabolism. To investigate changes in the metabolism during replicative ageing, we developed an multi-scale mathematical model using budding yeast as a model organism. The model consists of three interconnected modules: a Boolean model of the signalling network, an enzyme-constrained flux balance model of the central carbon metabolism and a dynamic model of growth and protein damage accumulation with discrete cell divisions. The model can explain known features of replicative ageing, like average lifespan and increase in generation time during successive division, in yeast wildtype cells by a decreasing pool of functional enzymes and an increasing energy demand for maintenance. We further used the model to identify three consecutive metabolic phases, that a cell can undergo during its life, and their influence on the replicative potential, and proposed an intervention span for lifespan control.
Collapse
Affiliation(s)
- Barbara Schnitzer
- Department of Mathematical Sciences, Chalmers University of Technology, Gothenburg, Sweden
- Department of Mathematical Sciences, University of Gothenburg, Gothenburg, Sweden
| | - Linnea Österberg
- Department of Mathematical Sciences, Chalmers University of Technology, Gothenburg, Sweden
- Department of Mathematical Sciences, University of Gothenburg, Gothenburg, Sweden
- Department of Biology and Biological Engineering, Chalmers University of Technology, Gothenburg, Sweden
| | - Iro Skopa
- Department of Mathematical Sciences, Chalmers University of Technology, Gothenburg, Sweden
- Department of Mathematical Sciences, University of Gothenburg, Gothenburg, Sweden
| | - Marija Cvijovic
- Department of Mathematical Sciences, Chalmers University of Technology, Gothenburg, Sweden
- Department of Mathematical Sciences, University of Gothenburg, Gothenburg, Sweden
- * E-mail:
| |
Collapse
|
6
|
Werle SD, Ikonomi N, Schwab JD, Kraus JM, Weidner FM, Lenhard Rudolph K, Pfister AS, Schuler R, Kühl M, Kestler HA. Identification of dynamic driver sets controlling phenotypical landscapes. Comput Struct Biotechnol J 2022; 20:1603-1617. [PMID: 35465155 PMCID: PMC9010550 DOI: 10.1016/j.csbj.2022.03.034] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2022] [Revised: 03/30/2022] [Accepted: 03/30/2022] [Indexed: 11/03/2022] Open
Abstract
Controlling phenotypical landscapes is of vital interest to modern biology. This task becomes highly demanding because cellular decisions involve complex networks engaging in crosstalk interactions. Previous work on control theory indicates that small sets of compounds can control single phenotypes. However, a dynamic approach is missing to determine the drivers of the whole network dynamics. By analyzing 35 biologically motivated Boolean networks, we developed a method to identify small sets of compounds sufficient to decide on the entire phenotypical landscape. These compounds do not strictly prefer highly related compounds and show a smaller impact on the stability of the attractor landscape. The dynamic driver sets include many intervention targets and cellular reprogramming drivers in human networks. Finally, by using a new comprehensive model of colorectal cancer, we provide a complete workflow on how to implement our approach to shift from in silico to in vitro guided experiments.
Collapse
|
7
|
Taou N, Lones M. Optimising Boolean Synthetic Regulatory Networks to Control Cell States. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2021; 18:2649-2658. [PMID: 32078555 DOI: 10.1109/tcbb.2020.2973636] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Controlling the dynamics of gene regulatory networks is a challenging problem. In recent years, a number of control methods have been proposed, but most of these approaches do not address the problem of how they could be implemented in practice. In this paper, we consider the idea of using a synthetic regulatory network as a closed-loop controller that can control and respond to the dynamics of a cell's native regulatory network in situ. We explore this idea using a computational model in which both native and synthetic regulatory networks are represented by Boolean networks. We then use an evolutionary algorithm to optimise both the structure and parameters of the synthetic Boolean network. To test this approach, we look at whether controllers can be optimised to target specific steady states in five different Boolean regulatory circuit models. Our results show that in most cases the controllers are able to drive the dynamics of the target system to a specified steady state, often using few interventions, and further experiments using random Boolean networks show that the approach scales well to larger controlled networks.
Collapse
|
8
|
Weidner FM, Schwab JD, Werle SD, Ikonomi N, Lausser L, Kestler HA. Capturing dynamic relevance in Boolean networks using graph theoretical measures. Bioinformatics 2021; 37:3530-3537. [PMID: 33983406 PMCID: PMC8545349 DOI: 10.1093/bioinformatics/btab277] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2020] [Revised: 03/19/2021] [Accepted: 04/22/2021] [Indexed: 11/14/2022] Open
Abstract
Motivation Interaction graphs are able to describe regulatory dependencies between compounds without capturing dynamics. In contrast, mathematical models that are based on interaction graphs allow to investigate the dynamics of biological systems. However, since dynamic complexity of these models grows exponentially with their size, exhaustive analyses of the dynamics and consequently screening all possible interventions eventually becomes infeasible. Thus, we designed an approach to identify dynamically relevant compounds based on the static network topology. Results Here, we present a method only based on static properties to identify dynamically influencing nodes. Coupling vertex betweenness and determinative power, we could capture relevant nodes for changing dynamics with an accuracy of 75% in a set of 35 published logical models. Further analyses of the selected compounds’ connectivity unravelled a new class of not highly connected nodes with high impact on the networks’ dynamics, which we call gatekeepers. We validated our method’s working concept on logical models, which can be readily scaled up to complex interaction networks, where dynamic analyses are not even feasible. Availability and implementation Code is freely available at https://github.com/sysbio-bioinf/BNStatic. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Felix M Weidner
- Institute of Medical Systems Biology, Ulm University, Germany.,International Graduate School of Molecular Medicine, Ulm University, Germany
| | - Julian D Schwab
- Institute of Medical Systems Biology, Ulm University, Germany
| | - Silke D Werle
- Institute of Medical Systems Biology, Ulm University, Germany.,International Graduate School of Molecular Medicine, Ulm University, Germany
| | - Nensi Ikonomi
- Institute of Medical Systems Biology, Ulm University, Germany.,International Graduate School of Molecular Medicine, Ulm University, Germany
| | - Ludwig Lausser
- Institute of Medical Systems Biology, Ulm University, Germany
| | - Hans A Kestler
- Institute of Medical Systems Biology, Ulm University, Germany
| |
Collapse
|
9
|
Schwob MR, Zhan J, Dempsey A. Modeling Cell Communication with Time-Dependent Signaling Hypergraphs. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2021; 18:1151-1163. [PMID: 31449029 DOI: 10.1109/tcbb.2019.2937033] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Signaling pathways describe a group of molecules in a cell that collaborate to control one or more cell functions, such as cell division or cell death. The pathways communicate by sending signals between molecules, and this process is repeated until the terminal molecule is activated and the cell function is executed. Signaling pathways are often represented as directed graphs, which does not provide enough information when modeling cell functions and reactions. Recently, directed hypergraphs have been proposed to more accurately represent reactions such as protein activation and interaction. To further improve the representation of signaling pathways, time dependency must be considered to improve the representation of cell signaling at any given time. In this paper, the importance of time dependency in modeling signaling pathways is presented. An algorithm that finds the shortest a priori path using time-dependent hypergraphs to more robustly model signaling pathways is adopted. The shortest time-dependent hyperpaths representing signaling pathways are an improvement to the recent adoption of hypergraphs representing these pathways. The results display the improved representation of signaling pathways and motivate the adoption of time-dependent signaling hypergraphs.
Collapse
|
10
|
Gates AJ, Brattig Correia R, Wang X, Rocha LM. The effective graph reveals redundancy, canalization, and control pathways in biochemical regulation and signaling. Proc Natl Acad Sci U S A 2021; 118:e2022598118. [PMID: 33737396 PMCID: PMC8000424 DOI: 10.1073/pnas.2022598118] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
The ability to map causal interactions underlying genetic control and cellular signaling has led to increasingly accurate models of the complex biochemical networks that regulate cellular function. These network models provide deep insights into the organization, dynamics, and function of biochemical systems: for example, by revealing genetic control pathways involved in disease. However, the traditional representation of biochemical networks as binary interaction graphs fails to accurately represent an important dynamical feature of these multivariate systems: some pathways propagate control signals much more effectively than do others. Such heterogeneity of interactions reflects canalization-the system is robust to dynamical interventions in redundant pathways but responsive to interventions in effective pathways. Here, we introduce the effective graph, a weighted graph that captures the nonlinear logical redundancy present in biochemical network regulation, signaling, and control. Using 78 experimentally validated models derived from systems biology, we demonstrate that 1) redundant pathways are prevalent in biological models of biochemical regulation, 2) the effective graph provides a probabilistic but precise characterization of multivariate dynamics in a causal graph form, and 3) the effective graph provides an accurate explanation of how dynamical perturbation and control signals, such as those induced by cancer drug therapies, propagate in biochemical pathways. Overall, our results indicate that the effective graph provides an enriched description of the structure and dynamics of networked multivariate causal interactions. We demonstrate that it improves explainability, prediction, and control of complex dynamical systems in general and biochemical regulation in particular.
Collapse
Affiliation(s)
- Alexander J Gates
- Network Science Institute, Northeastern University, Boston, MA 02115;
| | - Rion Brattig Correia
- Instituto Gulbenkian de Ciência, 2780-156 Oeiras, Portugal
- Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Ministry of Education of Brazil, 70040-020 Brasília, DF, Brazil
| | - Xuan Wang
- Center for Social and Biomedical Complexity, Luddy School of Informatics, Computing & Engineering, Indiana University, Bloomington, IN 47408
| | - Luis M Rocha
- Instituto Gulbenkian de Ciência, 2780-156 Oeiras, Portugal;
- Center for Social and Biomedical Complexity, Luddy School of Informatics, Computing & Engineering, Indiana University, Bloomington, IN 47408
- Department of Systems Science and Industrial Engineering, Binghamton University, Binghamton, NY 13902
| |
Collapse
|
11
|
Baumgartner L, Wuertz-Kozak K, Le Maitre CL, Wignall F, Richardson SM, Hoyland J, Ruiz Wills C, González Ballester MA, Neidlin M, Alexopoulos LG, Noailly J. Multiscale Regulation of the Intervertebral Disc: Achievements in Experimental, In Silico, and Regenerative Research. Int J Mol Sci 2021; 22:E703. [PMID: 33445782 PMCID: PMC7828304 DOI: 10.3390/ijms22020703] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Revised: 12/22/2020] [Accepted: 12/24/2020] [Indexed: 12/17/2022] Open
Abstract
Intervertebral disc (IVD) degeneration is a major risk factor of low back pain. It is defined by a progressive loss of the IVD structure and functionality, leading to severe impairments with restricted treatment options due to the highly demanding mechanical exposure of the IVD. Degenerative changes in the IVD usually increase with age but at an accelerated rate in some individuals. To understand the initiation and progression of this disease, it is crucial to identify key top-down and bottom-up regulations' processes, across the cell, tissue, and organ levels, in health and disease. Owing to unremitting investigation of experimental research, the comprehension of detailed cell signaling pathways and their effect on matrix turnover significantly rose. Likewise, in silico research substantially contributed to a holistic understanding of spatiotemporal effects and complex, multifactorial interactions within the IVD. Together with important achievements in the research of biomaterials, manifold promising approaches for regenerative treatment options were presented over the last years. This review provides an integrative analysis of the current knowledge about (1) the multiscale function and regulation of the IVD in health and disease, (2) the possible regenerative strategies, and (3) the in silico models that shall eventually support the development of advanced therapies.
Collapse
Affiliation(s)
- Laura Baumgartner
- BCN MedTech, Department of Information and Communication Technologies, Universitat Pompeu Fabra, 08018 Barcelona, Spain; (L.B.); (C.R.W.); (M.A.G.B.)
| | - Karin Wuertz-Kozak
- Department of Biomedical Engineering, Rochester Institute of Technology (RIT), Rochester, NY 14623, USA;
- Schön Clinic Munich Harlaching, Spine Center, Academic Teaching Hospital and Spine Research Institute of the Paracelsus Medical University Salzburg (Austria), 81547 Munich, Germany
| | - Christine L. Le Maitre
- Biomolecular Sciences Research Centre, Sheffield Hallam University, Sheffield S1 1WB, UK;
| | - Francis Wignall
- Division of Cell Matrix Biology and Regenerative Medicine, School of Biological Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester Academic Health Sciences Centre, Oxford Road, Manchester M13 9PT, UK; (F.W.); (S.M.R.); (J.H.)
| | - Stephen M. Richardson
- Division of Cell Matrix Biology and Regenerative Medicine, School of Biological Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester Academic Health Sciences Centre, Oxford Road, Manchester M13 9PT, UK; (F.W.); (S.M.R.); (J.H.)
| | - Judith Hoyland
- Division of Cell Matrix Biology and Regenerative Medicine, School of Biological Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester Academic Health Sciences Centre, Oxford Road, Manchester M13 9PT, UK; (F.W.); (S.M.R.); (J.H.)
| | - Carlos Ruiz Wills
- BCN MedTech, Department of Information and Communication Technologies, Universitat Pompeu Fabra, 08018 Barcelona, Spain; (L.B.); (C.R.W.); (M.A.G.B.)
| | - Miguel A. González Ballester
- BCN MedTech, Department of Information and Communication Technologies, Universitat Pompeu Fabra, 08018 Barcelona, Spain; (L.B.); (C.R.W.); (M.A.G.B.)
- Catalan Institution for Research and Advanced Studies (ICREA), Pg. Lluis Companys 23, 08010 Barcelona, Spain
| | - Michael Neidlin
- Department of Mechanical Engineering, National Technical University of Athens, 15780 Athens, Greece; (M.N.); (L.G.A.)
| | - Leonidas G. Alexopoulos
- Department of Mechanical Engineering, National Technical University of Athens, 15780 Athens, Greece; (M.N.); (L.G.A.)
| | - Jérôme Noailly
- BCN MedTech, Department of Information and Communication Technologies, Universitat Pompeu Fabra, 08018 Barcelona, Spain; (L.B.); (C.R.W.); (M.A.G.B.)
| |
Collapse
|
12
|
Groote JF, Larsen KG. Symbolic Coloured SCC Decomposition. TOOLS AND ALGORITHMS FOR THE CONSTRUCTION AND ANALYSIS OF SYSTEMS 2021. [PMCID: PMC7984532 DOI: 10.1007/978-3-030-72013-1_4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Problems arising in many scientific disciplines are often modelled using edge-coloured directed graphs. These can be enormous in the number of both vertices and colours. Given such a graph, the original problem frequently translates to the detection of the graph’s strongly connected components, which is challenging at this scale. We propose a new, symbolic algorithm that computes all the monochromatic strongly connected components of an edge-coloured graph. In the worst case, the algorithm performs \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$O(p\cdot n\cdot \log n)$$\end{document}O(p·n·logn) symbolic steps, where p is the number of colours and n the number of vertices. We evaluate the algorithm using an experimental implementation based on Binary Decision Diagrams (BDDs) and large (up to \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2^{48}$$\end{document}248) coloured graphs produced by models appearing in systems biology.
Collapse
|
13
|
Abstract
Motivation A major challenge in molecular and cellular biology is to map out the regulatory networks of cells. As regulatory interactions can typically not be directly observed experimentally, various computational methods have been proposed to disentangling direct and indirect effects. Most of these rely on assumptions that are rarely met or cannot be adapted to a given context. Results We present a network inference method that is based on a simple response logic with minimal presumptions. It requires that we can experimentally observe whether or not some of the system’s components respond to perturbations of some other components, and then identifies the directed networks that most accurately account for the observed propagation of the signal. To cope with the intractable number of possible networks, we developed a logic programming approach that can infer networks of hundreds of nodes, while being robust to noisy, heterogeneous or missing data. This allows to directly integrate prior network knowledge and additional constraints such as sparsity. We systematically benchmark our method on KEGG pathways, and show that it outperforms existing approaches in DREAM3 and DREAM4 challenges. Applied to a novel perturbation dataset on PI3K and MAPK pathways in isogenic models of a colon cancer cell line, it generates plausible network hypotheses that explain distinct sensitivities toward various targeted inhibitors due to different PI3K mutants. Availability and implementation A Python/Answer Set Programming implementation can be accessed at github.com/GrossTor/response-logic. Data and analysis scripts are available at github.com/GrossTor/response-logic-projects. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Torsten Gross
- Institut für Pathologie, Charité-Universitätsmedizin, Berlin.,IRI Life Sciences, Humboldt Universität zu Berlin, Berlin.,Berlin Institute of Health, Berlin, Germany
| | | | - Yibing Yan
- Oncology Biomarker Development, Genentech Inc., South San Francisco, CA, USA
| | - Nils Blüthgen
- Institut für Pathologie, Charité-Universitätsmedizin, Berlin.,IRI Life Sciences, Humboldt Universität zu Berlin, Berlin.,Berlin Institute of Health, Berlin, Germany
| |
Collapse
|
14
|
Martínez-Méndez D, Villarreal C, Mendoza L, Huerta L. An Integrative Network Modeling Approach to T CD4 Cell Activation. Front Physiol 2020; 11:380. [PMID: 32425809 PMCID: PMC7212416 DOI: 10.3389/fphys.2020.00380] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2019] [Accepted: 03/30/2020] [Indexed: 12/15/2022] Open
Abstract
The adaptive immune response is initiated by the interaction of the T cell antigen receptor/CD3 complex (TCR) with a cognate peptide bound to a MHC molecule. This interaction, along with the activity of co-stimulatory molecules and cytokines in the microenvironment, enables cells to proliferate and produce soluble factors that stimulate other branches of the immune response for inactivation of infectious agents. The intracellular activation signals are reinforced, amplified and diversified by a complex network of biochemical interactions, and includes the activity of molecules that modulate the activation process and stimulate the metabolic changes necessary for fulfilling the cell energy demands. We present an approach to the analysis of the main early signaling events of T cell activation by proposing a concise 46-node hybrid Boolean model of the main steps of TCR and CD28 downstream signaling, encompassing the activity of the anergy factor Ndrg1, modulation of activation by CTLA-4, and the activity of the nutrient sensor AMPK as intrinsic players of the activation process. The model generates stable states that reflect the overcoming of activation signals and induction of anergy by the expression of Ndrg1 in the absence of co-stimulation. The model also includes the induction of CTLA-4 upon activation and its competition with CD28 for binding to the co-stimulatory CD80/86 molecules, leading to stable states that reflect the activation arrest. Furthermore, the model integrates the activity of AMPK to the general pathways driving differentiation to functional cell subsets (Th1, Th2, Th17, and Treg). Thus, the network topology incorporates basic mechanism associated to activation, regulation and induction of effector cell phenotypes. The model puts forth a conceptual framework for the integration of functionally relevant processes in the analysis of the T CD4 cell function.
Collapse
Affiliation(s)
- David Martínez-Méndez
- Instituto de Investigaciones Biomédicas, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - Carlos Villarreal
- Instituto de Física, Universidad Nacional Autónoma de México, Mexico City, Mexico.,Centro de Ciencias de la Complejidad, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - Luis Mendoza
- Instituto de Investigaciones Biomédicas, Universidad Nacional Autónoma de México, Mexico City, Mexico.,Centro de Ciencias de la Complejidad, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - Leonor Huerta
- Instituto de Investigaciones Biomédicas, Universidad Nacional Autónoma de México, Mexico City, Mexico
| |
Collapse
|
15
|
Réda C, Wilczyński B. Automated inference of gene regulatory networks using explicit regulatory modules. J Theor Biol 2020; 486:110091. [PMID: 31790679 DOI: 10.1016/j.jtbi.2019.110091] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2018] [Revised: 08/09/2019] [Accepted: 11/21/2019] [Indexed: 11/25/2022]
Abstract
Gene regulatory networks are a popular tool for modelling important biological phenomena, such as cell differentiation or oncogenesis. Efficient identification of the causal connections between genes, their products and regulating transcription factors, is key to understanding how defects in their function may trigger diseases. Modelling approaches should keep up with the ever more detailed descriptions of the biological phenomena at play, as provided by new experimental findings and technical improvements. In recent years, we have seen great improvements in mapping of specific binding sites of many transcription factors to distinct regulatory regions. Recent gene regulatory network models use binding measurements; but usually only to define gene-to-gene interactions, ignoring regulatory module structure. Moreover, current huge amount of transcriptomic data, and exploration of all possible cis-regulatory arrangements which can lead to the same transcriptomic response, makes manual model building both tedious and time-consuming. In our paper, we propose a method to specify possible regulatory connections in a given Boolean network, based on transcription factor binding evidence. This is implemented by an algorithm which expands a regular Boolean network model into a "cis-regulatory" Boolean network model. This expanded model explicitly defines regulatory regions as additional nodes in the network, and adds new, valuable biological insights to the system dynamics. The expanded model can automatically be compared with expression data. And, for each node, a regulatory function, consistent with the experimental data, can be found. The resulting models are usually more constrained (by biologically-motivated metadata), and can then be inspected in in silico simulations. The fully automated method for model identification has been implemented in Python, and the expansion algorithm in R. The method resorts to the Z3 Satisfiability Modulo Theories (SMT) solver, and is similar to the RE:IN application (Yordanov et al., 2016). It is available on https://github.com/regulomics/expansion-network.
Collapse
Affiliation(s)
- Clémence Réda
- École Normale Supérieure Paris-Saclay, 61 avenue du Président Wilson, 94230 Cachan, France.
| | - Bartek Wilczyński
- Faculty of Mathematics, Informatics and Mechanics, University of Warsaw, ulica Stefana Banacha 2, Warsaw 02-097, Poland
| |
Collapse
|
16
|
Gouveia F, Lynce I, Monteiro PT. Revision of Boolean Models of Regulatory Networks Using Stable State Observations. J Comput Biol 2020; 27:144-155. [PMID: 31794671 DOI: 10.1089/cmb.2019.0289] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
Models of biological regulatory networks are essential to understand cellular processes. However, the definition of such models is still mostly manually performed, and consequently prone to error. Moreover, as new experimental data are acquired, models need to be revised and updated. Here, we propose a model revision procedure and associated tool, capable of providing the set of minimal repairs to render a model consistent with a set of experimental observations. We consider four possible repair operations, using a lexicographic optimization criterion, giving preference to function repairs over topological ones. Also, we consider observations at stable state discarding the model dynamics. In this article, we extend our previous work to tackle the problem of repairing nodes with multiple reasons of inconsistency. We evaluate our tool on five publicly available logical models. We perform random changes considering several parameter configurations to assess the tool repairing capabilities. Whenever a model is repaired under the time limit, the tool successfully produces the optimal solutions to repair the model. Instances were generated without the previous limitation to validate this extended approach.
Collapse
Affiliation(s)
- Filipe Gouveia
- Department of Computer Science and Engineering, INESC-ID/Instituto Superior Técnico, Universidade de Lisboa, Lisboa, Portugal
| | - Inês Lynce
- Department of Computer Science and Engineering, INESC-ID/Instituto Superior Técnico, Universidade de Lisboa, Lisboa, Portugal
| | - Pedro T Monteiro
- Department of Computer Science and Engineering, INESC-ID/Instituto Superior Técnico, Universidade de Lisboa, Lisboa, Portugal
| |
Collapse
|
17
|
Buetti-Dinh A, Herold M, Christel S, El Hajjami M, Delogu F, Ilie O, Bellenberg S, Wilmes P, Poetsch A, Sand W, Vera M, Pivkin IV, Friedman R, Dopson M. Reverse engineering directed gene regulatory networks from transcriptomics and proteomics data of biomining bacterial communities with approximate Bayesian computation and steady-state signalling simulations. BMC Bioinformatics 2020; 21:23. [PMID: 31964336 PMCID: PMC6975020 DOI: 10.1186/s12859-019-3337-9] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2019] [Accepted: 12/30/2019] [Indexed: 12/12/2022] Open
Abstract
BACKGROUND Network inference is an important aim of systems biology. It enables the transformation of OMICs datasets into biological knowledge. It consists of reverse engineering gene regulatory networks from OMICs data, such as RNAseq or mass spectrometry-based proteomics data, through computational methods. This approach allows to identify signalling pathways involved in specific biological functions. The ability to infer causality in gene regulatory networks, in addition to correlation, is crucial for several modelling approaches and allows targeted control in biotechnology applications. METHODS We performed simulations according to the approximate Bayesian computation method, where the core model consisted of a steady-state simulation algorithm used to study gene regulatory networks in systems for which a limited level of details is available. The simulations outcome was compared to experimentally measured transcriptomics and proteomics data through approximate Bayesian computation. RESULTS The structure of small gene regulatory networks responsible for the regulation of biological functions involved in biomining were inferred from multi OMICs data of mixed bacterial cultures. Several causal inter- and intraspecies interactions were inferred between genes coding for proteins involved in the biomining process, such as heavy metal transport, DNA damage, replication and repair, and membrane biogenesis. The method also provided indications for the role of several uncharacterized proteins by the inferred connection in their network context. CONCLUSIONS The combination of fast algorithms with high-performance computing allowed the simulation of a multitude of gene regulatory networks and their comparison to experimentally measured OMICs data through approximate Bayesian computation, enabling the probabilistic inference of causality in gene regulatory networks of a multispecies bacterial system involved in biomining without need of single-cell or multiple perturbation experiments. This information can be used to influence biological functions and control specific processes in biotechnology applications.
Collapse
Affiliation(s)
- Antoine Buetti-Dinh
- Institute of Computational Science, Faculty of Informatics, Università della Svizzera Italiana, Via Giuseppe Buffi 13, Lugano, CH-6900 Switzerland
- Swiss Institute of Bioinformatics, Quartier Sorge – Batiment Genopode, Lausanne, CH-1015 Switzerland
- Department of Chemistry and Biomedical Sciences, Linnæus University, Hus Vita, Kalmar, SE-391 82 Sweden
- Linnæus University Centre for Biomaterials Chemistry, Linnæus University, Hus Vita, Kalmar, SE-391 82 Sweden
- Centre for Ecology and Evolution in Microbial Model Systems, Linnæus University, Hus Vita, Kalmar, SE-391 82 Sweden
| | - Malte Herold
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Stephan Christel
- Centre for Ecology and Evolution in Microbial Model Systems, Linnæus University, Hus Vita, Kalmar, SE-391 82 Sweden
| | | | - Francesco Delogu
- Faculty of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences, Oslo, Norway
| | - Olga Ilie
- Institute of Computational Science, Faculty of Informatics, Università della Svizzera Italiana, Via Giuseppe Buffi 13, Lugano, CH-6900 Switzerland
- Swiss Institute of Bioinformatics, Quartier Sorge – Batiment Genopode, Lausanne, CH-1015 Switzerland
| | - Sören Bellenberg
- Centre for Ecology and Evolution in Microbial Model Systems, Linnæus University, Hus Vita, Kalmar, SE-391 82 Sweden
| | - Paul Wilmes
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, Belvaux, Luxembourg
| | - Ansgar Poetsch
- Plant Biochemistry, Ruhr University Bochum, Bochum, Germany
- Center for Marine and Molecular Biotechnology, QNLM, Qingdao, China
- College of Marine Life Sciences, Ocean University of China, Qingdao, China
| | - Wolfgang Sand
- Faculty of Chemistry, Essen, Germany
- College of Environmental Science and Engineering, Donghua University, Shanghai, People’s Republic of China
- Mining Academy and Technical University Freiberg, Freiberg, Germany
| | - Mario Vera
- Institute for Biological and Medical Engineering. Schools of Engineering, Medicine & Biological Sciences, Pontificia Universidad Católica de Chile, Santiago, Chile
- Department of Hydraulic & Environmental Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile
| | - Igor V. Pivkin
- Institute of Computational Science, Faculty of Informatics, Università della Svizzera Italiana, Via Giuseppe Buffi 13, Lugano, CH-6900 Switzerland
- Swiss Institute of Bioinformatics, Quartier Sorge – Batiment Genopode, Lausanne, CH-1015 Switzerland
| | - Ran Friedman
- Department of Chemistry and Biomedical Sciences, Linnæus University, Hus Vita, Kalmar, SE-391 82 Sweden
- Linnæus University Centre for Biomaterials Chemistry, Linnæus University, Hus Vita, Kalmar, SE-391 82 Sweden
| | - Mark Dopson
- Centre for Ecology and Evolution in Microbial Model Systems, Linnæus University, Hus Vita, Kalmar, SE-391 82 Sweden
| |
Collapse
|
18
|
AEON: Attractor Bifurcation Analysis of Parametrised Boolean Networks. COMPUTER AIDED VERIFICATION 2020. [PMCID: PMC7363220 DOI: 10.1007/978-3-030-53288-8_28] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Abstract
Boolean networks (BNs) provide an effective modelling tool for various phenomena from science and engineering. Any long-term behaviour of a BN eventually converges to a so-called attractor. Depending on various logical parameters, the structure and quality of attractors can undergo a significant change, known as a bifurcation. We present a tool for analysing bifurcations in asynchronous parametrised Boolean networks. To fight the state-space and parameter-space explosion problem the tool uses a parallel semi-symbolic algorithm.
Collapse
|
19
|
Sordo Vieira L, Laubenbacher RC, Murrugarra D. Control of Intracellular Molecular Networks Using Algebraic Methods. Bull Math Biol 2019; 82:2. [PMID: 31919596 DOI: 10.1007/s11538-019-00679-w] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2019] [Accepted: 12/02/2019] [Indexed: 10/25/2022]
Abstract
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.
Collapse
Affiliation(s)
- Luis Sordo Vieira
- The Jackson Laboratory for Genomic Medicine, 10 Discovery Drive, Farmington, CT, 06032, USA
| | - Reinhard C Laubenbacher
- The Jackson Laboratory for Genomic Medicine, 10 Discovery Drive, Farmington, CT, 06032, USA.,Center for Quantitative Medicine, UConn Health, Farmington, CT, 06032, USA
| | - David Murrugarra
- Department of Mathematics, University of Kentucky, Lexington, KY, 40506, USA.
| |
Collapse
|
20
|
Deritei D, Rozum J, Ravasz Regan E, Albert R. A feedback loop of conditionally stable circuits drives the cell cycle from checkpoint to checkpoint. Sci Rep 2019; 9:16430. [PMID: 31712566 PMCID: PMC6848090 DOI: 10.1038/s41598-019-52725-1] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2019] [Accepted: 10/22/2019] [Indexed: 12/12/2022] Open
Abstract
We perform logic-based network analysis on a model of the mammalian cell cycle. This model is composed of a Restriction Switch driving cell cycle commitment and a Phase Switch driving mitotic entry and exit. By generalizing the concept of stable motif, i.e., a self-sustaining positive feedback loop that maintains an associated state, we introduce the concept of a conditionally stable motif, the stability of which is contingent on external conditions. We show that the stable motifs of the Phase Switch are contingent on the state of three nodes through which it receives input from the rest of the network. Biologically, these conditions correspond to cell cycle checkpoints. Holding these nodes locked (akin to a checkpoint-free cell) transforms the Phase Switch into an autonomous oscillator that robustly toggles through the cell cycle phases G1, G2 and mitosis. The conditionally stable motifs of the Phase Switch Oscillator are organized into an ordered sequence, such that they serially stabilize each other but also cause their own destabilization. Along the way they channel the dynamics of the module onto a narrow path in state space, lending robustness to the oscillation. Self-destabilizing conditionally stable motifs suggest a general negative feedback mechanism leading to sustained oscillations.
Collapse
Affiliation(s)
- Dávid Deritei
- Department of Physics, Pennsylvania State University, University Park, PA, United States of America
- Department of Network and Data Science, Central European University, Budapest, Hungary
| | - Jordan Rozum
- Department of Physics, Pennsylvania State University, University Park, PA, United States of America
| | - Erzsébet Ravasz Regan
- Biochemistry and Molecular Biology, The College of Wooster, Wooster, OH, United States of America
| | - Réka Albert
- Department of Physics, Pennsylvania State University, University Park, PA, United States of America.
| |
Collapse
|
21
|
Enciso J, Pelayo R, Villarreal C. From Discrete to Continuous Modeling of Lymphocyte Development and Plasticity in Chronic Diseases. Front Immunol 2019; 10:1927. [PMID: 31481957 PMCID: PMC6710364 DOI: 10.3389/fimmu.2019.01927] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2018] [Accepted: 07/30/2019] [Indexed: 12/12/2022] Open
Abstract
The molecular events leading to differentiation, development, and plasticity of lymphoid cells have been subject of intense research due to their key roles in multiple pathologies, such as lymphoproliferative disorders, tumor growth maintenance and chronic diseases. The emergent roles of lymphoid cells and the use of high-throughput technologies have led to an extensive accumulation of experimental data allowing the reconstruction of gene regulatory networks (GRN) by integrating biochemical signals provided by the microenvironment with transcriptional modules of lineage-specific genes. Computational modeling of GRN has been useful for the identification of molecular switches involved in lymphoid specification, prediction of microenvironment-dependent cell plasticity, and analyses of signaling events occurring downstream the activation of antigen recognition receptors. Among most common modeling strategies to analyze the dynamical behavior of GRN, discrete dynamic models are widely used for their capacity to capture molecular interactions when a limited knowledge of kinetic parameters is present. However, they are less powerful when modeling complex systems sensitive to biochemical gradients. To compensate it, discrete models may be transformed into regulatory networks that includes state variables and parameters varying within a continuous range. This approach is based on a system of differential equations dynamics with regulatory interactions described by fuzzy logic propositions. Here, we discuss the applicability of this method on modeling of development and plasticity processes of adaptive lymphocytes, and its potential implications in the study of pathological landscapes associated to chronic diseases.
Collapse
Affiliation(s)
- Jennifer Enciso
- Centro de Investigación Biomédica de Oriente, Instituto Mexicano del Seguro Social, Mexico City, Mexico
- Programa de Doctorado en Ciencias Biomédicas, Universidad Nacional Autónoma de México, Mexico City, Mexico
- Centro de Ciencias de la Complejidad, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - Rosana Pelayo
- Centro de Investigación Biomédica de Oriente, Instituto Mexicano del Seguro Social, Mexico City, Mexico
| | - Carlos Villarreal
- Centro de Ciencias de la Complejidad, Universidad Nacional Autónoma de México, Mexico City, Mexico
- Departamento de Física Cuántica y Fotónica, Instituto de Física, Universidad Nacional Autónoma de México, Mexico City, Mexico
| |
Collapse
|
22
|
Ruklisa D, Brazma A, Cerans K, Schlitt T, Viksna J. Dynamics of gene regulatory networks and their dependence on network topology and quantitative parameters - the case of phage λ. BMC Bioinformatics 2019; 20:296. [PMID: 31151381 PMCID: PMC6544977 DOI: 10.1186/s12859-019-2909-z] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2019] [Accepted: 05/21/2019] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Gene regulatory networks can be modelled in various ways depending on the level of detail required and biological questions addressed. One of the earliest formalisms used for modeling is a Boolean network, although these models cannot describe most temporal aspects of a biological system. Differential equation models have also been used to model gene regulatory networks, but these frameworks tend to be too detailed for large models and many quantitative parameters might not be deducible in practice. Hybrid models bridge the gap between these two model classes - these are useful when concentration changes are important while the information about precise concentrations and binding site affinities is partial. RESULTS In this paper we study the stable behaviours of phage λ via a hybrid system based model. We identify wild type and mutant behaviours that arise for various orderings of binding site affinities. We propose experiments for detecting these behaviours: we suggest several ways of altering binding affinities with either mutations or genome rearrangements to achieve modified behaviours. The feasibility of these experiments is assessed. The interplay between the qualitative aspects of a network, e.g. network topology, and quantitative parameters, e.g. growth and degradation rates of proteins, is demonstrated. We also provide a software for exploring all feasible states of a hybrid system model and identifying all attractors. CONCLUSIONS The behaviours of phage λ are determined mainly by the topology of this network and by the mutual order of binding affinities. Exact affinities and growth and degradation rates of proteins fine tune the system. We show that only two stable behaviours are possible for phage λ if the main constraints of λ switch are preserved - these behaviours correspond to lysis and lysogeny. We identify several variants of both lysis and lysogeny - one wild type and one modified behaviour for each. We elucidate the necessary constraints for binding site affinities to achieve both wild type lysis and lysogeny. Our software is applicable to a wide range of biological models described as a hybrid system.
Collapse
Affiliation(s)
- Dace Ruklisa
- Newnham College, University of Cambridge, Sidgwick Avenue, Cambridge, CB3 9DF, UK.
| | - Alvis Brazma
- European Molecular Biology Laboratory, European Bioinformatics Institute, EMBL-EBI, Wellcome Genome Campus, Hinxton, CB10 1SD, UK
| | - Karlis Cerans
- Institute of Mathematics and Computer Science, University of Latvia, Rainis blvd. 29, Riga, LV 1459, Latvia
| | - Thomas Schlitt
- Division of Molecular Neuroscience, Faculty of Psychology, University of Basel, Basel, 4055, Switzerland
| | - Juris Viksna
- Institute of Mathematics and Computer Science, University of Latvia, Rainis blvd. 29, Riga, LV 1459, Latvia
| |
Collapse
|
23
|
A Boolean network control algorithm guided by forward dynamic programming. PLoS One 2019; 14:e0215449. [PMID: 31048917 PMCID: PMC6497256 DOI: 10.1371/journal.pone.0215449] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2018] [Accepted: 04/02/2019] [Indexed: 11/19/2022] Open
Abstract
Control problem in a biological system is the problem of finding an interventional policy for changing the state of the biological system from an undesirable state, e.g. disease, into a desirable healthy state. Boolean networks are utilized as a mathematical model for gene regulatory networks. This paper provides an algorithm to solve the control problem in Boolean networks. The proposed algorithm is implemented and applied on two biological systems: T-cell receptor network and Drosophila melanogaster network. Results show that the proposed algorithm works faster in solving the control problem over these networks, while having similar accuracy, in comparison to previous exact methods. Source code and a simple web service of the proposed algorithm is available at http://goliaei.ir/net-control/www/.
Collapse
|
24
|
Rozum JC, Albert R. Identifying (un)controllable dynamical behavior in complex networks. PLoS Comput Biol 2018; 14:e1006630. [PMID: 30532150 PMCID: PMC6301693 DOI: 10.1371/journal.pcbi.1006630] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2018] [Revised: 12/20/2018] [Accepted: 11/11/2018] [Indexed: 12/13/2022] Open
Abstract
We present a technique applicable in any dynamical framework to identify control-robust subsets of an interacting system. These robust subsystems, which we call stable modules, are characterized by constraints on the variables that make up the subsystem. They are robust in the sense that if the defining constraints are satisfied at a given time, they remain satisfied for all later times, regardless of what happens in the rest of the system, and can only be broken if the constrained variables are externally manipulated. We identify stable modules as graph structures in an expanded network, which represents causal links between variable constraints. A stable module represents a system "decision point", or trap subspace. Using the expanded network, small stable modules can be composed sequentially to form larger stable modules that describe dynamics on the system level. Collections of large, mutually exclusive stable modules describe the system's repertoire of long-term behaviors. We implement this technique in a broad class of dynamical systems and illustrate its practical utility via examples and algorithmic analysis of two published biological network models. In the segment polarity gene network of Drosophila melanogaster, we obtain a state-space visualization that reproduces by novel means the four possible cell fates and predicts the outcome of cell transplant experiments. In the T-cell signaling network, we identify six signaling elements that determine the high-signal response and show that control of an element connected to them cannot disrupt this response.
Collapse
Affiliation(s)
- Jordan C Rozum
- Department of Physics, The Pennsylvania State University, University Park, PA, USA
| | - Réka Albert
- Department of Physics, The Pennsylvania State University, University Park, PA, USA
- Department of Biology, The Pennsylvania State University, University Park, PA, USA
| |
Collapse
|
25
|
Integrative workflows for network analysis. Essays Biochem 2018; 62:549-561. [DOI: 10.1042/ebc20180005] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2018] [Revised: 09/10/2018] [Accepted: 09/14/2018] [Indexed: 01/17/2023]
Abstract
Due to genetic heterogeneity across patients, the identification of effective disease signatures and therapeutic targets is challenging. Addressing this challenge, we have previously developed a network-based approach, which integrates heterogeneous sources of biological information to identify disease specific core-regulatory networks. In particular, our workflow uses a multi-objective optimization function to calculate a ranking score for network components (e.g. feedback/feedforward loops) based on network properties, biomedical and high-throughput expression data. High ranked network components are merged to identify the core-regulatory network(s) that is then subjected to dynamical analysis using stimulus–response and in silico perturbation experiments for the identification of disease gene signatures and therapeutic targets. In a case study, we implemented our workflow to identify bladder and breast cancer specific core-regulatory networks underlying epithelial–mesenchymal transition from the E2F1 molecular interaction map.
In this study, we review our workflow and described how it has developed over time to understand the mechanisms underlying disease progression and prediction of signatures for clinical decision making.
Collapse
|
26
|
Pauleve L. Reduction of Qualitative Models of Biological Networks for Transient Dynamics Analysis. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2018; 15:1167-1179. [PMID: 28885158 DOI: 10.1109/tcbb.2017.2749225] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Qualitative models of dynamics of signalling pathways and gene regulatory networks allow for the capturing of temporal properties of biological networks while requiring few parameters. However, these discrete models typically suffer from the so-called state space explosion problem which makes the formal assessment of their potential behaviors very challenging. In this paper, we describe a method to reduce a qualitative model for enhancing the tractability of analysis of transient reachability properties. The reduction does not change the dimension of the model, but instead limits its degree of freedom, therefore reducing the set of states and transitions to consider. We rely on a transition-centered specification of qualitative models by the mean of automata networks. Our framework encompasses the usual asynchronous Boolean and multi-valued network, as well as 1-bounded Petri nets. Applied to different large-scale biological networks from the litterature, we show that the reduction can lead to a drastic improvement for the scalability of verification methods.
Collapse
|
27
|
P_UNSAT approach of attractor calculation for Boolean gene regulatory networks. J Theor Biol 2018; 447:171-177. [PMID: 29605228 DOI: 10.1016/j.jtbi.2018.03.037] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2017] [Revised: 03/27/2018] [Accepted: 03/28/2018] [Indexed: 11/21/2022]
Abstract
Boolean network models provide an efficient way for studying gene regulatory networks. The main dynamics of a Boolean network is determined by its attractors. Attractor calculation plays a key role for analyzing Boolean gene regulatory networks. An approach of attractor calculation was proposed in this study, which combined the predecessor approach and the logic unsatisfiability approach to accelerate attractor calculation. The proposed algorithm is effective to calculate all attractors for large-scale Boolean gene regulatory networks even the networks with a relatively large average degree.
Collapse
|
28
|
Ross BC, Boguslav M, Weeks H, Costello JC. Simulating heterogeneous populations using Boolean models. BMC SYSTEMS BIOLOGY 2018; 12:64. [PMID: 29879983 PMCID: PMC5992775 DOI: 10.1186/s12918-018-0591-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/23/2018] [Accepted: 05/24/2018] [Indexed: 11/10/2022]
Abstract
Background Certain biological processes, such as the development of cancer and immune activation, can be controlled by rare cellular events that are difficult to capture computationally through simulations of individual cells. Information about such rare events can be gleaned from an attractor analysis, for which a variety of methods exist (in particular for Boolean models). However, explicitly simulating a defined mixed population of cells in a way that tracks even the rarest subpopulations remains an open challenge. Results Here we show that when cellular states are described using a Boolean network model, one can exactly simulate the dynamics of non-interacting, highly heterogeneous populations directly, without having to model the various subpopulations. This strategy captures even the rarest outcomes of the model with no sampling error. Our method can incorporate heterogeneity in both cell state and, by augmenting the model, the underlying rules of the network as well (e.g., introducing loss-of-function genetic alterations). We demonstrate our method by using it to simulate a heterogeneous population of Boolean networks modeling the T-cell receptor, spanning ∼ 1020 distinct cellular states and mutational profiles. Conclusions We have developed a method for using Boolean models to perform a population-level simulation, in which the population consists of non-interacting individuals existing in different states. This approach can be used even when there are far too many distinct subpopulations to model individually. Electronic supplementary material The online version of this article (10.1186/s12918-018-0591-9) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Brian C Ross
- Computational Bioscience Program, University of Colorado Anschutz Medical Campus, 12801 E. 17th Ave., Aurora, CO, 80045, USA.,Department of Pharmacology, University of Colorado Anschutz Medical Campus, 12800 E. 19th Ave., Aurora, CO, 80045, USA
| | - Mayla Boguslav
- Computational Bioscience Program, University of Colorado Anschutz Medical Campus, 12801 E. 17th Ave., Aurora, CO, 80045, USA.,Department of Pharmacology, University of Colorado Anschutz Medical Campus, 12800 E. 19th Ave., Aurora, CO, 80045, USA
| | - Holly Weeks
- Department of Biostatistics and Informatics, University of Colorado Anschutz Medical Campus, 13001 E. 17th Place, Aurora, CO, 80045, USA
| | - James C Costello
- Computational Bioscience Program, University of Colorado Anschutz Medical Campus, 12801 E. 17th Ave., Aurora, CO, 80045, USA. .,Department of Pharmacology, University of Colorado Anschutz Medical Campus, 12800 E. 19th Ave., Aurora, CO, 80045, USA.
| |
Collapse
|
29
|
Yang G, Gómez Tejeda Zañudo J, Albert R. Target Control in Logical Models Using the Domain of Influence of Nodes. Front Physiol 2018; 9:454. [PMID: 29867523 PMCID: PMC5951947 DOI: 10.3389/fphys.2018.00454] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2018] [Accepted: 04/13/2018] [Indexed: 11/13/2022] Open
Abstract
Dynamical models of biomolecular networks are successfully used to understand the mechanisms underlying complex diseases and to design therapeutic strategies. Network control and its special case of target control, is a promising avenue toward developing disease therapies. In target control it is assumed that a small subset of nodes is most relevant to the system's state and the goal is to drive the target nodes into their desired states. An example of target control would be driving a cell to commit to apoptosis (programmed cell death). From the experimental perspective, gene knockout, pharmacological inhibition of proteins, and providing sustained external signals are among practical intervention techniques. We identify methodologies to use the stabilizing effect of sustained interventions for target control in Boolean network models of biomolecular networks. Specifically, we define the domain of influence (DOI) of a node (in a certain state) to be the nodes (and their corresponding states) that will be ultimately stabilized by the sustained state of this node regardless of the initial state of the system. We also define the related concept of the logical domain of influence (LDOI) of a node, and develop an algorithm for its identification using an auxiliary network that incorporates the regulatory logic. This way a solution to the target control problem is a set of nodes whose DOI can cover the desired target node states. We perform greedy randomized adaptive search in node state space to find such solutions. We apply our strategy to in silico biological network models of real systems to demonstrate its effectiveness.
Collapse
Affiliation(s)
- Gang Yang
- Department of Physics, Pennsylvania State University, University Park, PA, United States
| | - Jorge Gómez Tejeda Zañudo
- Department of Physics, Pennsylvania State University, University Park, PA, United States.,Department of Medical Oncology, Dana-Farber Cancer Institute, Boston, MA, United States.,Eli and Edythe L. Broad Institute of MIT and Harvard, Cambridge, MA, United States
| | - Réka Albert
- Department of Physics, Pennsylvania State University, University Park, PA, United States.,Department of Biology, Pennsylvania State University, University Park, PA, United States
| |
Collapse
|
30
|
Melkman AA, Cheng X, Ching WK, Akutsu T. Identifying a Probabilistic Boolean Threshold Network From Samples. IEEE TRANSACTIONS ON NEURAL NETWORKS AND LEARNING SYSTEMS 2018; 29:869-881. [PMID: 28129190 DOI: 10.1109/tnnls.2017.2648039] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
This paper studies the problem of exactly identifying the structure of a probabilistic Boolean network (PBN) from a given set of samples, where PBNs are probabilistic extensions of Boolean networks. Cheng et al. studied the problem while focusing on PBNs consisting of pairs of AND/OR functions. This paper considers PBNs consisting of Boolean threshold functions while focusing on those threshold functions that have unit coefficients. The treatment of Boolean threshold functions, and triplets and -tuplets of such functions, necessitates a deepening of the theoretical analyses. It is shown that wide classes of PBNs with such threshold functions can be exactly identified from samples under reasonable constraints, which include: 1) PBNs in which any number of threshold functions can be assigned provided that all have the same number of input variables and 2) PBNs consisting of pairs of threshold functions with different numbers of input variables. It is also shown that the problem of deciding the equivalence of two Boolean threshold functions is solvable in pseudopolynomial time but remains co-NP complete.
Collapse
|
31
|
Thiele S, Heise S, Hessenkemper W, Bongartz H, Fensky M, Schaper F, Klamt S. Designing optimal experiments to discriminate interaction graph models. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2018; 16:925-935. [PMID: 29993657 DOI: 10.1109/tcbb.2018.2812184] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Modern methods for the inference of cellular networks from experimental data often express nondeterminism through an ensemble of candidate models. To discriminate among these candidates new experiments need to be carried out. Theoretically, the number of possible experiments is exponential in the number of possible perturbations. In praxis, experiments are expensive and there exist several limiting constraints. Limiting factors exist on the combinations of perturbations that are technically possible, which components can be measured, and on the number of affordable experiments. Further, not all experiments are equally well suited to discriminate model candidates. The goal of optimal experiment design is to determine those experiments that discriminate most of the candidates while minimizing the costs. We present an approach for experiment planning with interaction graph models and sign consistency methods. This new approach can be used in combination with methods for network inference and consistency checking. We applied our method to study the Erythropoietin signal transduction in human kidney cells HEK293. We first used simulated experiment data from an ODE model to demonstrate in silico that our experimental design results in the inference of the gold standard model. Finally, we used the approach to plan in vivo experiments that discriminate model candidates for the Erythropoietin signal transduction in this cell line.
Collapse
|
32
|
Bloomingdale P, Nguyen VA, Niu J, Mager DE. Boolean network modeling in systems pharmacology. J Pharmacokinet Pharmacodyn 2018; 45:159-180. [PMID: 29307099 PMCID: PMC6531050 DOI: 10.1007/s10928-017-9567-4] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2017] [Accepted: 12/29/2017] [Indexed: 01/01/2023]
Abstract
Quantitative systems pharmacology (QSP) is an emerging discipline that aims to discover how drugs modulate the dynamics of biological components in molecular and cellular networks and the impact of those perturbations on human pathophysiology. The integration of systems-based experimental and computational approaches is required to facilitate the advancement of this field. QSP models typically consist of a series of ordinary differential equations (ODE). However, this mathematical framework requires extensive knowledge of parameters pertaining to biological processes, which is often unavailable. An alternative framework that does not require knowledge of system-specific parameters, such as Boolean network modeling, could serve as an initial foundation prior to the development of an ODE-based model. Boolean network models have been shown to efficiently describe, in a qualitative manner, the complex behavior of signal transduction and gene/protein regulatory processes. In addition to providing a starting point prior to quantitative modeling, Boolean network models can also be utilized to discover novel therapeutic targets and combinatorial treatment strategies. Identifying drug targets using a network-based approach could supplement current drug discovery methodologies and help to fill the innovation gap across the pharmaceutical industry. In this review, we discuss the process of developing Boolean network models and the various analyses that can be performed to identify novel drug targets and combinatorial approaches. An example for each of these analyses is provided using a previously developed Boolean network of signaling pathways in multiple myeloma. Selected examples of Boolean network models of human (patho-)physiological systems are also reviewed in brief.
Collapse
Affiliation(s)
- Peter Bloomingdale
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, University at Buffalo, The State University of New York, 431 Kapoor Hall, Buffalo, NY, 14214, USA
| | - Van Anh Nguyen
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, University at Buffalo, The State University of New York, 431 Kapoor Hall, Buffalo, NY, 14214, USA
| | - Jin Niu
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, University at Buffalo, The State University of New York, 431 Kapoor Hall, Buffalo, NY, 14214, USA
| | - Donald E Mager
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, University at Buffalo, The State University of New York, 431 Kapoor Hall, Buffalo, NY, 14214, USA.
| |
Collapse
|
33
|
Zhao XM, Li S. HISP: a hybrid intelligent approach for identifying directed signaling pathways. J Mol Cell Biol 2018; 9:453-462. [DOI: 10.1093/jmcb/mjx054] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2017] [Accepted: 12/20/2017] [Indexed: 01/15/2023] Open
Affiliation(s)
- Xing-Ming Zhao
- Institute of Science and Technology for Brain-Inspired Intelligence, Fudan University, Shanghai 200433, China
| | - Shan Li
- Department of Mathematics, Shanghai University, Shanghai 200444, China
| |
Collapse
|
34
|
Gómez Tejeda Zañudo J, Scaltriti M, Albert R. A network modeling approach to elucidate drug resistance mechanisms and predict combinatorial drug treatments in breast cancer. CANCER CONVERGENCE 2017; 1:5. [PMID: 29623959 PMCID: PMC5876695 DOI: 10.1186/s41236-017-0007-6] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2017] [Accepted: 11/30/2017] [Indexed: 02/08/2023] Open
Abstract
Background Mechanistic models of within-cell signal transduction networks can explain how these networks integrate internal and external inputs to give rise to the appropriate cellular response. These models can be fruitfully used in cancer cells, whose aberrant decision-making regarding their survival or death, proliferation or quiescence can be connected to errors in the state of nodes or edges of the signal transduction network. Results Here we present a comprehensive network, and discrete dynamic model, of signal transduction in ER+ breast cancer based on the literature of ER+, HER2+, and PIK3CA-mutant breast cancers. The network model recapitulates known resistance mechanisms to PI3K inhibitors and suggests other possibilities for resistance. The model also reveals known and novel combinatorial interventions that are more effective than PI3K inhibition alone. Conclusions The use of a logic-based, discrete dynamic model enables the identification of results that are mainly due to the organization of the signaling network, and those that also depend on the kinetics of individual events. Network-based models such as this will play an increasing role in the rational design of high-order therapeutic combinations. Electronic supplementary material The online version of this article (10.1186/s41236-017-0007-6) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Jorge Gómez Tejeda Zañudo
- 1Department of Physics, The Pennsylvania State University, University Park, PA 16802-6300 USA.,2Department of Medical Oncology, Dana-Farber Cancer Institute, Boston, MA 02215 USA.,3Broad Institute of Harvard and Massachusetts Institute of Technology, 7 Cambridge Center, Cambridge, MA 02142 USA
| | - Maurizio Scaltriti
- 4Human Oncology & Pathogenesis Program, Memorial Sloan Kettering Cancer Center, 1275 York Avenue, New York, NY 10065 USA.,5Department of Pathology, Memorial Sloan Kettering Cancer Center, 1275 York Avenue, New York, NY 10065 USA
| | - Réka Albert
- 1Department of Physics, The Pennsylvania State University, University Park, PA 16802-6300 USA.,6Department of Biology, The Pennsylvania State University, University Park, PA 16802-6300 USA
| |
Collapse
|
35
|
Bakker E, Tian K, Mutti L, Demonacos C, Schwartz JM, Krstic-Demonacos M. Insight into glucocorticoid receptor signalling through interactome model analysis. PLoS Comput Biol 2017; 13:e1005825. [PMID: 29107989 PMCID: PMC5690696 DOI: 10.1371/journal.pcbi.1005825] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2017] [Revised: 11/16/2017] [Accepted: 10/16/2017] [Indexed: 12/27/2022] Open
Abstract
Glucocorticoid hormones (GCs) are used to treat a variety of diseases because of their potent anti-inflammatory effect and their ability to induce apoptosis in lymphoid malignancies through the glucocorticoid receptor (GR). Despite ongoing research, high glucocorticoid efficacy and widespread usage in medicine, resistance, disease relapse and toxicity remain factors that need addressing. Understanding the mechanisms of glucocorticoid signalling and how resistance may arise is highly important towards improving therapy. To gain insight into this we undertook a systems biology approach, aiming to generate a Boolean model of the glucocorticoid receptor protein interaction network that encapsulates functional relationships between the GR, its target genes or genes that target GR, and the interactions between the genes that interact with the GR. This model named GEB052 consists of 52 nodes representing genes or proteins, the model input (GC) and model outputs (cell death and inflammation), connected by 241 logical interactions of activation or inhibition. 323 changes in the relationships between model constituents following in silico knockouts were uncovered, and steady-state analysis followed by cell-based microarray genome-wide model validation led to an average of 57% correct predictions, which was taken further by assessment of model predictions against patient microarray data. Lastly, semi-quantitative model analysis via microarray data superimposed onto the model with a score flow algorithm has also been performed, which demonstrated significantly higher correct prediction ratios (average of 80%), and the model has been assessed as a predictive clinical tool using published patient microarray data. In summary we present an in silico simulation of the glucocorticoid receptor interaction network, linked to downstream biological processes that can be analysed to uncover relationships between GR and its interactants. Ultimately the model provides a platform for future development both by directing laboratory research and allowing for incorporation of further components, encapsulating more interactions/genes involved in glucocorticoid receptor signalling. Here we present modelling of the glucocorticoid receptor (GR) signalling network. The GR is the effector for a class of drugs known as corticosteroids, which are widely used in medicine for their anti-inflammatory effects and ability to induce apoptosis in leukaemic cells. However, side effects, treatment-related toxicity and glucocorticoid resistance remain and therefore increased understanding of the glucocorticoid receptor mechanism of action may improve therapeutic outcomes. The GEB052 model presented herein has been used to generate predictions for how the network is altered between glucocorticoid-sensitive and glucocorticoid-resistant scenarios, and these predictions have been verified using published gene expression data from established cell lines (for both qualitative and semi-quantitative analysis). The model has also been preliminarily assessed as a predictive clinical tool by correlating model predictions with clinical outcomes of thirteen leukaemia patients. Thus, the GEB052 model demonstrates successful modelling to understand GR function. GEB052 provides accurate predictions and has indicated potential routes through which glucocorticoid resistance may arise. The work presented herein thus demonstrates a proof-of-principle of this modelling approach to furthering GR research, and provides insight into potential mechanisms of corticosteroids resistance.
Collapse
Affiliation(s)
- Emyr Bakker
- Biomedical Research Centre, School of Environment and Life Sciences, University of Salford, Salford, United Kingdom
| | - Kun Tian
- Biomedical Research Centre, School of Environment and Life Sciences, University of Salford, Salford, United Kingdom
| | - Luciano Mutti
- Biomedical Research Centre, School of Environment and Life Sciences, University of Salford, Salford, United Kingdom
| | - Constantinos Demonacos
- Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom
| | - Jean-Marc Schwartz
- Faculty of Biology, Medicine and Health, University of Manchester, Manchester, United Kingdom
- * E-mail: (JMS); (MKD)
| | - Marija Krstic-Demonacos
- Biomedical Research Centre, School of Environment and Life Sciences, University of Salford, Salford, United Kingdom
- * E-mail: (JMS); (MKD)
| |
Collapse
|
36
|
Cardelli L, Tribastone M, Tschaikowski M, Vandin A. Maximal aggregation of polynomial dynamical systems. Proc Natl Acad Sci U S A 2017; 114:10029-10034. [PMID: 28878023 PMCID: PMC5617256 DOI: 10.1073/pnas.1702697114] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Ordinary differential equations (ODEs) with polynomial derivatives are a fundamental tool for understanding the dynamics of systems across many branches of science, but our ability to gain mechanistic insight and effectively conduct numerical evaluations is critically hindered when dealing with large models. Here we propose an aggregation technique that rests on two notions of equivalence relating ODE variables whenever they have the same solution (backward criterion) or if a self-consistent system can be written for describing the evolution of sums of variables in the same equivalence class (forward criterion). A key feature of our proposal is to encode a polynomial ODE system into a finitary structure akin to a formal chemical reaction network. This enables the development of a discrete algorithm to efficiently compute the largest equivalence, building on approaches rooted in computer science to minimize basic models of computation through iterative partition refinements. The physical interpretability of the aggregation is shown on polynomial ODE systems for biochemical reaction networks, gene regulatory networks, and evolutionary game theory.
Collapse
Affiliation(s)
- Luca Cardelli
- Microsoft Research, Cambridge CB1 2FB, United Kingdom
- Department of Computing, University of Oxford, Oxford OX1 3QD, United Kingdom
| | | | | | | |
Collapse
|
37
|
Ben Abdallah E, Folschette M, Roux O, Magnin M. ASP-based method for the enumeration of attractors in non-deterministic synchronous and asynchronous multi-valued networks. Algorithms Mol Biol 2017; 12:20. [PMID: 28814968 PMCID: PMC5557630 DOI: 10.1186/s13015-017-0111-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2016] [Accepted: 07/26/2017] [Indexed: 01/31/2023] Open
Abstract
BACKGROUND This paper addresses the problem of finding attractors in biological regulatory networks. We focus here on non-deterministic synchronous and asynchronous multi-valued networks, modeled using automata networks (AN). AN is a general and well-suited formalism to study complex interactions between different components (genes, proteins,...). An attractor is a minimal trap domain, that is, a part of the state-transition graph that cannot be escaped. Such structures are terminal components of the dynamics and take the form of steady states (singleton) or complex compositions of cycles (non-singleton). Studying the effect of a disease or a mutation on an organism requires finding the attractors in the model to understand the long-term behaviors. RESULTS We present a computational logical method based on answer set programming (ASP) to identify all attractors. Performed without any network reduction, the method can be applied on any dynamical semantics. In this paper, we present the two most widespread non-deterministic semantics: the asynchronous and the synchronous updating modes. The logical approach goes through a complete enumeration of the states of the network in order to find the attractors without the necessity to construct the whole state-transition graph. We realize extensive computational experiments which show good performance and fit the expected theoretical results in the literature. CONCLUSION The originality of our approach lies on the exhaustive enumeration of all possible (sets of) states verifying the properties of an attractor thanks to the use of ASP. Our method is applied to non-deterministic semantics in two different schemes (asynchronous and synchronous). The merits of our methods are illustrated by applying them to biological examples of various sizes and comparing the results with some existing approaches. It turns out that our approach succeeds to exhaustively enumerate on a desktop computer, in a large model (100 components), all existing attractors up to a given size (20 states). This size is only limited by memory and computation time.
Collapse
Affiliation(s)
- Emna Ben Abdallah
- École Centrale de Nantes, LS2N UMR CNRS 6004, 1 rue de la Noë, 44321 Nantes, France
| | - Maxime Folschette
- Université de Nantes, LS2N UMR CNRS 6004, 1 rue de la Noë, 44321 Nantes, France
- Université Nice Sophia Antipolis, I3S UMR CNRS 7271, 2000, route des Lucioles, 06900 Nice, France
| | - Olivier Roux
- École Centrale de Nantes, LS2N UMR CNRS 6004, 1 rue de la Noë, 44321 Nantes, France
| | - Morgan Magnin
- École Centrale de Nantes, LS2N UMR CNRS 6004, 1 rue de la Noë, 44321 Nantes, France
- National Institute of Informatics, 2-1-2, Hitotsubashi, Chiyoda-ku, Tokyo, 101-8430 Japan
| |
Collapse
|
38
|
Unraveling a tumor type-specific regulatory core underlying E2F1-mediated epithelial-mesenchymal transition to predict receptor protein signatures. Nat Commun 2017; 8:198. [PMID: 28775339 PMCID: PMC5543083 DOI: 10.1038/s41467-017-00268-2] [Citation(s) in RCA: 61] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2016] [Accepted: 06/15/2017] [Indexed: 12/18/2022] Open
Abstract
Cancer is a disease of subverted regulatory pathways. In this paper, we reconstruct the regulatory network around E2F, a family of transcription factors whose deregulation has been associated to cancer progression, chemoresistance, invasiveness, and metastasis. We integrate gene expression profiles of cancer cell lines from two E2F1-driven highly aggressive bladder and breast tumors, and use network analysis methods to identify the tumor type-specific core of the network. By combining logic-based network modeling, in vitro experimentation, and gene expression profiles from patient cohorts displaying tumor aggressiveness, we identify and experimentally validate distinctive, tumor type-specific signatures of receptor proteins associated to epithelial-mesenchymal transition in bladder and breast cancer. Our integrative network-based methodology, exemplified in the case of E2F1-induced aggressive tumors, has the potential to support the design of cohort- as well as tumor type-specific treatments and ultimately, to fight metastasis and therapy resistance.Deregulation of E2F family transcription factors is associated with cancer progression and metastasis. Here, the authors construct a map of the regulatory network around the E2F family, and using gene expression profiles, identify tumour type-specific regulatory cores and receptor expression signatures associated with epithelial-mesenchymal transition in bladder and breast cancer.
Collapse
|
39
|
Traynard P, Tobalina L, Eduati F, Calzone L, Saez-Rodriguez J. Logic Modeling in Quantitative Systems Pharmacology. CPT-PHARMACOMETRICS & SYSTEMS PHARMACOLOGY 2017; 6:499-511. [PMID: 28681552 PMCID: PMC5572374 DOI: 10.1002/psp4.12225] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/06/2017] [Revised: 06/01/2017] [Accepted: 06/15/2017] [Indexed: 12/12/2022]
Abstract
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).
Collapse
Affiliation(s)
- Pauline Traynard
- Institut Curie, PSL Research University, Mines Paris Tech, Inserm, U900, Paris, France
| | - Luis Tobalina
- RWTH Aachen University, Faculty of Medicine, Joint Research Centre for Computational Biomedicine, Aachen, Germany
| | - Federica Eduati
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, UK
| | - Laurence Calzone
- Institut Curie, PSL Research University, Mines Paris Tech, Inserm, U900, Paris, France
| | - Julio Saez-Rodriguez
- RWTH Aachen University, Faculty of Medicine, Joint Research Centre for Computational Biomedicine, Aachen, Germany.,European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, UK
| |
Collapse
|
40
|
Amstein L, Ackermann J, Scheidel J, Fulda S, Dikic I, Koch I. Manatee invariants reveal functional pathways in signaling networks. BMC SYSTEMS BIOLOGY 2017; 11:72. [PMID: 28754124 PMCID: PMC5534052 DOI: 10.1186/s12918-017-0448-7] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/19/2016] [Accepted: 07/19/2017] [Indexed: 11/26/2022]
Abstract
Background Signal transduction pathways are important cellular processes to maintain the cell’s integrity. Their imbalance can cause severe pathologies. As signal transduction pathways feature complex regulations, they form intertwined networks. Mathematical models aim to capture their regulatory logic and allow an unbiased analysis of robustness and vulnerability of the signaling network. Pathway detection is yet a challenge for the analysis of signaling networks in the field of systems biology. A rigorous mathematical formalism is lacking to identify all possible signal flows in a network model. Results In this paper, we introduce the concept of Manatee invariants for the analysis of signal transduction networks. We present an algorithm for the characterization of the combinatorial diversity of signal flows, e.g., from signal reception to cellular response. We demonstrate the concept for a small model of the TNFR1-mediated NF- κB signaling pathway. Manatee invariants reveal all possible signal flows in the network. Further, we show the application of Manatee invariants for in silico knockout experiments. Here, we illustrate the biological relevance of the concept. Conclusions The proposed mathematical framework reveals the entire variety of signal flows in models of signaling systems, including cyclic regulations. Thereby, Manatee invariants allow for the analysis of robustness and vulnerability of signaling networks. The application to further analyses such as for in silico knockout was shown. The new framework of Manatee invariants contributes to an advanced examination of signaling systems. Electronic supplementary material The online version of this article (doi:10.1186/s12918-017-0448-7) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Leonie Amstein
- Molecular Bioinformatics, Institute of Computer Science, Goethe-University Frankfurt am Main, Robert-Mayer-Straße 11-15, Frankfurt am Main, 60325, Germany
| | - Jörg Ackermann
- Molecular Bioinformatics, Institute of Computer Science, Goethe-University Frankfurt am Main, Robert-Mayer-Straße 11-15, Frankfurt am Main, 60325, Germany
| | - Jennifer Scheidel
- Molecular Bioinformatics, Institute of Computer Science, Goethe-University Frankfurt am Main, Robert-Mayer-Straße 11-15, Frankfurt am Main, 60325, Germany
| | - Simone Fulda
- Institute for Experimental Cancer Research in Pediatrics, Goethe-University Frankfurt am Main, Komturstraße 3a, Frankfurt am Main, 60528, Germany.,German Cancer Consortium (DKTK), Heidelberg, Germany.,German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Ivan Dikic
- Institute of Biochemistry II, Goethe-University Hospital Frankfurt am Main, Theodor-Stern-Kai 7, Frankfurt am Main, 60590, Germany.,Buchmann Institute for Molecular Live Sciences, Max-von-Laue-Straße 15, Frankfurt am Main, 60438, Germany
| | - Ina Koch
- Molecular Bioinformatics, Institute of Computer Science, Goethe-University Frankfurt am Main, Robert-Mayer-Straße 11-15, Frankfurt am Main, 60325, Germany.
| |
Collapse
|
41
|
Teku GN, Vihinen M. Simulation of the dynamics of primary immunodeficiencies in CD4+ T-cells. PLoS One 2017; 12:e0176500. [PMID: 28448599 PMCID: PMC5407609 DOI: 10.1371/journal.pone.0176500] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2017] [Accepted: 04/11/2017] [Indexed: 01/05/2023] Open
Abstract
Primary immunodeficiencies (PIDs) form a large and heterogeneous group of mainly rare disorders that affect the immune system. T-cell deficiencies account for about one-tenth of PIDs, most of them being monogenic. Apart from genetic and clinical information, lots of other data are available for PID proteins and genes, including functions and interactions. Thus, it is possible to perform systems biology studies on the effects of PIDs on T-cell physiology and response. To achieve this, we reconstructed a T-cell network model based on literature mining and TPPIN, a previously published core T-cell network, and performed semi-quantitative dynamic network simulations on both normal and T-cell PID failure modes. The results for several loss-of-function PID simulations correspond to results of previously reported molecular studies. The simulations for TCR PTPRC, LCK, ZAP70 and ITK indicate profound changes to numerous proteins in the network. Significant effects were observed also in the BCL10, CARD11, MALT1, NEMO, IKKB and MAP3K14 simulations. No major effects were observed for PIDs that are caused by constitutively active proteins. The T-cell model facilitates the understanding of the underlying dynamics of PID disease processes. The approach confirms previous knowledge about T-cell signaling network and indicates several new important proteins that may be of interest when developing novel diagnosis and therapies to treat immunological defects.
Collapse
Affiliation(s)
- Gabriel N. Teku
- Department of Experimental Medical Science, Lund University, Lund, Sweden
| | - Mauno Vihinen
- Department of Experimental Medical Science, Lund University, Lund, Sweden
- * E-mail:
| |
Collapse
|
42
|
Adlung L, Kar S, Wagner MC, She B, Chakraborty S, Bao J, Lattermann S, Boerries M, Busch H, Wuchter P, Ho AD, Timmer J, Schilling M, Höfer T, Klingmüller U. Protein abundance of AKT and ERK pathway components governs cell type-specific regulation of proliferation. Mol Syst Biol 2017; 13:904. [PMID: 28123004 PMCID: PMC5293153 DOI: 10.15252/msb.20167258] [Citation(s) in RCA: 57] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Signaling through the AKT and ERK pathways controls cell proliferation. However, the integrated regulation of this multistep process, involving signal processing, cell growth and cell cycle progression, is poorly understood. Here, we study different hematopoietic cell types, in which AKT and ERK signaling is triggered by erythropoietin (Epo). Although these cell types share the molecular network topology for pro‐proliferative Epo signaling, they exhibit distinct proliferative responses. Iterating quantitative experiments and mathematical modeling, we identify two molecular sources for cell type‐specific proliferation. First, cell type‐specific protein abundance patterns cause differential signal flow along the AKT and ERK pathways. Second, downstream regulators of both pathways have differential effects on proliferation, suggesting that protein synthesis is rate‐limiting for faster cycling cells while slower cell cycles are controlled at the G1‐S progression. The integrated mathematical model of Epo‐driven proliferation explains cell type‐specific effects of targeted AKT and ERK inhibitors and faithfully predicts, based on the protein abundance, anti‐proliferative effects of inhibitors in primary human erythroid progenitor cells. Our findings suggest that the effectiveness of targeted cancer therapy might become predictable from protein abundance.
Collapse
Affiliation(s)
- Lorenz Adlung
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Sandip Kar
- Division of Theoretical Systems Biology, German Cancer Research Center (DKFZ), Heidelberg, Germany.,BioQuant Center, University of Heidelberg, Heidelberg, Germany.,Department of Chemistry, Indian Institute of Technology, Mumbai, India
| | - Marie-Christine Wagner
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Bin She
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Sajib Chakraborty
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Jie Bao
- Systems Biology of the Cellular Microenvironment Group, IMMZ, ALU, Freiburg, Germany
| | - Susen Lattermann
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Melanie Boerries
- Systems Biology of the Cellular Microenvironment Group, IMMZ, ALU, Freiburg, Germany.,German Cancer Consortium (DKTK), Freiburg, Germany.,German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Hauke Busch
- Systems Biology of the Cellular Microenvironment Group, IMMZ, ALU, Freiburg, Germany.,German Cancer Consortium (DKTK), Freiburg, Germany.,German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Patrick Wuchter
- Department of Medicine V, University of Heidelberg, Heidelberg, Germany.,Institute for Transfusion Medicine and Immunology, University of Heidelberg, Mannheim, Germany
| | - Anthony D Ho
- Department of Medicine V, University of Heidelberg, Heidelberg, Germany
| | - Jens Timmer
- Center for Biological Signaling Studies (BIOSS), Institute of Physics, University of Freiburg, Freiburg, Germany
| | - Marcel Schilling
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - Thomas Höfer
- Division of Theoretical Systems Biology, German Cancer Research Center (DKFZ), Heidelberg, Germany .,BioQuant Center, University of Heidelberg, Heidelberg, Germany
| | - Ursula Klingmüller
- Division of Systems Biology of Signal Transduction, German Cancer Research Center (DKFZ), Heidelberg, Germany .,Translational Lung Research Center (TLRC), Member of the German Center for Lung Research (DZL), Heidelberg, Germany
| |
Collapse
|
43
|
Dinwoodie IH. COMPUTATIONAL METHODS FOR ASYNCHRONOUS BASINS. DISCRETE AND CONTINUOUS DYNAMICAL SYSTEMS. SERIES B 2016; 21:3391-3405. [PMID: 28154501 PMCID: PMC5283826 DOI: 10.3934/dcdsb.2016103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
For a Boolean network we consider asynchronous updates and define the exclusive asynchronous basin of attraction for any steady state or cyclic attractor. An algorithm based on commutative algebra is presented to compute the exclusive basin. Finally its use for targeting desirable attractors by selective intervention on network nodes is illustrated with two examples, one cell signalling network and one sensor network measuring human mobility.
Collapse
|
44
|
He Z, Zhan M, Liu S, Fang Z, Yao C. An Algorithm for Finding the Singleton Attractors and Pre-Images in Strong-Inhibition Boolean Networks. PLoS One 2016; 11:e0166906. [PMID: 27861624 PMCID: PMC5115838 DOI: 10.1371/journal.pone.0166906] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2016] [Accepted: 11/04/2016] [Indexed: 11/18/2022] Open
Abstract
The detection of the singleton attractors is of great significance for the systematic study of genetic regulatory network. In this paper, we design an algorithm to compute the singleton attractors and pre-images of the strong-inhibition Boolean networks which is a biophysically plausible gene model. Our algorithm can not only identify accurately the singleton attractors, but also find easily the pre-images of the network. Based on extensive computational experiments, we show that the computational time of the algorithm is proportional to the number of the singleton attractors, which indicates the algorithm has much advantage in finding the singleton attractors for the networks with high average degree and less inhibitory interactions. Our algorithm may shed light on understanding the function and structure of the strong-inhibition Boolean networks.
Collapse
Affiliation(s)
- Zhiwei He
- Department of Mathematics, Shaoxing University, Shaoxing, China
| | - Meng Zhan
- State Key Laboratory of Advanced Electromagnetic Engineering and Technology, School of Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan, China
| | - Shuai Liu
- College of Science, Northwest A&F University, Yangling, China
| | - Zebo Fang
- Department of Physics, Shaoxing University, Shaoxing, China
| | - Chenggui Yao
- Department of Mathematics, Shaoxing University, Shaoxing, China
- * E-mail:
| |
Collapse
|
45
|
Wu H, von Kamp A, Leoncikas V, Mori W, Sahin N, Gevorgyan A, Linley C, Grabowski M, Mannan AA, Stoy N, Stewart GR, Ward LT, Lewis DJM, Sroka J, Matsuno H, Klamt S, Westerhoff HV, McFadden J, Plant NJ, Kierzek AM. MUFINS: multi-formalism interaction network simulator. NPJ Syst Biol Appl 2016; 2:16032. [PMID: 28725480 PMCID: PMC5516860 DOI: 10.1038/npjsba.2016.32] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2015] [Revised: 07/27/2016] [Accepted: 08/29/2016] [Indexed: 12/19/2022] Open
Abstract
Systems Biology has established numerous approaches for mechanistic modeling of molecular networks in the cell and a legacy of models. The current frontier is the integration of models expressed in different formalisms to address the multi-scale biological system organization challenge. We present MUFINS (MUlti-Formalism Interaction Network Simulator) software, implementing a unique set of approaches for multi-formalism simulation of interaction networks. We extend the constraint-based modeling (CBM) framework by incorporation of linear inhibition constraints, enabling for the first time linear modeling of networks simultaneously describing gene regulation, signaling and whole-cell metabolism at steady state. We present a use case where a logical hypergraph model of a regulatory network is expressed by linear constraints and integrated with a Genome-Scale Metabolic Network (GSMN) of mouse macrophage. We experimentally validate predictions, demonstrating application of our software in an iterative cycle of hypothesis generation, validation and model refinement. MUFINS incorporates an extended version of our Quasi-Steady State Petri Net approach to integrate dynamic models with CBM, which we demonstrate through a dynamic model of cortisol signaling integrated with the human Recon2 GSMN and a model of nutrient dynamics in physiological compartments. Finally, we implement a number of methods for deriving metabolic states from ~omics data, including our new variant of the iMAT congruency approach. We compare our approach with iMAT through the analysis of 262 individual tumor transcriptomes, recovering features of metabolic reprogramming in cancer. The software provides graphics user interface with network visualization, which facilitates use by researchers who are not experienced in coding and mathematical modeling environments.
Collapse
Affiliation(s)
- Huihai Wu
- School of Biosciences and Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Axel von Kamp
- Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg, Germany
| | - Vytautas Leoncikas
- School of Biosciences and Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Wataru Mori
- Graduate School of Science and Engineering and Faculty of Science, Yamaguchi University, Yoshida, Yamaguchi, Japan
| | - Nilgun Sahin
- Molecular Cell Physiology, VU University Amsterdam, Amsterdam, The Netherlands
| | | | - Catherine Linley
- School of Biosciences and Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Marek Grabowski
- Institute of Informatics, University of Warsaw, Warsaw, Poland
| | - Ahmad A Mannan
- School of Biosciences and Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Nicholas Stoy
- School of Biosciences and Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Graham R Stewart
- School of Biosciences and Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Lara T Ward
- Oncology DMPK, AstraZeneca, Alderley Park, Cheshire, UK
| | - David J M Lewis
- School of Biosciences and Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Jacek Sroka
- Institute of Informatics, University of Warsaw, Warsaw, Poland
| | - Hiroshi Matsuno
- Graduate School of Science and Engineering and Faculty of Science, Yamaguchi University, Yoshida, Yamaguchi, Japan
| | - Steffen Klamt
- Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg, Germany
| | - Hans V Westerhoff
- Molecular Cell Physiology, VU University Amsterdam, Amsterdam, The Netherlands
- Manchester Centre for Integrative Systems Biology, University of Manchester, Manchester, UK
- Synthetic Systems Biology, Netherlands Institute for Systems Biology, University of Amsterdam, Amsterdam, The Netherlands
| | - Johnjoe McFadden
- School of Biosciences and Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Nicholas J Plant
- School of Biosciences and Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Andrzej M Kierzek
- School of Biosciences and Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
- Simcyp Limited (a Certara Company), Sheffield, UK
| |
Collapse
|
46
|
Ostrowski M, Paulevé L, Schaub T, Siegel A, Guziolowski C. Boolean network identification from perturbation time series data combining dynamics abstraction and logic programming. Biosystems 2016; 149:139-153. [DOI: 10.1016/j.biosystems.2016.07.009] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2015] [Revised: 07/08/2016] [Accepted: 07/19/2016] [Indexed: 10/21/2022]
|
47
|
Model-Based Characterization of Inflammatory Gene Expression Patterns of Activated Macrophages. PLoS Comput Biol 2016; 12:e1005018. [PMID: 27464342 PMCID: PMC4963125 DOI: 10.1371/journal.pcbi.1005018] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2016] [Accepted: 06/08/2016] [Indexed: 12/14/2022] Open
Abstract
Macrophages are cells with remarkable plasticity. They integrate signals from their microenvironment leading to context-dependent polarization into classically (M1) or alternatively (M2) activated macrophages, representing two extremes of a broad spectrum of divergent phenotypes. Thereby, macrophages deliver protective and pro-regenerative signals towards injured tissue but, depending on the eliciting damage, may also be responsible for the generation and aggravation of tissue injury. Although incompletely understood, there is emerging evidence that macrophage polarization is critical for these antagonistic roles. To identify activation-specific expression patterns of chemokines and cytokines that may confer these distinct effects a systems biology approach was applied. A comprehensive literature-based Boolean model was developed to describe the M1 (LPS-activated) and M2 (IL-4/13-activated) polarization types. The model was validated using high-throughput transcript expression data from murine bone marrow derived macrophages. By dynamic modeling of gene expression, the chronology of pathway activation and autocrine signaling was estimated. Our results provide a deepened understanding of the physiological balance leading to M1/M2 activation, indicating the relevance of co-regulatory signals at the level of Akt1 or Akt2 that may be important for directing macrophage polarization. Macrophages are essential cells of the immune system and indispensable for a defense against bacterial infection. They reside as resting, immune modulatory cells in several tissues of the human body where they continuously sense inputs from their local environment. They react to stimuli such as toxins, injury or bacterial products in a process termed macrophage activation or polarization. For example, the bacterial component lipopolysaccharide induces so-called classical activation of macrophages into the M1 phenotype that secretes a number of inflammatory cytokines and chemokines leading to killing of bacteria and resolution of inflammation. Another prominent phenotype of macrophages is the M2 polarization state that is associated with wound healing and tissue regeneration. Unbalanced activation of macrophages is implicated in a number of diseases. An improved knowledge and extensive characterization of these macrophages as well as the factors determining their phenotypes will improve the understanding of the role of macrophages in disease progression.
Collapse
|
48
|
Node-independent elementary signaling modes: A measure of redundancy in Boolean signaling transduction networks. ACTA ACUST UNITED AC 2016. [DOI: 10.1017/nws.2016.4] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
AbstractThe redundancy of a system denotes the amount of duplicate components or mechanisms in it. For a network, especially one in which mass or information is being transferred from an origin to a destination, redundancy is related to the robustness of the system. Existing network measures of redundancy rely on local connectivity (e.g. clustering coefficients) or the existence of multiple paths. As in many systems there are functional dependencies between components and paths, a measure that not only characterizes the topology of a network, but also takes into account these functional dependencies, becomes most desirable.We propose a network redundancy measure in a prototypical model that contains functionally dependent directed paths: a Boolean model of a signal transduction network. The functional dependencies are made explicit by using an expanded network and the concept of elementary signaling modes (ESMs). We define the redundancy of a Boolean signal transduction network as the maximum number of node-independent ESMs and develop a methodology for identifying all maximal node-independent ESM combinations. We apply our measure to a number of signal transduction network models and show that it successfully distills known properties of the systems and offers new functional insights. The concept can be easily extended to similar related forms, e.g. edge-independent ESMs.
Collapse
|
49
|
Anastasio TJ. Temporal-logic analysis of microglial phenotypic conversion with exposure to amyloid-β. MOLECULAR BIOSYSTEMS 2016; 11:434-53. [PMID: 25406664 DOI: 10.1039/c4mb00457d] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Alzheimer Disease (AD) remains a leading killer with no adequate treatment. Ongoing research increasingly implicates the brain's immune system as a critical contributor to AD pathogenesis, but the complexity of the immune contribution poses a barrier to understanding. Here I use temporal logic to analyze a computational specification of the immune component of AD. Temporal logic is an extension of logic to propositions expressed in terms of time. It has traditionally been used to analyze computational specifications of complex engineered systems but applications to complex biological systems are now appearing. The inflammatory component of AD involves the responses of microglia to the peptide amyloid-β (Aβ), which is an inflammatory stimulus and a likely causative AD agent. Temporal-logic analysis of the model provides explanations for the puzzling findings that Aβ induces an anti-inflammatory and well as a pro-inflammatory response, and that Aβ is phagocytized by microglia in young but not in old animals. To potentially explain the first puzzle, the model suggests that interferon-γ acts as an "autocrine bridge" over which an Aβ-induced increase in pro-inflammatory cytokines leads to an increase in anti-inflammatory mediators also. To potentially explain the second puzzle, the model identifies a potential instability in signaling via insulin-like growth factor 1 that could explain the failure of old microglia to phagocytize Aβ. The model predicts that augmentation of insulin-like growth factor 1 signaling, and activation of protein kinase C in particular, could move old microglia from a neurotoxic back toward a more neuroprotective and phagocytic phenotype.
Collapse
Affiliation(s)
- Thomas J Anastasio
- Computational Neurobiology Laboratory, Department of Molecular and Integrative Physiology, and Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL, USA.
| |
Collapse
|
50
|
Paroni A, Graudenzi A, Caravagna G, Damiani C, Mauri G, Antoniotti M. CABeRNET: a Cytoscape app for augmented Boolean models of gene regulatory NETworks. BMC Bioinformatics 2016; 17:64. [PMID: 26846964 PMCID: PMC4743236 DOI: 10.1186/s12859-016-0914-z] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2015] [Accepted: 01/27/2016] [Indexed: 12/26/2022] Open
Abstract
BACKGROUND Dynamical models of gene regulatory networks (GRNs) are highly effective in describing complex biological phenomena and processes, such as cell differentiation and cancer development. Yet, the topological and functional characterization of real GRNs is often still partial and an exhaustive picture of their functioning is missing. RESULTS We here introduce CABERNET, a Cytoscape app for the generation, simulation and analysis of Boolean models of GRNs, specifically focused on their augmentation when a only partial topological and functional characterization of the network is available. By generating large ensembles of networks in which user-defined entities and relations are added to the original core, CABERNET allows to formulate hypotheses on the missing portions of real networks, as well to investigate their generic properties, in the spirit of complexity science. CONCLUSIONS CABERNET offers a series of innovative simulation and modeling functions and tools, including (but not being limited to) the dynamical characterization of the gene activation patterns ruling cell types and differentiation fates, and sophisticated robustness assessments, as in the case of gene knockouts. The integration within the widely used Cytoscape framework for the visualization and analysis of biological networks, makes CABERNET a new essential instrument for both the bioinformatician and the computational biologist, as well as a computational support for the experimentalist. An example application concerning the analysis of an augmented T-helper cell GRN is provided.
Collapse
Affiliation(s)
- Andrea Paroni
- Department of Informatics, Systems and Communication, University of Milan-Bicocca, Viale Sarca 336, Milan, I-20126, Italy
| | - Alex Graudenzi
- Department of Informatics, Systems and Communication, University of Milan-Bicocca, Viale Sarca 336, Milan, I-20126, Italy. .,Institute of Molecular Bioimaging and Physiology of the Italian National Research Council (IBFM-CNR), Via F.lli Cervi, 93, Segrate, I-20090, (MI), Italy.
| | - Giulio Caravagna
- Department of Informatics, Systems and Communication, University of Milan-Bicocca, Viale Sarca 336, Milan, I-20126, Italy. .,School of Informatics, University of Edinburgh, 10 Crichton St, Edinburgh, EH8 9AB, UK.
| | - Chiara Damiani
- Department of Informatics, Systems and Communication, University of Milan-Bicocca, Viale Sarca 336, Milan, I-20126, Italy. .,SYSBIO Centre of Systems Biology, Piazza della Scienza 2, 20126 Milano, Italy, Viale Sarca 336, Milan, I-20126, Italy.
| | - Giancarlo Mauri
- Department of Informatics, Systems and Communication, University of Milan-Bicocca, Viale Sarca 336, Milan, I-20126, Italy. .,Institute of Molecular Bioimaging and Physiology of the Italian National Research Council (IBFM-CNR), Via F.lli Cervi, 93, Segrate, I-20090, (MI), Italy. .,SYSBIO Centre of Systems Biology, Piazza della Scienza 2, 20126 Milano, Italy, Viale Sarca 336, Milan, I-20126, Italy.
| | - Marco Antoniotti
- Department of Informatics, Systems and Communication, University of Milan-Bicocca, Viale Sarca 336, Milan, I-20126, Italy. .,Milan Center for Neuroscience, University of Milan-Bicocca, Milan, Italy.
| |
Collapse
|