1
|
Marku M, Pancaldi V. From time-series transcriptomics to gene regulatory networks: A review on inference methods. PLoS Comput Biol 2023; 19:e1011254. [PMID: 37561790 PMCID: PMC10414591 DOI: 10.1371/journal.pcbi.1011254] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/12/2023] Open
Abstract
Inference of gene regulatory networks has been an active area of research for around 20 years, leading to the development of sophisticated inference algorithms based on a variety of assumptions and approaches. With the ever increasing demand for more accurate and powerful models, the inference problem remains of broad scientific interest. The abstract representation of biological systems through gene regulatory networks represents a powerful method to study such systems, encoding different amounts and types of information. In this review, we summarize the different types of inference algorithms specifically based on time-series transcriptomics, giving an overview of the main applications of gene regulatory networks in computational biology. This review is intended to give an updated reference of regulatory networks inference tools to biologists and researchers new to the topic and guide them in selecting the appropriate inference method that best fits their questions, aims, and experimental data.
Collapse
Affiliation(s)
- Malvina Marku
- CRCT, Université de Toulouse, Inserm, CNRS, Université Toulouse III-Paul Sabatier, Centre de Recherches en Cancérologie de Toulouse, Toulouse, France
| | - Vera Pancaldi
- CRCT, Université de Toulouse, Inserm, CNRS, Université Toulouse III-Paul Sabatier, Centre de Recherches en Cancérologie de Toulouse, Toulouse, France
- Barcelona Supercomputing Center, Barcelona, Spain
| |
Collapse
|
2
|
Marazzi L, Shah M, Balakrishnan S, Patil A, Vera-Licona P. NETISCE: a network-based tool for cell fate reprogramming. NPJ Syst Biol Appl 2022; 8:21. [PMID: 35725577 PMCID: PMC9209484 DOI: 10.1038/s41540-022-00231-y] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2021] [Accepted: 05/31/2022] [Indexed: 11/17/2022] Open
Abstract
The search for effective therapeutic targets in fields like regenerative medicine and cancer research has generated interest in cell fate reprogramming. This cellular reprogramming paradigm can drive cells to a desired target state from any initial state. However, methods for identifying reprogramming targets remain limited for biological systems that lack large sets of experimental data or a dynamical characterization. We present NETISCE, a novel computational tool for identifying cell fate reprogramming targets in static networks. In combination with machine learning algorithms, NETISCE estimates the attractor landscape and predicts reprogramming targets using signal flow analysis and feedback vertex set control, respectively. Through validations in studies of cell fate reprogramming from developmental, stem cell, and cancer biology, we show that NETISCE can predict previously identified cell fate reprogramming targets and identify potentially novel combinations of targets. NETISCE extends cell fate reprogramming studies to larger-scale biological networks without the need for full model parameterization and can be implemented by experimental and computational biologists to identify parts of a biological system relevant to the desired reprogramming task.
Collapse
Affiliation(s)
- Lauren Marazzi
- Center for Quantitative Medicine, University of Connecticut School of Medicine, Farmington, CT, 06030, USA
| | - Milan Shah
- Center for Quantitative Medicine, University of Connecticut School of Medicine, Farmington, CT, 06030, USA
| | - Shreedula Balakrishnan
- Center for Quantitative Medicine, University of Connecticut School of Medicine, Farmington, CT, 06030, USA
| | - Ananya Patil
- Center for Quantitative Medicine, University of Connecticut School of Medicine, Farmington, CT, 06030, USA
| | - Paola Vera-Licona
- Center for Quantitative Medicine, University of Connecticut School of Medicine, Farmington, CT, 06030, USA. .,Department of Cell Biology, University of Connecticut School of Medicine, Farmington, CT, 06030, USA. .,Center for Cell Analysis and Modeling, University of Connecticut School of Medicine, Farmington, CT, 06030, USA. .,Institute for Systems Genomics, University of Connecticut School of Medicine, Farmington, CT, 06030, USA.
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
Gondal MN, Butt RN, Shah OS, Sultan MU, Mustafa G, Nasir Z, Hussain R, Khawar H, Qazi R, Tariq M, Faisal A, Chaudhary SU. A Personalized Therapeutics Approach Using an In Silico Drosophila Patient Model Reveals Optimal Chemo- and Targeted Therapy Combinations for Colorectal Cancer. Front Oncol 2021; 11:692592. [PMID: 34336681 PMCID: PMC8323493 DOI: 10.3389/fonc.2021.692592] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Accepted: 06/30/2021] [Indexed: 12/18/2022] Open
Abstract
In silico models of biomolecular regulation in cancer, annotated with patient-specific gene expression data, can aid in the development of novel personalized cancer therapeutic strategies. Drosophila melanogaster is a well-established animal model that is increasingly being employed to evaluate such preclinical personalized cancer therapies. Here, we report five Boolean network models of biomolecular regulation in cells lining the Drosophila midgut epithelium and annotate them with colorectal cancer patient-specific mutation data to develop an in silico Drosophila Patient Model (DPM). We employed cell-type-specific RNA-seq gene expression data from the FlyGut-seq database to annotate and then validate these networks. Next, we developed three literature-based colorectal cancer case studies to evaluate cell fate outcomes from the model. Results obtained from analyses of the proposed DPM help: (i) elucidate cell fate evolution in colorectal tumorigenesis, (ii) validate cytotoxicity of nine FDA-approved CRC drugs, and (iii) devise optimal personalized treatment combinations. The personalized network models helped identify synergistic combinations of paclitaxel-regorafenib, paclitaxel-bortezomib, docetaxel-bortezomib, and paclitaxel-imatinib for treating different colorectal cancer patients. Follow-on therapeutic screening of six colorectal cancer patients from cBioPortal using this drug combination demonstrated a 100% increase in apoptosis and a 100% decrease in proliferation. In conclusion, this work outlines a novel roadmap for decoding colorectal tumorigenesis along with the development of personalized combinatorial therapeutics for preclinical translational studies.
Collapse
Affiliation(s)
- Mahnoor Naseer Gondal
- Biomedical Informatics Research Laboratory, Department of Biology, Lahore University of Management Sciences, Lahore, Pakistan
| | - Rida Nasir Butt
- Biomedical Informatics Research Laboratory, Department of Biology, Lahore University of Management Sciences, Lahore, Pakistan
| | - Osama Shiraz Shah
- Biomedical Informatics Research Laboratory, Department of Biology, Lahore University of Management Sciences, Lahore, Pakistan
| | - Muhammad Umer Sultan
- Biomedical Informatics Research Laboratory, Department of Biology, Lahore University of Management Sciences, Lahore, Pakistan
| | - Ghulam Mustafa
- Biomedical Informatics Research Laboratory, Department of Biology, Lahore University of Management Sciences, Lahore, Pakistan
| | - Zainab Nasir
- Biomedical Informatics Research Laboratory, Department of Biology, Lahore University of Management Sciences, Lahore, Pakistan
| | - Risham Hussain
- Biomedical Informatics Research Laboratory, Department of Biology, Lahore University of Management Sciences, Lahore, Pakistan
| | - Huma Khawar
- Biomedical Informatics Research Laboratory, Department of Biology, Lahore University of Management Sciences, Lahore, Pakistan
| | - Romena Qazi
- Department of Pathology, Shaukat Khanum Memorial Cancer Hospital and Research Centre, Lahore, Pakistan
| | - Muhammad Tariq
- Epigenetics Laboratory, Department of Biology, Lahore University of Management Sciences, Lahore, Pakistan
| | - Amir Faisal
- Cancer Therapeutics Laboratory, Department of Biology, Lahore University of Management Sciences, Lahore, Pakistan
| | - Safee Ullah Chaudhary
- Biomedical Informatics Research Laboratory, Department of Biology, Lahore University of Management Sciences, Lahore, Pakistan
| |
Collapse
|
6
|
AEON: Attractor Bifurcation Analysis of Parametrised Boolean Networks. COMPUTER AIDED VERIFICATION 2020. [PMCID: PMC7363220 DOI: 10.1007/978-3-030-53288-8_28] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
Abstract
Boolean networks (BNs) provide an effective modelling tool for various phenomena from science and engineering. Any long-term behaviour of a BN eventually converges to a so-called attractor. Depending on various logical parameters, the structure and quality of attractors can undergo a significant change, known as a bifurcation. We present a tool for analysing bifurcations in asynchronous parametrised Boolean networks. To fight the state-space and parameter-space explosion problem the tool uses a parallel semi-symbolic algorithm.
Collapse
|
7
|
Hou X, Li M, Jia C, Zhang X, Wang Y. Attractor - a new turning point in drug discovery. DRUG DESIGN DEVELOPMENT AND THERAPY 2019; 13:2957-2968. [PMID: 31686779 PMCID: PMC6709805 DOI: 10.2147/dddt.s216397] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/23/2019] [Accepted: 07/28/2019] [Indexed: 11/23/2022]
Abstract
Drug discovery for complex diseases can be viewed as a challenging problem in which the influence of compounds on dynamic features of disease system should be considered, especially the strategies escaping from the disease attractors. Moreover, escaping from the disease-related attractors has been proved to be a cue for the treatment of the complex diseases. The drug discovery methodology based on the attractor theory indicates new solutions for target identification, drug discovery and drug combination design. The methodology is based on the holism level of the organism and the features of system dynamics, so it has advantages for the classification of complex diseases and drug discovery. Currently, research results of this method have increased, which expand the insight scope for drug discovery. This article introduces the major drug discovery methods in the history of pharmacy development and their characteristics, so as to illustrate the reasons and inevitability of the appearance of attractor method, its position in the history of pharmacy development, and its advantages for drug discovery and design, thereby to prove that the attractor method can indeed become the next major drug development method. In addition, it provides a comprehensive description about the concept of attractor, the pipeline of attractor analysis, the common methods of each process and its research progress, so as to provide a macroscopic framework and optional methods and tools for the follow-up researchers.
Collapse
Affiliation(s)
- Xucan Hou
- Department of Traditional Chinese Medicine Information Fusion and Utilization, Beijing University of Chinese Medicine, Beijing, People's Republic of China
| | - Meng Li
- Department of Traditional Chinese Medicine Information Fusion and Utilization, Beijing University of Chinese Medicine, Beijing, People's Republic of China
| | - Congmin Jia
- Department of Traditional Chinese Medicine Information Fusion and Utilization, Beijing University of Chinese Medicine, Beijing, People's Republic of China
| | - Xianbao Zhang
- Department of Traditional Chinese Medicine Information Fusion and Utilization, Beijing University of Chinese Medicine, Beijing, People's Republic of China
| | - Yun Wang
- Department of Traditional Chinese Medicine Information Fusion and Utilization, Beijing University of Chinese Medicine, Beijing, People's Republic of China
| |
Collapse
|
8
|
Menini L, Possieri C, Tornambè A. Boolean network analysis through the joint use of linear algebra and algebraic geometry. J Theor Biol 2019; 472:46-53. [PMID: 30991072 DOI: 10.1016/j.jtbi.2019.04.012] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2018] [Revised: 03/06/2019] [Accepted: 04/11/2019] [Indexed: 11/26/2022]
Abstract
Among the various phenomena that can be modeled by Boolean networks, i.e., discrete-time dynamical systems with binary state variables, gene regulatory interactions are especially well known. Therefore, the analysis of Boolean networks is critical, e.g., to identify genetic pathways and to predict the effects of mutations on the cell functionality. Two methodologies (i.e., the semi-tensor product and the Gröbner bases over finite fields) have recently been proposed to tackle the problem of determining cycles and attractors (with the corresponding basin of attraction) for such systems. Here, it is shown that, by suitably coupling methodologies taken from these two fields (i.e., linear algebra and algebraic geometry), it is not only possible to determine cycles and attractors, but also to find closed-form solutions of the Boolean network. Such a goal is pursued by finding an immersion that recasts the Boolean dynamics in a linear form and by computing the closed-form solution of the latter system. The effectiveness of this technique is demonstrated by fully computing the solutions of the Boolean network modeling the differentiation of the Th-lymphocyte, a type of white blood cells involved in the human adaptive immune system.
Collapse
Affiliation(s)
- Laura Menini
- Dipartimento di Ingegneria Industriale, Università di Roma Tor Vergata, Roma 00133, Italy.
| | - Corrado Possieri
- Dipartimento di Elettronica e Telecomunicazioni, Politecnico di Torino, Torino 10129, Italy.
| | - Antonio Tornambè
- Dipartimento di Ingegneria Civile e Ingegneria Informatica, Università di Roma Tor Vergata, Roma 00133, Italy.
| |
Collapse
|
9
|
Mendes ND, Henriques R, Remy E, Carneiro J, Monteiro PT, Chaouiya C. Estimating Attractor Reachability in Asynchronous Logical Models. Front Physiol 2018; 9:1161. [PMID: 30245634 PMCID: PMC6137237 DOI: 10.3389/fphys.2018.01161] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2018] [Accepted: 08/02/2018] [Indexed: 12/12/2022] Open
Abstract
Logical models are well-suited to capture salient dynamical properties of regulatory networks. For networks controlling cell fate decisions, cell fates are associated with model attractors (stable states or cyclic attractors) whose identification and reachability properties are particularly relevant. While synchronous updates assume unlikely instantaneous or identical rates associated with component changes, the consideration of asynchronous updates is more realistic but, for large models, may hinder the analysis of the resulting non-deterministic concurrent dynamics. This complexity hampers the study of asymptotical behaviors, and most existing approaches suffer from efficiency bottlenecks, being generally unable to handle cyclical attractors and quantify attractor reachability. Here, we propose two algorithms providing probability estimates of attractor reachability in asynchronous dynamics. The first algorithm, named Firefront, exhaustively explores the state space from an initial state, and provides quasi-exact evaluations of the reachability probabilities of model attractors. The algorithm progresses in breadth, propagating the probabilities of each encountered state to its successors. Second, Avatar is an adapted Monte Carlo approach, better suited for models with large and intertwined transient and terminal cycles. Avatar iteratively explores the state space by randomly selecting trajectories and by using these random walks to estimate the likelihood of reaching an attractor. Unlike Monte Carlo simulations, Avatar is equipped to avoid getting trapped in transient cycles and to identify cyclic attractors. Firefront and Avatar are validated and compared to related methods, using as test cases logical models of synthetic and biological networks. Both algorithms are implemented as new functionalities of GINsim 3.0, a well-established software tool for logical modeling, providing executable GUI, Java API, and scripting facilities.
Collapse
Affiliation(s)
| | - Rui Henriques
- Department of Computer Science and Engineering, Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal.,Instituto de Engenharia de Sistemas e Computadores Investigação e Desenvolvimento, Lisbon, Portugal
| | - Elisabeth Remy
- Aix Marseille University, CNRS, Centrale Marseille, I2M UMR 7373, Marseille, France
| | | | - Pedro T Monteiro
- Department of Computer Science and Engineering, Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal.,Instituto de Engenharia de Sistemas e Computadores Investigação e Desenvolvimento, Lisbon, Portugal
| | | |
Collapse
|