1
|
Musilova J, Vafek Z, Puniya BL, Zimmer R, Helikar T, Sedlar K. Augusta: From RNA-Seq to gene regulatory networks and Boolean models. Comput Struct Biotechnol J 2024; 23:783-790. [PMID: 38312198 PMCID: PMC10837063 DOI: 10.1016/j.csbj.2024.01.013] [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: 10/26/2023] [Revised: 01/17/2024] [Accepted: 01/19/2024] [Indexed: 02/06/2024] Open
Abstract
Computational models of gene regulations help to understand regulatory mechanisms and are extensively used in a wide range of areas, e.g., biotechnology or medicine, with significant benefits. Unfortunately, there are only a few computational gene regulatory models of whole genomes allowing static and dynamic analysis due to the lack of sophisticated tools for their reconstruction. Here, we describe Augusta, an open-source Python package for Gene Regulatory Network (GRN) and Boolean Network (BN) inference from the high-throughput gene expression data. Augusta can reconstruct genome-wide models suitable for static and dynamic analyses. Augusta uses a unique approach where the first estimation of a GRN inferred from expression data is further refined by predicting transcription factor binding motifs in promoters of regulated genes and by incorporating verified interactions obtained from databases. Moreover, a refined GRN is transformed into a draft BN by searching in the curated model database and setting logical rules to incoming edges of target genes, which can be further manually edited as the model is provided in the SBML file format. The approach is applicable even if information about the organism under study is not available in the databases, which is typically the case for non-model organisms including most microbes. Augusta can be operated from the command line and, thus, is easy to use for automated prediction of models for various genomes. The Augusta package is freely available at github.com/JanaMus/Augusta. Documentation and tutorials are available at augusta.readthedocs.io.
Collapse
Affiliation(s)
- Jana Musilova
- Department of Biomedical Engineering, Faculty of Electrical Engineering and Communication, Brno University of Technology, Brno 61600, Czech Republic
- Department of Biochemistry, University of Nebraska-Lincoln, Lincoln 68588, NE, USA
| | - Zdenek Vafek
- Department of Biochemistry, University of Nebraska-Lincoln, Lincoln 68588, NE, USA
- Institute of Forensic Engineering, Brno University of Technology, Brno 61200, Czech Republic
| | - Bhanwar Lal Puniya
- Department of Biochemistry, University of Nebraska-Lincoln, Lincoln 68588, NE, USA
| | - Ralf Zimmer
- Department of Informatics, Ludwig-Maximilians-Universität München, Munich 80539, Germany
| | - Tomas Helikar
- Department of Biochemistry, University of Nebraska-Lincoln, Lincoln 68588, NE, USA
| | - Karel Sedlar
- Department of Biomedical Engineering, Faculty of Electrical Engineering and Communication, Brno University of Technology, Brno 61600, Czech Republic
- Department of Informatics, Ludwig-Maximilians-Universität München, Munich 80539, Germany
| |
Collapse
|
2
|
Andreu-Vilarroig C, Villanueva RJ, González-Parra G. Mathematical modeling for estimating influenza vaccine efficacy: A case study of the Valencian Community, Spain. Infect Dis Model 2024; 9:744-762. [PMID: 38689854 PMCID: PMC11058883 DOI: 10.1016/j.idm.2024.04.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2024] [Revised: 04/02/2024] [Accepted: 04/10/2024] [Indexed: 05/02/2024] Open
Abstract
Vaccine efficacy and its quantification is a crucial concept for the proper design of public health vaccination policies. In this work we proposed a mathematical model to estimate the efficacy of the influenza vaccine in a real-word scenario. In particular, our model is a SEIR-type epidemiological model, which distinguishes vaccinated and unvaccinated populations. Mathematically, its dynamics is governed by a nonlinear system of ordinary differential equations, where the non-linearity arises from the effective contacts between susceptible and infected individuals. Two key aspects of this study is that we use a vaccine distribution over time that is based on real data specific to the elderly people in the Valencian Community and the calibration process takes into account that over one influenza season a specific proportion of the population becomes infected with influenza. To consider the effectiveness of the vaccine, the model incorporates a parameter, the vaccine attenuation factor, which is related with the vaccine efficacy against the influenza virus. With this framework, in order to calibrate the model parameters and to obtain an influenza vaccine efficacy estimation, we considered the 2016-2017 influenza season in the Valencian Community, Spain, using the influenza reported cases of vaccinated and unvaccinated. In order to ensure the model identifiability, we choose to deterministically calibrate the parameters for different scenarios and we find the one with the minimum error in order to determine the vaccine efficacy. The calibration results suggest that the influenza vaccine developed for 2016-2017 influenza season has an efficacy of approximately 76.7%, and that the risk of becoming infected is five times higher for an unvaccinated individual in comparison with a vaccinated one. This estimation partially agrees with some previous studies related to the influenza vaccine. This study presents a new integrated mathematical approach to study the influenza vaccine efficacy and gives further insight into this important public health topic.
Collapse
Affiliation(s)
- Carlos Andreu-Vilarroig
- Instituto de Matemática Multidisciplinar, Universitat Politècnica de València, Valencia, Spain
| | - Rafael J. Villanueva
- Instituto de Matemática Multidisciplinar, Universitat Politècnica de València, Valencia, Spain
| | - Gilberto González-Parra
- Instituto de Matemática Multidisciplinar, Universitat Politècnica de València, Valencia, Spain
- Department of Mathematics, New Mexico Tech, Socorro, NM, USA
| |
Collapse
|
3
|
Bhat SS, Kulkarni SR, Uttarkar A, Niranjan V. Computational Insights into Papaveroline as an In Silico Drug Candidate for Alzheimer's Disease via Fyn Tyrosine Kinase Inhibition. Mol Biotechnol 2024:10.1007/s12033-024-01236-0. [PMID: 39004678 DOI: 10.1007/s12033-024-01236-0] [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: 12/06/2023] [Accepted: 06/25/2024] [Indexed: 07/16/2024]
Abstract
Alzheimer's disease (AD) poses a significant global health challenge, necessitating the exploration of novel therapeutic strategies. Fyn Tyrosine Kinase has emerged as a key player in AD pathogenesis, making it an attractive target for drug development. This study focuses on investigating the potential of Papaveroline as a drug candidate for AD by targeting Fyn Tyrosine Kinase. The research employed high-throughput virtual screening and QSAR analysis were conducted to identify compounds with optimal drug-like properties, emphasizing adherence to ADMET parameters for further evaluation. Molecular dynamics simulations to analyze the binding interactions between Papaveroline and Staurosporine with Fyn Tyrosine Kinase over a 200-ns period. The study revealed detailed insights into the binding mechanisms and stability of the Papaveroline-Fyn complex, showcasing the compound's potential as an inhibitor of Fyn Tyrosine Kinase. Comparative analysis with natural compounds and a reference compound highlighted Papaveroline's unique characteristics and promising therapeutic implications for AD treatment. Overall, the findings underscore Papaveroline's potential as a valuable drug candidate for targeting Fyn Tyrosine Kinase in AD therapy, offering new avenues for drug discovery in neurodegenerative diseases. This study contributes to advancing our understanding of molecular interactions in AD pathogenesis and paves the way for further research and development in this critical area.
Collapse
Affiliation(s)
- Shreya Satyanarayan Bhat
- Department of Biotechnology, R V College of Engineering (Affiliated to Visvesvaraya Technological University, Belagavi 590018), Bangalore, 560059, India
| | - Spoorthi R Kulkarni
- Department of Biotechnology, R V College of Engineering (Affiliated to Visvesvaraya Technological University, Belagavi 590018), Bangalore, 560059, India
| | - Akshay Uttarkar
- Department of Biotechnology, R V College of Engineering (Affiliated to Visvesvaraya Technological University, Belagavi 590018), Bangalore, 560059, India
| | - Vidya Niranjan
- Department of Biotechnology, R V College of Engineering (Affiliated to Visvesvaraya Technological University, Belagavi 590018), Bangalore, 560059, India.
| |
Collapse
|
4
|
Ramos RH, de Oliveira Lage Ferreira C, Simao A. Human protein-protein interaction networks: A topological comparison review. Heliyon 2024; 10:e27278. [PMID: 38562502 PMCID: PMC10982977 DOI: 10.1016/j.heliyon.2024.e27278] [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: 10/23/2023] [Revised: 02/26/2024] [Accepted: 02/27/2024] [Indexed: 04/04/2024] Open
Abstract
Protein-Protein Interaction Networks aim to model the interactome, providing a powerful tool for understanding the complex relationships governing cellular processes. These networks have numerous applications, including functional enrichment, discovering cancer driver genes, identifying drug targets, and more. Various databases make protein-protein networks available for many species, including Homo sapiens. This work topologically compares four Homo sapiens networks using a coarse-to-fine approach, comparing global characteristics, sub-network topology, specific nodes centrality, and interaction significance. Results show that the four human protein networks share many common protein-encoding genes and some global measures, but significantly differ in the interactions and neighbourhood. Small sub-networks from cancer pathways performed better than the whole networks, indicating an improved topological consistency in functional pathways. The centrality analysis shows that the same genes play different roles in different networks. We discuss how studies and analyses that rely on protein-protein networks for humans should consider their similarities and distinctions.
Collapse
Affiliation(s)
- Rodrigo Henrique Ramos
- University of São Paulo, São Carlos, SP, Brazil
- Federal Institute of São Paulo, São Carlos, SP, Brazil
| | | | | |
Collapse
|
5
|
Arjmand B, Khodadoost M, Jahani Sherafat S, Rezaei Tavirani M, Ahmadi N, Rezaei Tavirani S. Introducing critical proteins related to liver ischemia/reperfusion injury. GASTROENTEROLOGY AND HEPATOLOGY FROM BED TO BENCH 2024; 17:87-92. [PMID: 38737933 PMCID: PMC11080694 DOI: 10.22037/ghfbb.v17i1.2555] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/02/2023] [Accepted: 08/07/2023] [Indexed: 05/14/2024]
Abstract
Aim The current study aimed to introduce the key proteins involved in liver ischemia/reperfusion (I/R) injury through protein-protein interaction (PPI) analysis. Background Liver transplantation (LT) is a well-known treatment for liver diseases that threaten patients with mortality. LT is a complex operation, and several risks, including liver I/R injury, affect its success. Improving LT requires detection of its molecular mechanism. Experiments have revealed that high throughput methods such as proteomics in combination with bioinformatics are useful tools for analyzing the molecular mechanism of disease. Methods The differentially expressed proteins (DEPs) involved in liver I/R injury were extracted from the literature. The queried DEPs plus the first 100 neighbors were included in a network through STRING database using Cytoscape software. Degree, betweenness centrality, closeness centrality, and stress were considered to determine the central nodes. The queried DEPs were assessed by action map analysis using the CluePedia application of Cytoscape software. The key proteins were identified by comparing network analysis and action map evaluation results. Results Six proteins, namely ALB, INS, GAPDH, CAT, IL6, and TNF, among the added first neighbors were determined as the central first neighbors. MPO, CRP, MMP9, and HMOX1 were selected as central DEPs among the queried proteins. Action map analysis confirmed the PPI findings. The final evaluation revealed that MMP9 in combination with CRP and HMOX1 plays a critical role in liver I/R injury. Conclusion The significant role of MMP9 in liver I/R injury was detected in this study. Two central proteins (CRP and HMOX1) were shown to have a regulatory effect on MMP9; CRP activated MMP9, while HMXO1 downregulated it.
Collapse
Affiliation(s)
- Babak Arjmand
- Cell Therapy and Regenerative Medicine Research Center, Endocrinology and Metabolism Molecular-Cellular Sciences Institute, Tehran University of Medical Sciences, Tehran, Iran
| | - Mahmood Khodadoost
- School of Traditional Medicine, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Somayeh Jahani Sherafat
- Laser Application in Medical Sciences Research Center, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Mostafa Rezaei Tavirani
- Proteomics Research Center, Faculty of Paramedical Sciences, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Nayebali Ahmadi
- Proteomics Research Center, Faculty of Paramedical Sciences, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| | - Sina Rezaei Tavirani
- Proteomics Research Center, Faculty of Paramedical Sciences, Shahid Beheshti University of Medical Sciences, Tehran, Iran
| |
Collapse
|
6
|
Fang Y, Jiang Y, Wei L, Ma Q, Ren Z, Yuan Q, Wei DQ. DeepProSite: structure-aware protein binding site prediction using ESMFold and pretrained language model. Bioinformatics 2023; 39:btad718. [PMID: 38015872 PMCID: PMC10723037 DOI: 10.1093/bioinformatics/btad718] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2023] [Revised: 11/04/2023] [Accepted: 11/27/2023] [Indexed: 11/30/2023] Open
Abstract
MOTIVATION Identifying the functional sites of a protein, such as the binding sites of proteins, peptides, or other biological components, is crucial for understanding related biological processes and drug design. However, existing sequence-based methods have limited predictive accuracy, as they only consider sequence-adjacent contextual features and lack structural information. RESULTS In this study, DeepProSite is presented as a new framework for identifying protein binding site that utilizes protein structure and sequence information. DeepProSite first generates protein structures from ESMFold and sequence representations from pretrained language models. It then uses Graph Transformer and formulates binding site predictions as graph node classifications. In predicting protein-protein/peptide binding sites, DeepProSite outperforms state-of-the-art sequence- and structure-based methods on most metrics. Moreover, DeepProSite maintains its performance when predicting unbound structures, in contrast to competing structure-based prediction methods. DeepProSite is also extended to the prediction of binding sites for nucleic acids and other ligands, verifying its generalization capability. Finally, an online server for predicting multiple types of residue is established as the implementation of the proposed DeepProSite. AVAILABILITY AND IMPLEMENTATION The datasets and source codes can be accessed at https://github.com/WeiLab-Biology/DeepProSite. The proposed DeepProSite can be accessed at https://inner.wei-group.net/DeepProSite/.
Collapse
Affiliation(s)
- Yitian Fang
- State Key Laboratory of Microbial Metabolism, Shanghai-Islamabad-Belgrade Joint Innovation Center on Antibacterial Resistances, Joint International Research Laboratory of Metabolic & Developmental Sciences and School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200040, China
- Peng Cheng Laboratory, Shenzhen 518055, China
| | - Yi Jiang
- Department of Biomedical Informatics, College of Medicine, The Ohio State University, Columbus, OH 43210, USA
| | - Leyi Wei
- School of Software, Shandong University, Jinan, Shandong 250100, China
| | - Qin Ma
- Department of Biomedical Informatics, College of Medicine, The Ohio State University, Columbus, OH 43210, USA
| | | | - Qianmu Yuan
- School of Computer Science and Engineering, Sun Yat-sen University, Guangzhou 510000, China
| | - Dong-Qing Wei
- State Key Laboratory of Microbial Metabolism, Shanghai-Islamabad-Belgrade Joint Innovation Center on Antibacterial Resistances, Joint International Research Laboratory of Metabolic & Developmental Sciences and School of Life Sciences and Biotechnology, Shanghai Jiao Tong University, Shanghai 200040, China
- Peng Cheng Laboratory, Shenzhen 518055, China
| |
Collapse
|
7
|
Ayoup MS, Wahby Y, Abdel-Hamid H, Abu-Serie MM, Ramadan S, Barakat A, Teleb M, Ismail MMF. Reinvestigation of Passerini and Ugi scaffolds as multistep apoptotic inducers via dual modulation of caspase 3/7 and P53-MDM2 signaling for halting breast cancer. RSC Adv 2023; 13:27722-27737. [PMID: 37736568 PMCID: PMC10509784 DOI: 10.1039/d3ra04029a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2023] [Accepted: 08/28/2023] [Indexed: 09/23/2023] Open
Abstract
Selective induction of breast cancer apoptosis is viewed as the mainstay of various ongoing oncology drug discovery programs. Passerini scaffolds have been recently exploited as selective apoptosis inducers via a caspase 3/7 dependent pathway. Herein, the optimized Passerini caspase activators were manipulated to synergistically induce P53-dependent apoptosis via modulating the closely related P53-MDM2 signaling axis. The adopted design rationale and synthetic routes relied on mimicking the general thematic features of lead MDM2 inhibitors incorporating multiple aromatic rings. Accordingly, the cyclization of representative Passerini derivatives and related Ugi compounds into the corresponding diphenylimidazolidine and spiro derivative was performed, resembling the nutlin-based and spiro MDM-2 inhibitors, respectively. The study was also extended to explore the apoptotic induction capacity of the scaffold after simplification and modifications. MTT assay on MCF-7 and MDA-MB231 breast cancer cells compared to normal fibroblasts (WI-38) revealed their promising cytotoxic activities. The flexible Ugi derivatives 3 and 4, cyclic analog 8, Passerini adduct 12, and the thiosemicarbazide derivative 17 were identified as the study hits regarding cytotoxic potency and selectivity, being over 10-folds more potent (IC50 = 0.065-0.096 μM) and safer (SI = 4.4-18.7) than doxorubicin (IC50 = 0.478 μM, SI = 0.569) on MCF-7 cells. They promoted apoptosis induction via caspase 3/7 activation (3.1-4.1 folds) and P53 induction (up to 4 folds). Further apoptosis studies revealed that these compounds enhanced gene expression of BAX by 2 folds and suppressed Bcl-2 expression by 4.29-7.75 folds in the treated MCF-7 cells. Docking simulations displayed their plausible binding modes with the molecular targets and highlighted their structural determinants of activities for further optimization studies. Finally, in silico prediction of the entire library was computationally performed, showing that most of them could be envisioned as drug-like candidates.
Collapse
Affiliation(s)
- Mohammed Salah Ayoup
- Chemistry Department, Faculty of Science, Alexandria University P. O. Box 426 Alexandria 21321 Egypt
| | - Yasmin Wahby
- Chemistry Department, Faculty of Science, Alexandria University P. O. Box 426 Alexandria 21321 Egypt
| | - Hamida Abdel-Hamid
- Chemistry Department, Faculty of Science, Alexandria University P. O. Box 426 Alexandria 21321 Egypt
| | - Marwa M Abu-Serie
- Medical Biotechnology Department, Genetic Engineering and Biotechnology Research Institute, City of Scientific Research and Technological Applications (SRTA-City) Egypt
| | - Sherif Ramadan
- Chemistry Department, Michigan State University East Lansing MI 48824 USA
- Department of Chemistry, Benha University Benha Egypt
| | - Assem Barakat
- Department of Chemistry, College of Science, King Saud University P. O. Box 2455 Riyadh 11451 Saudi Arabia
| | - Mohamed Teleb
- Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Alexandria University Alexandria 21521 Egypt
| | - Magda M F Ismail
- Department of Pharmaceutical Medicinal Chemistry, Faculty of Pharmacy (Girls), Al-Azhar University Cairo 11754 Egypt
| |
Collapse
|
8
|
Vera-Siguenza E, Escribano-Gonzalez C, Serrano-Gonzalo I, Eskla KL, Spill F, Tennant D. Mathematical reconstruction of the metabolic network in an in-vitro multiple myeloma model. PLoS Comput Biol 2023; 19:e1011374. [PMID: 37713666 PMCID: PMC10503963 DOI: 10.1371/journal.pcbi.1011374] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2022] [Accepted: 07/19/2023] [Indexed: 09/17/2023] Open
Abstract
It is increasingly apparent that cancer cells, in addition to remodelling their metabolism to survive and proliferate, adapt and manipulate the metabolism of other cells. This property may be a telling sign that pre-clinical tumour metabolism studies exclusively utilising in-vitro mono-culture models could prove to be limited for uncovering novel metabolic targets able to translate into clinical therapies. Although this is increasingly recognised, and work towards addressing the issue is becoming routinary much remains poorly understood. For instance, knowledge regarding the biochemical mechanisms through which cancer cells manipulate non-cancerous cell metabolism, and the subsequent impact on their survival and proliferation remains limited. Additionally, the variations in these processes across different cancer types and progression stages, and their implications for therapy, also remain largely unexplored. This study employs an interdisciplinary approach that leverages the predictive power of mathematical modelling to enrich experimental findings. We develop a functional multicellular in-silico model that facilitates the qualitative and quantitative analysis of the metabolic network spawned by an in-vitro co-culture model of bone marrow mesenchymal stem- and myeloma cell lines. To procure this model, we devised a bespoke human genome constraint-based reconstruction workflow that combines aspects from the legacy mCADRE & Metabotools algorithms, the novel redHuman algorithm, along with 13C-metabolic flux analysis. Our workflow transforms the latest human metabolic network matrix (Recon3D) into two cell-specific models coupled with a metabolic network spanning a shared growth medium. When cross-validating our in-silico model against the in-vitro model, we found that the in-silico model successfully reproduces vital metabolic behaviours of its in-vitro counterpart; results include cell growth predictions, respiration rates, as well as support for observations which suggest cross-shuttling of redox-active metabolites between cells.
Collapse
Affiliation(s)
- Elias Vera-Siguenza
- Institute of Metabolism and Systems Research, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
- Watson School of Mathematics, University of Birmingham, Birmingham, United Kingdom
| | - Cristina Escribano-Gonzalez
- Institute of Metabolism and Systems Research, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
| | - Irene Serrano-Gonzalo
- Instituto de Investigación Sanitaria Aragón, Fundación Española para el Estudio y Terapéutica de la enfermedad de Gaucher y otras Lisosomales, Zaragoza, España
| | - Kattri-Liis Eskla
- Institute of Metabolism and Systems Research, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
- Department of Physiology, Institute of Biomedicine and Translational Medicine, University of Tartu, Tartu, Estonia
| | - Fabian Spill
- Watson School of Mathematics, University of Birmingham, Birmingham, United Kingdom
| | - Daniel Tennant
- Institute of Metabolism and Systems Research, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
| |
Collapse
|
9
|
Sarkar MS, Mia MM, Amin MA, Hossain MS, Islam MZ. Bioinformatics and network biology approach to identifying type 2 diabetes genes and pathways that influence the progression of breast cancer. Heliyon 2023; 9:e16151. [PMID: 37234659 PMCID: PMC10205526 DOI: 10.1016/j.heliyon.2023.e16151] [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/30/2022] [Revised: 04/28/2023] [Accepted: 05/07/2023] [Indexed: 05/28/2023] Open
Abstract
Breast cancer is the second most prevalent malignancy affecting women. Postmenopausal women breast tumor is one of the top causes of death in women, accounting for 23% of cancer cases. Type 2 diabetes, a worldwide pandemic, has been connected to a heightened risk of several malignancies, although its association with breast cancer is still uncertain. In comparison to non-diabetic women, women with T2DM had a 23% elevated likelihood of developing breast cancer. It is difficult to determine causative or genetic susceptibility that connect T2DM and breast cancer. We created a large-scale network-based quantitative approach employing unbiased methods to discover abnormally amplified genes in both T2DM and breast cancer, to solve these issues. We performed transcriptome analysis to uncover identical genetic biomarkers and pathways to clarify the connection between T2DM and breast cancer patients. In this study, two RNA-seq datasets (GSE103001 and GSE86468) from the Gene Expression Omnibus (GEO) are used to identify mutually differentially expressed genes (DEGs) for breast cancer and T2DM, as well as common pathways and prospective medicines. Firstly, 45 shared genes (30 upregulated and 15 downregulated) between T2D and breast cancer were detected. We employed gene ontology and pathway enrichment to characterize prevalent DEGs' molecular processes and signal transduction pathways and observed that T2DM has certain connections to the progression of breast cancer. Using several computational and statistical approaches, we created a protein-protein interactions (PPI) network and revealed hub genes. These hub genes can be potential biomarkers, which may also lead to new therapeutic strategies for investigated diseases. We conducted TF-gene interactions, gene-microRNA interactions, protein-drug interactions, and gene-disease associations to find potential connections between T2DM and breast cancer pathologies. We assume that the potential drugs that emerged from this study could be useful therapeutic values. Researchers, doctors, biotechnologists, and many others may benefit from this research.
Collapse
Affiliation(s)
- Md Sumon Sarkar
- Department of Pharmacy, Islamic University, Kushtia-7003, Bangladesh
| | - Md Misor Mia
- Department of Pharmacy, Islamic University, Kushtia-7003, Bangladesh
| | - Md Al Amin
- Department of Computer Science & Engineering, Prime University, Dhaka-1216, Bangladesh
| | - Md Sojib Hossain
- Department of Mathematics, Govt. Bangla College, Dhaka-1216, Bangladesh
| | - Md Zahidul Islam
- Department of Information & Communication Technology, Islamic University, Kushtia-7003, Bangladesh
| |
Collapse
|
10
|
Evolutionarily conserved gene expression patterns for affective disorders revealed using cross-species brain transcriptomic analyses in humans, rats and zebrafish. Sci Rep 2022; 12:20836. [PMID: 36460699 PMCID: PMC9718822 DOI: 10.1038/s41598-022-22688-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2022] [Accepted: 10/18/2022] [Indexed: 12/03/2022] Open
Abstract
Widespread, debilitating and often treatment-resistant, depression and other stress-related neuropsychiatric disorders represent an urgent unmet biomedical and societal problem. Although animal models of these disorders are commonly used to study stress pathogenesis, they are often difficult to translate across species into valuable and meaningful clinically relevant data. To address this problem, here we utilized several cross-species/cross-taxon approaches to identify potential evolutionarily conserved differentially expressed genes and their sets. We also assessed enrichment of these genes for transcription factors DNA-binding sites down- and up- stream from their genetic sequences. For this, we compared our own RNA-seq brain transcriptomic data obtained from chronically stressed rats and zebrafish with publicly available human transcriptomic data for patients with major depression and their respective healthy control groups. Utilizing these data from the three species, we next analyzed their differential gene expression, gene set enrichment and protein-protein interaction networks, combined with validated tools for data pooling. This approach allowed us to identify several key brain proteins (GRIA1, DLG1, CDH1, THRB, PLCG2, NGEF, IKZF1 and FEZF2) as promising, evolutionarily conserved and shared affective 'hub' protein targets, as well as to propose a novel gene set that may be used to further study affective pathogenesis. Overall, these approaches may advance cross-species brain transcriptomic analyses, and call for further cross-species studies into putative shared molecular mechanisms of affective pathogenesis.
Collapse
|
11
|
Wang C, Lye X, Kaalia R, Kumar P, Rajapakse JC. Deep learning and multi-omics approach to predict drug responses in cancer. BMC Bioinformatics 2022; 22:632. [PMID: 36443676 PMCID: PMC9703655 DOI: 10.1186/s12859-022-04964-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2022] [Accepted: 09/25/2022] [Indexed: 11/29/2022] Open
Abstract
BACKGROUND Cancers are genetically heterogeneous, so anticancer drugs show varying degrees of effectiveness on patients due to their differing genetic profiles. Knowing patient's responses to numerous cancer drugs are needed for personalized treatment for cancer. By using molecular profiles of cancer cell lines available from Cancer Cell Line Encyclopedia (CCLE) and anticancer drug responses available in the Genomics of Drug Sensitivity in Cancer (GDSC), we will build computational models to predict anticancer drug responses from molecular features. RESULTS We propose a novel deep neural network model that integrates multi-omics data available as gene expressions, copy number variations, gene mutations, reverse phase protein array expressions, and metabolomics expressions, in order to predict cellular responses to known anti-cancer drugs. We employ a novel graph embedding layer that incorporates interactome data as prior information for prediction. Moreover, we propose a novel attention layer that effectively combines different omics features, taking their interactions into account. The network outperformed feedforward neural networks and reported 0.90 for [Formula: see text] values for prediction of drug responses from cancer cell lines data available in CCLE and GDSC. CONCLUSION The outstanding results of our experiments demonstrate that the proposed method is capable of capturing the interactions of genes and proteins, and integrating multi-omics features effectively. Furthermore, both the results of ablation studies and the investigations of the attention layer imply that gene mutation has a greater influence on the prediction of drug responses than other omics data types. Therefore, we conclude that our approach can not only predict the anti-cancer drug response precisely but also provides insights into reaction mechanisms of cancer cell lines and drugs as well.
Collapse
Affiliation(s)
- Conghao Wang
- grid.59025.3b0000 0001 2224 0361School of Computer Science and Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore, 639798 Singapore
| | - Xintong Lye
- grid.59025.3b0000 0001 2224 0361School of Computer Science and Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore, 639798 Singapore
| | - Rama Kaalia
- grid.59025.3b0000 0001 2224 0361School of Computer Science and Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore, 639798 Singapore
| | - Parvin Kumar
- grid.59025.3b0000 0001 2224 0361School of Computer Science and Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore, 639798 Singapore
| | - Jagath C. Rajapakse
- grid.59025.3b0000 0001 2224 0361School of Computer Science and Engineering, Nanyang Technological University, 50 Nanyang Avenue, Singapore, 639798 Singapore
| |
Collapse
|
12
|
Filo M, Kumar S, Khammash M. A hierarchy of biomolecular proportional-integral-derivative feedback controllers for robust perfect adaptation and dynamic performance. Nat Commun 2022; 13:2119. [PMID: 35440114 PMCID: PMC9018779 DOI: 10.1038/s41467-022-29640-7] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2021] [Accepted: 03/25/2022] [Indexed: 01/05/2023] Open
Abstract
Proportional-Integral-Derivative (PID) feedback controllers are the most widely used controllers in industry. Recently, the design of molecular PID-controllers has been identified as an important goal for synthetic biology and the field of cybergenetics. In this paper, we consider the realization of PID-controllers via biomolecular reactions. We propose an array of topologies offering a compromise between simplicity and high performance. We first demonstrate that different biomolecular PI-controllers exhibit different performance-enhancing capabilities. Next, we introduce several derivative controllers based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop, which can be augmented by PI-components to yield PID-controllers. We demonstrate that PID-controllers can enhance stability and dynamic performance, and can also reduce stochastic noise. Finally, we provide an experimental demonstration using a hybrid setup where in silico PID-controllers regulate a genetic circuit in single yeast cells.
Collapse
Affiliation(s)
- Maurice Filo
- Department of Biosystems Science and Engineering, ETH Zürich, Mattenstrasse 26, 4058, Basel, Switzerland
| | - Sant Kumar
- Department of Biosystems Science and Engineering, ETH Zürich, Mattenstrasse 26, 4058, Basel, Switzerland
| | - Mustafa Khammash
- Department of Biosystems Science and Engineering, ETH Zürich, Mattenstrasse 26, 4058, Basel, Switzerland.
| |
Collapse
|
13
|
Röhl A, Baek SH, Kachroo P, Morrow JD, Tantisira K, Silverman EK, Weiss ST, Sharma A, Glass K, DeMeo DL. Protein interaction networks provide insight into fetal origins of chronic obstructive pulmonary disease. Respir Res 2022; 23:69. [PMID: 35331221 PMCID: PMC8944072 DOI: 10.1186/s12931-022-01963-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2021] [Accepted: 02/08/2022] [Indexed: 12/31/2022] Open
Abstract
BACKGROUND Chronic obstructive pulmonary disease (COPD) is a leading cause of death in adults that may have origins in early lung development. It is a complex disease, influenced by multiple factors including genetic variants and environmental factors. Maternal smoking during pregnancy may influence the risk for diseases during adulthood, potentially through epigenetic modifications including methylation. METHODS In this work, we explore the fetal origins of COPD by utilizing lung DNA methylation marks associated with in utero smoke (IUS) exposure, and evaluate the network relationships between methylomic and transcriptomic signatures associated with adult lung tissue from former smokers with and without COPD. To identify potential pathobiological mechanisms that may link fetal lung, smoke exposure and adult lung disease, we study the interactions (physical and functional) of identified genes using protein-protein interaction networks. RESULTS We build IUS-exposure and COPD modules, which identify connected subnetworks linking fetal lung smoke exposure to adult COPD. Studying the relationships and connectivity among the different modules for fetal smoke exposure and adult COPD, we identify enriched pathways, including the AGE-RAGE and focal adhesion pathways. CONCLUSIONS The modules identified in our analysis add new and potentially important insights to understanding the early life molecular perturbations related to the pathogenesis of COPD. We identify AGE-RAGE and focal adhesion as two biologically plausible pathways that may reveal lung developmental contributions to COPD. We were not only able to identify meaningful modules but were also able to study interconnections between smoke exposure and lung disease, augmenting our knowledge about the fetal origins of COPD.
Collapse
Affiliation(s)
- Annika Röhl
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, 02115, USA.
| | - Seung Han Baek
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, 02115, USA
| | - Priyadarshini Kachroo
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, 02115, USA
| | - Jarrett D Morrow
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, 02115, USA
| | - Kelan Tantisira
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, 02115, USA
- Division of Pediatric Respiratory Medicine, University of California San Diego, San Diego, USA
| | - Edwin K Silverman
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, 02115, USA
- Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, MA, USA
| | - Scott T Weiss
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, 02115, USA
| | - Amitabh Sharma
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, 02115, USA
- Center for Complex Network Research, Northeastern University, Boston, MA, USA
| | - Kimberly Glass
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, 02115, USA
| | - Dawn L DeMeo
- Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA, 02115, USA
- Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, MA, USA
| |
Collapse
|
14
|
Yuan Y, Qu A. High-Order Joint Embedding for Multi-Level Link Prediction. J Am Stat Assoc 2022. [DOI: 10.1080/01621459.2021.2005608] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Yubai Yuan
- Department of Statistics, University of California, Irvine, Irvine, CA
| | - Annie Qu
- Department of Statistics, University of California, Irvine, Irvine, CA
| |
Collapse
|
15
|
Wet-dry cycles protect surface-colonizing bacteria from major antibiotic classes. THE ISME JOURNAL 2022; 16:91-100. [PMID: 34253853 PMCID: PMC8692528 DOI: 10.1038/s41396-021-01051-4] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/21/2021] [Revised: 06/24/2021] [Accepted: 06/25/2021] [Indexed: 02/06/2023]
Abstract
Diverse antibiotic compounds are abundant in microbial habitats undergoing recurrent wet-dry cycles, such as soil, root and leaf surfaces, and the built environment. These antibiotics play a central role in microbial warfare and competition, thus affecting population dynamics and the composition of natural microbial communities. Yet, the impact of wet-dry cycles on bacterial response to antibiotics has been scarcely explored. Using the bacterium E. coli as a model organism, we show through a combination of experiments and computational modeling, that wet-dry cycles protect bacteria from beta-lactams. This is due to the combined effect of several mechanisms including tolerance induced by high salt concentrations and slow cell-growth, which are inherently associated with microscopic surface wetness-a hydration state typical to 'dry' periods. Moreover, we find evidence for a cross-protection effect, where lethal doses of antibiotic considerably increase bacterial survival during the dry periods. This work focuses on beta-lactams, yet similar protection was observed for additional major antibiotic classes. Our findings shed new light on how we understand bacterial response to antibiotics, with broad implications for population dynamics, interspecies interactions, and the evolution of antibiotic resistance in vast terrestrial microbial habitats.
Collapse
|
16
|
Mokhtar MM, El Allali A, Hegazy MEF, Atia MAM. PlantPathMarks (PPMdb): an interactive hub for pathways-based markers in plant genomes. Sci Rep 2021; 11:21300. [PMID: 34716373 PMCID: PMC8556342 DOI: 10.1038/s41598-021-00504-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2021] [Accepted: 09/06/2021] [Indexed: 11/12/2022] Open
Abstract
Over the past decade, the problem of finding an efficient gene-targeting marker set or signature for plant trait characterization has remained challenging. Many databases focusing on pathway mining have been released with one major deficiency, as they lack to develop marker sets that target only genes controlling a specific pathway or certain biological process. Herein, we present the PlantPathMarks database (PPMdb) as a comprehensive, web-based, user-friendly, and interactive hub for pathway-based markers in plant genomes. Based on our newly developed pathway gene set mining approach, two novel pathway-based marker systems called pathway gene-targeted markers (PGTMs) and pathway microsatellite-targeted markers (PMTMs) were developed as a novel class of annotation-based markers. In the PPMdb database, 2,690,742 pathway-based markers reflecting 9,894 marker panels were developed across 82 plant genomes. The markers include 691,555 PGTMs and 1,999,187 PMTMs. Across these genomes, 165,378 enzyme-coding genes were mapped against 126 KEGG reference pathway maps. PPMdb is furnished with three interactive visualization tools (Map Browse, JBrowse and Species Comparison) to visualize, map, and compare the developed markers over their KEGG reference pathway maps. All the stored marker panels can be freely downloaded. PPMdb promises to create a radical shift in the paradigm of the area of molecular marker research. The use of PPMdb as a mega-tool represents an impediment for non-bioinformatician plant scientists and breeders. PPMdb is freely available at http://ppmdb.easyomics.org.
Collapse
Affiliation(s)
- Morad M Mokhtar
- African Genome Center, Mohammed VI Polytechnic University, Ben Guerir, Morocco
| | - Achraf El Allali
- African Genome Center, Mohammed VI Polytechnic University, Ben Guerir, Morocco.
| | | | - Mohamed A M Atia
- Molecular Genetics and Genome Mapping Laboratory, Genome Mapping Department, Agricultural Genetic Engineering Research Institute (AGERI), Agriculture Research Center (ARC), Giza, 12619, Egypt.
| |
Collapse
|
17
|
Li S, Zhang Z, Li X, Tan Y, Wang L, Chen Z. An iteration model for identifying essential proteins by combining comprehensive PPI network with biological information. BMC Bioinformatics 2021; 22:430. [PMID: 34496745 PMCID: PMC8425031 DOI: 10.1186/s12859-021-04300-7] [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: 08/04/2020] [Accepted: 07/08/2021] [Indexed: 11/10/2022] Open
Abstract
Background Essential proteins have great impacts on cell survival and development, and played important roles in disease analysis and new drug design. However, since it is inefficient and costly to identify essential proteins by using biological experiments, then there is an urgent need for automated and accurate detection methods. In recent years, the recognition of essential proteins in protein interaction networks (PPI) has become a research hotspot, and many computational models for predicting essential proteins have been proposed successively. Results In order to achieve higher prediction performance, in this paper, a new prediction model called TGSO is proposed. In TGSO, a protein aggregation degree network is constructed first by adopting the node density measurement method for complex networks. And simultaneously, a protein co-expression interactive network is constructed by combining the gene expression information with the network connectivity, and a protein co-localization interaction network is constructed based on the subcellular localization data. And then, through integrating these three kinds of newly constructed networks, a comprehensive protein–protein interaction network will be obtained. Finally, based on the homology information, scores can be calculated out iteratively for different proteins, which can be utilized to estimate the importance of proteins effectively. Moreover, in order to evaluate the identification performance of TGSO, we have compared TGSO with 13 different latest competitive methods based on three kinds of yeast databases. And experimental results show that TGSO can achieve identification accuracies of 94%, 82% and 72% out of the top 1%, 5% and 10% candidate proteins respectively, which are to some degree superior to these state-of-the-art competitive models. Conclusions We constructed a comprehensive interactive network based on multi-source data to reduce the noise and errors in the initial PPI, and combined with iterative methods to improve the accuracy of necessary protein prediction, and means that TGSO may be conducive to the future development of essential protein recognition as well.
Collapse
Affiliation(s)
- Shiyuan Li
- College of Computer Engineering and Applied Mathematics, Changsha University, Changsha, 410022, China.,Hunan Province Key Laboratory of Industrial Internet Technology and Security, Changsha University, Changsha, 410022, China
| | - Zhen Zhang
- College of Electronic Information and Electrical Engineering, Changsha University, Changsha, 410022, China
| | - Xueyong Li
- College of Computer Engineering and Applied Mathematics, Changsha University, Changsha, 410022, China.,Hunan Province Key Laboratory of Industrial Internet Technology and Security, Changsha University, Changsha, 410022, China
| | - Yihong Tan
- College of Computer Engineering and Applied Mathematics, Changsha University, Changsha, 410022, China. .,Hunan Province Key Laboratory of Industrial Internet Technology and Security, Changsha University, Changsha, 410022, China.
| | - Lei Wang
- College of Computer Engineering and Applied Mathematics, Changsha University, Changsha, 410022, China.,Hunan Province Key Laboratory of Industrial Internet Technology and Security, Changsha University, Changsha, 410022, China
| | - Zhiping Chen
- College of Computer Engineering and Applied Mathematics, Changsha University, Changsha, 410022, China. .,Hunan Province Key Laboratory of Industrial Internet Technology and Security, Changsha University, Changsha, 410022, China.
| |
Collapse
|
18
|
Reconstruction and analysis of genome-scale metabolic model of weak Crabtree positive yeast Lachancea kluyveri. Sci Rep 2020; 10:16314. [PMID: 33004914 PMCID: PMC7530994 DOI: 10.1038/s41598-020-73253-3] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2020] [Accepted: 09/04/2020] [Indexed: 01/15/2023] Open
Abstract
Lachancea kluyveri, a weak Crabtree positive yeast, has been extensively studied for its unique URC pyrimidine catabolism pathway. It produces more biomass than Saccharomyces cerevisiae due to the underlying weak Crabtree effect and resorts to fermentation only in oxygen limiting conditions that renders it as a suitable industrial host. The yeast also produces ethyl acetate as a major overflow metabolite in aerobic conditions. Here, we report the first genome-scale metabolic model, iPN730, of L. kluyveri comprising of 1235 reactions, 1179 metabolites, and 730 genes distributed in 8 compartments. The in silico viability in different media conditions and the growth characteristics in various carbon sources show good agreement with experimental data. Dynamic flux balance analysis describes the growth dynamics, substrate utilization and product formation kinetics in various oxygen-limited conditions. We have also demonstrated the effect of switching carbon sources on the production of ethyl acetate under varying oxygen uptake rates. A phenotypic phase plane analysis described the energetic cost penalty of ethyl acetate and ethanol production on the specific growth rate of L. kluyveri. We generated the context specific models of L. kluyveri growing on uracil or ammonium salts as the sole nitrogen source. Differential flux calculated using flux variability analysis helped us in highlighting pathways like purine, histidine, riboflavin and pyrimidine metabolism associated with uracil degradation. The genome-scale metabolic construction of L. kluyveri will provide a better understanding of metabolism behind ethyl acetate production as well as uracil catabolism (pyrimidine degradation) pathway. iPN730 is an addition to genome-scale metabolic models of non-conventional yeasts that will facilitate system-wide omics analysis to understand fungal metabolic diversity.
Collapse
|
19
|
Correia K, Mahadevan R. Pan‐Genome‐Scale Network Reconstruction: Harnessing Phylogenomics Increases the Quantity and Quality of Metabolic Models. Biotechnol J 2020; 15:e1900519. [DOI: 10.1002/biot.201900519] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2019] [Revised: 07/22/2020] [Indexed: 12/31/2022]
Affiliation(s)
- Kevin Correia
- Department of Chemical Engineering and Applied Chemistry University of Toronto 200 College Street Toronto Ontario M5S 3E5 Canada
| | - Radhakrishnan Mahadevan
- Department of Chemical Engineering and Applied Chemistry University of Toronto 200 College Street Toronto Ontario M5S 3E5 Canada
- Institute of Biomedical Engineering University of Toronto 164 College Street Toronto Ontario M5S 3G9 Canada
| |
Collapse
|
20
|
Çubuk C, Hidalgo MR, Amadoz A, Rian K, Salavert F, Pujana MA, Mateo F, Herranz C, Carbonell-Caballero J, Dopazo J. Differential metabolic activity and discovery of therapeutic targets using summarized metabolic pathway models. NPJ Syst Biol Appl 2019; 5:7. [PMID: 30854222 PMCID: PMC6397295 DOI: 10.1038/s41540-019-0087-2] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2018] [Accepted: 02/13/2019] [Indexed: 12/13/2022] Open
Abstract
In spite of the increasing availability of genomic and transcriptomic data, there is still a gap between the detection of perturbations in gene expression and the understanding of their contribution to the molecular mechanisms that ultimately account for the phenotype studied. Alterations in the metabolism are behind the initiation and progression of many diseases, including cancer. The wealth of available knowledge on metabolic processes can therefore be used to derive mechanistic models that link gene expression perturbations to changes in metabolic activity that provide relevant clues on molecular mechanisms of disease and drug modes of action (MoA). In particular, pathway modules, which recapitulate the main aspects of metabolism, are especially suitable for this type of modeling. We present Metabolizer, a web-based application that offers an intuitive, easy-to-use interactive interface to analyze differences in pathway metabolic module activities that can also be used for class prediction and in silico prediction of knock-out (KO) effects. Moreover, Metabolizer can automatically predict the optimal KO intervention for restoring a diseased phenotype. We provide different types of validations of some of the predictions made by Metabolizer. Metabolizer is a web tool that allows understanding molecular mechanisms of disease or the MoA of drugs within the context of the metabolism by using gene expression measurements. In addition, this tool automatically suggests potential therapeutic targets for individualized therapeutic interventions.
Collapse
Affiliation(s)
- Cankut Çubuk
- Clinical Bioinformatics Area, Fundación Progreso y Salud, CDCA, Hospital Virgen del Rocío, 41013, Sevilla, Spain
| | - Marta R Hidalgo
- Bioinformatics and Biostatistics Unit, Centro de Investigación Príncipe Felipe (CIPF), 46012, Valencia, Spain
| | - Alicia Amadoz
- Department of Bioinformatics, Igenomix S.L, 46980, Valencia, Spain
| | - Kinza Rian
- Clinical Bioinformatics Area, Fundación Progreso y Salud, CDCA, Hospital Virgen del Rocío, 41013, Sevilla, Spain
| | - Francisco Salavert
- Bioinformatics in Rare Diseases (BiER), Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER), FPS, Hospital Virgen del Rocío, 41013, Sevilla, Spain
| | - Miguel A Pujana
- ProCURE. Catalan Institute of Oncology, Bellvitge Institute for Biomedical Research (IDIBELL), L'Hospitalet del Llobregat, 08908, Barcelona, Spain
| | - Francesca Mateo
- ProCURE. Catalan Institute of Oncology, Bellvitge Institute for Biomedical Research (IDIBELL), L'Hospitalet del Llobregat, 08908, Barcelona, Spain
| | - Carmen Herranz
- ProCURE. Catalan Institute of Oncology, Bellvitge Institute for Biomedical Research (IDIBELL), L'Hospitalet del Llobregat, 08908, Barcelona, Spain
| | - Jose Carbonell-Caballero
- Chromatin and Gene expression Lab, Gene Regulation, Stem Cells and Cancer Program, Centre de Regulació Genòmica (CRG), The Barcelona Institute of Science and Technology, PRBB, 08003, Barcelona, Spain
| | - Joaquín Dopazo
- Clinical Bioinformatics Area, Fundación Progreso y Salud, CDCA, Hospital Virgen del Rocío, 41013, Sevilla, Spain.
- Bioinformatics in Rare Diseases (BiER), Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER), FPS, Hospital Virgen del Rocío, 41013, Sevilla, Spain.
- Functional Genomics Node, INB-ELIXIR-es, FPS, Hospital Virgen del Rocío, 41013, Sevilla, Spain.
| |
Collapse
|
21
|
Wu S, Xia X, Wang D, Zhou Z, Wang G. Gene function and expression regulation of RuvRCAB in bacterial Cr(VI), As(III), Sb(III), and Cd(II) resistance. Appl Microbiol Biotechnol 2019; 103:2701-2713. [DOI: 10.1007/s00253-019-09666-6] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2018] [Revised: 01/23/2019] [Accepted: 01/24/2019] [Indexed: 12/24/2022]
|
22
|
Zou L, Chen J, Ruan S. Modeling and analyzing the transmission dynamics of visceral leishmaniasis. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2018; 14:1585-1604. [PMID: 29161877 DOI: 10.3934/mbe.2017082] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
In this paper, we develop a mathematical model to study the transmission dynamics of visceral leishmaniasis. Three populations: dogs, sandflies and humans, are considered in the model. Based on recent studies, we include vertical transmission of dogs in the spread of the disease. We also investigate the impact of asymptomatic humans and dogs as secondary reservoirs of the parasites. The basic reproduction number and sensitivity analysis show that the control of dog-sandfly transmission is more important for the elimination of the disease. Vaccination of susceptible dogs, treatment of infective dogs, as well as control of vertical transmission in dogs are effective prevention and control measures for visceral leishmaniasis.
Collapse
Affiliation(s)
- Lan Zou
- Department of Mathematics, Sichuan University, Chengdu, Sichuan 610064, China
| | - Jing Chen
- Department of Mathematics, University of Miami, Coral Gables, FL 33146, United States
| | - Shigui Ruan
- Department of Mathematics, University of Miami, Coral Gables, FL 33146, United States
| |
Collapse
|
23
|
Tebani A, Afonso C, Bekri S. Advances in metabolome information retrieval: turning chemistry into biology. Part II: biological information recovery. J Inherit Metab Dis 2018; 41:393-406. [PMID: 28842777 PMCID: PMC5959951 DOI: 10.1007/s10545-017-0080-0] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/13/2017] [Revised: 07/27/2017] [Accepted: 07/28/2017] [Indexed: 12/11/2022]
Abstract
This work reports the second part of a review intending to give the state of the art of major metabolic phenotyping strategies. It particularly deals with inherent advantages and limits regarding data analysis issues and biological information retrieval tools along with translational challenges. This Part starts with introducing the main data preprocessing strategies of the different metabolomics data. Then, it describes the main data analysis techniques including univariate and multivariate aspects. It also addresses the challenges related to metabolite annotation and characterization. Finally, functional analysis including pathway and network strategies are discussed. The last section of this review is devoted to practical considerations and current challenges and pathways to bring metabolomics into clinical environments.
Collapse
Affiliation(s)
- Abdellah Tebani
- Department of Metabolic Biochemistry, Rouen University Hospital, 76000, Rouen, France
- Normandie Université, UNIROUEN, CHU Rouen, IRIB, INSERM U1245, 76000, Rouen, France
- Normandie Université, UNIROUEN, INSA Rouen, CNRS, COBRA, 76000, Rouen, France
| | - Carlos Afonso
- Normandie Université, UNIROUEN, INSA Rouen, CNRS, COBRA, 76000, Rouen, France
| | - Soumeya Bekri
- Department of Metabolic Biochemistry, Rouen University Hospital, 76000, Rouen, France.
- Normandie Université, UNIROUEN, CHU Rouen, IRIB, INSERM U1245, 76000, Rouen, France.
| |
Collapse
|
24
|
Scharm M, Gebhardt T, Touré V, Bagnacani A, Salehzadeh-Yazdi A, Wolkenhauer O, Waltemath D. Evolution of computational models in BioModels Database and the Physiome Model Repository. BMC SYSTEMS BIOLOGY 2018; 12:53. [PMID: 29650016 PMCID: PMC5898004 DOI: 10.1186/s12918-018-0553-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/13/2017] [Accepted: 02/21/2018] [Indexed: 11/25/2022]
Abstract
BACKGROUND A useful model is one that is being (re)used. The development of a successful model does not finish with its publication. During reuse, models are being modified, i.e. expanded, corrected, and refined. Even small changes in the encoding of a model can, however, significantly affect its interpretation. Our motivation for the present study is to identify changes in models and make them transparent and traceable. METHODS We analysed 13734 models from BioModels Database and the Physiome Model Repository. For each model, we studied the frequencies and types of updates between its first and latest release. To demonstrate the impact of changes, we explored the history of a Repressilator model in BioModels Database. RESULTS We observed continuous updates in the majority of models. Surprisingly, even the early models are still being modified. We furthermore detected that many updates target annotations, which improves the information one can gain from models. To support the analysis of changes in model repositories we developed MoSt, an online tool for visualisations of changes in models. The scripts used to generate the data and figures for this study are available from GitHub https://github.com/binfalse/BiVeS-StatsGenerator and as a Docker image at https://hub.docker.com/r/binfalse/bives-statsgenerator/ . The website https://most.bio.informatik.uni-rostock.de/ provides interactive access to model versions and their evolutionary statistics. CONCLUSION The reuse of models is still impeded by a lack of trust and documentation. A detailed and transparent documentation of all aspects of the model, including its provenance, will improve this situation. Knowledge about a model's provenance can avoid the repetition of mistakes that others already faced. More insights are gained into how the system evolves from initial findings to a profound understanding. We argue that it is the responsibility of the maintainers of model repositories to offer transparent model provenance to their users.
Collapse
Affiliation(s)
- Martin Scharm
- Department of Systems Biology and Bioinformatics, University of Rostock, Rostock, 18051 Germany
| | - Tom Gebhardt
- Department of Systems Biology and Bioinformatics, University of Rostock, Rostock, 18051 Germany
| | - Vasundra Touré
- Department of Biology, Norwegian University of Science and Technology, Trondheim, 7491 Norway
| | - Andrea Bagnacani
- Department of Systems Biology and Bioinformatics, University of Rostock, Rostock, 18051 Germany
| | - Ali Salehzadeh-Yazdi
- Department of Systems Biology and Bioinformatics, University of Rostock, Rostock, 18051 Germany
| | - Olaf Wolkenhauer
- Department of Systems Biology and Bioinformatics, University of Rostock, Rostock, 18051 Germany
- Stellenbosch Institute for Advanced Study, Wallenberg Research Centre at Stellenbosch University, Stellenbosch, 7600 South Africa
| | - Dagmar Waltemath
- Department of Systems Biology and Bioinformatics, University of Rostock, Rostock, 18051 Germany
| |
Collapse
|
25
|
Oertel A, Matros A, Hartmann A, Arapitsas P, Dehmer KJ, Martens S, Mock HP. Metabolite profiling of red and blue potatoes revealed cultivar and tissue specific patterns for anthocyanins and other polyphenols. PLANTA 2017; 246:281-297. [PMID: 28664422 DOI: 10.1007/s00425-017-2718-4] [Citation(s) in RCA: 51] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/02/2017] [Accepted: 05/26/2017] [Indexed: 06/07/2023]
Abstract
Metabolite profiling of tuber flesh and peel for selected colored potato varieties revealed cultivar and tissue specific profiles of anthocyanins and other polyphenols with variations in composition and concentration. Starchy tubers of Solanum tuberosum are a staple crop and food in many countries. Among cultivated potato varieties a huge biodiversity exists, including an increasing number of red and purple colored cultivars. This coloration relates to the accumulation of anthocyanins and is supposed to offer nutritional benefits possibly associated with the antioxidative capacity of anthocyanins. However, the anthocyanin composition and its relation to the overall polyphenol constitution in colored potato tubers have not been investigated closely. This study focuses on the phytochemical characterization of the phenolic composition of a variety of colored potato tubers, both for peel and flesh tissues. First, liquid chromatography (LC) separation coupled to UV and mass spectrometry (MS) detection of polyphenolic compounds of potato tubers from 57 cultivars was used to assign groups of potato cultivars differing in their anthocyanin and polyphenol profiles. Tissues from 19 selected cultivars were then analyzed by LC separation coupled to multiple reaction monitoring (MRM) to detect quantitative differences in anthocyanin and polyphenol composition. The measured intensities of 21 anthocyanins present in the analyzed potato cultivars and tissues could be correlated with the specific tuber coloration. Besides secondary metabolites well-known for potato tubers, the metabolic profiling led to the detection of two anthocyanins not described for potato tuber previously, which we tentatively annotated as pelargonidin feruloyl-xylosyl-glucosyl-galactoside and cyanidin 3-p-coumaroylrutinoside-5-glucoside. We detected significant correlations between some of the measured metabolites, as for example the negative correlation between the main anthocyanins of red and blue potato cultivars. Mainly hydroxylation and methylation patterns of the B-ring of dihydroflavonols, leading to the formation of specific anthocyanidin backbones, can be assigned to a distinct coloring of the potato cultivars and tuber tissues. However, basically the same glycosylation and acylation reactions occur regardless of the main anthocyanidin precursor present in the respective red and blue/purple tissue. Thus, the different anthocyanin profiles in red and blue potato cultivars likely relate to superior regulation of the expression and activities of hydroxylases and methyltransferases rather than to differences for downstream glycosyl- and acyltransferases. In this regard, the characterized potato cultivars represent a valuable resource for the molecular analysis of the genetic background and the regulation of anthocyanin side chain modification.
Collapse
Affiliation(s)
- Anne Oertel
- TRANSMIT GmbH, Project Division: PlantMetaChem (PMC), Kerkrader Straße 3, 35394, Giessen, Germany
- Department of Physiology and Cell Biology, Leibniz Institute of Plant Genetics and Crop Plant Research (IPK-Gatersleben), Corrensstraße 3, 06466, Stadt Seeland OT Gatersleben, Germany
| | - Andrea Matros
- Department of Physiology and Cell Biology, Leibniz Institute of Plant Genetics and Crop Plant Research (IPK-Gatersleben), Corrensstraße 3, 06466, Stadt Seeland OT Gatersleben, Germany
| | - Anja Hartmann
- Department of Physiology and Cell Biology, Leibniz Institute of Plant Genetics and Crop Plant Research (IPK-Gatersleben), Corrensstraße 3, 06466, Stadt Seeland OT Gatersleben, Germany
| | - Panagiotis Arapitsas
- Department of Food Quality and Nutrition, Edmund Mach Foundation, Research and Innovation Centre, Via E. Mach 1, 38010, San Michele all'Adige (TN), Italy
| | - Klaus J Dehmer
- Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), Genebank Department/GLKS, Parkweg 3a, 18190, Gross Luesewitz, Germany
| | - Stefan Martens
- TRANSMIT GmbH, Project Division: PlantMetaChem (PMC), Kerkrader Straße 3, 35394, Giessen, Germany
- Department of Food Quality and Nutrition, Edmund Mach Foundation, Research and Innovation Centre, Via E. Mach 1, 38010, San Michele all'Adige (TN), Italy
| | - Hans-Peter Mock
- Department of Physiology and Cell Biology, Leibniz Institute of Plant Genetics and Crop Plant Research (IPK-Gatersleben), Corrensstraße 3, 06466, Stadt Seeland OT Gatersleben, Germany.
| |
Collapse
|
26
|
Ribeiro HAL, Maioli TU, de Freitas LM, Tieri P, Castiglione F. Modeling Immune Response to Leishmania Species Indicates Adenosine As an Important Inhibitor of Th-Cell Activation. Front Cell Infect Microbiol 2017; 7:309. [PMID: 28775959 PMCID: PMC5517480 DOI: 10.3389/fcimb.2017.00309] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2016] [Accepted: 06/22/2017] [Indexed: 12/20/2022] Open
Abstract
Infection by Leishmania protozoan parasites can cause a variety of disease outcomes in humans and other mammals, from single self-healing cutaneous lesions to a visceral dissemination of the parasite. The correlation between chronic lesions and ecto-nucleotidase enzymes activity on the surface of the parasite is addressed here using damage caused in epithelial cells by nitric oxide. In order to explore the role of purinergic metabolism in lesion formation and the outcome of the infection, we implemented a cellular automata/lattice gas model involving major immune characters (Th1 and Th2 cells, IFN-γ, IL-4, IL-12, adenosine-Ado-, NO) and parasite players for the dynamic analysis of the disease progress. The model were analyzed using partial ranking correlation coefficient (PRCC) to indicate the components that most influence the disease progression. Results show that low Ado inhibition rate over Th-cells is shared by L. major and L. braziliensis, while in L. amazonensis infection the Ado inhibition rate over Th-cells reaches 30%. IL-4 inhibition rate over Th-cell priming to Th1 independent of IL-12 are exclusive of L. major. The lesion size and progression showed agreement with published biological data and the model was able to simulate cutaneous leishmaniasis outcomes. The sensitivity analysis suggested that Ado inhibition rate over Th-cells followed by Leishmania survival probability were the most important characteristics of the process, with PRCC of 0.89 and 0.77 respectively. The simulations also showed a non-linear relationship between Ado inhibition rate over Th-cells and lesion size measured as number of dead epithelial cells. In conclusion, this model can be a useful tool for the quantitative understanding of the immune response in leishmaniasis.
Collapse
Affiliation(s)
- Henrique A. L. Ribeiro
- Consiglio Nazionale delle Ricerche, Istituto per le Applicazioni del CalcoloRome, Italy
- Departamento de Nutrição, Universidade Federal de Minas GeraisBelo Horizonte, Brazil
| | - Tatiani U. Maioli
- Departamento de Nutrição, Universidade Federal de Minas GeraisBelo Horizonte, Brazil
| | | | - Paolo Tieri
- Consiglio Nazionale delle Ricerche, Istituto per le Applicazioni del CalcoloRome, Italy
| | - Filippo Castiglione
- Consiglio Nazionale delle Ricerche, Istituto per le Applicazioni del CalcoloRome, Italy
| |
Collapse
|
27
|
van der Ark KCH, van Heck RGA, Martins Dos Santos VAP, Belzer C, de Vos WM. More than just a gut feeling: constraint-based genome-scale metabolic models for predicting functions of human intestinal microbes. MICROBIOME 2017; 5:78. [PMID: 28705224 PMCID: PMC5512848 DOI: 10.1186/s40168-017-0299-x] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/13/2016] [Accepted: 07/05/2017] [Indexed: 05/14/2023]
Abstract
The human gut is colonized with a myriad of microbes, with substantial interpersonal variation. This complex ecosystem is an integral part of the gastrointestinal tract and plays a major role in the maintenance of homeostasis. Its dysfunction has been correlated to a wide array of diseases, but the understanding of causal mechanisms is hampered by the limited amount of cultured microbes, poor understanding of phenotypes, and the limited knowledge about interspecies interactions. Genome-scale metabolic models (GEMs) have been used in many different fields, ranging from metabolic engineering to the prediction of interspecies interactions. We provide showcase examples for the application of GEMs for gut microbes and focus on (i) the prediction of minimal, synthetic, or defined media; (ii) the prediction of possible functions and phenotypes; and (iii) the prediction of interspecies interactions. All three applications are key in understanding the role of individual species in the gut ecosystem as well as the role of the microbiota as a whole. Using GEMs in the described fashions has led to designs of minimal growth media, an increased understanding of microbial phenotypes and their influence on the host immune system, and dietary interventions to improve human health. Ultimately, an increased understanding of the gut ecosystem will enable targeted interventions in gut microbial composition to restore homeostasis and appropriate host-microbe crosstalk.
Collapse
Affiliation(s)
- Kees C H van der Ark
- Laboratory of Microbiology, Wageningen University, Stippeneng 4, 6708 WE, Wageningen, The Netherlands
| | - Ruben G A van Heck
- Laboratory of Systems and Synthetic Biology, Wageningen University, Stippeneng 4, 6708 WE, Wageningen, The Netherlands
| | - Vitor A P Martins Dos Santos
- Laboratory of Systems and Synthetic Biology, Wageningen University, Stippeneng 4, 6708 WE, Wageningen, The Netherlands
- LifeGlimmer GmbH, Markelstrasse 38, 12163, Berlin, Germany
| | - Clara Belzer
- Laboratory of Microbiology, Wageningen University, Stippeneng 4, 6708 WE, Wageningen, The Netherlands
| | - Willem M de Vos
- Laboratory of Microbiology, Wageningen University, Stippeneng 4, 6708 WE, Wageningen, The Netherlands.
- RPU Immunobiology, Department of Bacteriology and Immunology, University of Helsinki, Haartmanikatu 4, 002940, Helsinki, Finland.
| |
Collapse
|
28
|
External Noise and External Signal Induced Transition of Gene Switch and Coherence Resonance in the Genetic Regulatory System. Acta Biotheor 2017; 65:135-150. [PMID: 28315023 DOI: 10.1007/s10441-017-9307-6] [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: 10/14/2015] [Accepted: 03/10/2017] [Indexed: 10/19/2022]
Abstract
The transition of gene switch induced by external noises (multiplicative external noise and additive external noise) and external signals is investigated in the genetic regulatory system. Results show that the state-to-state transition of gene switch as well as resonant behaviors, such as the explicit coherence resonance (ECR), implicit coherence resonance (ICR) and control parameter coherence biresonance (CPCBR), can appear when noises are injected into the genetic regulatory system. The ECR is increased with the increase of the control parameter value when starting from the supercritical Hopf bifurcation parameter point, and there exists a critical control parameter value for the occurrence of ECR. However, the ICR is decreased as the control parameter value is increased when starting from the subcritical Hopf bifurcation point. In particular, the coherence of ECR is higher and more sensitive to noise than that of ICR. When an external signal is introduced into the system, the enhancement or suppression of the CPCBR and the number of peaks strongly depend on the frequency and amplitude of the external signal. Furthermore, the gene regulation system can selectively enhance or decrease the noise-induced oscillation signals at preferred frequency and amplitude of an external signal.
Collapse
|
29
|
Penas DR, González P, Egea JA, Doallo R, Banga JR. Parameter estimation in large-scale systems biology models: a parallel and self-adaptive cooperative strategy. BMC Bioinformatics 2017; 18:52. [PMID: 28109249 PMCID: PMC5251293 DOI: 10.1186/s12859-016-1452-4] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2016] [Accepted: 12/24/2016] [Indexed: 12/02/2022] Open
Abstract
Background The development of large-scale kinetic models is one of the current key issues in computational systems biology and bioinformatics. Here we consider the problem of parameter estimation in nonlinear dynamic models. Global optimization methods can be used to solve this type of problems but the associated computational cost is very large. Moreover, many of these methods need the tuning of a number of adjustable search parameters, requiring a number of initial exploratory runs and therefore further increasing the computation times. Here we present a novel parallel method, self-adaptive cooperative enhanced scatter search (saCeSS), to accelerate the solution of this class of problems. The method is based on the scatter search optimization metaheuristic and incorporates several key new mechanisms: (i) asynchronous cooperation between parallel processes, (ii) coarse and fine-grained parallelism, and (iii) self-tuning strategies. Results The performance and robustness of saCeSS is illustrated by solving a set of challenging parameter estimation problems, including medium and large-scale kinetic models of the bacterium E. coli, bakerés yeast S. cerevisiae, the vinegar fly D. melanogaster, Chinese Hamster Ovary cells, and a generic signal transduction network. The results consistently show that saCeSS is a robust and efficient method, allowing very significant reduction of computation times with respect to several previous state of the art methods (from days to minutes, in several cases) even when only a small number of processors is used. Conclusions The new parallel cooperative method presented here allows the solution of medium and large scale parameter estimation problems in reasonable computation times and with small hardware requirements. Further, the method includes self-tuning mechanisms which facilitate its use by non-experts. We believe that this new method can play a key role in the development of large-scale and even whole-cell dynamic models. Electronic supplementary material The online version of this article (doi:10.1186/s12859-016-1452-4) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- David R Penas
- BioProcess Engineering Group, IIM-CSIC, Eduardo Cabello 6, Vigo, 36208, Spain
| | - Patricia González
- Computer Architecture Group, Universidade da Coruña, Campus de Elviña s/n, Coruña, 15071 A, Spain
| | - Jose A Egea
- Department of Applied Mathematics and Statistics, Universidad Politécnica de Cartagena, c/ Dr. Fleming s/n, Cartagena, 30202, Spain
| | - Ramón Doallo
- Computer Architecture Group, Universidade da Coruña, Campus de Elviña s/n, Coruña, 15071 A, Spain
| | - Julio R Banga
- BioProcess Engineering Group, IIM-CSIC, Eduardo Cabello 6, Vigo, 36208, Spain.
| |
Collapse
|
30
|
Arceo CPP, Jose EC, Lao AR, Mendoza ER. Reaction networks and kinetics of biochemical systems. Math Biosci 2017; 283:13-29. [DOI: 10.1016/j.mbs.2016.10.004] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2016] [Revised: 09/29/2016] [Accepted: 10/15/2016] [Indexed: 12/12/2022]
|
31
|
Touré V, Mazein A, Waltemath D, Balaur I, Saqi M, Henkel R, Pellet J, Auffray C. STON: exploring biological pathways using the SBGN standard and graph databases. BMC Bioinformatics 2016; 17:494. [PMID: 27919219 PMCID: PMC5139139 DOI: 10.1186/s12859-016-1394-x] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2016] [Accepted: 11/29/2016] [Indexed: 01/16/2023] Open
Abstract
Background When modeling in Systems Biology and Systems Medicine, the data is often extensive, complex and heterogeneous. Graphs are a natural way of representing biological networks. Graph databases enable efficient storage and processing of the encoded biological relationships. They furthermore support queries on the structure of biological networks. Results We present the Java-based framework STON (SBGN TO Neo4j). STON imports and translates metabolic, signalling and gene regulatory pathways represented in the Systems Biology Graphical Notation into a graph-oriented format compatible with the Neo4j graph database. Conclusion STON exploits the power of graph databases to store and query complex biological pathways. This advances the possibility of: i) identifying subnetworks in a given pathway; ii) linking networks across different levels of granularity to address difficulties related to incomplete knowledge representation at single level; and iii) identifying common patterns between pathways in the database. Electronic supplementary material The online version of this article (doi:10.1186/s12859-016-1394-x) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Vasundra Touré
- Department of Systems Biology and Bioinformatics, University of Rostock, Rostock, 18051, Germany. .,European Institute for Systems Biology and Medicine (EISBM), CIRI UMR 5308, CNRS-ENS-UCBL-INSERM, Université de Lyon, 50 Avenue Tony Garnier, Lyon, 69007, France.
| | - Alexander Mazein
- European Institute for Systems Biology and Medicine (EISBM), CIRI UMR 5308, CNRS-ENS-UCBL-INSERM, Université de Lyon, 50 Avenue Tony Garnier, Lyon, 69007, France
| | - Dagmar Waltemath
- Department of Systems Biology and Bioinformatics, University of Rostock, Rostock, 18051, Germany
| | - Irina Balaur
- European Institute for Systems Biology and Medicine (EISBM), CIRI UMR 5308, CNRS-ENS-UCBL-INSERM, Université de Lyon, 50 Avenue Tony Garnier, Lyon, 69007, France
| | - Mansoor Saqi
- European Institute for Systems Biology and Medicine (EISBM), CIRI UMR 5308, CNRS-ENS-UCBL-INSERM, Université de Lyon, 50 Avenue Tony Garnier, Lyon, 69007, France
| | - Ron Henkel
- Scientific Databases and Visualization, Heidelberg Institute for Theoretical Studies, Heidelberg, Germany.,Department of Business Information Systems, University of Rostock, Rostock, 18051, Germany
| | - Johann Pellet
- European Institute for Systems Biology and Medicine (EISBM), CIRI UMR 5308, CNRS-ENS-UCBL-INSERM, Université de Lyon, 50 Avenue Tony Garnier, Lyon, 69007, France
| | - Charles Auffray
- European Institute for Systems Biology and Medicine (EISBM), CIRI UMR 5308, CNRS-ENS-UCBL-INSERM, Université de Lyon, 50 Avenue Tony Garnier, Lyon, 69007, France
| |
Collapse
|
32
|
Big data mining powers fungal research: recent advances in fission yeast systems biology approaches. Curr Genet 2016; 63:427-433. [PMID: 27730285 DOI: 10.1007/s00294-016-0657-4] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2016] [Revised: 10/04/2016] [Accepted: 10/05/2016] [Indexed: 01/05/2023]
Abstract
Biology research has entered into big data era. Systems biology approaches therefore become the powerful tools to obtain the whole landscape of how cell separate, grow, and resist the stresses. Fission yeast Schizosaccharomyces pombe is wonderful unicellular eukaryote model, especially studying its division and metabolism can facilitate to understanding the molecular mechanism of cancer and discovering anticancer agents. In this perspective, we discuss the recent advanced fission yeast systems biology tools, mainly focus on metabolomics profiling and metabolic modeling, protein-protein interactome and genetic interaction network, DNA sequencing and applications, and high-throughput phenotypic screening. We therefore hope this review can be useful for interested fungal researchers as well as bioformaticians.
Collapse
|
33
|
Sánchez BJ, Nielsen J. Genome scale models of yeast: towards standardized evaluation and consistent omic integration. Integr Biol (Camb) 2016; 7:846-58. [PMID: 26079294 DOI: 10.1039/c5ib00083a] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Genome scale models (GEMs) have enabled remarkable advances in systems biology, acting as functional databases of metabolism, and as scaffolds for the contextualization of high-throughput data. In the case of Saccharomyces cerevisiae (budding yeast), several GEMs have been published and are currently used for metabolic engineering and elucidating biological interactions. Here we review the history of yeast's GEMs, focusing on recent developments. We study how these models are typically evaluated, using both descriptive and predictive metrics. Additionally, we analyze the different ways in which all levels of omics data (from gene expression to flux) have been integrated in yeast GEMs. Relevant conclusions and current challenges for both GEM evaluation and omic integration are highlighted.
Collapse
Affiliation(s)
- Benjamín J Sánchez
- Department of Biology and Biological Engineering, Chalmers University of Technology, SE41296 Gothenburg, Sweden.
| | | |
Collapse
|
34
|
Sergushichev AA, Loboda AA, Jha AK, Vincent EE, Driggers EM, Jones RG, Pearce EJ, Artyomov MN. GAM: a web-service for integrated transcriptional and metabolic network analysis. Nucleic Acids Res 2016; 44:W194-200. [PMID: 27098040 PMCID: PMC4987878 DOI: 10.1093/nar/gkw266] [Citation(s) in RCA: 65] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2016] [Accepted: 04/04/2016] [Indexed: 01/14/2023] Open
Abstract
Novel techniques for high-throughput steady-state metabolomic profiling yield information about changes of nearly thousands of metabolites. Such metabolomic profiles, when analyzed together with transcriptional profiles, can reveal novel insights about underlying biological processes. While a number of conceptual approaches have been developed for data integration, easily accessible tools for integrated analysis of mammalian steady-state metabolomic and transcriptional data are lacking. Here we present GAM (‘genes and metabolites’): a web-service for integrated network analysis of transcriptional and steady-state metabolomic data focused on identification of the most changing metabolic subnetworks between two conditions of interest. In the web-service, we have pre-assembled metabolic networks for humans, mice, Arabidopsis and yeast and adapted exact solvers for an optimal subgraph search to work in the context of these metabolic networks. The output is the most regulated metabolic subnetwork of size controlled by false discovery rate parameters. The subnetworks are then visualized online and also can be downloaded in Cytoscape format for subsequent processing. The web-service is available at: https://artyomovlab.wustl.edu/shiny/gam/
Collapse
Affiliation(s)
- Alexey A Sergushichev
- Computer Technologies Department, ITMO University, Saint Petersburg, 197101, Russia Department of Pathology & Immunology, Washington University in St. Louis, St.Louis, MO 63110, USA
| | - Alexander A Loboda
- Computer Technologies Department, ITMO University, Saint Petersburg, 197101, Russia
| | | | - Emma E Vincent
- Goodman Cancer Research Centre, McGill University, Montreal, QC H3A 1A3, Canada Department of Physiology, McGill University, Montreal, QC H3G 1Y6, Canada
| | | | - Russell G Jones
- Goodman Cancer Research Centre, McGill University, Montreal, QC H3A 1A3, Canada Department of Physiology, McGill University, Montreal, QC H3G 1Y6, Canada
| | - Edward J Pearce
- Department of Immunometabolism, Max Planck Institute of Immunobiology and Epigenetics, Freiburg, D-79108, Germany
| | - Maxim N Artyomov
- Computer Technologies Department, ITMO University, Saint Petersburg, 197101, Russia Department of Pathology & Immunology, Washington University in St. Louis, St.Louis, MO 63110, USA
| |
Collapse
|
35
|
Lee S. Systems Biology - A Pivotal Research Methodology for Understanding the Mechanisms of Traditional Medicine. J Pharmacopuncture 2015; 18:11-8. [PMID: 26388998 PMCID: PMC4573803 DOI: 10.3831/kpi.2015.18.020] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2015] [Accepted: 08/31/2015] [Indexed: 12/14/2022] Open
Abstract
Objectives: Systems biology is a novel subject in the field of life science that aims at a systems’ level understanding of biological systems. Because of the significant progress in high-throughput technologies and molecular biology, systems biology occupies an important place in research during the post-genome era. Methods: The characteristics of systems biology and its applicability to traditional medicine research have been discussed from three points of view: data and databases, network analysis and inference, and modeling and systems prediction. Results: The existing databases are mostly associated with medicinal herbs and their activities, but new databases reflecting clinical situations and platforms to extract, visualize and analyze data easily need to be constructed. Network pharmacology is a key element of systems biology, so addressing the multi-component, multi-target aspect of pharmacology is important. Studies of network pharmacology highlight the drug target network and network target. Mathematical modeling and simulation are just in their infancy, but mathematical modeling of dynamic biological processes is a central aspect of systems biology. Computational simulations allow structured systems and their functional properties to be understood and the effects of herbal medicines in clinical situations to be predicted. Conclusion: Systems biology based on a holistic approach is a pivotal research methodology for understanding the mechanisms of traditional medicine. If systems biology is to be incorporated into traditional medicine, computational technologies and holistic insights need to be integrated.
Collapse
Affiliation(s)
- Soojin Lee
- Department of Physiology, College of Korean Medicine, Sangji University, Wonju, Korea
| |
Collapse
|
36
|
Abstract
A mathematical model which predicts the intraerythrocytic stages of Plasmodium falciparum infection was developed using data from malaria-infected mice. Variables selected accounted for levels of healthy red blood cells, merozoite (Plasmodium asexual phase) infected red blood cells, gametocyte (Plasmodium sexual phase) infected red blood cells and a phenomenological variable which accounts for the mean activity of the immune system of the host. The model built was able to reproduce the behavior of three different scenarios of malaria. It predicts the later dynamics of malaria-infected humans well after the first peak of parasitemia, the qualitative response of malaria-infected monkeys to vaccination and the changes observed in malaria-infected mice when they are treated with antimalarial drugs. The mathematical model was used to identify new targets to be focused on drug design. Optimization methodologies were applied to identify five targets for minimizing the parasite load; four of the targets thus identified have never before been taken into account in drug design. The potential targets include: 1) increasing the death rate of the gametocytes, 2) decreasing the invasion rate of the red blood cells by the merozoites, 3) increasing the transformation of merozoites into gametocytes, 4) decreasing the activation of the immune system by the gametocytes, and finally 5) a combination of the previous target with decreasing the recycling rate of the red blood cells. The first target is already used in current therapies, whereas the remainders are proposals for potential new targets. Furthermore, the combined target (the simultaneous decrease of the activation of IS by gRBC and the decrease of the influence of IS on the recycling of hRBC) is interesting, since this combination does not affect the parasite directly. Thus, it is not expected to generate selective pressure on the parasites, which means that it would not produce resistance in Plasmodium.
Collapse
Affiliation(s)
- Guido Santos
- Departamento de Bioquímica y Biología Molecular, Universidad de La Laguna, San Cristóbal de La Laguna, Tenerife, Spain
| | - Néstor V. Torres
- Departamento de Bioquímica y Biología Molecular, Universidad de La Laguna, San Cristóbal de La Laguna, Tenerife, Spain
- Instituto Universitario de Enfermedades Tropicales y Salud Pública de Canarias, Universidad de La Laguna, San Cristóbal de La Laguna. Tenerife, Spain
- * E-mail:
| |
Collapse
|
37
|
Stauch A, Duerr HP, Dujardin JC, Vanaerschot M, Sundar S, Eichner M. Treatment of visceral leishmaniasis: model-based analyses on the spread of antimony-resistant L. donovani in Bihar, India. PLoS Negl Trop Dis 2012; 6:e1973. [PMID: 23285309 PMCID: PMC3527335 DOI: 10.1371/journal.pntd.0001973] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2012] [Accepted: 11/03/2012] [Indexed: 11/29/2022] Open
Abstract
Background Pentavalent antimonials have been the mainstay of antileishmanial therapy for decades, but increasing failure rates under antimonial treatment have challenged further use of these drugs in the Indian subcontinent. Experimental evidence has suggested that parasites which are resistant against antimonials have superior survival skills than sensitive ones even in the absence of antimonial treatment. Methods and Findings We use simulation studies based on a mathematical L. donovani transmission model to identify parameters which can explain why treatment failure rates under antimonial treatment increased up to 65% in Bihar between 1980 and 1997. Model analyses suggest that resistance to treatment alone cannot explain the observed treatment failure rates. We explore two hypotheses referring to an increased fitness of antimony-resistant parasites: the additional fitness is (i) disease-related, by causing more clinical cases (higher pathogenicity) or more severe disease (higher virulence), or (ii) is transmission-related, by increasing the transmissibility from sand flies to humans or vice versa. Conclusions Both hypotheses can potentially explain the Bihar observations. However, increased transmissibility as an explanation appears more plausible because it can occur in the background of asymptomatically transmitted infection whereas disease-related factors would most probably be observable. Irrespective of the cause of fitness, parasites with a higher fitness will finally replace sensitive parasites, even if antimonials are replaced by another drug. The protozoan flagellate Leishmania donovani causes the neglected, life-threatening disease visceral leishmaniasis. Parasites are transmitted from man to man by the bite of the sand fly Phlebotomus argentipes, the vector of the disease. Pentavalent antimonials have been the mainstay of antileishmanial therapy for decades but rapidly increasing failure rates up to 65% observed between 1980 and 1997 in the state of Bihar, India, have challenged further use of these drugs. Comparative in vitro and in vivo experiments indicate that antimony-resistant parasites have a higher fitness than antimony-sensitive ones even in the absence of antimonial treatment. Simulation studies based on a previously published mathematical L. donovani transmission model suggest that resistance to antimonial treatment alone cannot explain the Bihar observations but that resistance together with higher fitness offers the potential to explain the data. After an antimony-resistant parasite with higher fitness has emerged, it will finally replace the antimony-sensitive ones, even in complete absence of antimonial treatment.
Collapse
Affiliation(s)
- Anette Stauch
- Department of Medical Biometry, University of Tuebingen, Tuebingen, Germany.
| | | | | | | | | | | |
Collapse
|