1
|
Murrugarra D, Veliz-Cuba A, Dimitrova E, Kadelka C, Wheeler M, Laubenbacher R. Modular Control of Biological Networks. ARXIV 2024:arXiv:2401.12477v2. [PMID: 38344220 PMCID: PMC10854280] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/17/2024]
Abstract
The concept of control is central to understanding and applications of biological network models. Some of their key structural features relate to control functions, through gene regulation, signaling, or metabolic mechanisms, and computational models need to encode these. Applications of models often focus on model-based control, such as in biomedicine or metabolic engineering. This paper presents an approach to model-based control that exploits two common features of biological networks, namely their modular structure and canalizing features of their regulatory mechanisms. The paper focuses on intracellular regulatory networks, represented by Boolean network models. A main result of this paper is that control strategies can be identified by focusing on one module at a time. This paper also presents a criterion based on canalizing features of the regulatory rules to identify modules that do not contribute to network control and can be excluded. For even moderately sized networks, finding global control inputs is computationally very challenging. The modular approach presented here leads to a highly efficient approach to solving this problem. This approach is applied to a published Boolean network model of blood cancer large granular lymphocyte (T-LGL) leukemia to identify a minimal control set that achieves a desired control objective.
Collapse
Affiliation(s)
- David Murrugarra
- Department of Mathematics, University of Kentucky, Lexington, KY 40506, USA
| | - Alan Veliz-Cuba
- Department of Mathematics, University of Dayton, Dayton, OH 45469, USA
| | - Elena Dimitrova
- Mathematics Department, California Polytechnic State University, San Luis Obispo, CA 93407, USA
| | - Claus Kadelka
- Department of Mathematics, Iowa State University, Ames, IA 50011, USA
| | - Matthew Wheeler
- Department of Medicine, University of Florida, Gainesville, FL 32610, USA
| | | |
Collapse
|
2
|
Memariani M, Memariani H. Antifungal properties of cathelicidin LL-37: current knowledge and future research directions. World J Microbiol Biotechnol 2023; 40:34. [PMID: 38057654 DOI: 10.1007/s11274-023-03852-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2023] [Accepted: 11/22/2023] [Indexed: 12/08/2023]
Abstract
The threat of fungal diseases is substantially underestimated worldwide, but they have serious consequences for humans, animals, and plants. Given the limited number of existing antifungal drugs together with the emergence of drug-resistant strains, many researchers have actively sought alternatives or adjuvants to antimycotics. The best way to tackle these issues is to unearth potential antifungal agents with new modes of action. Antimicrobial peptides are being hailed as a promising source of novel antimicrobials since they exhibit rapid and broad-spectrum microbicidal activities with a reduced likelihood of developing drug resistance. Recent years have witnessed an explosion in knowledge on microbicidal activity of LL-37, the sole human cathelicidin. Herein, we provide a summary of the current understanding about antifungal properties of LL-37, with particular emphasis on its molecular mechanisms. We further illustrate fruitful areas for future research. LL-37 is able to inhibit the growth of clinically and agronomically relevant fungi including Aspergillus, Candida, Colletotrichum, Fusarium, Malassezia, Pythium, and Trichophyton. Destruction of the cell wall integrity, membrane permeabilization, induction of oxidative stress, disruption of endoplasmic reticulum homeostasis, formation of autophagy-like structures, alterations in expression of numerous fungal genes, and inhibition of cell cycle progression are the key mechanisms underlying antifungal effects of LL-37. Burgeoning evidence also suggests that LL-37 may act as a potential anti-virulence peptide. It is hoped that this review will not only motivate researchers to conduct more detailed studies in this field, but also inspire further innovations in the design of LL-37-based drugs for the treatment of fungal infections.
Collapse
Affiliation(s)
- Mojtaba Memariani
- Department of Medical Microbiology, Tehran University of Medical Sciences, Tehran, Iran
| | - Hamed Memariani
- Department of Medical Microbiology, Tehran University of Medical Sciences, Tehran, Iran.
| |
Collapse
|
3
|
Kadelka C, Wheeler M, Veliz-Cuba A, Murrugarra D, Laubenbacher R. Modularity of biological systems: a link between structure and function. J R Soc Interface 2023; 20:20230505. [PMID: 37876275 PMCID: PMC10598444 DOI: 10.1098/rsif.2023.0505] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2023] [Accepted: 10/05/2023] [Indexed: 10/26/2023] Open
Abstract
This paper addresses two topics in systems biology, the hypothesis that biological systems are modular and the problem of relating structure and function of biological systems. The focus here is on gene regulatory networks, represented by Boolean network models, a commonly used tool. Most of the research on gene regulatory network modularity has focused on network structure, typically represented through either directed or undirected graphs. But since gene regulation is a highly dynamic process as it determines the function of cells over time, it is natural to consider functional modularity as well. One of the main results is that the structural decomposition of a network into modules induces an analogous decomposition of the dynamic structure, exhibiting a strong relationship between network structure and function. An extensive simulation study provides evidence for the hypothesis that modularity might have evolved to increase phenotypic complexity while maintaining maximal dynamic robustness to external perturbations.
Collapse
Affiliation(s)
- Claus Kadelka
- Department of Mathematics, Iowa State University, Ames, IA, USA
| | - Matthew Wheeler
- Department of Medicine, University of Florida, Gainesville, FL, USA
| | - Alan Veliz-Cuba
- Department of Mathematics, University of Dayton, Dayton, OH, USA
| | - David Murrugarra
- Department of Mathematics, University of Kentucky, Lexington, KY, USA
| | | |
Collapse
|
4
|
Kadelka C, Wheeler M, Veliz-Cuba A, Murrugarra D, Laubenbacher R. Modularity of biological systems: a link between structure and function. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.09.11.557227. [PMID: 37745485 PMCID: PMC10515856 DOI: 10.1101/2023.09.11.557227] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/26/2023]
Abstract
This paper addresses two topics in systems biology, the hypothesis that biological systems are modular and the problem of relating structure and function of biological systems. The focus here is on gene regulatory networks, represented by Boolean network models, a commonly used tool. Most of the research on gene regulatory network modularity has focused on network structure, typically represented through either directed or undirected graphs. But since gene regulation is a highly dynamic process as it determines the function of cells over time, it is natural to consider functional modularity as well. One of the main results is that the structural decomposition of a network into modules induces an analogous decomposition of the dynamic structure, exhibiting a strong relationship between network structure and function. An extensive simulation study provides evidence for the hypothesis that modularity might have evolved to increase phenotypic complexity while maintaining maximal dynamic robustness to external perturbations.
Collapse
Affiliation(s)
- Claus Kadelka
- Department of Mathematics, Iowa State University, Ames, IA 50011, United States
| | - Matthew Wheeler
- Department of Medicine, University of Florida, Gainesville, FL, United States
| | - Alan Veliz-Cuba
- Department of Mathematics, University of Dayton, Dayton, OH, United States
| | - David Murrugarra
- Department of Mathematics, University of Kentucky, Lexington, KY, United States
| | | |
Collapse
|
5
|
Murrugarra D, Veliz-Cuba A. An information theoretic approach for the inference of Boolean networks and functions from data: BoCSE. PATTERNS (NEW YORK, N.Y.) 2022; 3:100617. [PMID: 36419439 PMCID: PMC9676479 DOI: 10.1016/j.patter.2022.100617] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Abstract
Building predictive models from data is an important and challenging task in many fields including biology, medicine, engineering, and economy. In this issue, Sun et al.1 present a method for the inference of Boolean networks along with practical applications.
Collapse
|
6
|
Gupta S, Panda PK, Luo W, Hashimoto RF, Ahuja R. Network analysis reveals that the tumor suppressor lncRNA GAS5 acts as a double-edged sword in response to DNA damage in gastric cancer. Sci Rep 2022; 12:18312. [PMID: 36316351 PMCID: PMC9622883 DOI: 10.1038/s41598-022-21492-x] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2022] [Accepted: 09/28/2022] [Indexed: 11/14/2022] Open
Abstract
The lncRNA GAS5 acts as a tumor suppressor and is downregulated in gastric cancer (GC). In contrast, E2F1, an important transcription factor and tumor promoter, directly inhibits miR-34c expression in GC cell lines. Furthermore, in the corresponding GC cell lines, lncRNA GAS5 directly targets E2F1. However, lncRNA GAS5 and miR-34c remain to be studied in conjunction with GC. Here, we present a dynamic Boolean network to classify gene regulation between these two non-coding RNAs (ncRNAs) in GC. This is the first study to show that lncRNA GAS5 can positively regulate miR-34c in GC through a previously unknown molecular pathway coupling lncRNA/miRNA. We compared our network to several in-vivo/in-vitro experiments and obtained an excellent agreement. We revealed that lncRNA GAS5 regulates miR-34c by targeting E2F1. Additionally, we found that lncRNA GAS5, independently of p53, inhibits GC proliferation through the ATM/p38 MAPK signaling pathway. Accordingly, our results support that E2F1 is an engaging target of drug development in tumor growth and aggressive proliferation of GC, and favorable results can be achieved through tumor suppressor lncRNA GAS5/miR-34c axis in GC. Thus, our findings unlock a new avenue for GC treatment in response to DNA damage by these ncRNAs.
Collapse
Affiliation(s)
- Shantanu Gupta
- grid.11899.380000 0004 1937 0722Departamento de Ciência da Computação, Instituto de Matemática e Estatística, Universidade de São Paulo, Rua do Matão 1010, São Paulo, SP 05508-090 Brasil
| | - Pritam Kumar Panda
- grid.8993.b0000 0004 1936 9457Condensed Matter Theory Group, Materials Theory Division, Department of Physics and Astronomy, Uppsala University, Box 516, 751 20 Uppsala, Sweden
| | - Wei Luo
- grid.8993.b0000 0004 1936 9457Condensed Matter Theory Group, Materials Theory Division, Department of Physics and Astronomy, Uppsala University, Box 516, 751 20 Uppsala, Sweden
| | - Ronaldo F. Hashimoto
- grid.11899.380000 0004 1937 0722Departamento de Ciência da Computação, Instituto de Matemática e Estatística, Universidade de São Paulo, Rua do Matão 1010, São Paulo, SP 05508-090 Brasil
| | - Rajeev Ahuja
- grid.8993.b0000 0004 1936 9457Condensed Matter Theory Group, Materials Theory Division, Department of Physics and Astronomy, Uppsala University, Box 516, 751 20 Uppsala, Sweden ,grid.462391.b0000 0004 1769 8011Department of Physics, Indian Institute of Technology Ropar, Rupnagar, Punjab 140001 India
| |
Collapse
|
7
|
Plaugher D, Aguilar B, Murrugarra D. Uncovering potential interventions for pancreatic cancer patients via mathematical modeling. J Theor Biol 2022; 548:111197. [PMID: 35752283 DOI: 10.1016/j.jtbi.2022.111197] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Revised: 05/20/2022] [Accepted: 06/08/2022] [Indexed: 12/24/2022]
Abstract
Pancreatic Ductal Adenocarcinoma (PDAC) is widely known for its poor prognosis because it is often diagnosed when the cancer is in a later stage. We built a Boolean model to analyze the microenvironment of pancreatic cancer in order to better understand the interplay between pancreatic cancer, stellate cells, and their signaling cytokines. Specifically, we have used our model to study the impact of inducing four common mutations: KRAS, TP53, SMAD4, and CDKN2A. After implementing the various mutation combinations, we used our stochastic simulator to derive aggressiveness scores based on simulated attractor probabilities and long-term trajectory approximations. These aggression scores were then corroborated with clinical data. Moreover, we found sets of control targets that are effective among common mutations. These control sets contain nodes within both the pancreatic cancer cell and the pancreatic stellate cell, including PIP3, RAF, PIK3 and BAX in pancreatic cancer cell as well as ERK and PIK3 in the pancreatic stellate cell. Many of these nodes were found to be differentially expressed among pancreatic cancer patients in the TCGA database. Furthermore, literature suggests that many of these nodes can be targeted by drugs currently in circulation. The results herein help provide a proof of concept in the path towards personalized medicine through a means of mathematical systems biology. All data and code used for running simulations, statistical analysis, and plotting is available on a GitHub repository athttps://github.com/drplaugher/PCC_Mutations.
Collapse
Affiliation(s)
- Daniel Plaugher
- Department of Mathematics, University of Kentucky, Lexington, KY, USA.
| | | | - David Murrugarra
- Department of Mathematics, University of Kentucky, Lexington, KY, USA.
| |
Collapse
|
8
|
Mazaya M, Kwon YK. In Silico Pleiotropy Analysis in KEGG Signaling Networks Using a Boolean Network Model. Biomolecules 2022; 12:biom12081139. [PMID: 36009032 PMCID: PMC9406064 DOI: 10.3390/biom12081139] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Revised: 08/10/2022] [Accepted: 08/15/2022] [Indexed: 11/16/2022] Open
Abstract
Pleiotropy, which refers to the ability of different mutations on the same gene to cause different pathological effects in human genetic diseases, is important in understanding system-level biological diseases. Although some biological experiments have been proposed, still little is known about pleiotropy on gene–gene dynamics, since most previous studies have been based on correlation analysis. Therefore, a new perspective is needed to investigate pleiotropy in terms of gene–gene dynamical characteristics. To quantify pleiotropy in terms of network dynamics, we propose a measure called in silico Pleiotropic Scores (sPS), which represents how much a gene is affected against a pair of different types of mutations on a Boolean network model. We found that our model can identify more candidate pleiotropic genes that are not known to be pleiotropic than the experimental database. In addition, we found that many types of functionally important genes tend to have higher sPS values than other genes; in other words, they are more pleiotropic. We investigated the relations of sPS with the structural properties in the signaling network and found that there are highly positive relations to degree, feedback loops, and centrality measures. This implies that the structural characteristics are principles to identify new pleiotropic genes. Finally, we found some biological evidence showing that sPS analysis is relevant to the real pleiotropic data and can be considered a novel candidate for pleiotropic gene research. Taken together, our results can be used to understand the dynamics pleiotropic characteristics in complex biological systems in terms of gene–phenotype relations.
Collapse
Affiliation(s)
- Maulida Mazaya
- Research Center for Computing, National Research and Innovation Agency (BRIN), Cibinong Science Center, Jl. Raya Jakarta-Bogor KM 46, Cibinong 16911, West Java, Indonesia
| | - Yung-Keun Kwon
- School of IT Convergence, University of Ulsan, 93 Daehak-ro, Nam-gu, Ulsan 44610, Korea
- Correspondence:
| |
Collapse
|
9
|
Werle SD, Ikonomi N, Schwab JD, Kraus JM, Weidner FM, Rudolph KL, Pfister AS, Schuler R, Kühl M, Kestler HA. Identification of dynamic driver sets controlling phenotypical landscapes. Comput Struct Biotechnol J 2022; 20:1603-1617. [PMID: 35465155 PMCID: PMC9010550 DOI: 10.1016/j.csbj.2022.03.034] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2022] [Revised: 03/30/2022] [Accepted: 03/30/2022] [Indexed: 11/03/2022] Open
Abstract
Controlling phenotypical landscapes is of vital interest to modern biology. This task becomes highly demanding because cellular decisions involve complex networks engaging in crosstalk interactions. Previous work on control theory indicates that small sets of compounds can control single phenotypes. However, a dynamic approach is missing to determine the drivers of the whole network dynamics. By analyzing 35 biologically motivated Boolean networks, we developed a method to identify small sets of compounds sufficient to decide on the entire phenotypical landscape. These compounds do not strictly prefer highly related compounds and show a smaller impact on the stability of the attractor landscape. The dynamic driver sets include many intervention targets and cellular reprogramming drivers in human networks. Finally, by using a new comprehensive model of colorectal cancer, we provide a complete workflow on how to implement our approach to shift from in silico to in vitro guided experiments.
Collapse
Affiliation(s)
- Silke D Werle
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - Nensi Ikonomi
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - Julian D Schwab
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - Johann M Kraus
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - Felix M Weidner
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - K Lenhard Rudolph
- Leibniz Institute of Aging - Fritz Lipman Institute, 07745 Jena, Thuringia, Germany
| | - Astrid S Pfister
- Institute of Biochemistry and Molecular Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - Rainer Schuler
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - Michael Kühl
- Institute of Biochemistry and Molecular Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| | - Hans A Kestler
- Institute of Medical Systems Biology, Ulm University, 89081 Ulm, Baden-Wuerttemberg, Germany
| |
Collapse
|
10
|
Dynamical Analysis of a Boolean Network Model of the Oncogene Role of lncRNA ANRIL and lncRNA UFC1 in Non-Small Cell Lung Cancer. Biomolecules 2022; 12:biom12030420. [PMID: 35327612 PMCID: PMC8946683 DOI: 10.3390/biom12030420] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2022] [Revised: 03/01/2022] [Accepted: 03/01/2022] [Indexed: 12/14/2022] Open
Abstract
Long non-coding RNA (lncRNA) such as ANRIL and UFC1 have been verified as oncogenic genes in non-small cell lung cancer (NSCLC). It is well known that the tumor suppressor microRNA-34a (miR-34a) is downregulated in NSCLC. Furthermore, miR-34a induces senescence and apoptosis in breast, glioma, cervical cancer including NSCLC by targeting Myc. Recent evidence suggests that these two lncRNAs act as a miR-34a sponge in corresponding cancers. However, the biological functions between these two non-coding RNAs (ncRNAs) have not yet been studied in NSCLC. Therefore, we present a Boolean model to analyze the gene regulation between these two ncRNAs in NSCLC. We compared our model to several experimental studies involving gain- or loss-of-function genes in NSCLC cells and achieved an excellent agreement. Additionally, we predict three positive circuits involving miR-34a/E2F1/ANRIL, miR-34a/E2F1/UFC1, and miR-34a/Myc/ANRIL. Our circuit- perturbation analysis shows that these circuits are important for regulating cell-fate decisions such as senescence and apoptosis. Thus, our Boolean network permits an explicit cell-fate mechanism associated with NSCLC. Therefore, our results support that ANRIL and/or UFC1 is an attractive target for drug development in tumor growth and aggressive proliferation of NSCLC, and that a valuable outcome can be achieved through the miRNA-34a/Myc pathway.
Collapse
|
11
|
Rozum JC, Deritei D, Park KH, Gómez Tejeda Zañudo J, Albert R. pystablemotifs: Python library for attractor identification and control in Boolean networks. Bioinformatics 2022; 38:1465-1466. [PMID: 34875008 DOI: 10.1093/bioinformatics/btab825] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2021] [Revised: 10/22/2021] [Accepted: 12/02/2021] [Indexed: 01/05/2023] Open
Abstract
SUMMARY pystablemotifs is a Python 3 library for analyzing Boolean networks. Its non-heuristic and exhaustive attractor identification algorithm was previously presented in Rozum et al. (2021). Here, we illustrate its performance improvements over similar methods and discuss how it uses outputs of the attractor identification process to drive a system to one of its attractors from any initial state. We implement six attractor control algorithms, five of which are new in this work. By design, these algorithms can return different control strategies, allowing for synergistic use. We also give a brief overview of the other tools implemented in pystablemotifs. AVAILABILITY AND IMPLEMENTATION The source code is on GitHub at https://github.com/jcrozum/pystablemotifs/. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Jordan C Rozum
- Department of Physics, The Pennsylvania State University, University Park, PA 16802, USA
| | - Dávid Deritei
- Department of Molecular Biology, Semmelweis University, Budapest 1085, Hungary.,Channing Division of Network Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA 02115, USA
| | - Kyu Hyong Park
- Department of Physics, The Pennsylvania State University, University Park, PA 16802, USA
| | - Jorge Gómez Tejeda Zañudo
- Eli and Edythe L. Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA.,Department of Medical Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA 02115, 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
|