1
|
Wei L, Mast FD, Aitchison JD, Kaushansky A. Systems-level reconstruction of kinase phosphosignaling networks regulating endothelial barrier integrity using temporal data. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.08.01.606198. [PMID: 39149238 PMCID: PMC11326140 DOI: 10.1101/2024.08.01.606198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/17/2024]
Abstract
Phosphosignaling networks control cellular processes. We built kinase-mediated regulatory networks elicited by thrombin stimulation of brain endothelial cells using two computational strategies: Temporal Pathway Synthesizer (TPS), which uses phosphoproetiomics data as input, and Temporally REsolved KInase Network Generation (TREKING), which uses kinase inhibitor screens. TPS and TREKING predicted overlapping barrier-regulatory kinases connected with unique network topology. Each strategy effectively describes regulatory signaling networks and is broadly applicable across biological systems.
Collapse
Affiliation(s)
- Ling Wei
- Center for Global Infectious Disease Research, Seattle Children’s Research Institute, Seattle, WA 98109, United States
| | - Fred D. Mast
- Center for Global Infectious Disease Research, Seattle Children’s Research Institute, Seattle, WA 98109, United States
- Department of Pediatrics, University of Washington, Seattle, WA 98105, United States
| | - John D. Aitchison
- Center for Global Infectious Disease Research, Seattle Children’s Research Institute, Seattle, WA 98109, United States
- Department of Pediatrics, University of Washington, Seattle, WA 98105, United States
- Department of Biochemistry, University of Washington, Seattle, WA 98105, United States
| | - Alexis Kaushansky
- Center for Global Infectious Disease Research, Seattle Children’s Research Institute, Seattle, WA 98109, United States
- Department of Pediatrics, University of Washington, Seattle, WA 98105, United States
- Department of Global Health, University of Washington, Seattle, WA 98105, United States
| |
Collapse
|
2
|
Zhang D, Gao S, Liu ZP, Gao R. LogicGep: Boolean networks inference using symbolic regression from time-series transcriptomic profiling data. Brief Bioinform 2024; 25:bbae286. [PMID: 38886006 PMCID: PMC11182660 DOI: 10.1093/bib/bbae286] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2024] [Revised: 05/09/2024] [Accepted: 06/06/2024] [Indexed: 06/20/2024] Open
Abstract
Reconstructing the topology of gene regulatory network from gene expression data has been extensively studied. With the abundance functional transcriptomic data available, it is now feasible to systematically decipher regulatory interaction dynamics in a logic form such as a Boolean network (BN) framework, which qualitatively indicates how multiple regulators aggregated to affect a common target gene. However, inferring both the network topology and gene interaction dynamics simultaneously is still a challenging problem since gene expression data are typically noisy and data discretization is prone to information loss. We propose a new method for BN inference from time-series transcriptional profiles, called LogicGep. LogicGep formulates the identification of Boolean functions as a symbolic regression problem that learns the Boolean function expression and solve it efficiently through multi-objective optimization using an improved gene expression programming algorithm. To avoid overly emphasizing dynamic characteristics at the expense of topology structure ones, as traditional methods often do, a set of promising Boolean formulas for each target gene is evolved firstly, and a feed-forward neural network trained with continuous expression data is subsequently employed to pick out the final solution. We validated the efficacy of LogicGep using multiple datasets including both synthetic and real-world experimental data. The results elucidate that LogicGep adeptly infers accurate BN models, outperforming other representative BN inference algorithms in both network topology reconstruction and the identification of Boolean functions. Moreover, the execution of LogicGep is hundreds of times faster than other methods, especially in the case of large network inference.
Collapse
Affiliation(s)
- Dezhen Zhang
- Center of Intelligent Medicine, School of Control Science and Engineering, Shandong University, Jinan, Shandong 250061, China
| | - Shuhua Gao
- Center of Intelligent Medicine, School of Control Science and Engineering, Shandong University, Jinan, Shandong 250061, China
| | - Zhi-Ping Liu
- Center of Intelligent Medicine, School of Control Science and Engineering, Shandong University, Jinan, Shandong 250061, China
| | - Rui Gao
- Center of Intelligent Medicine, School of Control Science and Engineering, Shandong University, Jinan, Shandong 250061, China
| |
Collapse
|
3
|
Meimetis N, Lauffenburger DA, Nilsson A. Inference of drug off-target effects on cellular signaling using interactome-based deep learning. iScience 2024; 27:109509. [PMID: 38591003 PMCID: PMC11000001 DOI: 10.1016/j.isci.2024.109509] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2023] [Revised: 02/04/2024] [Accepted: 03/13/2024] [Indexed: 04/10/2024] Open
Abstract
Many diseases emerge from dysregulated cellular signaling, and drugs are often designed to target specific signaling proteins. Off-target effects are, however, common and may ultimately result in failed clinical trials. Here we develop a computer model of the cell's transcriptional response to drugs for improved understanding of their mechanisms of action. The model is based on ensembles of artificial neural networks and simultaneously infers drug-target interactions and their downstream effects on intracellular signaling. With this, it predicts transcription factors' activities, while recovering known drug-target interactions and inferring many new ones, which we validate with an independent dataset. As a case study, we analyze the effects of the drug Lestaurtinib on downstream signaling. Alongside its intended target, FLT3, the model predicts an inhibition of CDK2 that enhances the downregulation of the cell cycle-critical transcription factor FOXM1. Our approach can therefore enhance our understanding of drug signaling for therapeutic design.
Collapse
Affiliation(s)
- Nikolaos Meimetis
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Douglas A. Lauffenburger
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Avlant Nilsson
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
- Department of Cell and Molecular Biology, SciLifeLab, Karolinska Institutet, Stockholm, Sweden
- Department of Biology and Biological Engineering, Chalmers University of Technology, Gothenburg, SE 41296, Sweden
| |
Collapse
|
4
|
Patidar K, Versypt ANF. Logic-Based Modeling of Inflammatory Macrophage Crosstalk with Glomerular Endothelial Cells in Diabetic Kidney Disease. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.04.04.535594. [PMID: 37066138 PMCID: PMC10104015 DOI: 10.1101/2023.04.04.535594] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Diabetic kidney disease is a complication in 1 out of 3 patients with diabetes. Aberrant glucose metabolism in diabetes leads to an immune response causing inflammation and to structural and functional damage in the glomerular cells of the kidney. Complex cellular signaling lies at the core of metabolic and functional derangement. Unfortunately, the mechanism underlying the role of inflammation in glomerular endothelial cell dysfunction during diabetic kidney disease is not fully understood. Computational models in systems biology allow the integration of experimental evidence and cellular signaling networks to understand mechanisms involved in disease progression. We built a logic-based ordinary differential equations model to study macrophage-dependent inflammation in glomerular endothelial cells during diabetic kidney disease progression. We studied the crosstalk between macrophages and glomerular endothelial cells in the kidney using a protein signaling network stimulated with glucose and lipopolysaccharide. The network and model were built using the open-source software package Netflux. This modeling approach overcomes the complexity of studying network models and the need for extensive mechanistic details. The model simulations were fitted and validated against available biochemical data from in vitro experiments. The model identified mechanisms responsible for dysregulated signaling in macrophages and glomerular endothelial cells during diabetic kidney disease. In addition, we investigated the influence of signaling interactions and species that on glomerular endothelial cell morphology through selective knockdown and downregulation. We found that partial knockdown of VEGF receptor 1, PLC-γ, adherens junction proteins, and calcium partially recovered the endothelial cell fenestration size. Our model findings contribute to understanding signaling and molecular perturbations that affect the glomerular endothelial cells in the early stage of diabetic kidney disease.
Collapse
|
5
|
Esser-Skala W, Fortelny N. Reliable interpretability of biology-inspired deep neural networks. NPJ Syst Biol Appl 2023; 9:50. [PMID: 37816807 PMCID: PMC10564878 DOI: 10.1038/s41540-023-00310-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2023] [Accepted: 09/15/2023] [Indexed: 10/12/2023] Open
Abstract
Deep neural networks display impressive performance but suffer from limited interpretability. Biology-inspired deep learning, where the architecture of the computational graph is based on biological knowledge, enables unique interpretability where real-world concepts are encoded in hidden nodes, which can be ranked by importance and thereby interpreted. In such models trained on single-cell transcriptomes, we previously demonstrated that node-level interpretations lack robustness upon repeated training and are influenced by biases in biological knowledge. Similar studies are missing for related models. Here, we test and extend our methodology for reliable interpretability in P-NET, a biology-inspired model trained on patient mutation data. We observe variability of interpretations and susceptibility to knowledge biases, and identify the network properties that drive interpretation biases. We further present an approach to control the robustness and biases of interpretations, which leads to more specific interpretations. In summary, our study reveals the broad importance of methods to ensure robust and bias-aware interpretability in biology-inspired deep learning.
Collapse
Affiliation(s)
- Wolfgang Esser-Skala
- Computational Systems Biology Group, Department of Biosciences and Medical Biology, University of Salzburg, Hellbrunner Straße 34, 5020, Salzburg, Austria
| | - Nikolaus Fortelny
- Computational Systems Biology Group, Department of Biosciences and Medical Biology, University of Salzburg, Hellbrunner Straße 34, 5020, Salzburg, Austria.
| |
Collapse
|
6
|
Glazer BJ, Lifferth JT, Lopez CF. Automatic mechanistic inference from large families of Boolean models generated by Monte Carlo tree search. Front Cell Dev Biol 2023; 11:1198359. [PMID: 37691824 PMCID: PMC10485623 DOI: 10.3389/fcell.2023.1198359] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2023] [Accepted: 08/07/2023] [Indexed: 09/12/2023] Open
Abstract
Many important processes in biology, such as signaling and gene regulation, can be described using logic models. These logic models are typically built to behaviorally emulate experimentally observed phenotypes, which are assumed to be steady states of a biological system. Most models are built by hand and therefore researchers are only able to consider one or perhaps a few potential mechanisms. We present a method to automatically synthesize Boolean logic models with a specified set of steady states. Our method, called MC-Boomer, is based on Monte Carlo Tree Search an efficient, parallel search method using reinforcement learning. Our approach enables users to constrain the model search space using prior knowledge or biochemical interaction databases, thus leading to generation of biologically plausible mechanistic hypotheses. Our approach can generate very large numbers of data-consistent models. To help develop mechanistic insight from these models, we developed analytical tools for multi-model inference and model selection. These tools reveal the key sets of interactions that govern the behavior of the models. We demonstrate that MC-Boomer works well at reconstructing randomly generated models. Then, using single time point measurements and reasonable biological constraints, our method generates hundreds of thousands of candidate models that match experimentally validated in-vivo behaviors of the Drosophila segment polarity network. Finally we outline how our multi-model analysis procedures elucidate potentially novel biological mechanisms and provide opportunities for model-driven experimental validation.
Collapse
Affiliation(s)
- Bryan J. Glazer
- Department of Biomedical Informatics, Vanderbilt University, Nashville, TN, United States
| | - Jonathan T. Lifferth
- Department of Human Genetics, Vanderbilt University, Nashville, TN, United States
| | - Carlos F. Lopez
- Department of Biochemistry, Vanderbilt University, Nashville, TN, United States
- Altos Labs, Redwood City, CA, United States
| |
Collapse
|
7
|
Bachman JA, Gyori BM, Sorger PK. Automated assembly of molecular mechanisms at scale from text mining and curated databases. Mol Syst Biol 2023; 19:e11325. [PMID: 36938926 PMCID: PMC10167483 DOI: 10.15252/msb.202211325] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2022] [Revised: 02/24/2023] [Accepted: 02/27/2023] [Indexed: 03/21/2023] Open
Abstract
The analysis of omic data depends on machine-readable information about protein interactions, modifications, and activities as found in protein interaction networks, databases of post-translational modifications, and curated models of gene and protein function. These resources typically depend heavily on human curation. Natural language processing systems that read the primary literature have the potential to substantially extend knowledge resources while reducing the burden on human curators. However, machine-reading systems are limited by high error rates and commonly generate fragmentary and redundant information. Here, we describe an approach to precisely assemble molecular mechanisms at scale using multiple natural language processing systems and the Integrated Network and Dynamical Reasoning Assembler (INDRA). INDRA identifies full and partial overlaps in information extracted from published papers and pathway databases, uses predictive models to improve the reliability of machine reading, and thereby assembles individual pieces of information into non-redundant and broadly usable mechanistic knowledge. Using INDRA to create high-quality corpora of causal knowledge we show it is possible to extend protein-protein interaction databases and explain co-dependencies in the Cancer Dependency Map.
Collapse
Affiliation(s)
- John A Bachman
- Laboratory of Systems Pharmacology, Harvard Medical School, Boston, MA, USA
| | - Benjamin M Gyori
- Laboratory of Systems Pharmacology, Harvard Medical School, Boston, MA, USA
| | - Peter K Sorger
- Laboratory of Systems Pharmacology, Harvard Medical School, Boston, MA, USA.,Department of Systems Biology, Harvard Medical School, Boston, MA, USA
| |
Collapse
|
8
|
Lo Surdo P, Iannuccelli M, Contino S, Castagnoli L, Licata L, Cesareni G, Perfetto L. SIGNOR 3.0, the SIGnaling network open resource 3.0: 2022 update. Nucleic Acids Res 2022; 51:D631-D637. [PMID: 36243968 PMCID: PMC9825604 DOI: 10.1093/nar/gkac883] [Citation(s) in RCA: 34] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Revised: 09/21/2022] [Accepted: 09/30/2022] [Indexed: 01/29/2023] Open
Abstract
The SIGnaling Network Open Resource (SIGNOR 3.0, https://signor.uniroma2.it) is a public repository that captures causal information and represents it according to an 'activity-flow' model. SIGNOR provides freely-accessible static maps of causal interactions that can be tailored, pruned and refined to build dynamic and predictive models. Each signaling relationship is annotated with an effect (up/down-regulation) and with the mechanism (e.g. binding, phosphorylation, transcriptional activation, etc.) causing the regulation of the target entity. Since its latest release, SIGNOR has undergone a significant upgrade including: (i) a new website that offers an improved user experience and novel advanced search and graph tools; (ii) a significant content growth adding up to a total of approx. 33,000 manually-annotated causal relationships between more than 8900 biological entities; (iii) an increase in the number of manually annotated pathways, currently including pathways deregulated by SARS-CoV-2 infection or involved in neurodevelopment synaptic transmission and metabolism, among others; (iv) additional features such as new model to represent metabolic reactions and a new confidence score assigned to each interaction.
Collapse
Affiliation(s)
| | - Marta Iannuccelli
- Department of Biology, University of Rome ‘Tor Vergata’, Rome 00133, Italy
| | - Silvia Contino
- Department of Biology, University of Rome ‘Tor Vergata’, Rome 00133, Italy
| | - Luisa Castagnoli
- Department of Biology, University of Rome ‘Tor Vergata’, Rome 00133, Italy
| | | | | | - Livia Perfetto
- To whom correspondence should be addressed. Tel: +39 0672594315;
| |
Collapse
|
9
|
Thuillier K, Baroukh C, Bockmayr A, Cottret L, Paulevé L, Siegel A. MERRIN: MEtabolic regulation rule INference from time series data. Bioinformatics 2022; 38:ii127-ii133. [PMID: 36124795 DOI: 10.1093/bioinformatics/btac479] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
MOTIVATION Many techniques have been developed to infer Boolean regulations from a prior knowledge network (PKN) and experimental data. Existing methods are able to reverse-engineer Boolean regulations for transcriptional and signaling networks, but they fail to infer regulations that control metabolic networks. RESULTS We present a novel approach to infer Boolean rules for metabolic regulation from time-series data and a PKN. Our method is based on a combination of answer set programming and linear programming. By solving both combinatorial and linear arithmetic constraints, we generate candidate Boolean regulations that can reproduce the given data when coupled to the metabolic network. We evaluate our approach on a core regulated metabolic network and show how the quality of the predictions depends on the available kinetic, fluxomics or transcriptomics time-series data. AVAILABILITY AND IMPLEMENTATION Software available at https://github.com/bioasp/merrin. SUPPLEMENTARY INFORMATION Supplementary data are available at https://doi.org/10.5281/zenodo.6670164.
Collapse
Affiliation(s)
- Kerian Thuillier
- INRIA, CNRS, IRISA, University of Rennes, Rennes F-35000, France
| | - Caroline Baroukh
- LIPME, INRAE, CNRS, Université de Toulouse, Castanet-Tolosan F-31326, France
| | - Alexander Bockmayr
- Institute of Mathematics, Freie Universität Berlin, Berlin D-14195, Germany
| | - Ludovic Cottret
- LIPME, INRAE, CNRS, Université de Toulouse, Castanet-Tolosan F-31326, France
| | - Loïc Paulevé
- Univ. Bordeaux, Bordeaux INP, CNRS, LaBRI, UMR5800, F-33400 Talence, France
| | - Anne Siegel
- INRIA, CNRS, IRISA, University of Rennes, Rennes F-35000, France
| |
Collapse
|
10
|
Müssel C, Ikonomi N, Werle SD, Weidner FM, Maucher M, Schwab JD, Kestler HA. CANTATA - prediction of missing links in Boolean networks using genetic programming. Bioinformatics 2022; 38:4893-4900. [PMID: 36094334 PMCID: PMC9620829 DOI: 10.1093/bioinformatics/btac623] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Revised: 08/25/2022] [Accepted: 09/09/2022] [Indexed: 11/27/2022] Open
Abstract
Motivation Biological processes are complex systems with distinct behaviour. Despite the growing amount of available data, knowledge is sparse and often insufficient to investigate the complex regulatory behaviour of these systems. Moreover, different cellular phenotypes are possible under varying conditions. Mathematical models attempt to unravel these mechanisms by investigating the dynamics of regulatory networks. Therefore, a major challenge is to combine regulations and phenotypical information as well as the underlying mechanisms. To predict regulatory links in these models, we established an approach called CANTATA to support the integration of information into regulatory networks and retrieve potential underlying regulations. This is achieved by optimizing both static and dynamic properties of these networks. Results Initial results show that the algorithm predicts missing interactions by recapitulating the known phenotypes while preserving the original topology and optimizing the robustness of the model. The resulting models allow for hypothesizing about the biological impact of certain regulatory dependencies. Availability and implementation Source code of the application, example files and results are available at https://github.com/sysbio-bioinf/Cantata. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Christoph Müssel
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - Nensi Ikonomi
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - Silke D Werle
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - Felix M Weidner
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - Markus Maucher
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - Julian D Schwab
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - Hans A Kestler
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| |
Collapse
|
11
|
Lee J, Choi SR, Cho KH. Network Dynamics Caused by Genomic Alteration Determine the Therapeutic Response to FGFR Inhibitors for Lung Cancer. Biomolecules 2022; 12:biom12091197. [PMID: 36139037 PMCID: PMC9496101 DOI: 10.3390/biom12091197] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2022] [Revised: 08/13/2022] [Accepted: 08/25/2022] [Indexed: 11/16/2022] Open
Abstract
Recently, FGFR inhibitors have been highlighted as promising targeted drugs due to the high prevalence of FGFR1 amplification in cancer patients. Although various potential biomarkers for FGFR inhibitors have been suggested, their functional effects have been shown to be limited due to the complexity of the cancer signaling network and the heterogenous genomic conditions of patients. To overcome such limitations, we have reconstructed a lung cancer network model by integrating a cell line genomic database and analyzing the model in order to understand the underlying mechanism of heterogeneous drug responses. Here, we identify novel genomic context-specific candidates that can increase the efficacy of FGFR inhibitors. Furthermore, we suggest optimal targets that can induce more effective therapeutic responses than that of FGFR inhibitors in each of the FGFR-resistant lung cancer cells through computational simulations at a system level. Our findings provide new insights into the regulatory mechanism of differential responses to FGFR inhibitors for optimal therapeutic strategies in lung cancer.
Collapse
Affiliation(s)
| | | | - Kwang-Hyun Cho
- Correspondence: ; Tel.: +82-42-350-4325; Fax: +82-42-350-4310
| |
Collapse
|
12
|
ELIMINATOR: essentiality analysis using multisystem networks and integer programming. BMC Bioinformatics 2022; 23:324. [PMID: 35933325 PMCID: PMC9357337 DOI: 10.1186/s12859-022-04855-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2022] [Accepted: 07/21/2022] [Indexed: 11/28/2022] Open
Abstract
A gene is considered as essential when it is indispensable for cells to grow and replicate in a certain environment. However, gene essentiality is not a structural property but rather a contextual one, which depends on the specific biological conditions affecting the cell. This circumstantial essentiality of genes is what brings the attention of scientist since we can identify genes essential for cancer cells but not essential for healthy cells. This same contextuality makes their identification extremely challenging. Huge experimental efforts such as Project Achilles where the essentiality of thousands of genes is measured together with a plethora of molecular data (transcriptomics, copy number, mutations, etc.) in over one thousand cell lines can shed light on the causality behind the essentiality of a gene in a given environment. Here, we present an in-silico method for the identification of patient-specific essential genes using constraint-based modelling (CBM). Our method expands the ideas behind traditional CBM to accommodate multisystem networks. In essence, it first calculates the minimum number of lowly expressed genes required to be activated by the cell to sustain life as defined by a set of requirements; and second, it performs an exhaustive in-silico gene knockout to find those that lead to the need of activating additional lowly expressed genes. We validated the proposed methodology using a set of 452 cancer cell lines derived from the Cancer Cell Line Encyclopedia where an exhaustive experimental large-scale gene knockout study using CRISPR (Achilles Project) evaluates the impact of each removal. We also show that the integration of different essentiality predictions per gene, what we called Essentiality Congruity Score, reduces the number of false positives. Finally, we explored our method in a breast cancer patient dataset, and our results showed high concordance with previous publications. These findings suggest that identifying genes whose activity is fundamental to sustain cellular life in a patient-specific manner is feasible using in-silico methods. The patient-level gene essentiality predictions can pave the way for precision medicine by identifying potential drug targets whose deletion can induce death in tumour cells.
Collapse
|
13
|
Emadi A, Ozen M, Abdi A. A hybrid model to study how late long-term potentiation is affected by faulty molecules in an intraneuronal signaling network regulating transcription factor CREB. INTEGRATIVE BIOLOGY : QUANTITATIVE BIOSCIENCES FROM NANO TO MACRO 2022; 14:111-125. [PMID: 35901510 DOI: 10.1093/intbio/zyac011] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/24/2022] [Revised: 06/01/2022] [Accepted: 07/04/2022] [Indexed: 11/12/2022]
Abstract
Systems biology analysis of intracellular signaling networks has tremendously expanded our understanding of normal and diseased cell behaviors and has revealed paths to finding proper therapeutic molecular targets. When it comes to neurons in the human brain, analysis of intraneuronal signaling networks provides invaluable information on learning, memory and cognition-related disorders, as well as potential therapeutic targets. However, neurons in the human brain form a highly complex neural network that, among its many roles, is also responsible for learning, memory formation and cognition. Given the impairment of these processes in mental and psychiatric disorders, one can envision that analyzing interneuronal processes, together with analyzing intraneuronal signaling networks, can result in a better understanding of the pathology and, subsequently, more effective target discovery. In this paper, a hybrid model is introduced, composed of the long-term potentiation (LTP) interneuronal process and an intraneuronal signaling network regulating CREB. LTP refers to an increased synaptic strength over a long period of time among neurons, typically induced upon occurring an activity that generates high-frequency stimulations (HFS) in the brain, and CREB is a transcription factor known to be highly involved in important functions of the cognitive and executive human brain such as learning and memory. The hybrid LTP-signaling model is analyzed using a proposed molecular fault diagnosis method. It allows to study the importance of various signaling molecules according to how much they affect an intercellular phenomenon when they are faulty, i.e. dysfunctional. This paper is intended to suggest another angle for understanding the pathology and therapeutic target discovery by classifying and ranking various intraneuronal signaling molecules based on how much their faulty behaviors affect an interneuronal process. Possible relations between the introduced hybrid analysis and the previous purely intracellular analysis are investigated in the paper as well.
Collapse
Affiliation(s)
- Ali Emadi
- Center for Wireless Information Processing, Dept. Electrical & Computer Engineering, New Jersey Institute of Technology (NJIT), Newark, NJ 07102, USA
| | - Mustafa Ozen
- Center for Wireless Information Processing, Dept. Electrical & Computer Engineering, New Jersey Institute of Technology (NJIT), Newark, NJ 07102, USA.,Dept. of Biochemistry, Vanderbilt University, Nashville, TN 37240, USA
| | - Ali Abdi
- Center for Wireless Information Processing, Dept. Electrical & Computer Engineering, New Jersey Institute of Technology (NJIT), Newark, NJ 07102, USA.,Dept. of Biological Sciences, New Jersey Institute of Technology (NJIT), Newark, NJ 07102, USA
| |
Collapse
|
14
|
Bosdriesz E, Fernandes Neto JM, Sieber A, Bernards R, Blüthgen N, Wessels LF. Identifying mutant-specific multi-drug combinations using Comparative Network Reconstruction. iScience 2022; 25:104760. [PMID: 35992065 PMCID: PMC9385552 DOI: 10.1016/j.isci.2022.104760] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2022] [Revised: 05/30/2022] [Accepted: 07/11/2022] [Indexed: 10/28/2022] Open
|
15
|
Bhavani GS, Palanisamy A. SNAIL driven by a feed forward loop motif promotes TGF βinduced epithelial to mesenchymal transition. Biomed Phys Eng Express 2022; 8. [PMID: 35700712 DOI: 10.1088/2057-1976/ac7896] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2022] [Accepted: 06/14/2022] [Indexed: 11/12/2022]
Abstract
Epithelial to Mesenchymal Transition (EMT) plays an important role in tissue regeneration, embryonic development, and cancer metastasis. Several signaling pathways are known to regulate EMT, among which the modulation of TGFβ(Transforming Growth Factor-β) induced EMT is crucial in several cancer types. Several mathematical models were built to explore the role of core regulatory circuit of ZEB/miR-200, SNAIL/miR-34 double negative feedback loops in modulating TGFβinduced EMT. Different emergent behavior including tristability, irreversible switching, existence of hybrid EMT states were inferred though these models. Some studies have explored the role of TGFβreceptor activation, SMADs nucleocytoplasmic shuttling and complex formation. Recent experiments have revealed that MDM2 along with SMAD complex regulates SNAIL expression driven EMT. Encouraged by this, in the present study we developed a mathematical model for p53/MDM2 dependent TGFβinduced EMT regulation. Inclusion of p53 brings in an additional mechanistic perspective in exploring the EM transition. The network formulated comprises a C1FFL moderating SNAIL expression involving MDM2 and SMAD complex, which functions as a noise filter and persistent detector. The C1FFL was also observed to operate as a coincidence detector driving the SNAIL dependent downstream signaling into phenotypic switching decision. Systems modelling and analysis of the devised network, displayed interesting dynamic behavior, systems response to various inputs stimulus, providing a better understanding of p53/MDM2 dependent TGF-βinduced Epithelial to Mesenchymal Transition.
Collapse
|
16
|
A probabilistic Boolean model on hair follicle cell fate regulation by TGF-β. Biophys J 2022; 121:2638-2652. [PMID: 35714600 DOI: 10.1016/j.bpj.2022.05.035] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2021] [Revised: 05/20/2022] [Accepted: 05/23/2022] [Indexed: 11/24/2022] Open
Abstract
Hair follicles (HFs) are mini skin organs that undergo cyclic growth. Various signals regulate HF cell fate decisions jointly. Recent experimental results suggest that transforming growth factor beta (TGF-β) exhibits a dual role in HF cell fate regulation that can be either anti- or pro-apoptosis. To understand the underlying mechanisms of HF cell fate control, we develop a novel probabilistic Boolean network (pBN) model on the HF epithelial cell gene regulation dynamics. First, the model is derived from literature, then refined using single-cell RNA sequencing data. Using the model, we both explore the mechanisms underlying HF cell fate decisions and make predictions that could potentially guide future experiments: 1) we propose that a threshold-like switch in the TGF-β strength may necessitate the dual roles of TGF-β in either activating apoptosis or cell proliferation, in cooperation with Bmp and tumor necrosis factor (TNF) and at different stages of a follicle growth cycle; 2) our model shows concordance with the high-activator-low-inhibitor theory of anagen initiation; 3) we predict that TNF may be more effective in catagen initiation than TGF-β, and they may cooperate in a two-step fashion; 4) finally, predictions of gene knockout and overexpression reveal the roles in HF cell fate regulations of each gene. Attractor and motif analysis from the associated Boolean networks reveal the relations between the topological structure of the gene regulation network and the cell fate regulation mechanism. A discrete spatial model equipped with the pBN illustrates how TGF-β and TNF cooperate in initiating and driving the apoptosis wave during catagen.
Collapse
|
17
|
Ozen M, S.Emamian E, Abdi A. Exploring extreme signaling failures in intracellular molecular networks. Comput Biol Med 2022; 148:105692. [DOI: 10.1016/j.compbiomed.2022.105692] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2022] [Revised: 05/29/2022] [Accepted: 05/31/2022] [Indexed: 11/27/2022]
|
18
|
Post JN, Loerakker S, Merks R, Carlier A. Implementing computational modeling in tissue engineering: where disciplines meet. Tissue Eng Part A 2022; 28:542-554. [PMID: 35345902 DOI: 10.1089/ten.tea.2021.0215] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
In recent years, the mathematical and computational sciences have developed novel methodologies and insights that can aid in designing advanced bioreactors, microfluidic set-ups or organ-on-chip devices, in optimizing culture conditions, or predicting long-term behavior of engineered tissues in vivo. In this review, we introduce the concept of computational models and how they can be integrated in an interdisciplinary workflow for Tissue Engineering and Regenerative Medicine (TERM). We specifically aim this review of general concepts and examples at experimental scientists with little or no computational modeling experience. We also describe the contribution of computational models in understanding TERM processes and in advancing the TERM field by providing novel insights.
Collapse
Affiliation(s)
- Janine Nicole Post
- University of Twente, 3230, Tissue Regeneration, Enschede, Overijssel, Netherlands;
| | - Sandra Loerakker
- Eindhoven University of Technology, 3169, Department of Biomedical Engineering, Eindhoven, Noord-Brabant, Netherlands.,Eindhoven University of Technology, 3169, Institute for Complex Molecular Systems, Eindhoven, Noord-Brabant, Netherlands;
| | - Roeland Merks
- Leiden University, 4496, Institute for Biology Leiden and Mathematical Institute, Leiden, Zuid-Holland, Netherlands;
| | - Aurélie Carlier
- Maastricht University, 5211, MERLN Institute for Technology-Inspired Regenerative Medicine, Universiteitssingel 40, 6229 ER Maastricht, Maastricht, Netherlands, 6200 MD;
| |
Collapse
|
19
|
Montagud A, Béal J, Tobalina L, Traynard P, Subramanian V, Szalai B, Alföldi R, Puskás L, Valencia A, Barillot E, Saez-Rodriguez J, Calzone L. Patient-specific Boolean models of signalling networks guide personalised treatments. eLife 2022; 11:e72626. [PMID: 35164900 PMCID: PMC9018074 DOI: 10.7554/elife.72626] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2021] [Accepted: 01/27/2022] [Indexed: 11/22/2022] Open
Abstract
Prostate cancer is the second most occurring cancer in men worldwide. To better understand the mechanisms of tumorigenesis and possible treatment responses, we developed a mathematical model of prostate cancer which considers the major signalling pathways known to be deregulated. We personalised this Boolean model to molecular data to reflect the heterogeneity and specific response to perturbations of cancer patients. A total of 488 prostate samples were used to build patient-specific models and compared to available clinical data. Additionally, eight prostate cell line-specific models were built to validate our approach with dose-response data of several drugs. The effects of single and combined drugs were tested in these models under different growth conditions. We identified 15 actionable points of interventions in one cell line-specific model whose inactivation hinders tumorigenesis. To validate these results, we tested nine small molecule inhibitors of five of those putative targets and found a dose-dependent effect on four of them, notably those targeting HSP90 and PI3K. These results highlight the predictive power of our personalised Boolean models and illustrate how they can be used for precision oncology.
Collapse
Affiliation(s)
- Arnau Montagud
- Institut Curie, PSL Research UniversityParisFrance
- INSERM, U900ParisFrance
- MINES ParisTech, PSL Research University, CBIO-Centre for Computational BiologyParisFrance
- Barcelona Supercomputing Center (BSC), Plaça Eusebi Güell, 1-3BarcelonaSpain
| | - Jonas Béal
- Institut Curie, PSL Research UniversityParisFrance
- INSERM, U900ParisFrance
- MINES ParisTech, PSL Research University, CBIO-Centre for Computational BiologyParisFrance
| | - Luis Tobalina
- Faculty of Medicine, Joint Research Centre for Computational Biomedicine (JRC-COMBINE), RWTH Aachen UniversityAachenGermany
| | - Pauline Traynard
- Institut Curie, PSL Research UniversityParisFrance
- INSERM, U900ParisFrance
- MINES ParisTech, PSL Research University, CBIO-Centre for Computational BiologyParisFrance
| | - Vigneshwari Subramanian
- Faculty of Medicine, Joint Research Centre for Computational Biomedicine (JRC-COMBINE), RWTH Aachen UniversityAachenGermany
| | - Bence Szalai
- Faculty of Medicine, Joint Research Centre for Computational Biomedicine (JRC-COMBINE), RWTH Aachen UniversityAachenGermany
- Semmelweis University, Faculty of Medicine, Department of PhysiologyBudapestHungary
| | | | | | - Alfonso Valencia
- Barcelona Supercomputing Center (BSC), Plaça Eusebi Güell, 1-3BarcelonaSpain
- ICREA, Pg. Lluís Companys 23BarcelonaSpain
| | - Emmanuel Barillot
- Institut Curie, PSL Research UniversityParisFrance
- INSERM, U900ParisFrance
- MINES ParisTech, PSL Research University, CBIO-Centre for Computational BiologyParisFrance
| | - Julio Saez-Rodriguez
- Faculty of Medicine, Joint Research Centre for Computational Biomedicine (JRC-COMBINE), RWTH Aachen UniversityAachenGermany
- Faculty of Medicine and Heidelberg University Hospital, Institute of Computational Biomedicine, Heidelberg UniversityHeidelbergGermany
| | - Laurence Calzone
- Institut Curie, PSL Research UniversityParisFrance
- INSERM, U900ParisFrance
- MINES ParisTech, PSL Research University, CBIO-Centre for Computational BiologyParisFrance
| |
Collapse
|
20
|
Putnins M, Campagne O, Mager DE, Androulakis IP. From data to QSP models: a pipeline for using Boolean networks for hypothesis inference and dynamic model building. J Pharmacokinet Pharmacodyn 2022; 49:101-115. [PMID: 34988912 PMCID: PMC9876619 DOI: 10.1007/s10928-021-09797-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2021] [Accepted: 10/27/2021] [Indexed: 01/27/2023]
Abstract
Quantitative Systems Pharmacology (QSP) models capture the physiological underpinnings driving the response to a drug and express those in a semi-mechanistic way, often involving ordinary differential equations (ODEs). The process of developing a QSP model generally starts with the definition of a set of reasonable hypotheses that would support a mechanistic interpretation of the expected response which are used to form a network of interacting elements. This is a hypothesis-driven and knowledge-driven approach, relying on prior information about the structure of the network. However, with recent advances in our ability to generate large datasets rapidly, often in a hypothesis-neutral manner, the opportunity emerges to explore data-driven approaches to establish the network topologies and models in a robust, repeatable manner. In this paper, we explore the possibility of developing complex network representations of physiological responses to pharmaceuticals using a logic-based analysis of available data and then convert the logic relations to dynamic ODE-based models. We discuss an integrated pipeline for converting data to QSP models. This pipeline includes using k-means clustering to binarize continuous data, inferring likely network relationships using a Best-Fit Extension method to create a Boolean network, and finally converting the Boolean network to a continuous ODE model. We utilized an existing QSP model for the dual-affinity re-targeting antibody flotetuzumab to demonstrate the robustness of the process. Key output variables from the QSP model were used to generate a continuous data set for use in the pipeline. This dataset was used to reconstruct a possible model. This reconstruction had no false-positive relationships, and the output of each of the species was similar to that of the original QSP model. This demonstrates the ability to accurately infer relationships in a hypothesis-neutral manner without prior knowledge of a system using this pipeline.
Collapse
Affiliation(s)
- M. Putnins
- Biomedical Engineering Department, Rutgers University, Piscataway, USA
| | - O. Campagne
- Department of Pharmaceutical Sciences, University at Buffalo, State University of New York, Buffalo, USA
| | - D. E. Mager
- Department of Pharmaceutical Sciences, University at Buffalo, State University of New York, Buffalo, USA
| | - I. P. Androulakis
- Biomedical Engineering Department, Rutgers University, Piscataway, USA,Chemical & Biochemical Engineering Department, Rutgers University, Piscataway, USA
| |
Collapse
|
21
|
Boolean function metrics can assist modelers to check and choose logical rules. J Theor Biol 2022; 538:111025. [DOI: 10.1016/j.jtbi.2022.111025] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2021] [Revised: 12/07/2021] [Accepted: 01/10/2022] [Indexed: 12/25/2022]
|
22
|
Jablonski KP, Pirkl M, Ćevid D, Bühlmann P, Beerenwinkel N. Identifying cancer pathway dysregulations using differential causal effects. Bioinformatics 2021; 38:1550-1559. [PMID: 34927666 PMCID: PMC8896597 DOI: 10.1093/bioinformatics/btab847] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2021] [Revised: 11/05/2021] [Accepted: 12/14/2021] [Indexed: 02/03/2023] Open
Abstract
MOTIVATION Signaling pathways control cellular behavior. Dysregulated pathways, for example, due to mutations that cause genes and proteins to be expressed abnormally, can lead to diseases, such as cancer. RESULTS We introduce a novel computational approach, called Differential Causal Effects (dce), which compares normal to cancerous cells using the statistical framework of causality. The method allows to detect individual edges in a signaling pathway that are dysregulated in cancer cells, while accounting for confounding. Hence, technical artifacts have less influence on the results and dce is more likely to detect the true biological signals. We extend the approach to handle unobserved dense confounding, where each latent variable, such as, for example, batch effects or cell cycle states, affects many covariates. We show that dce outperforms competing methods on synthetic datasets and on CRISPR knockout screens. We validate its latent confounding adjustment properties on a GTEx (Genotype-Tissue Expression) dataset. Finally, in an exploratory analysis on breast cancer data from TCGA (The Cancer Genome Atlas), we recover known and discover new genes involved in breast cancer progression. AVAILABILITY AND IMPLEMENTATION The method dce is freely available as an R package on Bioconductor (https://bioconductor.org/packages/release/bioc/html/dce.html) as well as on https://github.com/cbg-ethz/dce. The GitHub repository also contains the Snakemake workflows needed to reproduce all results presented here. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
| | | | - Domagoj Ćevid
- Seminar for Statistics, ETH Zürich, 8092 Zürich, Switzerland
| | - Peter Bühlmann
- Seminar for Statistics, ETH Zürich, 8092 Zürich, Switzerland
| | | |
Collapse
|
23
|
Rocca A, Kholodenko BN. Can Systems Biology Advance Clinical Precision Oncology? Cancers (Basel) 2021; 13:6312. [PMID: 34944932 PMCID: PMC8699328 DOI: 10.3390/cancers13246312] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2021] [Accepted: 12/10/2021] [Indexed: 12/13/2022] Open
Abstract
Precision oncology is perceived as a way forward to treat individual cancer patients. However, knowing particular cancer mutations is not enough for optimal therapeutic treatment, because cancer genotype-phenotype relationships are nonlinear and dynamic. Systems biology studies the biological processes at the systems' level, using an array of techniques, ranging from statistical methods to network reconstruction and analysis, to mathematical modeling. Its goal is to reconstruct the complex and often counterintuitive dynamic behavior of biological systems and quantitatively predict their responses to environmental perturbations. In this paper, we review the impact of systems biology on precision oncology. We show examples of how the analysis of signal transduction networks allows to dissect resistance to targeted therapies and inform the choice of combinations of targeted drugs based on tumor molecular alterations. Patient-specific biomarkers based on dynamical models of signaling networks can have a greater prognostic value than conventional biomarkers. These examples support systems biology models as valuable tools to advance clinical and translational oncological research.
Collapse
Affiliation(s)
- Andrea Rocca
- Hygiene and Public Health, Local Health Unit of Romagna, 47121 Forlì, Italy
| | - Boris N. Kholodenko
- Systems Biology Ireland, School of Medicine, University College Dublin, Belfield, D04 V1W8 Dublin, Ireland
- Conway Institute of Biomolecular and Biomedical Research, University College Dublin, Belfield, D04 V1W8 Dublin, Ireland
- Department of Pharmacology, Yale University School of Medicine, New Haven, CT 06520, USA
| |
Collapse
|
24
|
Choi SR, Hwang CY, Lee J, Cho KH. Network Analysis Identifies Regulators of Basal-Like Breast Cancer Reprogramming and Endocrine Therapy Vulnerability. Cancer Res 2021; 82:320-333. [PMID: 34845001 DOI: 10.1158/0008-5472.can-21-0621] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2021] [Revised: 07/13/2021] [Accepted: 11/15/2021] [Indexed: 11/16/2022]
Abstract
Basal-like breast cancer is the most aggressive breast cancer subtype with the worst prognosis. Despite its high recurrence rate, chemotherapy is the only treatment for basal-like breast cancer, which lacks expression of hormone receptors. In contrast, luminal A tumors express ERα and can undergo endocrine therapy for treatment. Previous studies have tried to develop effective treatments for basal-like patients using various therapeutics but failed due to the complex and dynamic nature of the disease. In this study, we performed a transcriptomic analysis of patients with breast cancer to construct a simplified but essential molecular regulatory network model. Network control analysis identified potential targets and elucidated the underlying mechanisms of reprogramming basal-like cancer cells into luminal A cells. Inhibition of BCL11A and HDAC1/2 effectively drove basal-like cells to transition to luminal A cells and increased ERα expression, leading to increased tamoxifen sensitivity. High expression of BCL11A and HDAC1/2 correlated with poor prognosis in patients with breast cancer. These findings identify mechanisms regulating breast cancer phenotypes and suggest the potential to reprogram basal-like breast cancer cells to enhance their targetability. SIGNIFICANCE: A network model enables investigation of mechanisms regulating the basal-to-luminal transition in breast cancer, identifying BCL11A and HDAC1/2 as optimal targets that can induce basal-like breast cancer reprogramming and endocrine therapy sensitivity.
Collapse
Affiliation(s)
- Sea R Choi
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, Republic of Korea
| | - Chae Young Hwang
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, Republic of Korea
| | - Jonghoon Lee
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, Republic of Korea
| | - Kwang-Hyun Cho
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, Republic of Korea.
| |
Collapse
|
25
|
Computational Phosphorylation Network Reconstruction: An Update on Methods and Resources. Methods Mol Biol 2021. [PMID: 34270057 DOI: 10.1007/978-1-0716-1625-3_15] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2023]
Abstract
Most proteins undergo some form of modification after translation, and phosphorylation is one of the most relevant and ubiquitous post-translational modifications. The succession of protein phosphorylation and dephosphorylation catalyzed by protein kinase and phosphatase, respectively, constitutes a key mechanism of molecular information flow in cellular systems. The protein interactions of kinases, phosphatases, and their regulatory subunits and substrates are the main part of phosphorylation networks. To elucidate the landscape of phosphorylation events has been a central goal pursued by both experimental and computational approaches. Substrate specificity (e.g., sequence, structure) or the phosphoproteome has been utilized in an array of different statistical learning methods to infer phosphorylation networks. In this chapter, different computational phosphorylation network inference-related methods and resources are summarized and discussed.
Collapse
|
26
|
Bernardo-Faura M, Rinas M, Wirbel J, Pertsovskaya I, Pliaka V, Messinis DE, Vila G, Sakellaropoulos T, Faigle W, Stridh P, Behrens JR, Olsson T, Martin R, Paul F, Alexopoulos LG, Villoslada P, Saez-Rodriguez J. Prediction of combination therapies based on topological modeling of the immune signaling network in multiple sclerosis. Genome Med 2021; 13:117. [PMID: 34271980 PMCID: PMC8284018 DOI: 10.1186/s13073-021-00925-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2019] [Accepted: 06/14/2021] [Indexed: 11/21/2022] Open
Abstract
Background Multiple sclerosis (MS) is a major health problem, leading to a significant disability and patient suffering. Although chronic activation of the immune system is a hallmark of the disease, its pathogenesis is poorly understood, while current treatments only ameliorate the disease and may produce severe side effects. Methods Here, we applied a network-based modeling approach based on phosphoproteomic data to uncover the differential activation in signaling wiring between healthy donors, untreated patients, and those under different treatments. Based in the patient-specific networks, we aimed to create a new approach to identify drug combinations that revert signaling to a healthy-like state. We performed ex vivo multiplexed phosphoproteomic assays upon perturbations with multiple drugs and ligands in primary immune cells from 169 subjects (MS patients, n=129 and matched healthy controls, n=40). Patients were either untreated or treated with fingolimod, natalizumab, interferon-β, glatiramer acetate, or the experimental therapy epigallocatechin gallate (EGCG). We generated for each donor a dynamic logic model by fitting a bespoke literature-derived network of MS-related pathways to the perturbation data. Last, we developed an approach based on network topology to identify deregulated interactions whose activity could be reverted to a “healthy-like” status by combination therapy. The experimental autoimmune encephalomyelitis (EAE) mouse model of MS was used to validate the prediction of combination therapies. Results Analysis of the models uncovered features of healthy-, disease-, and drug-specific signaling networks. We predicted several combinations with approved MS drugs that could revert signaling to a healthy-like state. Specifically, TGF-β activated kinase 1 (TAK1) kinase, involved in Transforming growth factor β-1 proprotein (TGF-β), Toll-like receptor, B cell receptor, and response to inflammation pathways, was found to be highly deregulated and co-druggable with all MS drugs studied. One of these predicted combinations, fingolimod with a TAK1 inhibitor, was validated in an animal model of MS. Conclusions Our approach based on donor-specific signaling networks enables prediction of targets for combination therapy for MS and other complex diseases. Supplementary Information The online version contains supplementary material available at 10.1186/s13073-021-00925-8.
Collapse
Affiliation(s)
- Marti Bernardo-Faura
- European Molecular Biology Laboratory, European Bioinformatics Institute, Cambridge, UK.,Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus UAB, Bellaterra, Barcelona, Spain
| | - Melanie Rinas
- Joint Research Center for Computational Biomedicine (JRC-COMBINE), Faculty of Medicine, RWTH-Aachen University, Aachen, Germany
| | - Jakob Wirbel
- European Molecular Biology Laboratory, European Bioinformatics Institute, Cambridge, UK.,Joint Research Center for Computational Biomedicine (JRC-COMBINE), Faculty of Medicine, RWTH-Aachen University, Aachen, Germany
| | - Inna Pertsovskaya
- Institut d' Investigacions Biomèdiques August Pi Sunyer (IDIBAPS), Barcelona, Spain
| | - Vicky Pliaka
- School of Mechanical Engineering, National Technical University of Athens, Zografou, Greece
| | | | - Gemma Vila
- Institut d' Investigacions Biomèdiques August Pi Sunyer (IDIBAPS), Barcelona, Spain
| | | | | | - Pernilla Stridh
- Department of Clinical Neuroscience, Karolinska Institutet, Stockholm, Sweden
| | - Janina R Behrens
- NeuroCure Clinical Research Center and Department of Neurology, Charité University Medicine Berlin, Berlin, Germany
| | - Tomas Olsson
- Department of Clinical Neuroscience, Karolinska Institutet, Stockholm, Sweden
| | | | - Friedemann Paul
- NeuroCure Clinical Research Center and Department of Neurology, Charité University Medicine Berlin, Berlin, Germany
| | - Leonidas G Alexopoulos
- School of Mechanical Engineering, National Technical University of Athens, Zografou, Greece. .,ProtATonce Ltd., Athens, Greece.
| | - Pablo Villoslada
- Institut d' Investigacions Biomèdiques August Pi Sunyer (IDIBAPS), Barcelona, Spain.
| | - Julio Saez-Rodriguez
- European Molecular Biology Laboratory, European Bioinformatics Institute, Cambridge, UK. .,Joint Research Center for Computational Biomedicine (JRC-COMBINE), Faculty of Medicine, RWTH-Aachen University, Aachen, Germany. .,Institute for Computational Biomedicine, Heidelberg University Hospital and Faculty of Medicine, Heidelberg University, Bioquant, Heidelberg, Germany.
| |
Collapse
|
27
|
Zhang L, Liu G, Kong M, Li T, Wu D, Zhou X, Yang C, Xia L, Yang Z, Chen L. Revealing dynamic regulations and the related key proteins of myeloma-initiating cells by integrating experimental data into a systems biological model. Bioinformatics 2021; 37:1554-1561. [PMID: 31350562 DOI: 10.1093/bioinformatics/btz542] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2018] [Revised: 06/17/2019] [Accepted: 07/19/2019] [Indexed: 12/24/2022] Open
Abstract
MOTIVATION The growth and survival of myeloma cells are greatly affected by their surrounding microenvironment. To understand the molecular mechanism and the impact of stiffness on the fate of myeloma-initiating cells (MICs), we develop a systems biological model to reveal the dynamic regulations by integrating reverse-phase protein array data and the stiffness-associated pathway. RESULTS We not only develop a stiffness-associated signaling pathway to describe the dynamic regulations of the MICs, but also clearly identify three critical proteins governing the MIC proliferation and death, including FAK, mTORC1 and NFκB, which are validated to be related with multiple myeloma by our immunohistochemistry experiment, computation and manually reviewed evidences. Moreover, we demonstrate that the systematic model performs better than widely used parameter estimation algorithms for the complicated signaling pathway. AVAILABILITY AND IMPLEMENTATION We can not only use the systems biological model to infer the stiffness-associated genetic signaling pathway and locate the critical proteins, but also investigate the important pathways, proteins or genes for other type of the cancer. Thus, it holds universal scientific significance. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Le Zhang
- College of Computer Science.,Medical Big Data Center, Sichuan University, Chengdu 610065, China.,Chongqqing Zhongdi Medical Information Technology Co., Ltd, Chongqing 401320, China
| | - Guangdi Liu
- College of Computer and Information Science, Southwest University, Chongqing 400715, China.,Library of Chengdu University, Chengdu University, Chengdu 610106, China
| | - Meijing Kong
- College of Computer and Information Science, Southwest University, Chongqing 400715, China
| | - Tingting Li
- College of Mathematics and Statistics, Southwest University, Chongqing 400715, China
| | - Dan Wu
- Department of Radiology, Wake Forest University School of Medicine, Winston-Salem, NC 27157, USA
| | - Xiaobo Zhou
- Department of Radiology, Wake Forest University School of Medicine, Winston-Salem, NC 27157, USA
| | - Chuanwei Yang
- Systems Biology, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
| | - Lei Xia
- Cancer Center, Research Institute of Surgery, Daping Hospital, Third Military Medical University, Chongqing 400042, China
| | - Zhenzhou Yang
- Cancer Center, Research Institute of Surgery, Daping Hospital, Third Military Medical University, Chongqing 400042, China
| | - Luonan Chen
- Key Laboratory of Systems Biology, CAS Center for Excellence in Molecular Cell Science, Institute of Biochemistry and Cell Biology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai 200031, China.,Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming 650223, China.,Shanghai Research Center for Brain Science and Brain-Inspired Intelligence, Shanghai 201210, China
| |
Collapse
|
28
|
Pugliese GM, Latini S, Massacci G, Perfetto L, Sacco F. Combining Mass Spectrometry-Based Phosphoproteomics with a Network-Based Approach to Reveal FLT3-Dependent Mechanisms of Chemoresistance. Proteomes 2021; 9:19. [PMID: 33925552 PMCID: PMC8167576 DOI: 10.3390/proteomes9020019] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Revised: 04/22/2021] [Accepted: 04/24/2021] [Indexed: 12/19/2022] Open
Abstract
FLT3 mutations are the most frequently identified genetic alterations in acute myeloid leukemia (AML) and are associated with poor clinical outcome, relapse and chemotherapeutic resistance. Elucidating the molecular mechanisms underlying FLT3-dependent pathogenesis and drug resistance is a crucial goal of biomedical research. Given the complexity and intricacy of protein signaling networks, deciphering the molecular basis of FLT3-driven drug resistance requires a systems approach. Here we discuss how the recent advances in mass spectrometry (MS)-based (phospho) proteomics and multiparametric analysis accompanied by emerging computational approaches offer a platform to obtain and systematically analyze cell-specific signaling networks and to identify new potential therapeutic targets.
Collapse
Affiliation(s)
- Giusj Monia Pugliese
- Department of Biology, University of Rome Tor Vergata, Via delle Ricerca Scientifica 1, 00133 Rome, Italy; (G.M.P.); (S.L.); (G.M.)
| | - Sara Latini
- Department of Biology, University of Rome Tor Vergata, Via delle Ricerca Scientifica 1, 00133 Rome, Italy; (G.M.P.); (S.L.); (G.M.)
| | - Giorgia Massacci
- Department of Biology, University of Rome Tor Vergata, Via delle Ricerca Scientifica 1, 00133 Rome, Italy; (G.M.P.); (S.L.); (G.M.)
| | - Livia Perfetto
- Fondazione Human Technopole, Department of Biology, Via Cristina Belgioioso 171, 20157 Milan, Italy;
| | - Francesca Sacco
- Department of Biology, University of Rome Tor Vergata, Via delle Ricerca Scientifica 1, 00133 Rome, Italy; (G.M.P.); (S.L.); (G.M.)
| |
Collapse
|
29
|
Gene regulatory networks exhibit several kinds of memory: quantification of memory in biological and random transcriptional networks. iScience 2021; 24:102131. [PMID: 33748699 PMCID: PMC7970124 DOI: 10.1016/j.isci.2021.102131] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2020] [Revised: 12/09/2020] [Accepted: 01/26/2021] [Indexed: 02/08/2023] Open
Abstract
Gene regulatory networks (GRNs) process important information in developmental biology and biomedicine. A key knowledge gap concerns how their responses change over time. Hypothesizing long-term changes of dynamics induced by transient prior events, we created a computational framework for defining and identifying diverse types of memory in candidate GRNs. We show that GRNs from a wide range of model systems are predicted to possess several types of memory, including Pavlovian conditioning. Associative memory offers an alternative strategy for the biomedical use of powerful drugs with undesirable side effects, and a novel approach to understanding the variability and time-dependent changes of drug action. We find evidence of natural selection favoring GRN memory. Vertebrate GRNs overall exhibit more memory than invertebrate GRNs, and memory is most prevalent in differentiated metazoan cell networks compared with undifferentiated cells. Timed stimuli are a powerful alternative for biomedical control of complex in vivo dynamics without genomic editing or transgenes. Gene regulatory networks' dynamics are modified by transient stimuli GRNs have several different types of memory, including associative conditioning Evolution favored GRN memory, and differentiated cells have the most memory capacity Training GRNs offers a novel biomedical strategy not dependent on genetic rewiring
Collapse
|
30
|
Integrating Patient-Specific Information into Logic Models of Complex Diseases: Application to Acute Myeloid Leukemia. J Pers Med 2021; 11:jpm11020117. [PMID: 33578936 PMCID: PMC7916657 DOI: 10.3390/jpm11020117] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2020] [Revised: 02/05/2021] [Accepted: 02/05/2021] [Indexed: 12/12/2022] Open
Abstract
High throughput technologies such as deep sequencing and proteomics are increasingly becoming mainstream in clinical practice and support diagnosis and patient stratification. Developing computational models that recapitulate cell physiology and its perturbations in disease is a required step to help with the interpretation of results of high content experiments and to devise personalized treatments. As complete cell-models are difficult to achieve, given limited experimental information and insurmountable computational problems, approximate approaches should be considered. We present here a general approach to modeling complex diseases by embedding patient-specific genomics data into actionable logic models that take into account prior knowledge. We apply the strategy to acute myeloid leukemia (AML) and assemble a network of logical relationships linking most of the genes that are found frequently mutated in AML patients. We derive Boolean models from this network and we show that by priming the model with genomic data we can infer relevant patient-specific clinical features. Here we propose that the integration of literature-derived causal networks with patient-specific data should be explored to help bedside decisions.
Collapse
|
31
|
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: 25] [Impact Index Per Article: 8.3] [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
|
32
|
Howell RSM, Klemm C, Thorpe PH, Csikász-Nagy A. Unifying the mechanism of mitotic exit control in a spatiotemporal logical model. PLoS Biol 2020; 18:e3000917. [PMID: 33180788 PMCID: PMC7685450 DOI: 10.1371/journal.pbio.3000917] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2020] [Revised: 11/24/2020] [Accepted: 10/09/2020] [Indexed: 11/18/2022] Open
Abstract
The transition from mitosis into the first gap phase of the cell cycle in budding yeast is controlled by the Mitotic Exit Network (MEN). The network interprets spatiotemporal cues about the progression of mitosis and ensures that release of Cdc14 phosphatase occurs only after completion of key mitotic events. The MEN has been studied intensively; however, a unified understanding of how localisation and protein activity function together as a system is lacking. In this paper, we present a compartmental, logical model of the MEN that is capable of representing spatial aspects of regulation in parallel to control of enzymatic activity. We show that our model is capable of correctly predicting the phenotype of the majority of mutants we tested, including mutants that cause proteins to mislocalise. We use a continuous time implementation of the model to demonstrate that Cdc14 Early Anaphase Release (FEAR) ensures robust timing of anaphase, and we verify our findings in living cells. Furthermore, we show that our model can represent measured cell-cell variation in Spindle Position Checkpoint (SPoC) mutants. This work suggests a general approach to incorporate spatial effects into logical models. We anticipate that the model itself will be an important resource to experimental researchers, providing a rigorous platform to test hypotheses about regulation of mitotic exit.
Collapse
Affiliation(s)
- Rowan S M Howell
- The Francis Crick Institute, London, United Kingdom.,Randall Centre for Cell and Molecular Biophysics, King's College London, London, United Kingdom
| | - Cinzia Klemm
- School of Biological and Chemical Sciences, Queen Mary University, London, United Kingdom
| | - Peter H Thorpe
- School of Biological and Chemical Sciences, Queen Mary University, London, United Kingdom
| | - Attila Csikász-Nagy
- Randall Centre for Cell and Molecular Biophysics, King's College London, London, United Kingdom.,Faculty of Information Technology and Bionics, Pázmány Péter Catholic University, Budapest, Hungary
| |
Collapse
|
33
|
Malekpour SA, Alizad-Rahvar AR, Sadeghi M. LogicNet: probabilistic continuous logics in reconstructing gene regulatory networks. BMC Bioinformatics 2020; 21:318. [PMID: 32690031 PMCID: PMC7372900 DOI: 10.1186/s12859-020-03651-x] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2019] [Accepted: 07/10/2020] [Indexed: 11/10/2022] Open
Abstract
Background Gene Regulatory Networks (GRNs) have been previously studied by using Boolean/multi-state logics. While the gene expression values are usually scaled into the range [0, 1], these GRN inference methods apply a threshold to discretize the data, resulting in missing information. Most of studies apply fuzzy logics to infer the logical gene-gene interactions from continuous data. However, all these approaches require an a priori known network structure. Results Here, by introducing a new probabilistic logic for continuous data, we propose a novel logic-based approach (called the LogicNet) for the simultaneous reconstruction of the GRN structure and identification of the logics among the regulatory genes, from the continuous gene expression data. In contrast to the previous approaches, the LogicNet does not require an a priori known network structure to infer the logics. The proposed probabilistic logic is superior to the existing fuzzy logics and is more relevant to the biological contexts than the fuzzy logics. The performance of the LogicNet is superior to that of several Mutual Information-based and regression-based tools for reconstructing GRNs. Conclusions The LogicNet reconstructs GRNs and logic functions without requiring prior knowledge of the network structure. Moreover, in another application, the LogicNet can be applied for logic function detection from the known regulatory genes-target interactions. We also conclude that computational modeling of the logical interactions among the regulatory genes significantly improves the GRN reconstruction accuracy.
Collapse
Affiliation(s)
- Seyed Amir Malekpour
- School of Biological Sciences, Institute for Research in Fundamental Sciences (IPM), Tehran, Iran.
| | - Amir Reza Alizad-Rahvar
- School of Biological Sciences, Institute for Research in Fundamental Sciences (IPM), Tehran, Iran
| | - Mehdi Sadeghi
- National Institute of Genetic Engineering and Biotechnology, Tehran, Iran
| |
Collapse
|
34
|
Banzhaf M, Resendis-Antonio O, Zepeda-Mendoza ML. Uncovering the Dynamic Mechanisms of the Pseudomonas Aeruginosa Quorum Sensing and Virulence Networks Using Boolean Modelling. IEEE Trans Nanobioscience 2020; 19:394-402. [DOI: 10.1109/tnb.2020.2977820] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
35
|
Abstract
Making decisions on how best to treat cancer patients requires the integration of different data sets, including genomic profiles, tumour histopathology, radiological images, proteomic analysis and more. This wealth of biological information calls for novel strategies to integrate such information in a meaningful, predictive and experimentally verifiable way. In this Perspective we explain how executable computational models meet this need. Such models provide a means for comprehensive data integration, can be experimentally validated, are readily interpreted both biologically and clinically, and have the potential to predict effective therapies for different cancer types and subtypes. We explain what executable models are and how they can be used to represent the dynamic biological behaviours inherent in cancer, and demonstrate how such models, when coupled with automated reasoning, facilitate our understanding of the mechanisms by which oncogenic signalling pathways regulate tumours. We explore how executable models have impacted the field of cancer research and argue that extending them to represent a tumour in a specific patient (that is, an avatar) will pave the way for improved personalized treatments and precision medicine. Finally, we highlight some of the ongoing challenges in developing executable models and stress that effective cross-disciplinary efforts are key to forward progress in the field.
Collapse
Affiliation(s)
- Matthew A Clarke
- Department of Biochemistry, University of Cambridge, Cambridge, UK
| | - Jasmin Fisher
- Department of Biochemistry, University of Cambridge, Cambridge, UK.
- UCL Cancer Institute, University College London, London, UK.
| |
Collapse
|
36
|
Hastings JF, O'Donnell YEI, Fey D, Croucher DR. Applications of personalised signalling network models in precision oncology. Pharmacol Ther 2020; 212:107555. [PMID: 32320730 DOI: 10.1016/j.pharmthera.2020.107555] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Accepted: 04/07/2020] [Indexed: 02/07/2023]
Abstract
As our ability to provide in-depth, patient-specific characterisation of the molecular alterations within tumours rapidly improves, it is becoming apparent that new approaches will be required to leverage the power of this data and derive the full benefit for each individual patient. Systems biology approaches are beginning to emerge within this field as a potential method of incorporating large volumes of network level data and distilling a coherent, clinically-relevant prediction of drug response. However, the initial promise of this developing field is yet to be realised. Here we argue that in order to develop these precise models of individual drug response and tailor treatment accordingly, we will need to develop mathematical models capable of capturing both the dynamic nature of drug-response signalling networks and key patient-specific information such as mutation status or expression profiles. We also review the modelling approaches commonly utilised within this field, and outline recent examples of their use in furthering the application of systems biology for a precision medicine approach to cancer treatment.
Collapse
Affiliation(s)
- Jordan F Hastings
- The Kinghorn Cancer Centre, Garvan Institute of Medical Research, Sydney, Australia
| | | | - Dirk Fey
- Systems Biology Ireland, University College Dublin, Belfield, Dublin 4, Ireland; School of Medicine and Medical Science, University College Dublin, Belfield, Dublin 4, Ireland
| | - David R Croucher
- The Kinghorn Cancer Centre, Garvan Institute of Medical Research, Sydney, Australia; School of Medicine and Medical Science, University College Dublin, Belfield, Dublin 4, Ireland; St Vincent's Hospital Clinical School, University of New South Wales, Sydney, NSW 2052, Australia.
| |
Collapse
|
37
|
Feedback analysis identifies a combination target for overcoming adaptive resistance to targeted cancer therapy. Oncogene 2020; 39:3803-3820. [PMID: 32157217 DOI: 10.1038/s41388-020-1255-y] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2019] [Revised: 02/26/2020] [Accepted: 02/28/2020] [Indexed: 12/19/2022]
Abstract
Targeted drugs aim to treat cancer by directly inhibiting oncogene activity or oncogenic pathways, but drug resistance frequently emerges. Due to the intricate dynamics of cancer signaling networks, which contain complex feedback regulations, cancer cells can rewire these networks to adapt to and counter the cytotoxic effects of a drug, thereby limiting the efficacy of targeted therapies. To identify a combinatorial drug target that can overcome such a limitation, we developed a Boolean network simulation and analysis framework and applied this approach to a large-scale signaling network of colorectal cancer with integrated genomic information. We discovered Src as a critical combination drug target that can overcome the adaptive resistance to the targeted inhibition of mitogen-activated protein kinase pathway by blocking the essential feedback regulation responsible for resistance. The proposed framework is generic and can be widely used to identify drug targets that can overcome adaptive resistance to targeted therapies.
Collapse
|
38
|
Mukherjee S, Nazemi M, Jonkers I, Geris L. Use of Computational Modeling to Study Joint Degeneration: A Review. Front Bioeng Biotechnol 2020; 8:93. [PMID: 32185167 PMCID: PMC7058554 DOI: 10.3389/fbioe.2020.00093] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2019] [Accepted: 01/31/2020] [Indexed: 12/13/2022] Open
Abstract
Osteoarthritis (OA), a degenerative joint disease, is the most common chronic condition of the joints, which cannot be prevented effectively. Computational modeling of joint degradation allows to estimate the patient-specific progression of OA, which can aid clinicians to estimate the most suitable time window for surgical intervention in osteoarthritic patients. This paper gives an overview of the different approaches used to model different aspects of joint degeneration, thereby focusing mostly on the knee joint. The paper starts by discussing how OA affects the different components of the joint and how these are accounted for in the models. Subsequently, it discusses the different modeling approaches that can be used to answer questions related to OA etiology, progression and treatment. These models are ordered based on their underlying assumptions and technologies: musculoskeletal models, Finite Element models, (gene) regulatory models, multiscale models and data-driven models (artificial intelligence/machine learning). Finally, it is concluded that in the future, efforts should be made to integrate the different modeling techniques into a more robust computational framework that should not only be efficient to predict OA progression but also easily allow a patient’s individualized risk assessment as screening tool for use in clinical practice.
Collapse
Affiliation(s)
- Satanik Mukherjee
- Prometheus, Division of Skeletal Tissue Engineering, KU Leuven, Leuven, Belgium.,Biomechanics Section, KU Leuven, Leuven, Belgium
| | - Majid Nazemi
- GIGA in silico Medicine, University of Liège, Liège, Belgium
| | - Ilse Jonkers
- Human Movement Biomechanics Research Group, Department of Movement Sciences, KU Leuven, Leuven, Belgium
| | - Liesbet Geris
- Prometheus, Division of Skeletal Tissue Engineering, KU Leuven, Leuven, Belgium.,Biomechanics Section, KU Leuven, Leuven, Belgium.,GIGA in silico Medicine, University of Liège, Liège, Belgium
| |
Collapse
|
39
|
Eduati F, Jaaks P, Wappler J, Cramer T, Merten CA, Garnett MJ, Saez‐Rodriguez J. Patient-specific logic models of signaling pathways from screenings on cancer biopsies to prioritize personalized combination therapies. Mol Syst Biol 2020; 16:e8664. [PMID: 32073727 PMCID: PMC7029724 DOI: 10.15252/msb.20188664] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2018] [Revised: 01/27/2020] [Accepted: 01/28/2020] [Indexed: 12/22/2022] Open
Abstract
Mechanistic modeling of signaling pathways mediating patient-specific response to therapy can help to unveil resistance mechanisms and improve therapeutic strategies. Yet, creating such models for patients, in particular for solid malignancies, is challenging. A major hurdle to build these models is the limited material available that precludes the generation of large-scale perturbation data. Here, we present an approach that couples ex vivo high-throughput screenings of cancer biopsies using microfluidics with logic-based modeling to generate patient-specific dynamic models of extrinsic and intrinsic apoptosis signaling pathways. We used the resulting models to investigate heterogeneity in pancreatic cancer patients, showing dissimilarities especially in the PI3K-Akt pathway. Variation in model parameters reflected well the different tumor stages. Finally, we used our dynamic models to efficaciously predict new personalized combinatorial treatments. Our results suggest that our combination of microfluidic experiments and mathematical model can be a novel tool toward cancer precision medicine.
Collapse
Affiliation(s)
- Federica Eduati
- European Molecular Biology Laboratory (EMBL)Genome Biology UnitHeidelbergGermany
- European Molecular Biology LaboratoryEuropean Bioinformatics Institute (EMBL‐EBI)HinxtonUK
- Joint Research Centre for Computational Biomedicine (JRC‐COMBINE)Faculty of MedicineRWTH Aachen UniversityAachenGermany
- Department of Biomedical EngineeringEindhoven University of TechnologyEindhovenThe Netherlands
| | | | - Jessica Wappler
- Department SurgeryMolecular Tumor BiologyRWTH University HospitalAachenGermany
| | - Thorsten Cramer
- Department SurgeryMolecular Tumor BiologyRWTH University HospitalAachenGermany
- ESCAM – European Surgery Center Aachen MaastrichtAachenGermany
- ESCAM – European Surgery Center Aachen MaastrichtMaastrichtThe Netherlands
| | - Christoph A Merten
- European Molecular Biology Laboratory (EMBL)Genome Biology UnitHeidelbergGermany
| | | | - Julio Saez‐Rodriguez
- European Molecular Biology LaboratoryEuropean Bioinformatics Institute (EMBL‐EBI)HinxtonUK
- Joint Research Centre for Computational Biomedicine (JRC‐COMBINE)Faculty of MedicineRWTH Aachen UniversityAachenGermany
- Institute for Computational BiomedicineFaculty of MedicineBIOQUANT‐CenterHeidelberg UniversityHeidelbergGermany
| |
Collapse
|
40
|
Bruner A, Sharan R. A Robustness Analysis of Dynamic Boolean Models of Cellular Circuits. J Comput Biol 2020; 27:133-143. [PMID: 31770006 DOI: 10.1089/cmb.2019.0290] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
With ever growing amounts of omics data, the next challenge in biological research is the interpretation of these data to gain mechanistic insights about cellular function. Dynamic models of cellular circuits that capture the activity levels of proteins and other molecules over time offer great expressive power by allowing the simulation of the effects of specific internal or external perturbations on the workings of the cell. However, the study of such models is at its infancy and no large-scale analysis of the robustness of real models to changing conditions has been conducted to date. Here we provide a computational framework to study the robustness of such models using a combination of stochastic simulations and integer linear programming techniques. We apply our framework to a large collection of cellular circuits and benchmark the results against randomized models. We find that the steady states of real circuits tend to be more robust in multiple aspects compared with their randomized counterparts.
Collapse
Affiliation(s)
- Ariel Bruner
- Blavatnik School of Computer Science, Tel Aviv University, Tel Aviv, Israel
| | - Roded Sharan
- Blavatnik School of Computer Science, Tel Aviv University, Tel Aviv, Israel
| |
Collapse
|
41
|
Romers J, Thieme S, Münzner U, Krantz M. A scalable method for parameter-free simulation and validation of mechanistic cellular signal transduction network models. NPJ Syst Biol Appl 2020; 6:2. [PMID: 31934349 PMCID: PMC6954118 DOI: 10.1038/s41540-019-0120-5] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2018] [Accepted: 11/20/2019] [Indexed: 11/09/2022] Open
Abstract
The metabolic modelling community has established the gold standard for bottom-up systems biology with reconstruction, validation and simulation of mechanistic genome-scale models. Similar methods have not been established for signal transduction networks, where the representation of complexes and internal states leads to scalability issues in both model formulation and execution. While rule- and agent-based methods allow efficient model definition and execution, respectively, model parametrisation introduces an additional layer of uncertainty due to the sparsity of reliably measured parameters. Here, we present a scalable method for parameter-free simulation of mechanistic signal transduction networks. It is based on rxncon and uses a bipartite Boolean logic with separate update rules for reactions and states. Using two generic update rules, we enable translation of any rxncon model into a unique Boolean model, which can be used for network validation and simulation-allowing the prediction of system-level function directly from molecular mechanistic data. Through scalable model definition and simulation, and the independence of quantitative parameters, it opens up for simulation and validation of mechanistic genome-scale models of signal transduction networks.
Collapse
Affiliation(s)
- Jesper Romers
- Institute for Biology, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Sebastian Thieme
- Institute for Biology, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Ulrike Münzner
- Institute for Biology, Humboldt-Universität zu Berlin, Berlin, Germany
- Bioinformatics Center, Institute for Chemical Research, Kyoto University, Uji, Japan
| | - Marcus Krantz
- Institute for Biology, Humboldt-Universität zu Berlin, Berlin, Germany
| |
Collapse
|
42
|
Bosdriesz E, Prahallad A, Klinger B, Sieber A, Bosma A, Bernards R, Blüthgen N, Wessels LFA. Comparative Network Reconstruction using mixed integer programming. Bioinformatics 2019; 34:i997-i1004. [PMID: 30423075 PMCID: PMC6129277 DOI: 10.1093/bioinformatics/bty616] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Motivation Signal-transduction networks are often aberrated in cancer cells, and new anti-cancer drugs that specifically target oncogenes involved in signaling show great clinical promise. However, the effectiveness of such targeted treatments is often hampered by innate or acquired resistance due to feedbacks, crosstalks or network adaptations in response to drug treatment. A quantitative understanding of these signaling networks and how they differ between cells with different oncogenic mutations or between sensitive and resistant cells can help in addressing this problem. Results Here, we present Comparative Network Reconstruction (CNR), a computational method to reconstruct signaling networks based on possibly incomplete perturbation data, and to identify which edges differ quantitatively between two or more signaling networks. Prior knowledge about network topology is not required but can straightforwardly be incorporated. We extensively tested our approach using simulated data and applied it to perturbation data from a BRAF mutant, PTPN11 KO cell line that developed resistance to BRAF inhibition. Comparing the reconstructed networks of sensitive and resistant cells suggests that the resistance mechanism involves re-establishing wild-type MAPK signaling, possibly through an alternative RAF-isoform. Availability and implementation CNR is available as a python module at https://github.com/NKI-CCB/cnr. Additionally, code to reproduce all figures is available at https://github.com/NKI-CCB/CNR-analyses. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Evert Bosdriesz
- Division of Molecular Carcinogenesis, The Oncode Institute, The Netherlands Cancer Institute, Amsterdam, The Netherlands
| | - Anirudh Prahallad
- Division of Molecular Carcinogenesis, The Oncode Institute, The Netherlands Cancer Institute, Amsterdam, The Netherlands
| | - Bertram Klinger
- Institute of Pathology, Charité Universitätsmedizin, Berlin, Germany.,IRI Life Sciences, Humboldt University of Berlin, Berlin, Germany
| | - Anja Sieber
- Institute of Pathology, Charité Universitätsmedizin, Berlin, Germany.,IRI Life Sciences, Humboldt University of Berlin, Berlin, Germany
| | - Astrid Bosma
- Division of Molecular Carcinogenesis, The Oncode Institute, The Netherlands Cancer Institute, Amsterdam, The Netherlands
| | - René Bernards
- Division of Molecular Carcinogenesis, The Oncode Institute, The Netherlands Cancer Institute, Amsterdam, The Netherlands
| | - Nils Blüthgen
- Institute of Pathology, Charité Universitätsmedizin, Berlin, Germany.,IRI Life Sciences, Humboldt University of Berlin, Berlin, Germany.,Berlin Institute of Health, Berlin, Germany
| | - Lodewyk F A Wessels
- Division of Molecular Carcinogenesis, The Oncode Institute, The Netherlands Cancer Institute, Amsterdam, The Netherlands.,Faculty of EEMCS, Delft University of Technology, Delft, The Netherlands
| |
Collapse
|
43
|
Li Y, Jann T, Vera-Licona P. Benchmarking time-series data discretization on inference methods. Bioinformatics 2019; 35:3102-3109. [PMID: 30657860 DOI: 10.1093/bioinformatics/btz036] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2018] [Revised: 12/10/2018] [Accepted: 01/14/2019] [Indexed: 12/15/2022] Open
Abstract
SUMMARY The rapid development in quantitatively measuring DNA, RNA and protein has generated a great interest in the development of reverse-engineering methods, that is, data-driven approaches to infer the network structure or dynamical model of the system. Many reverse-engineering methods require discrete quantitative data as input, while many experimental data are continuous. Some studies have started to reveal the impact that the choice of data discretization has on the performance of reverse-engineering methods. However, more comprehensive studies are still greatly needed to systematically and quantitatively understand the impact that discretization methods have on inference methods. Furthermore, there is an urgent need for systematic comparative methods that can help select between discretization methods. In this work, we consider four published intracellular networks inferred with their respective time-series datasets. We discretized the data using different discretization methods. Across all datasets, changing the data discretization to a more appropriate one improved the reverse-engineering methods' performance. We observed no universal best discretization method across different time-series datasets. Thus, we propose DiscreeTest, a two-step evaluation metric for ranking discretization methods for time-series data. The underlying assumption of DiscreeTest is that an optimal discretization method should preserve the dynamic patterns observed in the original data across all variables. We used the same datasets and networks to show that DiscreeTest is able to identify an appropriate discretization among several candidate methods. To our knowledge, this is the first time that a method for benchmarking and selecting an appropriate discretization method for time-series data has been proposed. AVAILABILITY AND IMPLEMENTATION All the datasets, reverse-engineering methods and source code used in this paper are available in Vera-Licona's lab Github repository: https://github.com/VeraLiconaResearchGroup/Benchmarking_TSDiscretizations. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Yuezhe Li
- R.D. Berlin Center for Cell Analysis and Modeling, University of Connecticut School of Medicine, Farmington, CT, USA
| | - Tiffany Jann
- Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, CA, USA
| | - Paola Vera-Licona
- Center for Quantitative Medicine, University of Connecticut School of Medicine, Farmington, CT, USA.,Department of Cell Biology, University of Connecticut School of Medicine, Farmington, CT, USA.,Department of Pediatrics, University of Connecticut School of Medicine, Farmington, CT, USA.,Institute for Systems Genomics, University of Connecticut School of Medicine, Farmington, CT, USA
| |
Collapse
|
44
|
Gabel M, Hohl T, Imle A, Fackler OT, Graw F. FAMoS: A Flexible and dynamic Algorithm for Model Selection to analyse complex systems dynamics. PLoS Comput Biol 2019; 15:e1007230. [PMID: 31419221 PMCID: PMC6697322 DOI: 10.1371/journal.pcbi.1007230] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2019] [Accepted: 06/30/2019] [Indexed: 01/12/2023] Open
Abstract
Most biological systems are difficult to analyse due to a multitude of interacting components and the concomitant lack of information about the essential dynamics. Finding appropriate models that provide a systematic description of such biological systems and that help to identify their relevant factors and processes can be challenging given the sheer number of possibilities. Model selection algorithms that evaluate the performance of a multitude of different models against experimental data provide a useful tool to identify appropriate model structures. However, many algorithms addressing the analysis of complex dynamical systems, as they are often used in biology, compare a preselected number of models or rely on exhaustive searches of the total model space which might be unfeasible dependent on the number of possibilities. Therefore, we developed an algorithm that is able to perform model selection on complex systems and searches large model spaces in a dynamical way. Our algorithm includes local and newly developed non-local search methods that can prevent the algorithm from ending up in local minima of the model space by accounting for structurally similar processes. We tested and validated the algorithm based on simulated data and showed its flexibility for handling different model structures. We also used the algorithm to analyse experimental data on the cell proliferation dynamics of CD4+ and CD8+ T cells that were cultured under different conditions. Our analyses indicated dynamical changes within the proliferation potential of cells that was reduced within tissue-like 3D ex vivo cultures compared to suspension. Due to the flexibility in handling various model structures, the algorithm is applicable to a large variety of different biological problems and represents a useful tool for the data-oriented evaluation of complex model spaces. Identifying the systematic interactions of multiple components within a complex biological system can be challenging due to the number of potential processes and the concomitant lack of information about the essential dynamics. Selection algorithms that allow an automated evaluation of a large number of different models provide a useful tool in identifying the systematic relationships between experimental data. However, many of the existing model selection algorithms are not able to address complex model structures, such as systems of differential equations, and partly rely on local or exhaustive search methods which are inappropriate for the analysis of various biological systems. Therefore, we developed a flexible model selection algorithm that performs a robust and dynamical search of large model spaces to identify complex systems dynamics and applied it to the analysis of T cell proliferation dynamics within different culture conditions. The algorithm, which is available as an R-package, provides an advanced tool for the analysis of complex systems behaviour and, due to its flexible structure, can be applied to a large variety of biological problems.
Collapse
Affiliation(s)
- Michael Gabel
- Center for Modelling and Simulation in the Biosciences, BioQuant-Center, Heidelberg University, Heidelberg, Germany
- * E-mail: (MG); (FG)
| | - Tobias Hohl
- Center for Modelling and Simulation in the Biosciences, BioQuant-Center, Heidelberg University, Heidelberg, Germany
| | - Andrea Imle
- Department of Infectious Diseases, Centre for Integrative Infectious Disease Research (CIID), Integrative Virology, University Hospital Heidelberg, Heidelberg, Germany
| | - Oliver T. Fackler
- Department of Infectious Diseases, Centre for Integrative Infectious Disease Research (CIID), Integrative Virology, University Hospital Heidelberg, Heidelberg, Germany
| | - Frederik Graw
- Center for Modelling and Simulation in the Biosciences, BioQuant-Center, Heidelberg University, Heidelberg, Germany
- * E-mail: (MG); (FG)
| |
Collapse
|
45
|
Silverbush D, Sharan R. A systematic approach to orient the human protein-protein interaction network. Nat Commun 2019; 10:3015. [PMID: 31289271 PMCID: PMC6617457 DOI: 10.1038/s41467-019-10887-6] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2018] [Accepted: 06/06/2019] [Indexed: 11/16/2022] Open
Abstract
The protein-protein interaction (PPI) network of an organism serves as a skeleton for its signaling circuitry, which mediates cellular response to environmental and genetic cues. Understanding this circuitry could improve the prediction of gene function and cellular behavior in response to diverse signals. To realize this potential, one has to comprehensively map PPIs and their directions of signal flow. While the quality and the volume of identified human PPIs improved dramatically over the last decade, the directions of these interactions are still mostly unknown, thus precluding subsequent prediction and modeling efforts. Here we present a systematic approach to orient the human PPI network using drug response and cancer genomic data. We provide a diffusion-based method for the orientation task that significantly outperforms existing methods. The oriented network leads to improved prioritization of cancer driver genes and drug targets compared to the state-of-the-art unoriented network.
Collapse
Affiliation(s)
- Dana Silverbush
- The Blavatnik School of Computer Science, Tel Aviv University, Tel Aviv, 69978, Israel
| | - Roded Sharan
- The Blavatnik School of Computer Science, Tel Aviv University, Tel Aviv, 69978, Israel.
| |
Collapse
|
46
|
Pacini C, Koziol MJ. Bioinformatics challenges and perspectives when studying the effect of epigenetic modifications on alternative splicing. Philos Trans R Soc Lond B Biol Sci 2019; 373:rstb.2017.0073. [PMID: 29685977 PMCID: PMC5915717 DOI: 10.1098/rstb.2017.0073] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/14/2017] [Indexed: 02/07/2023] Open
Abstract
It is widely known that epigenetic modifications are important in regulating transcription, but several have also been reported in alternative splicing. The regulation of pre-mRNA splicing is important to explain proteomic diversity and the misregulation of splicing has been implicated in many diseases. Here, we give a brief overview of the role of epigenetics in alternative splicing and disease. We then discuss the bioinformatics methods that can be used to model interactions between epigenetic marks and regulators of splicing. These models can be used to identify alternative splicing and epigenetic changes across different phenotypes. This article is part of a discussion meeting issue ‘Frontiers in epigenetic chemical biology’.
Collapse
Affiliation(s)
- Clare Pacini
- Wellcome Trust Cancer Research UK Gurdon Institute, University of Cambridge, Tennis Court Road, Cambridge, CB2 1QN, UK.,Department of Zoology, University of Cambridge, Downing Street, Cambridge, CB2 3EJ, UK
| | - Magdalena J Koziol
- Wellcome Trust Cancer Research UK Gurdon Institute, University of Cambridge, Tennis Court Road, Cambridge, CB2 1QN, UK .,Department of Zoology, University of Cambridge, Downing Street, Cambridge, CB2 3EJ, UK
| |
Collapse
|
47
|
Drug Effect Prediction by Integrating L1000 Genomic and Proteomic Big Data. Methods Mol Biol 2019. [PMID: 30848468 DOI: 10.1007/978-1-4939-9089-4_16] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
Abstract
The library of integrated Network-Based Cellular Signatures (LINCS) project aims to create a network-based understanding of biology by cataloging changes in gene expression and signal transduction. Gene expression and proteomic data in LINCS L1000 are cataloged for human cancer cells treated with compounds and genetic reagents. For understanding the related cell pathways and facilitating drug discovery, we developed binary linear programming (BLP) to infer cell-specific pathways and identify compounds' effects using L1000 gene expression and phosphoproteomics data. A generic pathway map for the MCF7 breast cancer cell line was built. Within them, BLP extracted the cell-specific pathways, which reliably predicted the compounds' effects. In this way, the potential drug effects are revealed by our models.
Collapse
|
48
|
Béal J, Montagud A, Traynard P, Barillot E, Calzone L. Personalization of Logical Models With Multi-Omics Data Allows Clinical Stratification of Patients. Front Physiol 2019; 9:1965. [PMID: 30733688 PMCID: PMC6353844 DOI: 10.3389/fphys.2018.01965] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2018] [Accepted: 12/31/2018] [Indexed: 12/26/2022] Open
Abstract
Logical models of cancer pathways are typically built by mining the literature for relevant experimental observations. They are usually generic as they apply for large cohorts of individuals. As a consequence, they generally do not capture the heterogeneity of patient tumors and their therapeutic responses. We present here a novel framework, referred to as PROFILE, to tailor logical models to a particular biological sample such as a patient tumor. This methodology permits to compare the model simulations to individual clinical data, i.e., survival time. Our approach focuses on integrating mutation data, copy number alterations (CNA), and expression data (transcriptomics or proteomics) to logical models. These data need first to be either binarized or set between 0 and 1, and can then be incorporated in the logical model by modifying the activity of the node, the initial conditions or the state transition rates. The use of MaBoSS, a tool based on Monte-Carlo kinetic algorithm to perform stochastic simulations on logical models results in model state probabilities, and allows for a semi-quantitative study of the model phenotypes and perturbations. As a proof of concept, we use a published generic model of cancer signaling pathways and molecular data from METABRIC breast cancer patients. For this example, we test several combinations of data incorporation and discuss that, with these data, the most comprehensive patient-specific cancer models are obtained by modifying the nodes' activity of the model with mutations, in combination or not with CNA data, and altering the transition rates with RNA expression. We conclude that these model simulations show good correlation with clinical data such as patients' Nottingham prognostic index (NPI) subgrouping and survival time. We observe that two highly relevant cancer phenotypes derived from personalized models, Proliferation and Apoptosis, are biologically consistent prognostic factors: patients with both high proliferation and low apoptosis have the worst survival rate, and conversely. Our approach aims to combine the mechanistic insights of logical modeling with multi-omics data integration to provide patient-relevant models. This work leads to the use of logical modeling for precision medicine and will eventually facilitate the choice of patient-specific drug treatments by physicians.
Collapse
Affiliation(s)
| | | | | | - Emmanuel Barillot
- Institut Curie, PSL Research University, Mines Paris Tech, Inserm, U900, Paris, France
| | - Laurence Calzone
- Institut Curie, PSL Research University, Mines Paris Tech, Inserm, U900, Paris, France
| |
Collapse
|
49
|
Lesage R, Kerkhofs J, Geris L. Computational Modeling and Reverse Engineering to Reveal Dominant Regulatory Interactions Controlling Osteochondral Differentiation: Potential for Regenerative Medicine. Front Bioeng Biotechnol 2018; 6:165. [PMID: 30483498 PMCID: PMC6243751 DOI: 10.3389/fbioe.2018.00165] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2018] [Accepted: 10/22/2018] [Indexed: 01/11/2023] Open
Abstract
The specialization of cartilage cells, or chondrogenic differentiation, is an intricate and meticulously regulated process that plays a vital role in both bone formation and cartilage regeneration. Understanding the molecular regulation of this process might help to identify key regulatory factors that can serve as potential therapeutic targets, or that might improve the development of qualitative and robust skeletal tissue engineering approaches. However, each gene involved in this process is influenced by a myriad of feedback mechanisms that keep its expression in a desirable range, making the prediction of what will happen if one of these genes defaults or is targeted with drugs, challenging. Computer modeling provides a tool to simulate this intricate interplay from a network perspective. This paper aims to give an overview of the current methodologies employed to analyze cell differentiation in the context of skeletal tissue engineering in general and osteochondral differentiation in particular. In network modeling, a network can either be derived from mechanisms and pathways that have been reported in the literature (knowledge-based approach) or it can be inferred directly from the data (data-driven approach). Combinatory approaches allow further optimization of the network. Once a network is established, several modeling technologies are available to interpret dynamically the relationships that have been put forward in the network graph (implication of the activation or inhibition of certain pathways on the evolution of the system over time) and to simulate the possible outcomes of the established network such as a given cell state. This review provides for each of the aforementioned steps (building, optimizing, and modeling the network) a brief theoretical perspective, followed by a concise overview of published works, focusing solely on applications related to cell fate decisions, cartilage differentiation and growth plate biology. Particular attention is paid to an in-house developed example of gene regulatory network modeling of growth plate chondrocyte differentiation as all the aforementioned steps can be illustrated. In summary, this paper discusses and explores a series of tools that form a first step toward a rigorous and systems-level modeling of osteochondral differentiation in the context of regenerative medicine.
Collapse
Affiliation(s)
- Raphaelle Lesage
- Prometheus, Division of Skeletal Tissue Engineering Leuven, KU Leuven, Leuven, Belgium.,Biomechanics Section, KU Leuven, Leuven, Belgium
| | - Johan Kerkhofs
- Prometheus, Division of Skeletal Tissue Engineering Leuven, KU Leuven, Leuven, Belgium.,Biomechanics Section, KU Leuven, Leuven, Belgium
| | - Liesbet Geris
- Prometheus, Division of Skeletal Tissue Engineering Leuven, KU Leuven, Leuven, Belgium.,Biomechanics Section, KU Leuven, Leuven, Belgium.,Biomechanics Research Unit, GIGA in silico Medicine, University of Liège, Liège, Belgium
| |
Collapse
|
50
|
Oftadeh MO, Marashi SA. Accounting for robustness in modeling signal transduction responses. J Recept Signal Transduct Res 2018; 38:442-447. [DOI: 10.1080/10799893.2019.1572762] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
| | - Sayed-Amir Marashi
- Department of Biotechnology, College of Science, University of Tehran, Tehran, Iran
| |
Collapse
|