1
|
Jayathilake PG, Victori P, Pavillet CE, Lee CH, Voukantsis D, Miar A, Arora A, Harris AL, Morten KJ, Buffa FM. Metabolic symbiosis between oxygenated and hypoxic tumour cells: An agent-based modelling study. PLoS Comput Biol 2024; 20:e1011944. [PMID: 38489376 DOI: 10.1371/journal.pcbi.1011944] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2023] [Revised: 03/27/2024] [Accepted: 02/24/2024] [Indexed: 03/17/2024] Open
Abstract
Deregulated metabolism is one of the hallmarks of cancer. It is well-known that tumour cells tend to metabolize glucose via glycolysis even when oxygen is available and mitochondrial respiration is functional. However, the lower energy efficiency of aerobic glycolysis with respect to mitochondrial respiration makes this behaviour, namely the Warburg effect, counter-intuitive, although it has now been recognized as source of anabolic precursors. On the other hand, there is evidence that oxygenated tumour cells could be fuelled by exogenous lactate produced from glycolysis. We employed a multi-scale approach that integrates multi-agent modelling, diffusion-reaction, stoichiometric equations, and Boolean networks to study metabolic cooperation between hypoxic and oxygenated cells exposed to varying oxygen, nutrient, and inhibitor concentrations. The results show that the cooperation reduces the depletion of environmental glucose, resulting in an overall advantage of using aerobic glycolysis. In addition, the oxygen level was found to be decreased by symbiosis, promoting a further shift towards anaerobic glycolysis. However, the oxygenated and hypoxic populations may gradually reach quasi-equilibrium. A sensitivity analysis using Latin hypercube sampling and partial rank correlation shows that the symbiotic dynamics depends on properties of the specific cell such as the minimum glucose level needed for glycolysis. Our results suggest that strategies that block glucose transporters may be more effective to reduce tumour growth than those blocking lactate intake transporters.
Collapse
Affiliation(s)
| | - Pedro Victori
- Department of Oncology, Medical Sciences Division, University of Oxford, Oxford, United Kingdom
| | - Clara E Pavillet
- Department of Oncology, Medical Sciences Division, University of Oxford, Oxford, United Kingdom
- MRC Weatherall Institute of Molecular Medicine, Radcliffe Department of Medicine, Medical Sciences Division, University of Oxford, Oxford, United Kingdom
- Department of Computing Sciences and Institute for Data Science and Analytics, Bocconi University, Milan, Italy
| | - Chang Heon Lee
- Department of Oncology, Medical Sciences Division, University of Oxford, Oxford, United Kingdom
| | - Dimitrios Voukantsis
- Department of Oncology, Medical Sciences Division, University of Oxford, Oxford, United Kingdom
| | - Ana Miar
- Department of Oncology, Medical Sciences Division, University of Oxford, Oxford, United Kingdom
| | - Anjali Arora
- Department of Oncology, Medical Sciences Division, University of Oxford, Oxford, United Kingdom
| | - Adrian L Harris
- Department of Oncology, Medical Sciences Division, University of Oxford, Oxford, United Kingdom
| | - Karl J Morten
- Nuffield Department of Women's and Reproductive Health, University of Oxford, Oxford, United Kingdom
| | - Francesca M Buffa
- Department of Oncology, Medical Sciences Division, University of Oxford, Oxford, United Kingdom
- Department of Computing Sciences and Institute for Data Science and Analytics, Bocconi University, Milan, Italy
| |
Collapse
|
2
|
Arfat Y, Zafar I, Sehgal SA, Ayaz M, Sajid M, Khan JM, Ahsan M, Rather MA, Khan AA, Alshehri JM, Akash S, Nepovimova E, Kuca K, Sharma R. In silico designing of multiepitope-based-peptide (MBP) vaccine against MAPK protein express for Alzheimer's disease in Zebrafish. Heliyon 2023; 9:e22204. [PMID: 38058625 PMCID: PMC10695983 DOI: 10.1016/j.heliyon.2023.e22204] [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: 05/17/2023] [Revised: 10/24/2023] [Accepted: 11/06/2023] [Indexed: 12/08/2023] Open
Abstract
Understanding the role of the mitogen-activated protein kinases (MAPKs) signalling pathway is essential in advancing treatments for neurodegenerative disorders like Alzheimer's. In this study, we investigate in-silico techniques involving computer-based methods to extract the MAPK1 sequence. Our applied methods enable us to analyze the protein's structure, evaluate its properties, establish its evolutionary relationships, and assess its prevalence in populations. We also predict epitopes, assess their ability to trigger immune responses, and check for allergenicity using advanced computational tools to understand their immunological properties comprehensively. We apply virtual screening, docking, and structure modelling to identify promising drug candidates, analyze their interactions, and enhance drug design processes. We identified a total of 30 cell-targeting molecules against the MAPK1 protein, where we selected top 10 CTL epitopes (PAGGGPNPG, GGGPNPGSG, SAPAGGGPN, AVSAPAGGG, AGGGPNPGS, ATAAVSAPA, TAAVSAPAG, ENIIGINDI, INDIIRTPT, and NDIIRTPTI) for further evaluation to determine their potential efficacy, safety, and suitability for vaccine design based on strong binding potential. The potential to cover a large portion of the world's population with these vaccines is substantial-88.5 % for one type and 99.99 % for another. In exploring the molecular docking analyses, we examined a library of compounds from the ZINC database. Among them, we identified twelve compounds with the lowest binding energy. Critical residues in the MAPK1 protein, such as VAL48, LYS63, CYS175, ASP176, LYS160, ALA61, LEU165, TYR45, SER162, ARG33, PRO365, PHE363, ILE40, ASN163, and GLU42, are pivotal for interactions with these compounds. Our result suggests that these compounds could influence the protein's behaviour. Moreover, our docking analyses revealed that the predicted peptides have a strong affinity for the MAPK1 protein. These peptides form stable complexes, indicating their potential as potent inhibitors. This study contributes to the identification of new drug compounds and the screening of their desired properties. These compounds could potentially help reduce the excessive activity of MAPK1, which is linked to Alzheimer's disease.
Collapse
Affiliation(s)
- Yasir Arfat
- Department of Biotechnology, Faculty of Life Sciences, University of Okara, Okara, 56300, Pakistan
| | - Imran Zafar
- Department of Bioinformatics and Computational Biology, Virtual University, Punjab, 54700, Pakistan
| | - Sheikh Arslan Sehgal
- Department of Bioinformatics, Institute of Biochemistry, Biotechnology and Bioinformatics, The Islamia University of Bahawalpur, Bahawalpur, Pakistan
| | - Mazhar Ayaz
- Department of Parasitology, Faculty of Veterinary Science, Cholistan University of Veterinary and Animal Sciences, Bahawalpur, Pakistan
| | - Muhammad Sajid
- Department of Biotechnology, Faculty of Life Sciences, University of Okara, Okara, 56300, Pakistan
| | - Jamal Muhammad Khan
- Department of Parasitology, Faculty of Veterinary Science, Cholistan University of Veterinary and Animal Sciences, Bahawalpur, Pakistan
| | - Muhammad Ahsan
- Department of Environmental Sciences, Institute of Environmental and Agricultural Sciences, University of Okara, Okara, 56300, Pakistan
| | - Mohd Ashraf Rather
- Division of Fish Genetics and Biotechnology, Faculty of Fisheries, Rangil- Gandarbal (SKAUST-K), India
| | - Azmat Ali Khan
- Pharmaceutical Biotechnology Laboratory, Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, Riyadh, 11451, Saudi Arabia
| | - Jamilah M. Alshehri
- Pharmaceutical Biotechnology Laboratory, Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, Riyadh, 11451, Saudi Arabia
| | - Shopnil Akash
- Department of Pharmacy, Faculty of Allied Health Sciences, Daffodil International, University, Dhaka, Bangladesh
| | - Eugenie Nepovimova
- Department of Chemistry, Faculty of Science, University of Hradec Králové, Hradec Králové, 50 003, Czech Republic
| | - Kamil Kuca
- Department of Chemistry, Faculty of Science, University of Hradec Králové, Hradec Králové, 50 003, Czech Republic
- Biomedical Research Center, University Hospital Hradec Kralove, 50005, Hradec Kralove, Czech Republic
| | - Rohit Sharma
- Department of Rasashastra and Bhaishajya Kalpana, Faculty of Ayurveda, Institute of Medical Sciences, Banaras Hindu University, Varanasi, India
| |
Collapse
|
3
|
Joo JI, Park H, Cho K. Normalizing Input-Output Relationships of Cancer Networks for Reversion Therapy. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2023; 10:e2207322. [PMID: 37269056 PMCID: PMC10460890 DOI: 10.1002/advs.202207322] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/12/2022] [Revised: 03/17/2023] [Indexed: 06/04/2023]
Abstract
Accumulated genetic alterations in cancer cells distort cellular stimulus-response (or input-output) relationships, resulting in uncontrolled proliferation. However, the complex molecular interaction network within a cell implicates a possibility of restoring such distorted input-output relationships by rewiring the signal flow through controlling hidden molecular switches. Here, a system framework of analyzing cellular input-output relationships in consideration of various genetic alterations and identifying possible molecular switches that can normalize the distorted relationships based on Boolean network modeling and dynamics analysis is presented. Such reversion is demonstrated by the analysis of a number of cancer molecular networks together with a focused case study on bladder cancer with in vitro experiments and patient survival data analysis. The origin of reversibility from an evolutionary point of view based on the redundancy and robustness intrinsically embedded in complex molecular regulatory networks is further discussed.
Collapse
Affiliation(s)
- Jae Il Joo
- Department of Bio and Brain EngineeringKorea Advanced Institute of Science and Technology (KAIST)Daejeon34141Republic of Korea
- Present address:
biorevert IncDaejeon34051Republic of Korea
| | - Hwa‐Jeong Park
- Department of Bio and Brain EngineeringKorea Advanced Institute of Science and Technology (KAIST)Daejeon34141Republic of Korea
- Present address:
Promega Corporationan affiliate of PromegaSouth Korea
| | - Kwang‐Hyun Cho
- Department of Bio and Brain EngineeringKorea Advanced Institute of Science and Technology (KAIST)Daejeon34141Republic of Korea
| |
Collapse
|
4
|
Argyris GA, Lluch Lafuente A, Tribastone M, Tschaikowski M, Vandin A. Reducing Boolean networks with backward equivalence. BMC Bioinformatics 2023; 24:212. [PMID: 37221494 DOI: 10.1186/s12859-023-05326-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2022] [Accepted: 05/05/2023] [Indexed: 05/25/2023] Open
Abstract
BACKGROUND Boolean Networks (BNs) are a popular dynamical model in biology where the state of each component is represented by a variable taking binary values that express, for instance, activation/deactivation or high/low concentrations. Unfortunately, these models suffer from the state space explosion, i.e., there are exponentially many states in the number of BN variables, which hampers their analysis. RESULTS We present Boolean Backward Equivalence (BBE), a novel reduction technique for BNs which collapses system variables that, if initialized with same value, maintain matching values in all states. A large-scale validation on 86 models from two online model repositories reveals that BBE is effective, since it is able to reduce more than 90% of the models. Furthermore, on such models we also show that BBE brings notable analysis speed-ups, both in terms of state space generation and steady-state analysis. In several cases, BBE allowed the analysis of models that were originally intractable due to the complexity. On two selected case studies, we show how one can tune the reduction power of BBE using model-specific information to preserve all dynamics of interest, and selectively exclude behavior that does not have biological relevance. CONCLUSIONS BBE complements existing reduction methods, preserving properties that other reduction methods fail to reproduce, and vice versa. BBE drops all and only the dynamics, including attractors, originating from states where BBE-equivalent variables have been initialized with different activation values The remaining part of the dynamics is preserved exactly, including the length of the preserved attractors, and their reachability from given initial conditions, without adding any spurious behaviours. Given that BBE is a model-to-model reduction technique, it can be combined with further reduction methods for BNs.
Collapse
Affiliation(s)
- Georgios A Argyris
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Lyngby, Denmark
| | - Alberto Lluch Lafuente
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Lyngby, Denmark
| | | | - Max Tschaikowski
- Department of Computer Science, University of Aalborg, Aalborg, Denmark
| | - Andrea Vandin
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Lyngby, Denmark.
- Department of Excellence EMbeDS and Institute of Economics, Sant'Anna School for Advanced Studies, Pisa, Italy.
| |
Collapse
|
5
|
Kim J, Hopper C, Cho KH. Statistical control of structural networks with limited interventions to minimize cellular phenotypic diversity represented by point attractors. Sci Rep 2023; 13:6275. [PMID: 37072458 PMCID: PMC10113376 DOI: 10.1038/s41598-023-33346-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2022] [Accepted: 04/12/2023] [Indexed: 05/03/2023] Open
Abstract
The underlying genetic networks of cells give rise to diverse behaviors known as phenotypes. Control of this cellular phenotypic diversity (CPD) may reveal key targets that govern differentiation during development or drug resistance in cancer. This work establishes an approach to control CPD that encompasses practical constraints, including model limitations, the number of simultaneous control targets, which targets are viable for control, and the granularity of control. Cellular networks are often limited to the structure of interactions, due to the practical difficulty of modeling interaction dynamics. However, these dynamics are essential to CPD. In response, our statistical control approach infers the CPD directly from the structure of a network, by considering an ensemble average function over all possible Boolean dynamics for each node in the network. These ensemble average functions are combined with an acyclic form of the network to infer the number of point attractors. Our approach is applied to several known biological models and shown to outperform existing approaches. Statistical control of CPD offers a new avenue to contend with systemic processes such as differentiation and cancer, despite practical limitations in the field.
Collapse
Affiliation(s)
- Jongwan Kim
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| | - Corbin Hopper
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| | - Kwang-Hyun Cho
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea.
| |
Collapse
|
6
|
An S, Jang SY, Park SM, Lee CK, Kim HM, Cho KH. Global stabilizing control of large-scale biomolecular regulatory networks. Bioinformatics 2023; 39:6998201. [PMID: 36688702 PMCID: PMC9891247 DOI: 10.1093/bioinformatics/btad045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2022] [Revised: 01/13/2023] [Accepted: 01/20/2023] [Indexed: 01/24/2023] Open
Abstract
MOTIVATION Cellular behavior is determined by complex non-linear interactions between numerous intracellular molecules that are often represented by Boolean network models. To achieve a desired cellular behavior with minimal intervention, we need to identify optimal control targets that can drive heterogeneous cellular states to the desired phenotypic cellular state with minimal node intervention. Previous attempts to realize such global stabilization were based solely on either network structure information or simple linear dynamics. Other attempts based on non-linear dynamics are not scalable. RESULTS Here, we investigate the underlying relationship between structurally identified control targets and optimal global stabilizing control targets based on non-linear dynamics. We discovered that optimal global stabilizing control targets can be identified by analyzing the dynamics between structurally identified control targets. Utilizing these findings, we developed a scalable global stabilizing control framework using both structural and dynamic information. Our framework narrows down the search space based on strongly connected components and feedback vertex sets then identifies global stabilizing control targets based on the canalization of Boolean network dynamics. We find that the proposed global stabilizing control is superior with respect to the number of control target nodes, scalability, and computational complexity. AVAILABILITY AND IMPLEMENTATION We provide a GitHub repository that contains the DCGS framework written in Python as well as biological random Boolean network datasets (https://github.com/sugyun/DCGS). SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
| | | | | | - Chun-Kyung Lee
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Hoon-Min Kim
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | | |
Collapse
|
7
|
Su C, Pang J. Target Control of Asynchronous Boolean Networks. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2023; 20:707-719. [PMID: 34882560 DOI: 10.1109/tcbb.2021.3133608] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
We study the target control of asynchronous Boolean networks, to identify interventions that can drive the dynamics of a given Boolean network from any initial state to the desired target attractor. Based on the application time, the control can be realised with three types of perturbations, including instantaneous, temporary and permanent perturbations. We develop efficient methods to compute the target control for a given target attractor with these three types of perturbations. We compare our methods with the stable motif-based control method on a variety of real-life biological networks to evaluate their performance. We show that our methods scale well for large Boolean networks and they are able to identify a rich set of solutions with a small number of perturbations.
Collapse
|
8
|
Temporary and permanent control of partially specified Boolean networks. Biosystems 2023; 223:104795. [PMID: 36377120 DOI: 10.1016/j.biosystems.2022.104795] [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: 02/09/2022] [Revised: 10/16/2022] [Accepted: 10/19/2022] [Indexed: 01/11/2023]
Abstract
Boolean networks (BNs) are a well-accepted modelling formalism in computational systems biology. Nevertheless, modellers often cannot identify only a single BN that matches the biological reality. The typical reasons for this is insufficient knowledge or a lack of experimental data. Formally, this uncertainty can be expressed using partially specified Boolean networks (PSBNs), which encode the wide range of network candidates into a single structure. In this paper, we target the control of PSBNs. The goal of BN control is to find perturbations which guarantee stabilisation of the system in the desired state. Specifically, we consider variable perturbations (gene knock-out and over-expression) with three types of application time-window: one-step, temporary, and permanent. While the control of fully specified BNs is a thoroughly explored topic, control of PSBNs introduces additional challenges that we address in this paper. In particular, the unspecified components of the model cause a significant amount of additional state space explosion. To address this issue, we propose a fully symbolic methodology that can represent the numerous system variants in a compact form. In fully specified models, the efficiency of a perturbation is characterised by the count of perturbed variables (the perturbation size). However, in the case of a PSBN, a perturbation might work only for a subset of concrete BN models. To that end, we introduce and quantify perturbation robustness. This metric characterises how efficient the given perturbation is with respect to the model uncertainty. Finally, we evaluate the novel control methods using non-trivial real-world PSBN models. We inspect the method's scalability and efficiency with respect to the size of the state space and the number of unspecified components. We also compare the robustness metrics for all three perturbation types. Our experiments support the hypothesis that one-step perturbations are significantly less robust than temporary and permanent ones.
Collapse
|
9
|
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
|
10
|
Beneš N, Brim L, Kadlecaj J, Pastva S, Šafránek D. Exploring attractor bifurcations in Boolean networks. BMC Bioinformatics 2022; 23:173. [PMID: 35546394 PMCID: PMC9092939 DOI: 10.1186/s12859-022-04708-9] [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: 11/23/2021] [Accepted: 04/19/2022] [Indexed: 11/10/2022] Open
Abstract
Background Boolean networks (BNs) provide an effective modelling formalism for various complex biochemical phenomena. Their long term behaviour is represented by attractors–subsets of the state space towards which the BN eventually converges. These are then typically linked to different biological phenotypes. Depending on various logical parameters, the structure and quality of attractors can undergo a significant change, known as a bifurcation. We present a methodology for analysing bifurcations in asynchronous parametrised Boolean networks. Results In this paper, we propose a computational framework employing advanced symbolic graph algorithms that enable the analysis of large networks with hundreds of Boolean variables. To visualise the results of this analysis, we developed a novel interactive presentation technique based on decision trees, allowing us to quickly uncover parameters crucial to the changes in the attractor landscape. As a whole, the methodology is implemented in our tool AEON. We evaluate the method’s applicability on a complex human cell signalling network describing the activity of type-1 interferons and related molecules interacting with SARS-COV-2 virion. In particular, the analysis focuses on explaining the potential suppressive role of the recently proposed drug molecule GRL0617 on replication of the virus. Conclusions The proposed method creates a working analogy to the concept of bifurcation analysis widely used in kinetic modelling to reveal the impact of parameters on the system’s stability. The important feature of our tool is its unique capability to work fast with large-scale networks with a relatively large extent of unknown information. The results obtained in the case study are in agreement with the recent biological findings.
Collapse
Affiliation(s)
- Nikola Beneš
- Faculty of Informatics, Masaryk University, Brno, Czechia.
| | - Luboš Brim
- Faculty of Informatics, Masaryk University, Brno, Czechia
| | - Jakub Kadlecaj
- Faculty of Informatics, Masaryk University, Brno, Czechia
| | - Samuel Pastva
- Faculty of Informatics, Masaryk University, Brno, Czechia
| | - David Šafránek
- Faculty of Informatics, Masaryk University, Brno, Czechia
| |
Collapse
|
11
|
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
|
12
|
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:72626. [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] [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)
| | - Jonas Béal
- Institut Curie, PSL Research University, Paris, France
| | - Luis Tobalina
- Faculty of Medicine, Joint Research Centre for Computational Biomedicine (JRC-COMBINE), RWTH Aachen University, Aachen, Germany
| | | | - Vigneshwari Subramanian
- Faculty of Medicine, Joint Research Centre for Computational Biomedicine (JRC-COMBINE), RWTH Aachen University, Aachen, Germany
| | - Bence Szalai
- Department of Physiology, Semmelweis University, Budapest, Hungary
| | | | | | | | | | - Julio Saez-Rodriguez
- Institute of Computational Biomedicine, Heidelberg University, Heidelberg, Germany
| | | |
Collapse
|
13
|
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]
|
14
|
Gondal MN, Chaudhary SU. Navigating Multi-Scale Cancer Systems Biology Towards Model-Driven Clinical Oncology and Its Applications in Personalized Therapeutics. Front Oncol 2021; 11:712505. [PMID: 34900668 PMCID: PMC8652070 DOI: 10.3389/fonc.2021.712505] [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: 05/20/2021] [Accepted: 10/26/2021] [Indexed: 12/19/2022] Open
Abstract
Rapid advancements in high-throughput omics technologies and experimental protocols have led to the generation of vast amounts of scale-specific biomolecular data on cancer that now populates several online databases and resources. Cancer systems biology models built using this data have the potential to provide specific insights into complex multifactorial aberrations underpinning tumor initiation, development, and metastasis. Furthermore, the annotation of these single- and multi-scale models with patient data can additionally assist in designing personalized therapeutic interventions as well as aid in clinical decision-making. Here, we have systematically reviewed the emergence and evolution of (i) repositories with scale-specific and multi-scale biomolecular cancer data, (ii) systems biology models developed using this data, (iii) associated simulation software for the development of personalized cancer therapeutics, and (iv) translational attempts to pipeline multi-scale panomics data for data-driven in silico clinical oncology. The review concludes that the absence of a generic, zero-code, panomics-based multi-scale modeling pipeline and associated software framework, impedes the development and seamless deployment of personalized in silico multi-scale models in clinical settings.
Collapse
Affiliation(s)
- Mahnoor Naseer Gondal
- Biomedical Informatics Research Laboratory, Department of Biology, Syed Babar Ali School of Science and Engineering, Lahore University of Management Sciences, Lahore, Pakistan
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, United States
| | - Safee Ullah Chaudhary
- Biomedical Informatics Research Laboratory, Department of Biology, Syed Babar Ali School of Science and Engineering, Lahore University of Management Sciences, Lahore, Pakistan
| |
Collapse
|
15
|
Selvaggio G, Chaouiya C, Janody F. In Silico Logical Modelling to Uncover Cooperative Interactions in Cancer. Int J Mol Sci 2021; 22:ijms22094897. [PMID: 34063110 PMCID: PMC8125147 DOI: 10.3390/ijms22094897] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2021] [Revised: 04/29/2021] [Accepted: 04/30/2021] [Indexed: 12/13/2022] Open
Abstract
The multistep development of cancer involves the cooperation between multiple molecular lesions, as well as complex interactions between cancer cells and the surrounding tumour microenvironment. The search for these synergistic interactions using experimental models made tremendous contributions to our understanding of oncogenesis. Yet, these approaches remain labour-intensive and challenging. To tackle such a hurdle, an integrative, multidisciplinary effort is required. In this article, we highlight the use of logical computational models, combined with experimental validations, as an effective approach to identify cooperative mechanisms and therapeutic strategies in the context of cancer biology. In silico models overcome limitations of reductionist approaches by capturing tumour complexity and by generating powerful testable hypotheses. We review representative examples of logical models reported in the literature and their validation. We then provide further analyses of our logical model of Epithelium to Mesenchymal Transition (EMT), searching for additional cooperative interactions involving inputs from the tumour microenvironment and gain of function mutations in NOTCH.
Collapse
Affiliation(s)
- Gianluca Selvaggio
- Fondazione the Microsoft Research—University of Trento Centre for Computational and Systems Biology (COSBI), Piazza Manifattura 1, 38068 Rovereto, Italy;
- Instituto Gulbenkian de Ciência, Rua da Quinta Grande 6, 2780-156 Oeiras, Portugal
| | - Claudine Chaouiya
- Instituto Gulbenkian de Ciência, Rua da Quinta Grande 6, 2780-156 Oeiras, Portugal
- CNRS, Centrale Marseille, I2M, Aix Marseille University, 13397 Marseille, France
- Correspondence: (C.C.); (F.J.)
| | - Florence Janody
- Instituto Gulbenkian de Ciência, Rua da Quinta Grande 6, 2780-156 Oeiras, Portugal
- i3S—Instituto de Investigação e Inovação em Saúde, Universidade do Porto, Rua Alfredo Allen, 208, 4200-135 Porto, Portugal
- IPATIMUP—Instituto de Patologia e Imunologia Molecular da Universidade do Porto, Rua Dr. Roberto Frias s/n, 4200-465 Porto, Portugal
- Correspondence: (C.C.); (F.J.)
| |
Collapse
|
16
|
Aghamiri SS, Singh V, Naldi A, Helikar T, Soliman S, Niarakis A. Automated inference of Boolean models from molecular interaction maps using CaSQ. Bioinformatics 2021; 36:4473-4482. [PMID: 32403123 PMCID: PMC7575051 DOI: 10.1093/bioinformatics/btaa484] [Citation(s) in RCA: 32] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 04/17/2020] [Accepted: 05/06/2020] [Indexed: 12/16/2022] Open
Abstract
Motivation Molecular interaction maps have emerged as a meaningful way of representing biological mechanisms in a comprehensive and systematic manner. However, their static nature provides limited insights to the emerging behaviour of the described biological system under different conditions. Computational modelling provides the means to study dynamic properties through in silico simulations and perturbations. We aim to bridge the gap between static and dynamic representations of biological systems with CaSQ, a software tool that infers Boolean rules based on the topology and semantics of molecular interaction maps built with CellDesigner. Results We developed CaSQ by defining conversion rules and logical formulas for inferred Boolean models according to the topology and the annotations of the starting molecular interaction maps. We used CaSQ to produce executable files of existing molecular maps that differ in size, complexity and the use of Systems Biology Graphical Notation (SBGN) standards. We also compared, where possible, the manually built logical models corresponding to a molecular map to the ones inferred by CaSQ. The tool is able to process large and complex maps built with CellDesigner (either following SBGN standards or not) and produce Boolean models in a standard output format, Systems Biology Marked Up Language-qualitative (SBML-qual), that can be further analyzed using popular modelling tools. References, annotations and layout of the CellDesigner molecular map are retained in the obtained model, facilitating interoperability and model reusability. Availability and implementation The present tool is available online: https://lifeware.inria.fr/∼soliman/post/casq/ and distributed as a Python package under the GNU GPLv3 license. The code can be accessed here: https://gitlab.inria.fr/soliman/casq. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Sara Sadat Aghamiri
- GenHotel, Département de Biologie, Univ. èvry, Université Paris-Saclay, Genopole, èvry 91025, France
| | - Vidisha Singh
- GenHotel, Département de Biologie, Univ. èvry, Université Paris-Saclay, Genopole, èvry 91025, France
| | - Aurélien Naldi
- Département de Biologie, Institut de Biologie de l'Ecole Normale Supérieure (IBENS), ècole Normale Supérieure, CNRS, INSERM, Université PSL, Paris 75005, France
| | - Tomáš Helikar
- Department of Biochemistry, University of Nebraska-Lincoln, Lincoln, NE 68588, USA
| | - Sylvain Soliman
- Lifeware Group, Inria Saclay-île de France, Palaiseau 91120, France
| | - Anna Niarakis
- GenHotel, Département de Biologie, Univ. èvry, Université Paris-Saclay, Genopole, èvry 91025, France
| |
Collapse
|
17
|
Kołat D, Kałuzińska Ż, Bednarek AK, Płuciennik E. Fragile Gene WWOX Guides TFAP2A/ TFAP2C-Dependent Actions Against Tumor Progression in Grade II Bladder Cancer. Front Oncol 2021; 11:621060. [PMID: 33718178 PMCID: PMC7947623 DOI: 10.3389/fonc.2021.621060] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2020] [Accepted: 01/18/2021] [Indexed: 02/05/2023] Open
Abstract
INTRODUCTION The presence of common fragile sites is associated with no-accidental chromosomal instability which occurs prior to carcinogenesis. The WWOX gene spans the second most active fragile site: FRA16D. Chromosomal breakage at this site is more common in bladder cancer patients who are tobacco smokers which suggests the importance of WWOX gene loss regarding bladder carcinogenesis. Tryptophan domains of WWOX are known to recognize motifs of other proteins such as AP-2α and AP-2γ allowing protein-protein interactions. While the roles of both AP-2 transcription factors are important for bladder carcinogenesis, their nature is different. Based on the literature, AP-2γ appears to be oncogenic, whereas AP-2α mainly exhibits tumor suppressor character. Presumably, the interaction between WWOX and both transcription factors regulates thousands of genes, hence the aim of the present study was to determine WWOX, AP-2α, and AP-2γ function in modulating biological processes of bladder cancer. METHODS RT-112 cell line (grade II bladder cancer) was subjected to two stable lentiviral transductions. Overall, this resulted in six variants to investigate distinct WWOX, AP-2α, or AP-2γ function as well as WWOX in collaboration with a particular transcription factor. Cellular models were examined with immunocytochemical staining and in terms of differences in biological processes using assays investigating cell viability, proliferation, apoptosis, adhesion, clonogenicity, migration, activity of metalloproteinases and 3D culture growth. RESULTS WWOX overexpression increased apoptosis but decreased cell viability, migration and large spatial colonies. AP-2α overexpression decreased tumor cell viability, migratory potential, matrix metalloproteinase-2 activity and clonogenicity. AP-2γ overexpression decreased matrix metalloproteinase-2 activity but increased wound healing, adhesion, clonogenicity and spatial colony formation. WWOX and AP-2α overexpression induced apoptosis but decreased cell viability, adhesion, matrix metalloproteinase-2 activity, overall number of cultured colonies and migration rate. WWOX and AP-2γ overexpression decreased tumor cell viability, proliferation potential, adhesion, clonogenicity and the ability to create spatial structures, but also increased apoptosis or migration rate. CONCLUSION Co-overexpression of WWOX with AP-2α or WWOX with AP-2γ resulted in a net anti-tumor effect. However, considering this research findings and the difference between AP-2α and AP-2γ, we suggest that this similarity is due to a divergent behavior of WWOX.
Collapse
|
18
|
Choo SM, Almomani LM, Cho KH. Boolean Feedforward Neural Network Modeling of Molecular Regulatory Networks for Cellular State Conversion. Front Physiol 2020; 11:594151. [PMID: 33335489 PMCID: PMC7736109 DOI: 10.3389/fphys.2020.594151] [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: 08/12/2020] [Accepted: 11/03/2020] [Indexed: 11/13/2022] Open
Abstract
The molecular regulatory network (MRN) within a cell determines cellular states and transitions between them. Thus, modeling of MRNs is crucial, but this usually requires extensive analysis of time-series measurements, which is extremely difficult to obtain from biological experiments. However, single-cell measurement data such as single-cell RNA-sequencing databases have recently provided a new insight into resolving this problem by ordering thousands of cells in pseudo-time according to their differential gene expressions. Neural network modeling can be employed by using temporal data as learning data. In contrast, Boolean network modeling of MRNs has a growing interest, as it is a parameter-free logical modeling and thereby robust to noisy data while still capturing essential dynamics of biological networks. In this study, we propose a Boolean feedforward neural network (FFN) modeling by combining neural network and Boolean network modeling approach to reconstruct a practical and useful MRN model from large temporal data. Furthermore, analyzing the reconstructed MRN model can enable us to identify control targets for potential cellular state conversion. Here, we show the usefulness of Boolean FFN modeling by demonstrating its applicability through a toy model and biological networks.
Collapse
Affiliation(s)
- Sang-Mok Choo
- Department of Mathematics, University of Ulsan, Ulsan, South Korea
| | | | - Kwang-Hyun Cho
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, South Korea
| |
Collapse
|
19
|
Manica M, Polig R, Purandare M, Mathis R, Hagleitner C, Martinez MR. FPGA Accelerated Analysis of Boolean Gene Regulatory Networks. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2020; 17:2141-2147. [PMID: 31494553 DOI: 10.1109/tcbb.2019.2936836] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Boolean models are a powerful abstraction for qualitative modeling of gene regulatory networks. With the recent availability of advanced high-throughput technologies, Boolean models have increasingly grown in size and complexity, posing a challenge for existing software simulation tools that have not scaled at the same speed. Field Programmable Gate Arrays (FPGAs) are powerful reconfigurable integrated circuits that can offer massive performance improvements. Due to their highly parallel nature, FPGAs are well suited to simulate complex molecular networks. We present here a new simulation framework for Boolean models, which first converts the model to Verilog, a standardized hardware description language, and then connects it to an execution core that runs on an FPGA coherently attached to a POWER8 processor. We report an order of magnitude speedup over a multi-threaded software simulation tool running on the same processor on a selection of Boolean models. Analysis on a T-cell large granular lymphocyte leukemia (T-LGL) demonstrates that our framework achieves consistent performance improvements resulting in new biological insights. In addition, we show that our solution allows to perform attractor detection at an unprecedented speed, exhibiting a speedup ranging from one to three orders of magnitude compared to alternative software solutions.
Collapse
|
20
|
Tsirvouli E, Touré V, Niederdorfer B, Vázquez M, Flobak Å, Kuiper M. A Middle-Out Modeling Strategy to Extend a Colon Cancer Logical Model Improves Drug Synergy Predictions in Epithelial-Derived Cancer Cell Lines. Front Mol Biosci 2020; 7:502573. [PMID: 33195403 PMCID: PMC7581946 DOI: 10.3389/fmolb.2020.502573] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2019] [Accepted: 09/22/2020] [Indexed: 11/23/2022] Open
Abstract
Cancer is a heterogeneous and complex disease and one of the leading causes of death worldwide. The high tumor heterogeneity between individuals affected by the same cancer type is accompanied by distinct molecular and phenotypic tumor profiles and variation in drug treatment response. In silico modeling of cancer as an aberrantly regulated system of interacting signaling molecules provides a basis to enhance our biological understanding of disease progression, and it offers the means to use computer simulations to test and optimize drug therapy designs on particular cancer types and subtypes. This sets the stage for precision medicine: the design of treatments tailored to individuals or groups of patients based on their tumor-specific molecular cancer profiles. Here, we show how a relatively large manually curated logical model can be efficiently enhanced further by including components highlighted by a multi-omics data analysis of data from Consensus Molecular Subtypes covering colorectal cancer. The model expansion was performed in a pathway-centric manner, following a partitioning of the model into functional subsystems, named modules. The resulting approach constitutes a middle-out modeling strategy enabling a data-driven expansion of a model from a generic and intermediate level of molecular detail to a model better covering relevant processes that are affected in specific cancer subtypes, comprising 183 biological entities and 603 interactions between them, partitioned in 25 functional modules of varying size and structure. We tested this model for its ability to correctly predict drug combination synergies, against a dataset of experimentally determined cell growth responses with 18 drugs in all combinations, on eight cancer cell lines. The results indicate that the extended model had an improved accuracy for drug synergy prediction for the majority of the experimentally tested cancer cell lines, although significant improvements of the model's predictive performance are still needed. Our study demonstrates how a tumor-data driven middle-out approach toward refining a logical model of a biological system can further customize a computer model to represent specific cancer cell lines and provide a basis for identifying synergistic effects of drugs targeting specific regulatory proteins. This approach bridges between preclinical cancer model data and clinical patient data and may thereby ultimately be of help to develop patient-specific in silico models that can steer treatment decisions in the clinic.
Collapse
Affiliation(s)
- Eirini Tsirvouli
- Department of Biology, Norwegian University of Science and Technology, Trondheim, Norway
| | - Vasundra Touré
- Department of Biology, Norwegian University of Science and Technology, Trondheim, Norway
| | - Barbara Niederdorfer
- Department of Clinical and Molecular Medicine, Norwegian University of Science and Technology, Trondheim, Norway
| | - Miguel Vázquez
- Department of Clinical and Molecular Medicine, Norwegian University of Science and Technology, Trondheim, Norway
| | - Åsmund Flobak
- Department of Clinical and Molecular Medicine, Norwegian University of Science and Technology, Trondheim, Norway
- The Cancer Clinic, St. Olav’s University Hospital, Trondheim, Norway
| | - Martin Kuiper
- Department of Biology, Norwegian University of Science and Technology, Trondheim, Norway
| |
Collapse
|
21
|
Hernandez C, Thomas-Chollier M, Naldi A, Thieffry D. Computational Verification of Large Logical Models-Application to the Prediction of T Cell Response to Checkpoint Inhibitors. Front Physiol 2020; 11:558606. [PMID: 33101049 PMCID: PMC7554341 DOI: 10.3389/fphys.2020.558606] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2020] [Accepted: 08/19/2020] [Indexed: 12/31/2022] Open
Abstract
At the crossroad between biology and mathematical modeling, computational systems biology can contribute to a mechanistic understanding of high-level biological phenomenon. But as knowledge accumulates, the size and complexity of mathematical models increase, calling for the development of efficient dynamical analysis methods. Here, we propose the use of two approaches for the development and analysis of complex cellular network models. A first approach, called "model verification" and inspired by unitary testing in software development, enables the formalization and automated verification of validation criteria for whole models or selected sub-parts. When combined with efficient analysis methods, this approach is suitable for continuous testing, thereby greatly facilitating model development. A second approach, called "value propagation," enables efficient analytical computation of the impact of specific environmental or genetic conditions on the dynamical behavior of some models. We apply these two approaches to the delineation and the analysis of a comprehensive model for T cell activation, taking into account CTLA4 and PD-1 checkpoint inhibitory pathways. While model verification greatly eases the delineation of logical rules complying with a set of dynamical specifications, propagation provides interesting insights into the different potential of CTLA4 and PD-1 immunotherapies. Both methods are implemented and made available in the all-inclusive CoLoMoTo Docker image, while the different steps of the model analysis are fully reported in two companion interactive jupyter notebooks, thereby ensuring the reproduction of our results.
Collapse
Affiliation(s)
- Céline Hernandez
- Institut de Biologie de l'ENS (IBENS), Département de Biologie, École Normale Supérieure, CNRS, INSERM, Université PSL, Paris, France
| | - Morgane Thomas-Chollier
- Institut de Biologie de l'ENS (IBENS), Département de Biologie, École Normale Supérieure, CNRS, INSERM, Université PSL, Paris, France.,Institut Universitaire de France, Paris, France
| | - Aurélien Naldi
- Institut de Biologie de l'ENS (IBENS), Département de Biologie, École Normale Supérieure, CNRS, INSERM, Université PSL, Paris, France
| | - Denis Thieffry
- Institut de Biologie de l'ENS (IBENS), Département de Biologie, École Normale Supérieure, CNRS, INSERM, Université PSL, Paris, France
| |
Collapse
|
22
|
Klarner H, Heinitz F, Nee S, Siebert H. Basins of Attraction, Commitment Sets, and Phenotypes of Boolean Networks. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2020; 17:1115-1124. [PMID: 30575543 DOI: 10.1109/tcbb.2018.2879097] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The attractors of Boolean networks and their basins have been shown to be highly relevant for model validation and predictive modeling, e.g., in systems biology. Yet, there are currently very few tools available that are able to compute and visualize not only attractors but also their basins. In the realm of asynchronous, non-deterministic modeling not only is the repertoire of software even more limited, but also the formal notions for basins of attraction are often lacking. In this setting, the difficulty both for theory and computation arises from the fact that states may be elements of several distinct basins. In this paper, we address this topic by partitioning the state space into sets that are committed to the same attractors. These commitment sets can easily be generalized to sets that are equivalent w.r.t. the long-term behaviors of pre-selected nodes which leads us to the notions of markers and phenotypes which we illustrate in a case study on bladder tumorigenesis. For every concept, we propose equivalent CTL model checking queries and an extension of the state of the art model checking software NuSMV is made available that is capable of computing the respective sets. All notions are fully integrated as three new modules in our Python package PyBoolNet, including functions for visualizing the basins, commitment sets, and phenotypes as quotient graphs and pie charts.
Collapse
|
23
|
Montagud A, Traynard P, Martignetti L, Bonnet E, Barillot E, Zinovyev A, Calzone L. Conceptual and computational framework for logical modelling of biological networks deregulated in diseases. Brief Bioinform 2020; 20:1238-1249. [PMID: 29237040 DOI: 10.1093/bib/bbx163] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2017] [Revised: 10/24/2017] [Indexed: 01/02/2023] Open
Abstract
Mathematical models can serve as a tool to formalize biological knowledge from diverse sources, to investigate biological questions in a formal way, to test experimental hypotheses, to predict the effect of perturbations and to identify underlying mechanisms. We present a pipeline of computational tools that performs a series of analyses to explore a logical model's properties. A logical model of initiation of the metastatic process in cancer is used as a transversal example. We start by analysing the structure of the interaction network constructed from the literature or existing databases. Next, we show how to translate this network into a mathematical object, specifically a logical model, and how robustness analyses can be applied to it. We explore the visualization of the stable states, defined as specific attractors of the model, and match them to cellular fates or biological read-outs. With the different tools we present here, we explain how to assign to each solution of the model a probability and how to identify genetic interactions using mutant phenotype probabilities. Finally, we connect the model to relevant experimental data: we present how some data analyses can direct the construction of the network, and how the solutions of a mathematical model can also be compared with experimental data, with a particular focus on high-throughput data in cancer biology. A step-by-step tutorial is provided as a Supplementary Material and all models, tools and scripts are provided on an accompanying website: https://github.com/sysbio-curie/Logical_modelling_pipeline.
Collapse
|
24
|
Schwab JD, Kühlwein SD, Ikonomi N, Kühl M, Kestler HA. Concepts in Boolean network modeling: What do they all mean? Comput Struct Biotechnol J 2020; 18:571-582. [PMID: 32257043 PMCID: PMC7096748 DOI: 10.1016/j.csbj.2020.03.001] [Citation(s) in RCA: 60] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2019] [Revised: 01/27/2020] [Accepted: 03/01/2020] [Indexed: 12/02/2022] Open
Abstract
Boolean network models are one of the simplest models to study complex dynamic behavior in biological systems. They can be applied to unravel the mechanisms regulating the properties of the system or to identify promising intervention targets. Since its introduction by Stuart Kauffman in 1969 for describing gene regulatory networks, various biologically based networks and tools for their analysis were developed. Here, we summarize and explain the concepts for Boolean network modeling. We also present application examples and guidelines to work with and analyze Boolean network models.
Collapse
Affiliation(s)
- Julian D Schwab
- Institute of Medical Systems Biology, Ulm University, Albert-Einstein-Allee 11, 89081 Ulm, Germany
| | - Silke D Kühlwein
- Institute of Medical Systems Biology, Ulm University, Albert-Einstein-Allee 11, 89081 Ulm, Germany
| | - Nensi Ikonomi
- Institute of Medical Systems Biology, Ulm University, Albert-Einstein-Allee 11, 89081 Ulm, Germany
| | - Michael Kühl
- Institute of Biochemistry and Molecular Biology, Ulm University, Albert-Einstein-Allee 11, 89081 Ulm, Germany
| | - Hans A Kestler
- Institute of Medical Systems Biology, Ulm University, Albert-Einstein-Allee 11, 89081 Ulm, Germany
| |
Collapse
|
25
|
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
|
26
|
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
|
27
|
Singh V, Kalliolias GD, Ostaszewski M, Veyssiere M, Pilalis E, Gawron P, Mazein A, Bonnet E, Petit-Teixeira E, Niarakis A. RA-map: building a state-of-the-art interactive knowledge base for rheumatoid arthritis. Database (Oxford) 2020; 2020:baaa017. [PMID: 32311035 PMCID: PMC7170216 DOI: 10.1093/database/baaa017] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2019] [Revised: 01/21/2020] [Accepted: 02/13/2020] [Indexed: 02/07/2023]
Abstract
Rheumatoid arthritis (RA) is a progressive, inflammatory autoimmune disease of unknown aetiology. The complex mechanism of aetiopathogenesis, progress and chronicity of the disease involves genetic, epigenetic and environmental factors. To understand the molecular mechanisms underlying disease phenotypes, one has to place implicated factors in their functional context. However, integration and organization of such data in a systematic manner remains a challenging task. Molecular maps are widely used in biology to provide a useful and intuitive way of depicting a variety of biological processes and disease mechanisms. Recent large-scale collaborative efforts such as the Disease Maps Project demonstrate the utility of such maps as versatile tools to organize and formalize disease-specific knowledge in a comprehensive way, both human and machine-readable. We present a systematic effort to construct a fully annotated, expert validated, state-of-the-art knowledge base for RA in the form of a molecular map. The RA map illustrates molecular and signalling pathways implicated in the disease. Signal transduction is depicted from receptors to the nucleus using the Systems Biology Graphical Notation (SBGN) standard representation. High-quality manual curation, use of only human-specific studies and focus on small-scale experiments aim to limit false positives in the map. The state-of-the-art molecular map for RA, using information from 353 peer-reviewed scientific publications, comprises 506 species, 446 reactions and 8 phenotypes. The species in the map are classified to 303 proteins, 61 complexes, 106 genes, 106 RNA entities, 2 ions and 7 simple molecules. The RA map is available online at ramap.elixir-luxembourg.org as an open-access knowledge base allowing for easy navigation and search of molecular pathways implicated in the disease. Furthermore, the RA map can serve as a template for omics data visualization.
Collapse
Affiliation(s)
- Vidisha Singh
- Laboratoire Européen de Recherche pour la Polyarthrite Rhumatoïde - Genhotel, Univ Evry, Université Paris-Saclay, 2, rue Gaston Crémieux, 91057 EVRY-GENOPOLE cedex, Evry, France
| | - George D Kalliolias
- Arthritis and Tissue Degeneration Program, Hospital for Special Surgery, 535 East 70th Street, New York, NY 10021, USA
- Weill Cornell Medical Center, Weill Department of Medicine, 525 East 68th Street, New York, NY 10065, USA
| | - Marek Ostaszewski
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, 6 Avenue du Swing, L-4367 Belvaux, Luxembourg
| | - Maëva Veyssiere
- Laboratoire Européen de Recherche pour la Polyarthrite Rhumatoïde - Genhotel, Univ Evry, Université Paris-Saclay, 2, rue Gaston Crémieux, 91057 EVRY-GENOPOLE cedex, Evry, France
| | - Eleftherios Pilalis
- eNIOS Applications P.C., R&D department, Alexandrou Pantou 25, 17671, Kallithea-Athens, Greece
| | - Piotr Gawron
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, 6 Avenue du Swing, L-4367 Belvaux, Luxembourg
| | - Alexander Mazein
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, 6 Avenue du Swing, L-4367 Belvaux, Luxembourg
| | - Eric Bonnet
- Centre National de Recherche en Génomique Humaine (CNRGH), CEA, 2 rue Gaston Crémieux, CP5706 91057 EVRY-GENOPOLE cedex, Evry, France
| | - Elisabeth Petit-Teixeira
- Laboratoire Européen de Recherche pour la Polyarthrite Rhumatoïde - Genhotel, Univ Evry, Université Paris-Saclay, 2, rue Gaston Crémieux, 91057 EVRY-GENOPOLE cedex, Evry, France
| | - Anna Niarakis
- Laboratoire Européen de Recherche pour la Polyarthrite Rhumatoïde - Genhotel, Univ Evry, Université Paris-Saclay, 2, rue Gaston Crémieux, 91057 EVRY-GENOPOLE cedex, Evry, France
| |
Collapse
|
28
|
Deritei D, Rozum J, Ravasz Regan E, Albert R. A feedback loop of conditionally stable circuits drives the cell cycle from checkpoint to checkpoint. Sci Rep 2019; 9:16430. [PMID: 31712566 PMCID: PMC6848090 DOI: 10.1038/s41598-019-52725-1] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2019] [Accepted: 10/22/2019] [Indexed: 12/12/2022] Open
Abstract
We perform logic-based network analysis on a model of the mammalian cell cycle. This model is composed of a Restriction Switch driving cell cycle commitment and a Phase Switch driving mitotic entry and exit. By generalizing the concept of stable motif, i.e., a self-sustaining positive feedback loop that maintains an associated state, we introduce the concept of a conditionally stable motif, the stability of which is contingent on external conditions. We show that the stable motifs of the Phase Switch are contingent on the state of three nodes through which it receives input from the rest of the network. Biologically, these conditions correspond to cell cycle checkpoints. Holding these nodes locked (akin to a checkpoint-free cell) transforms the Phase Switch into an autonomous oscillator that robustly toggles through the cell cycle phases G1, G2 and mitosis. The conditionally stable motifs of the Phase Switch Oscillator are organized into an ordered sequence, such that they serially stabilize each other but also cause their own destabilization. Along the way they channel the dynamics of the module onto a narrow path in state space, lending robustness to the oscillation. Self-destabilizing conditionally stable motifs suggest a general negative feedback mechanism leading to sustained oscillations.
Collapse
Affiliation(s)
- Dávid Deritei
- Department of Physics, Pennsylvania State University, University Park, PA, United States of America
- Department of Network and Data Science, Central European University, Budapest, Hungary
| | - Jordan Rozum
- Department of Physics, Pennsylvania State University, University Park, PA, United States of America
| | - Erzsébet Ravasz Regan
- Biochemistry and Molecular Biology, The College of Wooster, Wooster, OH, United States of America
| | - Réka Albert
- Department of Physics, Pennsylvania State University, University Park, PA, United States of America.
| |
Collapse
|
29
|
Kondratova M, Czerwinska U, Sompairac N, Amigorena SD, Soumelis V, Barillot E, Zinovyev A, Kuperstein I. A multiscale signalling network map of innate immune response in cancer reveals cell heterogeneity signatures. Nat Commun 2019; 10:4808. [PMID: 31641119 PMCID: PMC6805895 DOI: 10.1038/s41467-019-12270-x] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2018] [Accepted: 08/02/2019] [Indexed: 12/14/2022] Open
Abstract
The lack of integrated resources depicting the complexity of the innate immune response in cancer represents a bottleneck for high-throughput data interpretation. To address this challenge, we perform a systematic manual literature mining of molecular mechanisms governing the innate immune response in cancer and represent it as a signalling network map. The cell-type specific signalling maps of macrophages, dendritic cells, myeloid-derived suppressor cells and natural killers are constructed and integrated into a comprehensive meta map of the innate immune response in cancer. The meta-map contains 1466 chemical species as nodes connected by 1084 biochemical reactions, and it is supported by information from 820 articles. The resource helps to interpret single cell RNA-Seq data from macrophages and natural killer cells in metastatic melanoma that reveal different anti- or pro-tumor sub-populations within each cell type. Here, we report a new open source analytic platform that supports data visualisation and interpretation of tumour microenvironment activity in cancer. The complexity of the innate immune response to cancer makes interpretation of large data sets challenging. Here, the authors provide an integrated multi-scale map of signalling networks representing the different immune cells and their interactions and show its utility for data interpretation.
Collapse
Affiliation(s)
- Maria Kondratova
- Institut Curie, PSL Research University, Mines Paris Tech, Inserm, U900, 75005, Paris, France
| | - Urszula Czerwinska
- Institut Curie, PSL Research University, Mines Paris Tech, Inserm, U900, 75005, Paris, France.,Université Paris Descartes, Centre de Recherches Interdisciplinaires, Paris, France
| | - Nicolas Sompairac
- Institut Curie, PSL Research University, Mines Paris Tech, Inserm, U900, 75005, Paris, France.,Université Paris Descartes, Centre de Recherches Interdisciplinaires, Paris, France
| | | | - Vassili Soumelis
- Institut Curie, PSL Research University, Inserm, U932, 75005, Paris, France
| | - Emmanuel Barillot
- Institut Curie, PSL Research University, Mines Paris Tech, Inserm, U900, 75005, Paris, France
| | - Andrei Zinovyev
- Institut Curie, PSL Research University, Mines Paris Tech, Inserm, U900, 75005, Paris, France
| | - Inna Kuperstein
- Institut Curie, PSL Research University, Mines Paris Tech, Inserm, U900, 75005, Paris, France.
| |
Collapse
|
30
|
Choo SM, Park SM, Cho KH. Minimal intervening control of biomolecular networks leading to a desired cellular state. Sci Rep 2019; 9:13124. [PMID: 31511585 PMCID: PMC6739335 DOI: 10.1038/s41598-019-49571-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2019] [Accepted: 08/27/2019] [Indexed: 02/07/2023] Open
Abstract
A cell phenotype can be represented by an attractor state of the underlying molecular regulatory network, to which other network states eventually converge. Here, the set of states converging to each attractor is called its basin of attraction. A central question is how to drive a particular cell state toward a desired attractor with minimal interventions on the network system. We develop a general control framework of complex Boolean networks to provide an answer to this question by identifying control targets on which one-time temporary perturbation can induce a state transition to the boundary of a desired attractor basin. Examples are shown to illustrate the proposed control framework which is also applicable to other types of complex Boolean networks.
Collapse
Affiliation(s)
- Sang-Mok Choo
- Department of Mathematics, University of Ulsan, Ulsan, 44610, Republic of Korea
| | - Sang-Min Park
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| | - Kwang-Hyun Cho
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea.
| |
Collapse
|
31
|
Baudin A, Paul S, Su C, Pang J. Controlling large Boolean networks with single-step perturbations. Bioinformatics 2019; 35:i558-i567. [PMID: 31510648 PMCID: PMC6612811 DOI: 10.1093/bioinformatics/btz371] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
Motivation The control of Boolean networks has traditionally focussed on strategies where the perturbations are applied to the nodes of the network for an extended period of time. In this work, we study if and how a Boolean network can be controlled by perturbing a minimal set of nodes for a single-step and letting the system evolve afterwards according to its original dynamics. More precisely, given a Boolean network (BN), we compute a minimal subset Cmin of the nodes such that BN can be driven from any initial state in an attractor to another ‘desired’ attractor by perturbing some or all of the nodes of Cmin for a single-step. Such kind of control is attractive for biological systems because they are less time consuming than the traditional strategies for control while also being financially more viable. However, due to the phenomenon of state-space explosion, computing such a minimal subset is computationally inefficient and an approach that deals with the entire network in one-go, does not scale well for large networks. Results We develop a ‘divide-and-conquer’ approach by decomposing the network into smaller partitions, computing the minimal control on the projection of the attractors to these partitions and then composing the results to obtain Cmin for the whole network. We implement our method and test it on various real-life biological networks to demonstrate its applicability and efficiency. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Alexis Baudin
- Department of Computer Science, École Normale Supérieure Paris-Saclay, Cachan, France
| | - Soumya Paul
- Faculty of Science, Technology and Communication, University of Luxembourg, Esch-sur-Alzette, Luxembourg
| | - Cui Su
- Interdisciplinary Centre for Security, Reliability and Trust, University of Luxembourg, Luxembourg, Luxembourg
| | - Jun Pang
- Faculty of Science, Technology and Communication, University of Luxembourg, Esch-sur-Alzette, Luxembourg.,Interdisciplinary Centre for Security, Reliability and Trust, University of Luxembourg, Luxembourg, Luxembourg
| |
Collapse
|
32
|
Glucocorticoid Receptor Modulates EGFR Feedback upon Acquisition of Resistance to Monoclonal Antibodies. J Clin Med 2019; 8:jcm8050600. [PMID: 31052457 PMCID: PMC6572202 DOI: 10.3390/jcm8050600] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2019] [Revised: 04/26/2019] [Accepted: 04/28/2019] [Indexed: 12/20/2022] Open
Abstract
Evidences of a crosstalk between Epidermal Growth Factor Receptor (EGFR) and Glucocorticoid Receptor (GR) has been reported, ranging from the modulation of receptor levels or GR mediated transcriptional repression of EGFR target genes, with modifications of epigenetic markers. The present study focuses on the involvement of EGFR positive and negative feedback genes in the establishment of cetuximab (CTX) resistance in metastatic Colorectal Cancer (CRC) patients. We evaluated the expression profile of the EGFR ligands TGFA and HBEGF, along with the pro-inflammatory cytokines IL-1B and IL-8, which were previously reported to be negatively associated with monoclonal antibody response, both in mice and patient specimens. Among EGFR negative feedback loops, we focused on ERRFI1, DUSP1, LRIG3, and LRIG1. We observed that EGFR positive feedback genes are increased in CTX-resistant cells, whereas negative feedback genes are reduced. Next, we tested the expression of these genes in CTX-resistant cells upon GR modulation. We unveiled that GR activation leads to an increase in ERRFI1, DUSP1, and LRIG1, which were shown to restrict EGFR activity, along with a decrease in the EGFR activators (TGFA and IL-8). Finally, in a cohort of xenopatients, stratified for response to cetuximab, we observed an inverse association between the expression level of LRIG1 and CRC progression upon CTX treatment. Our model implies that combining GR modulation to EGFR inhibition may yield an effective treatment strategy in halting cancer progression.
Collapse
|
33
|
Rodríguez-Jorge O, Kempis-Calanis LA, Abou-Jaoudé W, Gutiérrez-Reyna DY, Hernandez C, Ramirez-Pliego O, Thomas-Chollier M, Spicuglia S, Santana MA, Thieffry D. Cooperation between T cell receptor and Toll-like receptor 5 signaling for CD4 + T cell activation. Sci Signal 2019; 12:12/577/eaar3641. [PMID: 30992399 DOI: 10.1126/scisignal.aar3641] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
CD4+ T cells recognize antigens through their T cell receptors (TCRs); however, additional signals involving costimulatory receptors, for example, CD28, are required for proper T cell activation. Alternative costimulatory receptors have been proposed, including members of the Toll-like receptor (TLR) family, such as TLR5 and TLR2. To understand the molecular mechanism underlying a potential costimulatory role for TLR5, we generated detailed molecular maps and logical models for the TCR and TLR5 signaling pathways and a merged model for cross-interactions between the two pathways. Furthermore, we validated the resulting model by analyzing how T cells responded to the activation of these pathways alone or in combination, in terms of the activation of the transcriptional regulators CREB, AP-1 (c-Jun), and NF-κB (p65). Our merged model accurately predicted the experimental results, showing that the activation of TLR5 can play a similar role to that of CD28 activation with respect to AP-1, CREB, and NF-κB activation, thereby providing insights regarding the cross-regulation of these pathways in CD4+ T cells.
Collapse
Affiliation(s)
- Otoniel Rodríguez-Jorge
- Centro de Investigación en Dinámica Celular, Instituto de Investigación en Ciencias Básicas y Aplicadas, Universidad Autónoma del Estado de Morelos, 62210 Cuernavaca, México.,Escuela de Estudios Superiores de Axochiapan, Universidad Autónoma del Estado de Morelos, 62951 Axochiapan, México
| | - Linda A Kempis-Calanis
- Centro de Investigación en Dinámica Celular, Instituto de Investigación en Ciencias Básicas y Aplicadas, Universidad Autónoma del Estado de Morelos, 62210 Cuernavaca, México
| | - Wassim Abou-Jaoudé
- Computational System Biology Team, Institut de Biologie de l'Ecole Normale Supérieure, CNRS UMR8197, INSERM U1024, École Normale Supérieure, Université PSL, 75005 Paris, France
| | - Darely Y Gutiérrez-Reyna
- Centro de Investigación en Dinámica Celular, Instituto de Investigación en Ciencias Básicas y Aplicadas, Universidad Autónoma del Estado de Morelos, 62210 Cuernavaca, México
| | - Céline Hernandez
- Computational System Biology Team, Institut de Biologie de l'Ecole Normale Supérieure, CNRS UMR8197, INSERM U1024, École Normale Supérieure, Université PSL, 75005 Paris, France
| | - Oscar Ramirez-Pliego
- Centro de Investigación en Dinámica Celular, Instituto de Investigación en Ciencias Básicas y Aplicadas, Universidad Autónoma del Estado de Morelos, 62210 Cuernavaca, México
| | - Morgane Thomas-Chollier
- Computational System Biology Team, Institut de Biologie de l'Ecole Normale Supérieure, CNRS UMR8197, INSERM U1024, École Normale Supérieure, Université PSL, 75005 Paris, France
| | | | - Maria A Santana
- Centro de Investigación en Dinámica Celular, Instituto de Investigación en Ciencias Básicas y Aplicadas, Universidad Autónoma del Estado de Morelos, 62210 Cuernavaca, México.
| | - Denis Thieffry
- Computational System Biology Team, Institut de Biologie de l'Ecole Normale Supérieure, CNRS UMR8197, INSERM U1024, École Normale Supérieure, Université PSL, 75005 Paris, France.
| |
Collapse
|
34
|
Ostaszewski M, Gebel S, Kuperstein I, Mazein A, Zinovyev A, Dogrusoz U, Hasenauer J, Fleming RMT, Le Novère N, Gawron P, Ligon T, Niarakis A, Nickerson D, Weindl D, Balling R, Barillot E, Auffray C, Schneider R. Community-driven roadmap for integrated disease maps. Brief Bioinform 2019; 20:659-670. [PMID: 29688273 PMCID: PMC6556900 DOI: 10.1093/bib/bby024] [Citation(s) in RCA: 40] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2018] [Revised: 03/02/2018] [Indexed: 01/07/2023] Open
Abstract
The Disease Maps Project builds on a network of scientific and clinical groups that exchange best practices, share information and develop systems biomedicine tools. The project aims for an integrated, highly curated and user-friendly platform for disease-related knowledge. The primary focus of disease maps is on interconnected signaling, metabolic and gene regulatory network pathways represented in standard formats. The involvement of domain experts ensures that the key disease hallmarks are covered and relevant, up-to-date knowledge is adequately represented. Expert-curated and computer readable, disease maps may serve as a compendium of knowledge, allow for data-supported hypothesis generation or serve as a scaffold for the generation of predictive mathematical models. This article summarizes the 2nd Disease Maps Community meeting, highlighting its important topics and outcomes. We outline milestones on the roadmap for the future development of disease maps, including creating and maintaining standardized disease maps; sharing parts of maps that encode common human disease mechanisms; providing technical solutions for complexity management of maps; and Web tools for in-depth exploration of such maps. A dedicated discussion was focused on mathematical modeling approaches, as one of the main goals of disease map development is the generation of mathematically interpretable representations to predict disease comorbidity or drug response and to suggest drug repositioning, altogether supporting clinical decisions.
Collapse
Affiliation(s)
- Marek Ostaszewski
- Luxembourg Centre for Systems Biomedicine, Universite du Luxembourg, 7 Avenue des Hauts-Fourneaux, L-4362 Esch-sur-Alzette, Luxembourg
| | - Stephan Gebel
- Luxembourg Centre for Systems Biomedicine, Universite du Luxembourg, 7 Avenue des Hauts-Fourneaux, L-4362 Esch-sur-Alzette, Luxembourg
| | - Inna Kuperstein
- Institut Curie, PSL Research University, INSERM U900, F-75005 Paris, France and CBIO-Centre for Computational Biology, MINES ParisTech, PSL Research University, F-75006 Paris, France
| | - Alexander Mazein
- European Institute for Systems Biology and Medicine, CIRI UMR5308, CNRS-ENS-UCBL-INSERM, Université de Lyon, 50 Avenue Tony Garnier, 69007 Lyon, France
| | - Andrei Zinovyev
- Institut Curie, PSL Research University, INSERM U900, F-75005 Paris, France and CBIO-Centre for Computational Biology, MINES ParisTech, PSL Research University, F-75006 Paris, France
| | - Ugur Dogrusoz
- Computer Engineering Department, Faculty of Engineering, Bilkent University, Ankara 06800, Turkey
| | - Jan Hasenauer
- Institute of Computational Biology, Helmholtz Zentrum München-German Research Center for Environmental Health, Ingolstädter Landstr. 1, 85764 Neuherberg, Germany
| | - Ronan M T Fleming
- Division of Systems Biomedicine and Pharmacology, Leiden Academic Centre for Drug Research, Faculty of Science, Leiden University, Leiden, Netherlands
| | - Nicolas Le Novère
- The Babraham Institute, Babraham Research Campus, Cambridge CB22 3AT, United Kingdom
| | - Piotr Gawron
- Luxembourg Centre for Systems Biomedicine, Universite du Luxembourg, 7 Avenue des Hauts-Fourneaux, L-4362 Esch-sur-Alzette, Luxembourg
| | - Thomas Ligon
- Faculty of Physics and Center for NanoScience (CeNS), Ludwig-Maximilians-Universität, 80539 München, Germany
| | - Anna Niarakis
- GenHotel EA3886, Univ Evry, Université Paris-Saclay, Evry 91025, France
| | - David Nickerson
- Auckland Bioengineering Institute, University of Auckland, Auckland, New Zealand
| | - Daniel Weindl
- Institute of Computational Biology, Helmholtz Zentrum München-German Research Center for Environmental Health, Ingolstädter Landstr. 1, 85764 Neuherberg, Germany
| | - Rudi Balling
- Luxembourg Centre for Systems Biomedicine, Universite du Luxembourg, 7 Avenue des Hauts-Fourneaux, L-4362 Esch-sur-Alzette, Luxembourg
| | - Emmanuel Barillot
- Institut Curie, PSL Research University, INSERM U900, F-75005 Paris, France and CBIO-Centre for Computational Biology, MINES ParisTech, PSL Research University, F-75006 Paris, France
| | - Charles Auffray
- European Institute for Systems Biology and Medicine, CIRI UMR5308, CNRS-ENS-UCBL-INSERM, Université de Lyon, 50 Avenue Tony Garnier, 69007 Lyon, France
| | - Reinhard Schneider
- Luxembourg Centre for Systems Biomedicine, Universite du Luxembourg, 7 Avenue des Hauts-Fourneaux, L-4362 Esch-sur-Alzette, Luxembourg
| |
Collapse
|
35
|
Sizek H, Hamel A, Deritei D, Campbell S, Ravasz Regan E. Boolean model of growth signaling, cell cycle and apoptosis predicts the molecular mechanism of aberrant cell cycle progression driven by hyperactive PI3K. PLoS Comput Biol 2019; 15:e1006402. [PMID: 30875364 PMCID: PMC6436762 DOI: 10.1371/journal.pcbi.1006402] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2018] [Revised: 03/27/2019] [Accepted: 02/12/2019] [Indexed: 02/07/2023] Open
Abstract
The PI3K/AKT signaling pathway plays a role in most cellular functions linked to cancer progression, including cell growth, proliferation, cell survival, tissue invasion and angiogenesis. It is generally recognized that hyperactive PI3K/AKT1 are oncogenic due to their boost to cell survival, cell cycle entry and growth-promoting metabolism. That said, the dynamics of PI3K and AKT1 during cell cycle progression are highly nonlinear. In addition to negative feedback that curtails their activity, protein expression of PI3K subunits has been shown to oscillate in dividing cells. The low-PI3K/low-AKT1 phase of these oscillations is required for cytokinesis, indicating that oncogenic PI3K may directly contribute to genome duplication. To explore this, we construct a Boolean model of growth factor signaling that can reproduce PI3K oscillations and link them to cell cycle progression and apoptosis. The resulting modular model reproduces hyperactive PI3K-driven cytokinesis failure and genome duplication and predicts the molecular drivers responsible for these failures by linking hyperactive PI3K to mis-regulation of Polo-like kinase 1 (Plk1) expression late in G2. To do this, our model captures the role of Plk1 in cell cycle progression and accurately reproduces multiple effects of its loss: G2 arrest, mitotic catastrophe, chromosome mis-segregation / aneuploidy due to premature anaphase, and cytokinesis failure leading to genome duplication, depending on the timing of Plk1 inhibition along the cell cycle. Finally, we offer testable predictions on the molecular drivers of PI3K oscillations, the timing of these oscillations with respect to division, and the role of altered Plk1 and FoxO activity in genome-level defects caused by hyperactive PI3K. Our model is an important starting point for the predictive modeling of cell fate decisions that include AKT1-driven senescence, as well as the non-intuitive effects of drugs that interfere with mitosis.
Collapse
Affiliation(s)
- Herbert Sizek
- Biochemistry and Molecular Biology, The College of Wooster, Wooster, OH, United States of America
| | - Andrew Hamel
- Biochemistry and Molecular Biology, The College of Wooster, Wooster, OH, United States of America
| | - Dávid Deritei
- Department of Physics, Pennsylvania State University, State College, PA, United States of America
- Department of Network and Data Science, Central European University, Budapest, Hungary
| | - Sarah Campbell
- Biochemistry and Molecular Biology, The College of Wooster, Wooster, OH, United States of America
| | - Erzsébet Ravasz Regan
- Biochemistry and Molecular Biology, The College of Wooster, Wooster, OH, United States of America
| |
Collapse
|
36
|
Voukantsis D, Kahn K, Hadley M, Wilson R, Buffa FM. Modeling genotypes in their microenvironment to predict single- and multi-cellular behavior. Gigascience 2019; 8:giz010. [PMID: 30715320 PMCID: PMC6423375 DOI: 10.1093/gigascience/giz010] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2018] [Revised: 11/13/2018] [Accepted: 01/16/2019] [Indexed: 12/21/2022] Open
Abstract
A cell's phenotype is the set of observable characteristics resulting from the interaction of the genotype with the surrounding environment, determining cell behavior. Deciphering genotype-phenotype relationships has been crucial to understanding normal and disease biology. Analysis of molecular pathways has provided an invaluable tool to such understanding; however, typically it does not consider the physical microenvironment, which is a key determinant of phenotype. In this study, we present a novel modeling framework that enables the study of the link between genotype, signaling networks, and cell behavior in a three-dimensional microenvironment. To achieve this, we bring together Agent-Based Modeling, a powerful computational modeling technique, and gene networks. This combination allows biological hypotheses to be tested in a controlled stepwise fashion, and it lends itself naturally to model a heterogeneous population of cells acting and evolving in a dynamic microenvironment, which is needed to predict the evolution of complex multi-cellular dynamics. Importantly, this enables modeling co-occurring intrinsic perturbations, such as mutations, and extrinsic perturbations, such as nutrient availability, and their interactions. Using cancer as a model system, we illustrate how this framework delivers a unique opportunity to identify determinants of single-cell behavior, while uncovering emerging properties of multi-cellular growth. This framework is freely available at http://www.microc.org.
Collapse
Affiliation(s)
- Dimitrios Voukantsis
- Computational Biology and Integrative Genomics, MRC/CRUK Oxford Institute, Departmemt of Oncology, University of Oxford, Old Road Campus, Oxford, Oxfordshire, OX3 7DQ, UK
| | - Kenneth Kahn
- Computational Biology and Integrative Genomics, MRC/CRUK Oxford Institute, Departmemt of Oncology, University of Oxford, Old Road Campus, Oxford, Oxfordshire, OX3 7DQ, UK
- Academic Information Technology Research Team, University of Oxford, 13 Bambury Road, Oxford, Oxfordshire, OX2 6NN, UK
| | - Martin Hadley
- Academic Information Technology Research Team, University of Oxford, 13 Bambury Road, Oxford, Oxfordshire, OX2 6NN, UK
| | - Rowan Wilson
- Academic Information Technology Research Team, University of Oxford, 13 Bambury Road, Oxford, Oxfordshire, OX2 6NN, UK
| | - Francesca M Buffa
- Computational Biology and Integrative Genomics, MRC/CRUK Oxford Institute, Departmemt of Oncology, University of Oxford, Old Road Campus, Oxford, Oxfordshire, OX3 7DQ, UK
| |
Collapse
|
37
|
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: 42] [Impact Index Per Article: 8.4] [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
|
38
|
Mizera A, Pang J, Qu H, Yuan Q. Taming Asynchrony for Attractor Detection in Large Boolean Networks. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2019; 16:31-42. [PMID: 29994682 DOI: 10.1109/tcbb.2018.2850901] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Boolean networks is a well-established formalism for modelling biological systems. A vital challenge for analyzing a Boolean network is to identify all the attractors. This becomes more challenging for large asynchronous Boolean networks, due to the asynchronous scheme. Existing methods are prohibited due to the well-known state-space explosion problem in large Boolean networks. In this paper, we tackle this challenge by proposing a SCC-based decomposition method. We prove the correctness of our proposed method and demonstrate its efficiency with two real-life biological networks.
Collapse
|
39
|
Rynkeviciene R, Simiene J, Strainiene E, Stankevicius V, Usinskiene J, Miseikyte Kaubriene E, Meskinyte I, Cicenas J, Suziedelis K. Non-Coding RNAs in Glioma. Cancers (Basel) 2018; 11:cancers11010017. [PMID: 30583549 PMCID: PMC6356972 DOI: 10.3390/cancers11010017] [Citation(s) in RCA: 90] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2018] [Revised: 12/18/2018] [Accepted: 12/19/2018] [Indexed: 12/12/2022] Open
Abstract
Glioma is the most aggressive brain tumor of the central nervous system. The ability of glioma cells to migrate, rapidly diffuse and invade normal adjacent tissue, their sustained proliferation, and heterogeneity contribute to an overall survival of approximately 15 months for most patients with high grade glioma. Numerous studies indicate that non-coding RNA species have critical functions across biological processes that regulate glioma initiation and progression. Recently, new data emerged, which shows that the cross-regulation between long non-coding RNAs and small non-coding RNAs contribute to phenotypic diversity of glioblastoma subclasses. In this paper, we review data of long non-coding RNA expression, which was evaluated in human glioma tissue samples during a five-year period. Thus, this review summarizes the following: (I) the role of non-coding RNAs in glioblastoma pathogenesis, (II) the potential application of non-coding RNA species in glioma-grading, (III) crosstalk between lncRNAs and miRNAs (IV) future perspectives of non-coding RNAs as biomarkers for glioma.
Collapse
Affiliation(s)
- Ryte Rynkeviciene
- Nacional Cancer Institute, Santariskiu str. 1, LT-08660 Vilnius, Lithuania.
| | - Julija Simiene
- Nacional Cancer Institute, Santariskiu str. 1, LT-08660 Vilnius, Lithuania.
- Institute of Biosciences, Life Sciences Center, Vilnius University, Sauletekio ave. 7, LT-08412 Vilnius, Lithuania.
| | - Egle Strainiene
- Nacional Cancer Institute, Santariskiu str. 1, LT-08660 Vilnius, Lithuania.
- Department of Chemistry and Bioengineering, Vilnius Gediminas Technical University, Sauletekio ave. 11, LT-10122 Vilnius, Lithuania.
| | - Vaidotas Stankevicius
- Nacional Cancer Institute, Santariskiu str. 1, LT-08660 Vilnius, Lithuania.
- Institute of Biotechnology, Vilnius University, LT-10257 Vilnius, Lithuania.
| | - Jurgita Usinskiene
- Nacional Cancer Institute, Santariskiu str. 1, LT-08660 Vilnius, Lithuania.
| | - Edita Miseikyte Kaubriene
- Nacional Cancer Institute, Santariskiu str. 1, LT-08660 Vilnius, Lithuania.
- Faculty of Medicine, Vilnius University, M.K. Cˇiurlionio 21, LT-03101 Vilnius, Lithuania.
| | - Ingrida Meskinyte
- Proteomics Center, Institute of Biochemistry, Vilnius University Life Sciences Center, Sauletekio al. 7, LT-10257 Vilnius, Lithuania.
- MAP Kinase Resource, Bioinformatics, Melchiorstrasse 9, 3027 Bern, Switzerland.
| | - Jonas Cicenas
- Proteomics Center, Institute of Biochemistry, Vilnius University Life Sciences Center, Sauletekio al. 7, LT-10257 Vilnius, Lithuania.
- MAP Kinase Resource, Bioinformatics, Melchiorstrasse 9, 3027 Bern, Switzerland.
- Energy and Biotechnology Engineering Institute, Aleksandro Stulginskio University, Studentų g. 11, LT-53361 Akademija, Lithuania.
| | - Kestutis Suziedelis
- Nacional Cancer Institute, Santariskiu str. 1, LT-08660 Vilnius, Lithuania.
- Institute of Biosciences, Life Sciences Center, Vilnius University, Sauletekio ave. 7, LT-08412 Vilnius, Lithuania.
| |
Collapse
|
40
|
Thobe K, Kuznia C, Sers C, Siebert H. Evaluating Uncertainty in Signaling Networks Using Logical Modeling. Front Physiol 2018; 9:1335. [PMID: 30364151 PMCID: PMC6191669 DOI: 10.3389/fphys.2018.01335] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2018] [Accepted: 09/04/2018] [Indexed: 01/01/2023] Open
Abstract
Systems biology studies the structure and dynamics of biological systems using mathematical approaches. Bottom-up approaches create models from prior knowledge but usually cannot cope with uncertainty, whereas top-down approaches infer models directly from data using statistical methods but mostly neglect valuable known information from former studies. Here, we want to present a workflow that includes prior knowledge while allowing for uncertainty in the modeling process. We build not one but all possible models that arise from the uncertainty using logical modeling and subsequently filter for those models in agreement with data in a top-down manner. This approach enables us to investigate new and more complex biological research questions, however, the encoding in such a framework is often not obvious and thus not easily accessible for researcher from life sciences. To mitigate this problem, we formulate a pipeline with specific templates to address some research questions common in signaling network analysis. To illustrate the potential of this approach, we applied the pipeline to growth factor signaling processes in two renal cancer cell lines. These two cell lines originate from similar tissue, but surprisingly showed a very different behavior toward the cancer drug Sorafenib. Thus our aim was to explore differences between these cell lines regarding three sources of uncertainty in one analysis: possible targets of Sorafenib, crosstalk between involved pathways, and the effect of a mutation in mammalian target of Rapamycin (mTOR) in one of the cell lines. We were able to show that the model pools from the cell lines are disjoint, thus the discrepancies in behavior originate from differences in the cellular wiring. Also the mutation in mTOR is not affecting its activity in the pathway. The results on Sorafenib, while not fully clarifying the mechanisms involved, illustrate the potential of this analysis for generating new hypotheses.
Collapse
Affiliation(s)
- Kirsten Thobe
- Group for Discrete Biomathematics, Department for Mathematics and Computer Science, Freie Universität Berlin, Berlin, Germany.,Group for Mathematical Modelling of Cellular Processes, Max-Delbrück Center for Molecular Medicine, Berlin, Germany
| | - Christina Kuznia
- Laboratory of Molecular Tumor Pathology, Institute of Pathology, Charité Universitätsmedizin Berlin, Berlin, Germany.,Laboratory of Bioorganic Synthesis, Department of Chemistry, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Christine Sers
- Laboratory of Molecular Tumor Pathology, Institute of Pathology, Charité Universitätsmedizin Berlin, Berlin, Germany
| | - Heike Siebert
- Group for Discrete Biomathematics, Department for Mathematics and Computer Science, Freie Universität Berlin, Berlin, Germany
| |
Collapse
|
41
|
Chen L, Kulasiri D, Samarasinghe S. A Novel Data-Driven Boolean Model for Genetic Regulatory Networks. Front Physiol 2018; 9:1328. [PMID: 30319440 PMCID: PMC6167558 DOI: 10.3389/fphys.2018.01328] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2017] [Accepted: 09/03/2018] [Indexed: 11/30/2022] Open
Abstract
A Boolean model is a simple, discrete and dynamic model without the need to consider the effects at the intermediate levels. However, little effort has been made into constructing activation, inhibition, and protein decay networks, which could indicate the direct roles of a gene (or its synthesized protein) as an activator or inhibitor of a target gene. Therefore, we propose to focus on the general Boolean functions at the subfunction level taking into account the effectiveness of protein decay, and further split the subfunctions into the activation and inhibition domains. As a consequence, we developed a novel data-driven Boolean model; namely, the Fundamental Boolean Model (FBM), to draw insights into gene activation, inhibition, and protein decay. This novel Boolean model provides an intuitive definition of activation and inhibition pathways and includes mechanisms to handle protein decay issues. To prove the concept of the novel model, we implemented a platform using R language, called FBNNet. Our experimental results show that the proposed FBM could explicitly display the internal connections of the mammalian cell cycle between genes separated into the connection types of activation, inhibition and protein decay. Moreover, the method we proposed to infer the gene regulatory networks for the novel Boolean model can be run in parallel and; hence, the computation cost is affordable. Finally, the novel Boolean model and related Fundamental Boolean Networks (FBNs) could show significant trajectories in genes to reveal how genes regulated each other over a given period. This new feature could facilitate further research on drug interventions to detect the side effects of a newly-proposed drug.
Collapse
Affiliation(s)
- Leshi Chen
- Computational Systems Biology Laboratory, Centre for Advanced Computational Solutions, Lincoln University, Lincoln, New Zealand
| | - Don Kulasiri
- Computational Systems Biology Laboratory, Centre for Advanced Computational Solutions, Lincoln University, Lincoln, New Zealand
| | - Sandhya Samarasinghe
- Integrated Systems Modelling Group, Centre for Advanced Computational Solutions, Lincoln University, Lincoln, New Zealand
| |
Collapse
|
42
|
Determining Relative Dynamic Stability of Cell States Using Boolean Network Model. Sci Rep 2018; 8:12077. [PMID: 30104572 PMCID: PMC6089891 DOI: 10.1038/s41598-018-30544-0] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2018] [Accepted: 08/02/2018] [Indexed: 01/05/2023] Open
Abstract
Cell state transition is at the core of biological processes in metazoan, which includes cell differentiation, epithelial-to-mesenchymal transition (EMT) and cell reprogramming. In these cases, it is important to understand the molecular mechanism of cellular stability and how the transitions happen between different cell states, which is controlled by a gene regulatory network (GRN) hard-wired in the genome. Here we use Boolean modeling of GRN to study the cell state transition of EMT and systematically compare four available methods to calculate the cellular stability of three cell states in EMT in both normal and genetically mutated cases. The results produced from four methods generally agree but do not totally agree with each other. We show that distribution of one-degree neighborhood of cell states, which are the nearest states by Hamming distance, causes the difference among the methods. From that, we propose a new method based on one-degree neighborhood, which is the simplest one and agrees with other methods to estimate the cellular stability in all scenarios of our EMT model. This new method will help the researchers in the field of cell differentiation and cell reprogramming to calculate cellular stability using Boolean model, and then rationally design their experimental protocols to manipulate the cell state transition.
Collapse
|
43
|
Yang JM, Lee CK, Cho KH. Global Stabilization of Boolean Networks to Control the Heterogeneity of Cellular Responses. Front Physiol 2018; 9:774. [PMID: 30072906 PMCID: PMC6060448 DOI: 10.3389/fphys.2018.00774] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2018] [Accepted: 06/04/2018] [Indexed: 12/13/2022] Open
Abstract
Boolean networks (BNs) have been widely used as a useful model for molecular regulatory networks in systems biology. In the state space of BNs, attractors represent particular cell phenotypes. For targeted therapy of cancer, there is a pressing need to control the heterogeneity of cellular responses to the targeted drug by reducing the number of attractors associated with the ill phenotypes of cancer cells. Here, we present a novel control scheme for global stabilization of BNs to a unique fixed point. Using a sufficient condition of global stabilization with respect to the adjacency matrix, we can determine a set of constant controls so that the controlled BN is steered toward an unspecified fixed point which can then be further transformed to a desired attractor by subsequent control. Our method is efficient in that it has polynomial complexity with respect to the number of state variables, while having exponential complexity with respect to in-degree of BNs. To demonstrate the applicability of the proposed control scheme, we conduct simulation studies using a regulation influence network describing the metastatic process of cells and the Mitogen-activated protein kinase (MAPK) signaling network that is crucial in cancer cell fate determination.
Collapse
Affiliation(s)
- Jung-Min Yang
- School of Electronics Engineering, Kyungpook National University, Daegu, South Korea
| | - Chun-Kyung Lee
- School of Electronics Engineering, Kyungpook National University, Daegu, South Korea
| | - Kwang-Hyun Cho
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology, Daejeon, South Korea
| |
Collapse
|
44
|
Fages F, Martinez T, Rosenblueth DA, Soliman S. Influence Networks Compared with Reaction Networks: Semantics, Expressivity and Attractors. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2018; 15:1138-1151. [PMID: 29994637 DOI: 10.1109/tcbb.2018.2805686] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Biochemical reaction networks are one of the most widely used formalisms in systems biology to describe the molecular mechanisms of high-level cell processes. However, modellers also reason with influence diagrams to represent the positive and negative influences between molecular species and may find an influence network useful in the process of building a reaction network. In this paper, we introduce a formalism of influence networks with forces, and equip it with a hierarchy of Boolean, Petri net, stochastic and differential semantics, similarly to reaction networks with rates. We show that the expressive power of influence networks is the same as that of reaction networks under the differential semantics, but weaker under the discrete semantics. Furthermore, the hierarchy of semantics leads us to consider a (positive) Boolean semantics that cannot test the absence of a species, that we compare with the (negative) Boolean semantics with test for absence of a species in gene regulatory networks à la Thomas. We study the monotonicity properties of the positive semantics and derive from them an algorithm to compute attractors in both the positive and negative Boolean semantics. We illustrate our results on models of the literature about the p53/Mdm2 DNA damage repair system, the circadian clock, and the influence of MAPK signaling on cell-fate decision in urinary bladder cancer.
Collapse
|
45
|
Pauleve L. Reduction of Qualitative Models of Biological Networks for Transient Dynamics Analysis. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2018; 15:1167-1179. [PMID: 28885158 DOI: 10.1109/tcbb.2017.2749225] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Qualitative models of dynamics of signalling pathways and gene regulatory networks allow for the capturing of temporal properties of biological networks while requiring few parameters. However, these discrete models typically suffer from the so-called state space explosion problem which makes the formal assessment of their potential behaviors very challenging. In this paper, we describe a method to reduce a qualitative model for enhancing the tractability of analysis of transient reachability properties. The reduction does not change the dimension of the model, but instead limits its degree of freedom, therefore reducing the set of states and transitions to consider. We rely on a transition-centered specification of qualitative models by the mean of automata networks. Our framework encompasses the usual asynchronous Boolean and multi-valued network, as well as 1-bounded Petri nets. Applied to different large-scale biological networks from the litterature, we show that the reduction can lead to a drastic improvement for the scalability of verification methods.
Collapse
|
46
|
Naldi A, Hernandez C, Abou-Jaoudé W, Monteiro PT, Chaouiya C, Thieffry D. Logical Modeling and Analysis of Cellular Regulatory Networks With GINsim 3.0. Front Physiol 2018; 9:646. [PMID: 29971008 PMCID: PMC6018412 DOI: 10.3389/fphys.2018.00646] [Citation(s) in RCA: 46] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2018] [Accepted: 05/11/2018] [Indexed: 11/13/2022] Open
Abstract
The logical formalism is well adapted to model large cellular networks, in particular when detailed kinetic data are scarce. This tutorial focuses on this well-established qualitative framework. Relying on GINsim (release 3.0), a software implementing this formalism, we guide the reader step by step toward the definition, the analysis and the simulation of a four-node model of the mammalian p53-Mdm2 network.
Collapse
Affiliation(s)
- Aurélien Naldi
- Computational Systems Biology Team, Institut de Biologie de l'Ecole Normale Supérieure (IBENS), École Normale Supérieure, Centre National de la Recherche Scientifique, Institut National de la Sante et de la Recherche Médicale, PSL Université, Paris, France
| | - Céline Hernandez
- Computational Systems Biology Team, Institut de Biologie de l'Ecole Normale Supérieure (IBENS), École Normale Supérieure, Centre National de la Recherche Scientifique, Institut National de la Sante et de la Recherche Médicale, PSL Université, Paris, France
| | - Wassim Abou-Jaoudé
- Computational Systems Biology Team, Institut de Biologie de l'Ecole Normale Supérieure (IBENS), École Normale Supérieure, Centre National de la Recherche Scientifique, Institut National de la Sante et de la Recherche Médicale, PSL Université, Paris, France
| | - Pedro T. Monteiro
- INESC-ID, Instituto Superior Técnico, University of Lisbon, Lisbon, Portugal
| | | | - Denis Thieffry
- Computational Systems Biology Team, Institut de Biologie de l'Ecole Normale Supérieure (IBENS), École Normale Supérieure, Centre National de la Recherche Scientifique, Institut National de la Sante et de la Recherche Médicale, PSL Université, Paris, France
| |
Collapse
|
47
|
P_UNSAT approach of attractor calculation for Boolean gene regulatory networks. J Theor Biol 2018; 447:171-177. [PMID: 29605228 DOI: 10.1016/j.jtbi.2018.03.037] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2017] [Revised: 03/27/2018] [Accepted: 03/28/2018] [Indexed: 11/21/2022]
Abstract
Boolean network models provide an efficient way for studying gene regulatory networks. The main dynamics of a Boolean network is determined by its attractors. Attractor calculation plays a key role for analyzing Boolean gene regulatory networks. An approach of attractor calculation was proposed in this study, which combined the predecessor approach and the logic unsatisfiability approach to accelerate attractor calculation. The proposed algorithm is effective to calculate all attractors for large-scale Boolean gene regulatory networks even the networks with a relatively large average degree.
Collapse
|
48
|
G T Zañudo J, Steinway SN, Albert R. Discrete dynamic network modeling of oncogenic signaling: Mechanistic insights for personalized treatment of cancer. ACTA ACUST UNITED AC 2018; 9:1-10. [PMID: 32954058 PMCID: PMC7487767 DOI: 10.1016/j.coisb.2018.02.002] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Targeted drugs disrupting proteins that are dysregulated in cancer have emerged as promising treatments because of their specificity to cancer cell aberrations and thus their improved side effect profile. However, their success remains limited, largely due to existing or emergent therapy resistance. We suggest that this is due to limited understanding of the entire relevant cellular landscape. A class of mathematical models called discrete dynamic network models can be used to understand the integrated effect of an individual tumor's aberrations. We review the recent literature on discrete dynamic models of cancer and highlight their predicted therapeutic strategies. We believe dynamic network modeling can be used to drive treatment decision-making in a personalized manner to direct improved treatments in cancer. Cancer is rooted in incorrect cellular decisions caused by genetic alterations. Dynamic models of signaling networks can map the relevant repertoire of alterations. Discrete dynamic network models can predict therapeutic interventions. Progress in personalized medicine needs integration of multiple data and model types.
Collapse
Affiliation(s)
- Jorge G T Zañudo
- Department of Physics, The Pennsylvania State University, University Park, PA 16802, USA.,Department of Medical Oncology, Dana-Farber Cancer Institute and Broad Institute of Harvard and MIT, Boston MA, USA
| | - Steven N Steinway
- Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, MD 21287, USA
| | - Réka Albert
- Department of Physics, The Pennsylvania State University, University Park, PA 16802, USA.,Department of Biology, The Pennsylvania State University, University Park, PA 16802, USA
| |
Collapse
|
49
|
Yang G, Gómez Tejeda Zañudo J, Albert R. Target Control in Logical Models Using the Domain of Influence of Nodes. Front Physiol 2018; 9:454. [PMID: 29867523 PMCID: PMC5951947 DOI: 10.3389/fphys.2018.00454] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2018] [Accepted: 04/13/2018] [Indexed: 11/13/2022] Open
Abstract
Dynamical models of biomolecular networks are successfully used to understand the mechanisms underlying complex diseases and to design therapeutic strategies. Network control and its special case of target control, is a promising avenue toward developing disease therapies. In target control it is assumed that a small subset of nodes is most relevant to the system's state and the goal is to drive the target nodes into their desired states. An example of target control would be driving a cell to commit to apoptosis (programmed cell death). From the experimental perspective, gene knockout, pharmacological inhibition of proteins, and providing sustained external signals are among practical intervention techniques. We identify methodologies to use the stabilizing effect of sustained interventions for target control in Boolean network models of biomolecular networks. Specifically, we define the domain of influence (DOI) of a node (in a certain state) to be the nodes (and their corresponding states) that will be ultimately stabilized by the sustained state of this node regardless of the initial state of the system. We also define the related concept of the logical domain of influence (LDOI) of a node, and develop an algorithm for its identification using an auxiliary network that incorporates the regulatory logic. This way a solution to the target control problem is a set of nodes whose DOI can cover the desired target node states. We perform greedy randomized adaptive search in node state space to find such solutions. We apply our strategy to in silico biological network models of real systems to demonstrate its effectiveness.
Collapse
Affiliation(s)
- Gang Yang
- Department of Physics, Pennsylvania State University, University Park, PA, United States
| | - Jorge Gómez Tejeda Zañudo
- Department of Physics, Pennsylvania State University, University Park, PA, United States.,Department of Medical Oncology, Dana-Farber Cancer Institute, Boston, MA, United States.,Eli and Edythe L. Broad Institute of MIT and Harvard, Cambridge, MA, United States
| | - Réka Albert
- Department of Physics, Pennsylvania State University, University Park, PA, United States.,Department of Biology, Pennsylvania State University, University Park, PA, United States
| |
Collapse
|
50
|
Monraz Gomez LC, Kondratova M, Ravel JM, Barillot E, Zinovyev A, Kuperstein I. Application of Atlas of Cancer Signalling Network in preclinical studies. Brief Bioinform 2018; 20:701-716. [DOI: 10.1093/bib/bby031] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2017] [Revised: 03/28/2018] [Indexed: 12/27/2022] Open
Affiliation(s)
- L Cristobal Monraz Gomez
- Institut Curie, PSL Research University, F-75005 Paris, France, INSERM, U900, F-75005 Paris, France and MINES ParisTech, PSL Research University, CBIO-Centre for Computational Biology, F-75006 Paris, France
| | - Maria Kondratova
- Institut Curie, PSL Research University, F-75005 Paris, France, INSERM, U900, F-75005 Paris, France and MINES ParisTech, PSL Research University, CBIO-Centre for Computational Biology, F-75006 Paris, France
| | - Jean-Marie Ravel
- Genetic Laboratory, Nancy's Regional University Hospital, Vandœuvre-lès-Nancy and INSERM UMR 954, Lorraine University, Vandœuvre-lès-Nancy
| | - Emmanuel Barillot
- Institut Curie, PSL Research University, F-75005 Paris, France, INSERM, U900, F-75005 Paris, France and MINES ParisTech, PSL Research University, CBIO-Centre for Computational Biology, F-75006 Paris, France
| | - Andrei Zinovyev
- Institut Curie, PSL Research University, F-75005 Paris, France, INSERM, U900, F-75005 Paris, France and MINES ParisTech, PSL Research University, CBIO-Centre for Computational Biology, F-75006 Paris, France
| | - Inna Kuperstein
- Institut Curie, PSL Research University, F-75005 Paris, France, INSERM, U900, F-75005 Paris, France and MINES ParisTech, PSL Research University, CBIO-Centre for Computational Biology, F-75006 Paris, France
| |
Collapse
|