1
|
Turanli B, Gulfidan G, Aydogan OO, Kula C, Selvaraj G, Arga KY. Genome-scale metabolic models in translational medicine: the current status and potential of machine learning in improving the effectiveness of the models. Mol Omics 2024; 20:234-247. [PMID: 38444371 DOI: 10.1039/d3mo00152k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
The genome-scale metabolic model (GEM) has emerged as one of the leading modeling approaches for systems-level metabolic studies and has been widely explored for a broad range of organisms and applications. Owing to the development of genome sequencing technologies and available biochemical data, it is possible to reconstruct GEMs for model and non-model microorganisms as well as for multicellular organisms such as humans and animal models. GEMs will evolve in parallel with the availability of biological data, new mathematical modeling techniques and the development of automated GEM reconstruction tools. The use of high-quality, context-specific GEMs, a subset of the original GEM in which inactive reactions are removed while maintaining metabolic functions in the extracted model, for model organisms along with machine learning (ML) techniques could increase their applications and effectiveness in translational research in the near future. Here, we briefly review the current state of GEMs, discuss the potential contributions of ML approaches for more efficient and frequent application of these models in translational research, and explore the extension of GEMs to integrative cellular models.
Collapse
Affiliation(s)
- Beste Turanli
- Marmara University, Faculty of Engineering, Department of Bioengineering, Istanbul, Turkey.
- Health Biotechnology Joint Research and Application Center of Excellence, Istanbul, Turkey
| | - Gizem Gulfidan
- Marmara University, Faculty of Engineering, Department of Bioengineering, Istanbul, Turkey.
| | - Ozge Onluturk Aydogan
- Marmara University, Faculty of Engineering, Department of Bioengineering, Istanbul, Turkey.
| | - Ceyda Kula
- Marmara University, Faculty of Engineering, Department of Bioengineering, Istanbul, Turkey.
- Health Biotechnology Joint Research and Application Center of Excellence, Istanbul, Turkey
| | - Gurudeeban Selvaraj
- Concordia University, Centre for Research in Molecular Modeling & Department of Chemistry and Biochemistry, Quebec, Canada
- Saveetha Institute of Medical and Technical Sciences (SIMATS), Saveetha Dental College and Hospital, Department of Biomaterials, Bioinformatics Unit, Chennai, India
| | - Kazim Yalcin Arga
- Marmara University, Faculty of Engineering, Department of Bioengineering, Istanbul, Turkey.
- Health Biotechnology Joint Research and Application Center of Excellence, Istanbul, Turkey
- Marmara University, Genetic and Metabolic Diseases Research and Investigation Center, Istanbul, Turkey
| |
Collapse
|
2
|
Chen K, Wang F. Cell-specific genome-scale metabolic modeling of SARS-CoV-2-infected lung to identify antiviral enzymes. FEBS Open Bio 2023; 13:2172-2186. [PMID: 37734920 PMCID: PMC10699103 DOI: 10.1002/2211-5463.13710] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Revised: 08/09/2023] [Accepted: 09/19/2023] [Indexed: 09/23/2023] Open
Abstract
Computational systems biology plays a key role in the discovery of suitable antiviral targets. We designed a cell-specific, constraint-based modeling technique for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected lungs. We used the gene sequence of the alpha variant of SARS-CoV-2 to build a viral biomass reaction (VBR). We also used the mass proportion of lipids between the viral biomass and its host cell to estimate the stoichiometric coefficients of viral lipids in the reaction. We then integrated the VBR, the gene expression of the alpha variant of SARS-CoV-2, and the generic human metabolic network Recon3D to reconstruct a cell-specific genome-scale metabolic model. An antiviral target discovery (AVTD) platform was introduced using this model to identify therapeutic drug targets for combating COVID-19. The AVTD platform not only identified antiviral genes for eliminating viral replication but also predicted side effects of treatments. Our computational results revealed that knocking out dihydroorotate dehydrogenase (DHODH) might reduce the synthesis rate of cytidine-5'-triphosphate and uridine-5'-triphosphate, which terminate the viral building blocks of DNA and RNA for SARS-CoV-2 replication. Our results also indicated that DHODH is a promising antiviral target that causes minor side effects, which is consistent with the results of recent reports. Moreover, we discovered that the genes that participate in the de novo biosynthesis of glycerophospholipids and ceramides become unidentifiable if the VBR does not involve the stoichiometry of lipids.
Collapse
Affiliation(s)
- Ke‐Lin Chen
- Department of Chemical EngineeringNational Chung Cheng UniversityChiayiTaiwan
| | - Feng‐Sheng Wang
- Department of Chemical EngineeringNational Chung Cheng UniversityChiayiTaiwan
| |
Collapse
|
3
|
Chu SW, Wang FS. Fuzzy optimization for identifying antiviral targets for treating SARS-CoV-2 infection in the heart. BMC Bioinformatics 2023; 24:364. [PMID: 37759157 PMCID: PMC10537911 DOI: 10.1186/s12859-023-05487-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2023] [Accepted: 09/18/2023] [Indexed: 09/29/2023] Open
Abstract
In this paper, a fuzzy hierarchical optimization framework is proposed for identifying potential antiviral targets for treating severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in the heart. The proposed framework comprises four objectives for evaluating the elimination of viral biomass growth and the minimization of side effects during treatment. In the application of the framework, Dulbecco's modified eagle medium (DMEM) and Ham's medium were used as uptake nutrients on an antiviral target discovery platform. The prediction results from the framework reveal that most of the antiviral enzymes in the aforementioned media are involved in fatty acid metabolism and amino acid metabolism. However, six enzymes involved in cholesterol biosynthesis in Ham's medium and three enzymes involved in glycolysis in DMEM are unable to eliminate the growth of the SARS-CoV-2 biomass. Three enzymes involved in glycolysis, namely BPGM, GAPDH, and ENO1, in DMEM combine with the supplemental uptake of L-cysteine to increase the cell viability grade and metabolic deviation grade. Moreover, six enzymes involved in cholesterol biosynthesis reduce and fail to reduce viral biomass growth in a culture medium if a cholesterol uptake reaction does not occur and occurs in this medium, respectively.
Collapse
Affiliation(s)
- Sz-Wei Chu
- Department of Chemical Engineering, National Chung Cheng University, Chiayi, 621301, Taiwan
| | - Feng-Sheng Wang
- Department of Chemical Engineering, National Chung Cheng University, Chiayi, 621301, Taiwan.
| |
Collapse
|
4
|
Leonidou N, Renz A, Mostolizadeh R, Dräger A. New workflow predicts drug targets against SARS-CoV-2 via metabolic changes in infected cells. PLoS Comput Biol 2023; 19:e1010903. [PMID: 36952396 PMCID: PMC10035753 DOI: 10.1371/journal.pcbi.1010903] [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/27/2022] [Accepted: 01/30/2023] [Indexed: 03/25/2023] Open
Abstract
COVID-19 is one of the deadliest respiratory diseases, and its emergence caught the pharmaceutical industry off guard. While vaccines have been rapidly developed, treatment options for infected people remain scarce, and COVID-19 poses a substantial global threat. This study presents a novel workflow to predict robust druggable targets against emerging RNA viruses using metabolic networks and information of the viral structure and its genome sequence. For this purpose, we implemented pymCADRE and PREDICATE to create tissue-specific metabolic models, construct viral biomass functions and predict host-based antiviral targets from more than one genome. We observed that pymCADRE reduces the computational time of flux variability analysis for internal optimizations. We applied these tools to create a new metabolic network of primary bronchial epithelial cells infected with SARS-CoV-2 and identified enzymatic reactions with inhibitory effects. The most promising reported targets were from the purine metabolism, while targeting the pyrimidine and carbohydrate metabolisms seemed to be promising approaches to enhance viral inhibition. Finally, we computationally tested the robustness of our targets in all known variants of concern, verifying our targets' inhibitory effects. Since laboratory tests are time-consuming and involve complex readouts to track processes, our workflow focuses on metabolic fluxes within infected cells and is applicable for rapid hypothesis-driven identification of potentially exploitable antivirals concerning various viruses and host cell types.
Collapse
Affiliation(s)
- Nantia Leonidou
- Computational Systems Biology of Infections and Antimicrobial-Resistant Pathogens, Institute for Bioinformatics and Medical Informatics (IBMI), Eberhard Karls University of Tübingen, Tübingen, Germany
- Department of Computer Science, Eberhard Karls University of Tübingen, Tübingen, Germany
- Cluster of Excellence ‘Controlling Microbes to Fight Infections’, Eberhard Karls University of Tübingen, Tübingen, Germany
| | - Alina Renz
- Computational Systems Biology of Infections and Antimicrobial-Resistant Pathogens, Institute for Bioinformatics and Medical Informatics (IBMI), Eberhard Karls University of Tübingen, Tübingen, Germany
- Department of Computer Science, Eberhard Karls University of Tübingen, Tübingen, Germany
- Cluster of Excellence ‘Controlling Microbes to Fight Infections’, Eberhard Karls University of Tübingen, Tübingen, Germany
| | - Reihaneh Mostolizadeh
- Computational Systems Biology of Infections and Antimicrobial-Resistant Pathogens, Institute for Bioinformatics and Medical Informatics (IBMI), Eberhard Karls University of Tübingen, Tübingen, Germany
- Department of Computer Science, Eberhard Karls University of Tübingen, Tübingen, Germany
- Cluster of Excellence ‘Controlling Microbes to Fight Infections’, Eberhard Karls University of Tübingen, Tübingen, Germany
- German Center for Infection Research (DZIF), partner site Tübingen, Germany
| | - Andreas Dräger
- Computational Systems Biology of Infections and Antimicrobial-Resistant Pathogens, Institute for Bioinformatics and Medical Informatics (IBMI), Eberhard Karls University of Tübingen, Tübingen, Germany
- Department of Computer Science, Eberhard Karls University of Tübingen, Tübingen, Germany
- Cluster of Excellence ‘Controlling Microbes to Fight Infections’, Eberhard Karls University of Tübingen, Tübingen, Germany
- German Center for Infection Research (DZIF), partner site Tübingen, Germany
| |
Collapse
|
5
|
Oarga A, Bannerman BP, Júlvez J. CONTRABASS: exploiting flux constraints in genome-scale models for the detection of vulnerabilities. Bioinformatics 2023; 39:7000333. [PMID: 36692133 PMCID: PMC9907045 DOI: 10.1093/bioinformatics/btad053] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2022] [Revised: 11/16/2022] [Accepted: 01/23/2023] [Indexed: 01/25/2023] Open
Abstract
MOTIVATION Despite the fact that antimicrobial resistance is an increasing health concern, the pace of production of new drugs is slow due to the high cost and uncertain success of the process. The development of high-throughput technologies has allowed the integration of biological data into detailed genome-scale models of multiple organisms. Such models can be exploited by means of computational methods to identify system vulnerabilities such as chokepoint reactions and essential reactions. These vulnerabilities are appealing drug targets that can lead to novel drug developments. However, the current approach to compute these vulnerabilities is only based on topological data and ignores the dynamic information of the model. This can lead to misidentified drug targets. RESULTS This work computes flux constraints that are consistent with a certain growth rate of the modelled organism, and integrates the computed flux constraints into the model to improve the detection of vulnerabilities. By exploiting these flux constraints, we are able to obtain a directionality of the reactions of metabolism consistent with a given growth rate of the model, and consequently, a more realistic detection of vulnerabilities can be performed. Several sets of reactions that are system vulnerabilities are defined and the relationships among them are studied. The approach for the detection of these vulnerabilities has been implemented in the Python tool CONTRABASS. Such tool, for which an online web server has also been implemented, computes flux constraints and generates a report with the detected vulnerabilities. AVAILABILITY AND IMPLEMENTATION CONTRABASS is available as an open source Python package at https://github.com/openCONTRABASS/CONTRABASS under GPL-3.0 License. An online web server is available at http://contrabass.unizar.es. SUPPLEMENTARY INFORMATION A glossary of terms are available at Bioinformatics online.
Collapse
Affiliation(s)
- Alexandru Oarga
- Department of Computer Science and Systems Engineering, University of Zaragoza, Zaragoza 50018, Spain
| | - Bridget P Bannerman
- Lucy Cavendish College, Biological Sciences, University of Cambridge, Cambridge CB3 0BU, UK.,Science Resources Foundation, Health Unit, London EC1V 2NX, UK
| | - Jorge Júlvez
- Department of Computer Science and Systems Engineering, University of Zaragoza, Zaragoza 50018, Spain
| |
Collapse
|
6
|
Bannerman BP, Oarga A, Júlvez J. Mycobacterial metabolic model development for drug target identification. GIGABYTE 2023; 2023:gigabyte80. [PMID: 37153490 PMCID: PMC10154535 DOI: 10.46471/gigabyte.80] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2022] [Accepted: 04/24/2023] [Indexed: 05/09/2023] Open
Abstract
Antibiotic resistance is increasing at an alarming rate, and three related mycobacteria are sources of widespread infections in humans. According to the World Health Organization, Mycobacterium leprae, which causes leprosy, is still endemic in tropical countries; Mycobacterium tuberculosis is the second leading infectious killer worldwide after COVID-19; and Mycobacteroides abscessus, a group of non-tuberculous mycobacteria, causes lung infections and other healthcare-associated infections in humans. Due to the rise in resistance to common antibacterial drugs, it is critical that we develop alternatives to traditional treatment procedures. Furthermore, an understanding of the biochemical mechanisms underlying pathogenic evolution is important for the treatment and management of these diseases. In this study, metabolic models have been developed for two bacterial pathogens, M. leprae and My. abscessus, and a new computational tool has been used to identify potential drug targets, which are referred to as bottleneck reactions. The genes, reactions, and pathways in each of these organisms have been highlighted; the potential drug targets can be further explored as broad-spectrum antibacterials and the unique drug targets for each pathogen are significant for precision medicine initiatives. The models and associated datasets described in this paper are available in GigaDB, Biomodels, and PatMeDB repositories.
Collapse
Affiliation(s)
- Bridget P. Bannerman
- Lucy Cavendish College, University of Cambridge, Lady Margaret Rd, Cambridge, CB3 0BU, UK
- Science Resources Foundation, 128 City Road, London, EC1V 2NX, UK
- Corresponding author. E-mail:
| | - Alexandru Oarga
- Department of Computer Science and Systems Engineering, University of Zaragoza, C/María de Luna n 1, 50018, Zaragoza, Spain
| | - Jorge Júlvez
- Department of Computer Science and Systems Engineering, University of Zaragoza, C/María de Luna n 1, 50018, Zaragoza, Spain
| |
Collapse
|
7
|
Multi-omics personalized network analyses highlight progressive disruption of central metabolism associated with COVID-19 severity. Cell Syst 2022; 13:665-681.e4. [PMID: 35933992 PMCID: PMC9263811 DOI: 10.1016/j.cels.2022.06.006] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2021] [Revised: 04/18/2022] [Accepted: 06/27/2022] [Indexed: 01/26/2023]
Abstract
The clinical outcome and disease severity in coronavirus disease 2019 (COVID-19) are heterogeneous, and the progression or fatality of the disease cannot be explained by a single factor like age or comorbidities. In this study, we used system-wide network-based system biology analysis using whole blood RNA sequencing, immunophenotyping by flow cytometry, plasma metabolomics, and single-cell-type metabolomics of monocytes to identify the potential determinants of COVID-19 severity at personalized and group levels. Digital cell quantification and immunophenotyping of the mononuclear phagocytes indicated a substantial role in coordinating the immune cells that mediate COVID-19 severity. Stratum-specific and personalized genome-scale metabolic modeling indicated monocarboxylate transporter family genes (e.g., SLC16A6), nucleoside transporter genes (e.g., SLC29A1), and metabolites such as α-ketoglutarate, succinate, malate, and butyrate could play a crucial role in COVID-19 severity. Metabolic perturbations targeting the central metabolic pathway (TCA cycle) can be an alternate treatment strategy in severe COVID-19.
Collapse
|
8
|
Wang FS, Chen KL, Chu SW. Human/SARS-CoV-2 Genome-Scale Metabolic Modeling to Discover Potential Antiviral Targets for COVID-19. J Taiwan Inst Chem Eng 2022; 133:104273. [PMID: 35186172 PMCID: PMC8843340 DOI: 10.1016/j.jtice.2022.104273] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Revised: 02/06/2022] [Accepted: 02/11/2022] [Indexed: 12/20/2022]
Affiliation(s)
- Feng-Sheng Wang
- Department of Chemical Engineering, National Chung Cheng University, Chiayi 621301, Taiwan
| | - Ke-Lin Chen
- Department of Chemical Engineering, National Chung Cheng University, Chiayi 621301, Taiwan
| | - Sz-Wei Chu
- Department of Chemical Engineering, National Chung Cheng University, Chiayi 621301, Taiwan
| |
Collapse
|