A metabolic network approach for the identification and prioritization of antimicrobial drug targets

Synergizing drug repurposing and target identification for neurodegenerative diseases

Despite dedicated research efforts, the absence of disease-curing remedies for neurodegenerative diseases (NDDs) continues to jeopardize human society and stands as a challenge. Drug repurposing is an attempt to find new functionality of existing drugs and take it as an opportunity to discourse the clinically unmet need to treat neurodegeneration. However, despite applying this approach to rediscover a drug, it can also be used to identify the target on which a drug could work. The primary objective of target identification is to unravel all the possibilities of detecting a new drug or repurposing an existing drug. Lately, scientists and researchers have been focusing on specific genes, a particular site in DNA, a protein, or a molecule that might be involved in the pathogenesis of the disease. However, the new era discusses directing the signaling mechanism involved in the disease progression, where receptors, ion channels, enzymes, and other carrier molecules play a huge role. This review aims to highlight how target identification can expedite the whole process of drug repurposing. Here, we first spot various target-identification methods and drug-repositioning studies, including drug-target and structure-based identification studies. Moreover, we emphasize various drug repurposing approaches in NDDs, namely, experimental-based, mechanism-based, and in silico approaches. Later, we draw attention to validation techniques and stress on drugs that are currently undergoing clinical trials in NDDs. Lastly, we underscore the future perspective of synergizing drug repurposing and target identification in NDDs and present an unresolved question to address the issue.

Cryosphere: a frozen home of microbes and a potential source for drug discovery

The world is concerned about the emergence of pathogens and the occurrence and spread of antibiotic resistance among pathogens. Drug development requires time to combat these issues. Consequently, drug development from natural sources is unavoidable. Cryosphere represents a gigantic source of microbes that could be the bioprospecting source of natural products with unique scaffolds as molecules or drug templates. This review focuses on the novel source of drug discovery and cryospheric environments as a potential source for microbial metabolites having potential medicinal applications. Furthermore, the problems encountered in discovering metabolites from cold-adapted microbes and their resolutions are discussed. By adopting modern practical approaches, the discovery of bioactive compounds might fulfill the demand for new drug development.

Metabolic Connectome and Its Role in the Prediction, Diagnosis, and Treatment of Complex Diseases

The interconnectivity of advanced biological systems is essential for their proper functioning. In modern connectomics, biological entities such as proteins, genes, RNA, DNA, and metabolites are often represented as nodes, while the physical, biochemical, or functional interactions between them are represented as edges. Among these entities, metabolites are particularly significant as they exhibit a closer relationship to an organism's phenotype compared to genes or proteins. Moreover, the metabolome has the ability to amplify small proteomic and transcriptomic changes, even those from minor genomic changes. Metabolic networks, which consist of complex systems comprising hundreds of metabolites and their interactions, play a critical role in biological research by mediating energy conversion and chemical reactions within cells. This review provides an introduction to common metabolic network models and their construction methods. It also explores the diverse applications of metabolic networks in elucidating disease mechanisms, predicting and diagnosing diseases, and facilitating drug development. Additionally, it discusses potential future directions for research in metabolic networks. Ultimately, this review serves as a valuable reference for researchers interested in metabolic network modeling, analysis, and their applications.

ICON-GEMs: integration of co-expression network in genome-scale metabolic models, shedding light through systems biology

BACKGROUND

Flux Balance Analysis (FBA) is a key metabolic modeling method used to simulate cellular metabolism under steady-state conditions. Its simplicity and versatility have led to various strategies incorporating transcriptomic and proteomic data into FBA, successfully predicting flux distribution and phenotypic results. However, despite these advances, the untapped potential lies in leveraging gene-related connections like co-expression patterns for valuable insights.

RESULTS

To fill this gap, we introduce ICON-GEMs, an innovative constraint-based model to incorporate gene co-expression network into the FBA model, facilitating more precise determination of flux distributions and functional pathways. In this study, transcriptomic data from both Escherichia coli and Saccharomyces cerevisiae were integrated into their respective genome-scale metabolic models. A comprehensive gene co-expression network was constructed as a global view of metabolic mechanism of the cell. By leveraging quadratic programming, we maximized the alignment between pairs of reaction fluxes and the correlation of their corresponding genes in the co-expression network. The outcomes notably demonstrated that ICON-GEMs outperformed existing methodologies in predictive accuracy. Flux variabilities over subsystems and functional modules also demonstrate promising results. Furthermore, a comparison involving different types of biological networks, including protein-protein interactions and random networks, reveals insights into the utilization of the co-expression network in genome-scale metabolic engineering.

CONCLUSION

ICON-GEMs introduce an innovative constrained model capable of simultaneous integration of gene co-expression networks, ready for board application across diverse transcriptomic data sets and multiple organisms. It is freely available as open-source at https://github.com/ThummaratPaklao/ICOM-GEMs.git .

A Genome-Scale Metabolic Model of Marine Heterotroph Vibrio splendidus Strain 1A01

While Vibrio splendidus is best known as an opportunistic pathogen in oysters, Vibrio splendidus strain 1A01 was first identified as an early colonizer of synthetic chitin particles incubated in seawater. To gain a better understanding of its metabolism, a genome-scale metabolic model (GSMM) of V. splendidus 1A01 was reconstructed. GSMMs enable us to simulate all metabolic reactions in a bacterial cell using flux balance analysis. A draft model was built using an automated pipeline from BioCyc. Manual curation was then performed based on experimental data, in part by gap-filling metabolic pathways and tailoring the model's biomass reaction to V. splendidus 1A01. The challenges of building a metabolic model for a marine microorganism like V. splendidus 1A01 are described. IMPORTANCE A genome-scale metabolic model of V. splendidus 1A01 was reconstructed in this work. We offer solutions to the technical problems associated with model reconstruction for a marine bacterial strain like V. splendidus 1A01, which arise largely from the high salt concentration found in both seawater and culture media that simulate seawater.

Rapid-SL identifies synthetic lethal sets with an arbitrary cardinality

The multidrug resistance of numerous pathogenic microorganisms is a serious challenge that raises global healthcare concerns. Multi-target medications and combinatorial therapeutics are much more effective than single-target drugs due to their synergistic impact on the systematic activities of microorganisms. Designing efficient combinatorial therapeutics can benefit from identification of synthetic lethals (SLs). An SL is a set of non-essential targets (i.e., reactions or genes) that prevent the proliferation of a microorganism when they are "knocked out" simultaneously. To facilitate the identification of SLs, we introduce Rapid-SL, a new multimodal implementation of the Fast-SL method, using the depth-first search algorithm. The advantages of Rapid-SL over Fast-SL include: (a) the enumeration of all SLs that have an arbitrary cardinality, (b) a shorter runtime due to search space reduction, (c) embarrassingly parallel computations, and (d) the targeted identification of SLs. Targeted identification is important because the enumeration of higher order SLs demands the examination of too many reaction sets. Accordingly, we present specific applications of Rapid-SL for the efficient targeted identification of SLs. In particular, we found up to 67% of all quadruple SLs by investigating about 1% of the search space. Furthermore, 307 sextuples, 476 septuples, and over 9000 octuples are found for Escherichia coli genome-scale model, iAF1260.

Bioinformatics strategies for studying the molecular mechanisms of fungal extracellular vesicles with a focus on infection and immune responses

Fungal extracellular vesicles (EVs) are released during pathogenesis and are found to be an opportunistic infection in most cases. EVs are immunocompetent with their host and have paved the way for new biomedical approaches to drug delivery and the treatment of complex diseases including cancer. With computing and processing advancements, the rise of bioinformatics tools for the evaluation of various parameters involved in fungal EVs has blossomed. In this review, we have complied and explored the bioinformatics tools to analyze the host-pathogen interaction, toxicity, omics and pathogenesis with an array of specific tools that have depicted the ability of EVs as vector/carrier for therapeutic agents and as a potential theme for immunotherapy. We have also discussed the generation and pathways involved in the production, transport, pathogenic action and immunological interactions of EVs in the host system. The incorporation of network pharmacology approaches has been discussed regarding fungal pathogens and their significance in drug discovery. To represent the overview, we have presented and demonstrated an in silico study model to portray the human Cryptococcal interactions.

A comprehensive review of Artificial Intelligence and Network based approaches to drug repurposing in Covid-19

Conventional drug discovery and development is tedious and time-taking process; because of which it has failed to keep the required pace to mitigate threats and cater demands of viral and re-occurring diseases, such as Covid-19. The main reasons of this delay in traditional drug development are: high attrition rates, extensive time requirements, and huge financial investment with significant risk. The effective solution to de novo drug discovery is drug repurposing. Previous studies have shown that the network-based approaches and analysis are versatile platform for repurposing as the network biology is used to model the interactions between variety of biological concepts. Herein, we provide a comprehensive background of machine learning and deep learning in drug repurposing while specifically focusing on the applications of network-based approach to drug repurposing in Covid-19, data sources, and tools used. Furthermore, use of network proximity, network diffusion, and AI on network-based drug repurposing for Covid-19 is well-explained. Finally, limitations of network-based approaches in general and specific to network are stated along with future recommendations for better network-based models.

Interrogation of Essentiality in the Reconstructed Haemophilus influenzae Metabolic Network Identifies Lipid Metabolism Antimicrobial Targets: Preclinical Evaluation of a FabH β-Ketoacyl-ACP Synthase Inhibitor

Expediting drug discovery to fight antibacterial resistance requires holistic approaches at system levels. In this study, we focused on the human-adapted pathogen Haemophilus influenzae, and by constructing a high-quality genome-scale metabolic model, we rationally identified new metabolic drug targets in this organism. Contextualization of available gene essentiality data within in silico predictions identified most genes involved in lipid metabolism as promising targets. We focused on the β-ketoacyl-acyl carrier protein synthase III FabH, responsible for catalyzing the first step in the FASII fatty acid synthesis pathway and feedback inhibition. Docking studies provided a plausible three-dimensional model of FabH in complex with the synthetic inhibitor 1-(5-(2-fluoro-5-(hydroxymethyl)phenyl)pyridin-2-yl)piperidine-4-acetic acid (FabHi). Validating our in silico predictions, FabHi reduced H. influenzae viability in a dose- and strain-dependent manner, and this inhibitory effect was independent of fabH gene expression levels. fabH allelic variation was observed among H. influenzae clinical isolates. Many of these polymorphisms, relevant for stabilization of the dimeric active form of FabH and/or activity, may modulate the inhibitory effect as part of a complex multifactorial process with the overall metabolic context emerging as a key factor tuning FabHi activity. Synergies with antibiotics were not observed and bacteria were not prone to develop resistance. Inhibitor administration during H. influenzae infection on a zebrafish septicemia infection model cleared bacteria without signs of host toxicity. Overall, we highlight the potential of H. influenzae metabolism as a source of drug targets, metabolic models as target-screening tools, and FASII targeting suitability to counteract this bacterial infection. IMPORTANCE Antimicrobial resistance drives the need of synergistically combined powerful computational tools and experimental work to accelerate target identification and drug development. Here, we present a high-quality metabolic model of H. influenzae and show its usefulness both as a computational framework for large experimental data set contextualization and as a tool to discover condition-independent drug targets. We focus on β-ketoacyl-acyl carrier protein synthase III FabH chemical inhibition by using a synthetic molecule with good synthetic and antimicrobial profiles that specifically binds to the active site. The mechanistic complexity of FabH inhibition may go beyond allelic variation, and the strain-dependent effect of the inhibitor tested supports the impact of metabolic context as a key factor driving bacterial cell behavior. Therefore, this study highlights the systematic metabolic evaluation of individual strains through computational frameworks to identify secondary metabolic hubs modulating drug response, which will facilitate establishing synergistic and/or more precise and robust antibacterial treatments.

Antimicrobial Compounds from Microorganisms

Antimicrobial resistance is an exigent public health concern owing to the emergence of novel strains of human resistant pathogens and the concurrent rise in multi-drug resistance. An influx of new antimicrobials is urgently required to improve the treatment outcomes of infectious diseases and save lives. Plant metabolites and bioactive compounds from chemical synthesis have found their efficacy to be dwindling, despite some of them being developed as drugs and used to treat human infections for several decades. Microorganisms are considered untapped reservoirs for promising biomolecules with varying structural and functional antimicrobial activity. The advent of cost-effective and convenient model organisms, state-of-the-art molecular biology, omics technology, and machine learning has enhanced the bioprospecting of novel antimicrobial drugs and the identification of new drug targets. This review summarizes antimicrobial compounds isolated from microorganisms and reports on the modern tools and strategies for exploiting promising antimicrobial drug candidates. The investigation identified a plethora of novel compounds from microbial sources with excellent antimicrobial activity against disease-causing human pathogens. Researchers could maximize the use of novel model systems and advanced biomolecular and computational tools in exploiting lead antimicrobials, consequently ameliorating antimicrobial resistance.

Bacterial metabolism and pathogenesis intimate intertwining: time for metabolic modelling to come into action

We take a snapshot of the recent understanding of bacterial metabolism and the bacterial-host metabolic interplay during infection, and highlight key outcomes and challenges for the practical implementation of bacterial metabolic modelling computational tools in the pathogenesis field.

Uncovering the Role of Metabolism in Oomycete-Host Interactions Using Genome-Scale Metabolic Models

Metabolism is the set of biochemical reactions of an organism that enables it to assimilate nutrients from its environment and to generate building blocks for growth and proliferation. It forms a complex network that is intertwined with the many molecular and cellular processes that take place within cells. Systems biology aims to capture the complexity of cells, organisms, or communities by reconstructing models based on information gathered by high-throughput analyses (omics data) and prior knowledge. One type of model is a genome-scale metabolic model (GEM) that allows studying the distributions of metabolic fluxes, i.e., the "mass-flow" through the network of biochemical reactions. GEMs are nowadays widely applied and have been reconstructed for various microbial pathogens, either in a free-living state or in interaction with their hosts, with the aim to gain insight into mechanisms of pathogenicity. In this review, we first introduce the principles of systems biology and GEMs. We then describe how metabolic modeling can contribute to unraveling microbial pathogenesis and host-pathogen interactions, with a specific focus on oomycete plant pathogens and in particular Phytophthora infestans. Subsequently, we review achievements obtained so far and identify and discuss potential pitfalls of current models. Finally, we propose a workflow for reconstructing high-quality GEMs and elaborate on the resources needed to advance a system biology approach aimed at untangling the intimate interactions between plants and pathogens.

A Systematic Strategy to Find Potential Therapeutic Targets for Pseudomonas aeruginosa Using Integrated Computational Models

Pseudomonas aeruginosa is an opportunistic human pathogen that has been a constant global health problem due to its ability to cause infection at different body sites and its resistance to a broad spectrum of clinically available antibiotics. The World Health Organization classified multidrug-resistant Pseudomonas aeruginosa among the top-ranked organisms that require urgent research and development of effective therapeutic options. Several approaches have been taken to achieve these goals, but they all depend on discovering potential drug targets. The large amount of data obtained from sequencing technologies has been used to create computational models of organisms, which provide a powerful tool for better understanding their biological behavior. In the present work, we applied a method to integrate transcriptome data with genome-scale metabolic networks of Pseudomonas aeruginosa. We submitted both metabolic and integrated models to dynamic simulations and compared their performance with published in vitro growth curves. In addition, we used these models to identify potential therapeutic targets and compared the results to analyze the assumption that computational models enriched with biological measurements can provide more selective and (or) specific predictions. Our results demonstrate that dynamic simulations from integrated models result in more accurate growth curves and flux distribution more coherent with biological observations. Moreover, identifying drug targets from integrated models is more selective as the predicted genes were a subset of those found in the metabolic models. Our analysis resulted in the identification of 26 non-host homologous targets. Among them, we highlighted five top-ranked genes based on lesser conservation with the human microbiome. Overall, some of the genes identified in this work have already been proposed by different approaches and (or) are already investigated as targets to antimicrobial compounds, reinforcing the benefit of using integrated models as a starting point to selecting biologically relevant therapeutic targets.

Drug repurposing: Iron in the fire for older drugs

Repositioning or "repurposing" of existing therapies for indications of alternative disease is an attractive approach that can generate lower costs and require a shorter approval time than developing a de novo drug. The development of experimental drugs is time-consuming, expensive, and limited to a fairly small number of targets. The incorporation of separate and complementary data should be used, as each type of data set exposes a specific feature of organism knowledge Drug repurposing opportunities are often focused on sporadic findings or on time-consuming pre-clinical drug tests which are often not guided by hypothesis. In comparison, repurposing in-silico drugs is a new, hypothesis-driven method that takes advantage of big-data use. Nonetheless, the widespread use of omics technology, enhanced data storage, data sense, machine learning algorithms, and computational modeling all give unparalleled knowledge of the methods of action of biological processes and drugs, providing wide availability, for both disease-related data and drug-related data. This review has taken an in-depth look at the current state, possibilities, and limitations of further progress in the field of drug repositioning.

Antibiotic resistance: Time of synthesis in a post-genomic age

Antibiotic resistance has been highlighted by international organizations, including World Health Organization, World Bank and United Nations, as one of the most relevant global health problems. Classical approaches to study this problem have focused in infected humans, mainly at hospitals. Nevertheless, antibiotic resistance can expand through different ecosystems and geographical allocations, hence constituting a One-Health, Global-Health problem, requiring specific integrative analytic tools. Antibiotic resistance evolution and transmission are multilayer, hierarchically organized processes with several elements (from genes to the whole microbiome) involved. However, their study has been traditionally gene-centric, each element independently studied. The development of robust-economically affordable whole genome sequencing approaches, as well as other -omic techniques as transcriptomics and proteomics, is changing this panorama. These technologies allow the description of a system, either a cell or a microbiome as a whole, overcoming the problems associated with gene-centric approaches. We are currently at the time of combining the information derived from -omic studies to have a more holistic view of the evolution and spread of antibiotic resistance. This synthesis process requires the accurate integration of -omic information into computational models that serve to analyse the causes and the consequences of acquiring AR, fed by curated databases capable of identifying the elements involved in the acquisition of resistance. In this review, we analyse the capacities and drawbacks of the tools that are currently in use for the global analysis of AR, aiming to identify the more useful targets for effective corrective interventions.

Exploring the druggable proteome of Candida species through comprehensive computational analysis

Candida albicans and non-albicans Candida spp. are major cause of systemic mycoses. Antifungal drugs such as azoles and polyenes are not efficient to successfully eradicate Candida infection owing to their fungistatic nature or low bioavailability. Here, we have adopted a comprehensive computational workflow for identification, prioritization and validation of targets from proteomes of Candida albicans and Candida tropicalis. The protocol involves identification of essential drug-target candidates using subtractive genomics, protein-protein interaction network properties and systems biology based methods. The essentiality of the novel metabolic and non-metabolic targets was established by performing in silico gene knockouts, under aerobic as well as anaerobic conditions, and in vitro drug inhibition assays respectively. Deletion of twelve genes that are involved in amino acid, secondary metabolite, and carbon metabolism showed zero growth in metabolic model under simulated conditions. The algorithm, used in this study, can be downloaded from http://pbit.bicnirrh.res.in/offline.php and executed locally.

A Pan-Genome Guided Metabolic Network Reconstruction of Five Propionibacterium Species Reveals Extensive Metabolic Diversity

Propionibacteria have been studied extensively since the early 1930s due to their relevance to industry and importance as human pathogens. Still, their unique metabolism is far from fully understood. This is partly due to their signature high GC content, which has previously hampered the acquisition of quality sequence data, the accurate annotation of the available genomes, and the functional characterization of genes. The recent completion of the genome sequences for several species has led researchers to reassess the taxonomical classification of the genus Propionibacterium, which has been divided into several new genres. Such data also enable a comparative genomic approach to annotation and provide a new opportunity to revisit our understanding of their metabolism. Using pan-genome analysis combined with the reconstruction of the first high-quality Propionibacterium genome-scale metabolic model and a pan-metabolic model of current and former members of the genus Propionibacterium, we demonstrate that despite sharing unique metabolic traits, these organisms have an unexpected diversity in central carbon metabolism and a hidden layer of metabolic complexity. This combined approach gave us new insights into the evolution of Propionibacterium metabolism and led us to propose a novel, putative ferredoxin-linked energy conservation strategy. The pan-genomic approach highlighted key differences in Propionibacterium metabolism that reflect adaptation to their environment. Results were mathematically captured in genome-scale metabolic reconstructions that can be used to further explore metabolism using metabolic modeling techniques. Overall, the data provide a platform to explore Propionibacterium metabolism and a tool for the rational design of strains.

Genome-Scale Metabolic Model of the Human Pathogen Candida albicans: A Promising Platform for Drug Target Prediction

Candida albicans is one of the most impactful fungal pathogens and the most common cause of invasive candidiasis, which is associated with very high mortality rates. With the rise in the frequency of multidrug-resistant clinical isolates, the identification of new drug targets and new drugs is crucial in overcoming the increase in therapeutic failure. In this study, the first validated genome-scale metabolic model for Candida albicans, iRV781, is presented. The model consists of 1221 reactions, 926 metabolites, 781 genes, and four compartments. This model was reconstructed using the open-source software tool merlin 4.0.2. It is provided in the well-established systems biology markup language (SBML) format, thus, being usable in most metabolic engineering platforms, such as OptFlux or COBRA. The model was validated, proving accurate when predicting the capability of utilizing different carbon and nitrogen sources when compared to experimental data. Finally, this genome-scale metabolic reconstruction was tested as a platform for the identification of drug targets, through the comparison between known drug targets and the prediction of gene essentiality in conditions mimicking the human host. Altogether, this model provides a promising platform for global elucidation of the metabolic potential of C. albicans, possibly guiding the identification of new drug targets to tackle human candidiasis.

Novel antimicrobial development using genome-scale metabolic model of Gram-negative pathogens: a review

Antimicrobial resistance (AMR) threatens the effective prevention and treatment of a wide range of infections. Governments around the world are beginning to devote effort for innovative treatment development to treat these resistant bacteria. Systems biology methods have been applied extensively to provide valuable insights into metabolic processes at system level. Genome-scale metabolic models serve as platforms for constraint-based computational techniques which aid in novel drug discovery. Tools for automated reconstruction of metabolic models have been developed to support system level metabolic analysis. We discuss features of such software platforms for potential users to best fit their purpose of research. In this work, we focus to review the development of genome-scale metabolic models of Gram-negative pathogens and also metabolic network approach for identification of antimicrobial drugs targets.

Novel therapeutic strategies for treating Pseudomonas aeruginosa infection

INTRODUCTION

Persistent infections caused by the superbug Pseudomonas aeruginosa and its resistance to multiple antimicrobial agents are huge threats to patients with cystic fibrosis as well as those with compromised immune systems. Multidrug-resistant P. aeruginosa has posed a major challenge to conventional antibiotics and therapeutic approaches, which show limited efficacy and cause serious side effects. The public demand for new antibiotics is enormous; yet, drug development pipelines have started to run dry with limited targets available for inventing new antibacterial drugs. Consequently, it is important to uncover potential therapeutic targets.

AREAS COVERED

The authors review the current state of drug development strategies that are promising in terms of the development of novel and potent drugs to treat P. aeruginosa infection.

EXPERT OPINION

The prevention of P. aeruginosa infection is increasingly challenging. Furthermore, targeting key virulence regulators has great potential for developing novel anti-P. aeruginosa drugs. Additional promising strategies include bacteriophage therapy, immunotherapies, and antimicrobial peptides. Additionally, the authors believe that in the coming years, the overall network of molecular regulatory mechanism of P. aeruginosa virulence will be fully elucidated, which will provide more novel and promising drug targets for treating P. aeruginosa infections.

Nanoinformatics and biomolecular nanomodeling: a novel move en route for effective cancer treatment

Empowering role of nanoinformatics in design and elucidation of nanoparticles for effective cancer treatment has made this field a fascinating area for researchers, inspiring them to enhance up the quality and efficacy of existing anticancer medicines. Theoretical and computational modeling is being seen as a forefront solution for problems related to surface chemistry, optimized geometry, or other properties in nanoparticle designing and drug delivery. The current review aims to acquaint with the insight story of the incubation of in silico tools and techniques in nanotechnology to develop better anticancer nanomedicines. The review also recapitulates the assets and liabilities of this field and present an outline of existing inventiveness and endeavors of nanoinformatics. We propose how nanoinformatics could hasten up the advancements in anticancer nanomedicines through use of computational tools, nanoparticles repositories & various modeling and simulation methods.

Integrating proteomic or transcriptomic data into metabolic models using linear bound flux balance analysis

Motivation

Transcriptomics and proteomics data have been integrated into constraint-based models to influence flux predictions. However, it has been reported recently for Escherichia coli and Saccharomyces cerevisiae, that model predictions from parsimonious flux balance analysis (pFBA), which does not use expression data, are as good or better than predictions from various algorithms that integrate transcriptomics or proteomics data into constraint-based models.

Results

In this paper, we describe a novel constraint-based method called Linear Bound Flux Balance Analysis (LBFBA), which uses expression data (either transcriptomic or proteomic) to predict metabolic fluxes. The method uses expression data to place soft constraints on individual fluxes, which can be violated. Parameters in the soft constraints are first estimated from a training expression and flux dataset before being used to predict fluxes from expression data in other conditions. We applied LBFBA to E.coli and S.cerevisiae datasets and found that LBFBA predictions were more accurate than pFBA predictions, with average normalized errors roughly half of those from pFBA. For the first time, we demonstrate a computational method that integrates expression data into constraint-based models and improves quantitative flux predictions over pFBA.

Availability and implementation

Code is available in the Supplementary data available at Bioinformatics online.

Supplementary information

Supplementary data are available at Bioinformatics online.

Metabolic Model of the Phytophthora infestans-Tomato Interaction Reveals Metabolic Switches during Host Colonization

The oomycete pathogen Phytophthora infestans causes potato and tomato late blight, a disease that is a serious threat to agriculture. P. infestans is a hemibiotrophic pathogen, and during infection, it scavenges nutrients from living host cells for its own proliferation. To date, the nutrient flux from host to pathogen during infection has hardly been studied, and the interlinked metabolisms of the pathogen and host remain poorly understood. Here, we reconstructed an integrated metabolic model of P. infestans and tomato (Solanum lycopersicum) by integrating two previously published models for both species. We used this integrated model to simulate metabolic fluxes from host to pathogen and explored the topology of the model to study the dependencies of the metabolism of P. infestans on that of tomato. This showed, for example, that P. infestans, a thiamine auxotroph, depends on certain metabolic reactions of the tomato thiamine biosynthesis. We also exploited dual-transcriptome data of a time course of a full late blight infection cycle on tomato leaves and integrated the expression of metabolic enzymes in the model. This revealed profound changes in pathogen-host metabolism during infection. As infection progresses, P. infestans performs less de novo synthesis of metabolites and scavenges more metabolites from tomato. This integrated metabolic model for the P. infestans-tomato interaction provides a framework to integrate data and generate hypotheses about in planta nutrition of P. infestans throughout its infection cycle.IMPORTANCE Late blight disease caused by the oomycete pathogen Phytophthora infestans leads to extensive yield losses in tomato and potato cultivation worldwide. To effectively control this pathogen, a thorough understanding of the mechanisms shaping the interaction with its hosts is paramount. While considerable work has focused on exploring host defense mechanisms and identifying P. infestans proteins contributing to virulence and pathogenicity, the nutritional strategies of the pathogen are mostly unresolved. Genome-scale metabolic models (GEMs) can be used to simulate metabolic fluxes and help in unravelling the complex nature of metabolism. We integrated a GEM of tomato with a GEM of P. infestans to simulate the metabolic fluxes that occur during infection. This yields insights into the nutrients that P. infestans obtains during different phases of the infection cycle and helps in generating hypotheses about nutrition in planta.

FindTargetsWEB: A User-Friendly Tool for Identification of Potential Therapeutic Targets in Metabolic Networks of Bacteria

Background: Healthcare-associated infections (HAIs) are a serious public health problem. They can be associated with morbidity and mortality and are responsible for the increase in patient hospitalization. Antimicrobial resistance among pathogens causing HAI has increased at alarming levels. In this paper, a robust method for analyzing genome-scale metabolic networks of bacteria is proposed in order to identify potential therapeutic targets, along with its corresponding web implementation, dubbed FindTargetsWEB. The proposed method assumes that every metabolic network presents fragile genes whose blockade will impair one or more metabolic functions, such as biomass accumulation. FindTargetsWEB automates the process of identification of such fragile genes using flux balance analysis (FBA), flux variability analysis (FVA), extended Systems Biology Markup Language (SBML) file parsing, and queries to three public repositories, i.e., KEGG, UniProt, and DrugBank. The web application was developed in Python using COBRApy and Django.

Results: The proposed method was demonstrated to be robust enough to process even non-curated, incomplete, or imprecise metabolic networks, in addition to integrated host-pathogen models. A list of potential therapeutic targets and their putative inhibitors was generated as a result of the analysis of Pseudomonas aeruginosa metabolic networks available in the literature and a curated version of the metabolic network of a multidrug-resistant P. aeruginosa strain belonging to a clone endemic in Brazil (P. aeruginosa ST277). Genome-scale metabolic networks of other gram-positive and gram-negative bacteria, such as Staphylococcus aureus, Klebsiella pneumoniae, and Haemophilus influenzae, were also analyzed using FindTargetsWEB. Multiple potential targets have been found using the proposed method in all metabolic networks, including some overlapping between two or more pathogens. Among the potential targets, several have been previously reported in the literature as targets for antimicrobial development, and many targets have approved drugs. Despite similarities in the metabolic network structure for closely related bacteria, we show that the method is able to selectively identify targets in pathogenic versus non-pathogenic organisms.

Conclusions: This new computational system can give insights into the identification of new candidate therapeutic targets for pathogenic bacteria and discovery of new antimicrobial drugs through genome-scale metabolic network analysis and heterogeneous data integration, even for non-curated or incomplete networks.

Antibiotic Discovery: Where Have We Come from, Where Do We Go?

Given the increase in antibiotic-resistant bacteria, alongside the alarmingly low rate of newly approved antibiotics for clinical usage, we are on the verge of not having effective treatments for many common infectious diseases. Historically, antibiotic discovery has been crucial in outpacing resistance and success is closely related to systematic procedures—platforms—that have catalyzed the antibiotic golden age, namely the Waksman platform, followed by the platforms of semi-synthesis and fully synthetic antibiotics. Said platforms resulted in the major antibiotic classes: aminoglycosides, amphenicols, ansamycins, beta-lactams, lipopeptides, diaminopyrimidines, fosfomycins, imidazoles, macrolides, oxazolidinones, streptogramins, polymyxins, sulphonamides, glycopeptides, quinolones and tetracyclines. During the genomics era came the target-based platform, mostly considered a failure due to limitations in translating drugs to the clinic. Therefore, cell-based platforms were re-instituted, and are still of the utmost importance in the fight against infectious diseases. Although the antibiotic pipeline is still lackluster, especially of new classes and novel mechanisms of action, in the post-genomic era, there is an increasingly large set of information available on microbial metabolism. The translation of such knowledge into novel platforms will hopefully result in the discovery of new and better therapeutics, which can sway the war on infectious diseases back in our favor.

Reconciling high-throughput gene essentiality data with metabolic network reconstructions

The identification of genes essential for bacterial growth and survival represents a promising strategy for the discovery of antimicrobial targets. Essential genes can be identified on a genome-scale using transposon mutagenesis approaches; however, variability between screens and challenges with interpretation of essentiality data hinder the identification of both condition-independent and condition-dependent essential genes. To illustrate the scope of these challenges, we perform a large-scale comparison of multiple published Pseudomonas aeruginosa gene essentiality datasets, revealing substantial differences between the screens. We then contextualize essentiality using genome-scale metabolic network reconstructions and demonstrate the utility of this approach in providing functional explanations for essentiality and reconciling differences between screens. Genome-scale metabolic network reconstructions also enable a high-throughput, quantitative analysis to assess the impact of media conditions on the identification of condition-independent essential genes. Our computational model-driven analysis provides mechanistic insight into essentiality and contributes novel insights for design of future gene essentiality screens and the identification of core metabolic processes.

With the rise of antibiotic resistance, there is a growing need to discover new therapeutic targets to treat bacterial infections. One attractive strategy is to target genes that are essential for growth and survival. Essential genes can be identified with transposon mutagenesis approaches; however, variability between screens and challenges with interpretation of essentiality data hinder the identification and analysis of essential genes. We performed a large-scale comparison of multiple gene essentiality screens of the microbial pathogen Pseudomonas aeruginosa. We implemented a computational model-driven approach to provide functional explanations for essentiality and reconcile differences between screens. The integration of computational modeling with high-throughput experimental screens may enable the identification of drug targets with high-confidence and provide greater understanding for the development of novel therapeutic strategies.

A review of network-based approaches to drug repositioning

Experimental drug development is time-consuming, expensive and limited to a relatively small number of targets. However, recent studies show that repositioning of existing drugs can function more efficiently than de novo experimental drug development to minimize costs and risks. Previous studies have proven that network analysis is a versatile platform for this purpose, as the biological networks are used to model interactions between many different biological concepts. The present study is an attempt to review network-based methods in predicting drug targets for drug repositioning. For each method, the preferred type of data set is described, and their advantages and limitations are discussed. For each method, we seek to provide a brief description, as well as an evaluation based on its performance metrics.We conclude that integrating distinct and complementary data should be used because each type of data set reveals a unique aspect of information about an organism. We also suggest that applying a standard set of evaluation metrics and data sets would be essential in this fast-growing research domain.

Metabolic Responses to Polymyxin Treatment in Acinetobacter baumannii ATCC 19606: Integrating Transcriptomics and Metabolomics with Genome-Scale Metabolic Modeling

Multidrug-resistant (MDR) Acinetobacter baumannii has emerged as a very problematic pathogen over the past decades, with a high incidence in nosocomial infections. Discovered in the late 1940s but abandoned in the 1970s, polymyxins (i.e., polymyxin B and colistin) have been revived as the last-line therapy against Gram-negative "superbugs," including MDR A. baumannii. Worryingly, resistance to polymyxins in A. baumannii has been increasingly reported, urging the development of novel antimicrobial therapies to rescue this last-line class of antibiotics. In the present study, we integrated genome-scale metabolic modeling with multiomics data to elucidate the mechanisms of cellular responses to colistin treatment in A. baumannii. A genome-scale metabolic model, iATCC19606, was constructed for strain ATCC 19606 based on the literature and genome annotation, containing 897 genes, 1,270 reactions, and 1,180 metabolites. After extensive curation, prediction of growth on 190 carbon sources using iATCC19606 achieved an overall accuracy of 84.3% compared to Biolog experimental results. Prediction of gene essentiality reached a high accuracy of 86.1% and 82.7% compared to two transposon mutant libraries of AB5075 and ATCC 17978, respectively. Further integrative modeling with our correlative transcriptomics and metabolomics data deciphered the complex regulation on metabolic responses to colistin treatment, including (i) upregulated fluxes through gluconeogenesis, the pentose phosphate pathway, and amino acid and nucleotide biosynthesis; (ii) downregulated TCA cycle and peptidoglycan and lipopolysaccharide biogenesis; and (iii) altered fluxes over respiratory chain. Our results elucidated the interplay of multiple metabolic pathways under colistin treatment in A. baumannii and provide key mechanistic insights into optimizing polymyxin combination therapy. IMPORTANCE Combating antimicrobial resistance has been highlighted as a critical global health priority. Due to the drying drug discovery pipeline, polymyxins have been employed as the last-line therapy against Gram-negative "superbugs"; however, the detailed mechanisms of antibacterial killing remain largely unclear, hampering the improvement of polymyxin therapy. Our integrative modeling using the constructed genome-scale metabolic model iATCC19606 and the correlative multiomics data provide the fundamental understanding of the complex metabolic responses to polymyxin treatment in A. baumannii at the systems level. The model iATCC19606 may have a significant potential in antimicrobial systems pharmacology research in A. baumannii.

Methods for automated genome-scale metabolic model reconstruction

In the era of next-generation sequencing and ubiquitous assembly and binning of metagenomes, new putative genome sequences are being produced from isolate and microbiome samples at ever-increasing rates. Genome-scale metabolic models have enormous utility for supporting the analysis and predictive characterization of these genomes based on sequence data. As a result, tools for rapid automated reconstruction of metabolic models are becoming critically important for supporting the analysis of new genome sequences. Many tools and algorithms have now emerged to support rapid model reconstruction and analysis. Here, we are comparing and contrasting the capabilities and output of a variety of these tools, including ModelSEED, Raven Toolbox, PathwayTools, SuBliMinal Toolbox and merlin.

Genome-wide characterization of Phytophthora infestans metabolism: a systems biology approach

Genome-scale metabolic models (GEMs) provide a functional view of the complex network of biochemical reactions in the living cell. Initially mainly applied to reconstruct the metabolism of model organisms, the availability of increasingly sophisticated reconstruction methods and more extensive biochemical databases now make it possible to reconstruct GEMs for less well-characterized organisms, and have the potential to unravel the metabolism in pathogen-host systems. Here, we present a GEM for the oomycete plant pathogen Phytophthora infestans as a first step towards an integrative model with its host. We predict the biochemical reactions in different cellular compartments and investigate the gene-protein-reaction associations in this model to obtain an impression of the biochemical capabilities of P. infestans. Furthermore, we generate life stage-specific models to place the transcriptomic changes of the genes encoding metabolic enzymes into a functional context. In sporangia and zoospores, there is an overall down-regulation, most strikingly reflected in the fatty acid biosynthesis pathway. To investigate the robustness of the GEM, we simulate gene deletions to predict which enzymes are essential for in vitro growth. This model is an essential first step towards an understanding of P. infestans and its interactions with plants as a system, which will help to formulate new hypotheses on infection mechanisms and disease prevention.

Biomedical applications of genome-scale metabolic network reconstructions of human pathogens

The growing global threat of antibiotic resistant human pathogens has coincided with improved methods for developing and using genome-scale metabolic network reconstructions. Consequently, there has been an increase in the number of high-quality reconstructions of relevant human and zoonotic pathogens. Novel biomedical applications of pathogen reconstructions focus on three key aspects of pathogen behavior: the evolution of antibiotic resistance, virulence factor production, and host-pathogen interactions. New methods using these reconstructions aim to improve understanding of microbe pathogenicity and guide the development of new therapeutic strategies. This review summarizes the latest ways that genome-scale metabolic network reconstructions have been used to study human pathogens and suggests future applications with the potential to mitigate infectious disease.

Genome-scale metabolic models as platforms for identification of novel genes as antimicrobial drug targets

The growing number of multidrug-resistant pathogenic bacteria is becoming a world leading challenge for the scientific community and for public health. However, advances in high-throughput technologies and whole-genome sequencing of bacterial pathogens make the construction of bacterial genome-scale metabolic models (GEMs) increasingly realistic. The use of GEMs as an alternative platforms will expedite identification of novel unconditionally essential genes and enzymes of target organisms with existing and forthcoming GEMs. This approach will follow the existing protocol for construction of high-quality GEMs, which could ultimately reduce the time, cost and labor-intensive processes involved in identification of novel antimicrobial drug targets in drug discovery pipelines. We discuss the current impact of existing GEMs of selected multidrug-resistant pathogenic bacteria for identification of novel antimicrobial drug targets and the challenges of closing the gap between genome-scale metabolic modeling and conventional experimental trial-and-error approaches in drug discovery pipelines.

Control of primary metabolism by a virulence regulatory network promotes robustness in a plant pathogen

Robustness is a key system-level property of living organisms to maintain their functions while tolerating perturbations. We investigate here how a regulatory network controlling multiple virulence factors impacts phenotypic robustness of a bacterial plant pathogen. We reconstruct a cell-scale model of Ralstonia solanacearum connecting a genome-scale metabolic network, a virulence macromolecule network, and a virulence regulatory network, which includes 63 regulatory components. We develop in silico methods to quantify phenotypic robustness under a broad set of conditions in high-throughput simulation analyses. This approach reveals that the virulence regulatory network exerts a control of the primary metabolism to promote robustness upon infection. The virulence regulatory network plugs into the primary metabolism mainly through the control of genes likely acquired via horizontal gene transfer, which results in a functional overlay with ancestral genes. These results support the view that robustness may be a selected trait that promotes pathogenic fitness upon infection.

Analysis of Piscirickettsia salmonis Metabolism Using Genome-Scale Reconstruction, Modeling, and Testing

Piscirickettsia salmonis is an intracellular bacterial fish pathogen that causes piscirickettsiosis, a disease with highly adverse impact in the Chilean salmon farming industry. The development of effective treatment and control methods for piscireckttsiosis is still a challenge. To meet it the number of studies on P. salmonis has grown in the last couple of years but many aspects of the pathogen's biology are still poorly understood. Studies on its metabolism are scarce and only recently a metabolic model for reference strain LF-89 was developed. We present a new genome-scale model for P. salmonis LF-89 with more than twice as many genes as in the previous model and incorporating specific elements of the fish pathogen metabolism. Comparative analysis with models of different bacterial pathogens revealed a lower flexibility in P. salmonis metabolic network. Through constraint-based analysis, we determined essential metabolites required for its growth and showed that it can benefit from different carbon sources tested experimentally in new defined media. We also built an additional model for strain A1-15972, and together with an analysis of P. salmonis pangenome, we identified metabolic features that differentiate two main species clades. Both models constitute a knowledge-base for P. salmonis metabolism and can be used to guide the efficient culture of the pathogen and the identification of specific drug targets.

Reconstruction of the metabolic network of Pseudomonas aeruginosa to interrogate virulence factor synthesis

Virulence-linked pathways in opportunistic pathogens are putative therapeutic targets that may be associated with less potential for resistance than targets in growth-essential pathways. However, efficacy of virulence-linked targets may be affected by the contribution of virulence-related genes to metabolism. We evaluate the complex interrelationships between growth and virulence-linked pathways using a genome-scale metabolic network reconstruction of Pseudomonas aeruginosa strain PA14 and an updated, expanded reconstruction of P. aeruginosa strain PAO1. The PA14 reconstruction accounts for the activity of 112 virulence-linked genes and virulence factor synthesis pathways that produce 17 unique compounds. We integrate eight published genome-scale mutant screens to validate gene essentiality predictions in rich media, contextualize intra-screen discrepancies and evaluate virulence-linked gene distribution across essentiality datasets. Computational screening further elucidates interconnectivity between inhibition of virulence factor synthesis and growth. Successful validation of selected gene perturbations using PA14 transposon mutants demonstrates the utility of model-driven screening of therapeutic targets.

Exploration of Nitrate Reductase Metabolic Pathway in Corynebacterium pseudotuberculosis

Based on the ability of nitrate reductase synthesis, Corynebacterium pseudotuberculosis is classified into two biovars: Ovis and Equi. Due to the presence of nitrate reductase, the Equi biovar can survive in absence of oxygen. On the other hand, Ovis biovar that does not have nitrate reductase is able to adapt to various ecological niches and can grow on certain carbon sources. Apart from these two biovars, some other strains are also able to carry out the reduction of nitrate. The enzymes that are involved in electron transport chain are also identified by in silico methods. Findings about pathogen metabolism can contribute to the identification of relationship between nitrate reductase and the C. pseudotuberculosis pathogenicity, virulence factors, and discovery of drug targets.

Multiscale Model of Mycobacterium tuberculosis Infection Maps Metabolite and Gene Perturbations to Granuloma Sterilization Predictions

Granulomas are a hallmark of tuberculosis. Inside granulomas, the pathogen Mycobacterium tuberculosis may enter a metabolically inactive state that is less susceptible to antibiotics. Understanding M. tuberculosis metabolism within granulomas could contribute to reducing the lengthy treatment required for tuberculosis and provide additional targets for new drugs. Two key adaptations of M. tuberculosis are a nonreplicating phenotype and accumulation of lipid inclusions in response to hypoxic conditions. To explore how these adaptations influence granuloma-scale outcomes in vivo, we present a multiscale in silico model of granuloma formation in tuberculosis. The model comprises host immunity, M. tuberculosis metabolism, M. tuberculosis growth adaptation to hypoxia, and nutrient diffusion. We calibrated our model to in vivo data from nonhuman primates and rabbits and apply the model to predict M. tuberculosis population dynamics and heterogeneity within granulomas. We found that bacterial populations are highly dynamic throughout infection in response to changing oxygen levels and host immunity pressures. Our results indicate that a nonreplicating phenotype, but not lipid inclusion formation, is important for long-term M. tuberculosis survival in granulomas. We used virtual M. tuberculosis knockouts to predict the impact of both metabolic enzyme inhibitors and metabolic pathways exploited to overcome inhibition. Results indicate that knockouts whose growth rates are below ∼66% of the wild-type growth rate in a culture medium featuring lipid as the only carbon source are unable to sustain infections in granulomas. By mapping metabolite- and gene-scale perturbations to granuloma-scale outcomes and predicting mechanisms of sterilization, our method provides a powerful tool for hypothesis testing and guiding experimental searches for novel antituberculosis interventions.

Reconstruction and applications of consensus yeast metabolic network based on RNA sequencing

One practical application of genome-scale metabolic reconstructions is to interrogate multispecies relationships. Here, we report a consensus metabolic model in four yeast species (Saccharomyces cerevisiae, S. paradoxus, S. mikatae, and S. bayanus) by integrating metabolic network simulations with RNA sequencing (RNA-seq) datasets. We generated high-resolution transcriptome maps of four yeast species through de novo assembly and genome-guided approaches. The transcriptomes were annotated and applied to build the consensus metabolic network, which was verified using independent RNA-seq experiments. The expression profiles reveal that the genes involved in amino acid and lipid metabolism are highly coexpressed. The diverse phenotypic characteristics, such as cellular growth and gene deletions, can be simulated using the metabolic model. We also explored the applications of the consensus model in metabolic engineering using yeast-specific reactions and biofuel production as examples. Similar strategies will benefit communities studying genome-scale metabolic networks of other organisms.

Integrated Modeling of Gene Regulatory and Metabolic Networks in Mycobacterium tuberculosis

Mycobacterium tuberculosis (MTB) is the causative bacterium of tuberculosis, a disease responsible for over a million deaths worldwide annually with a growing number of strains resistant to antibiotics. The development of better therapeutics would greatly benefit from improved understanding of the mechanisms associated with MTB responses to different genetic and environmental perturbations. Therefore, we expanded a genome-scale regulatory-metabolic model for MTB using the Probabilistic Regulation of Metabolism (PROM) framework. Our model, MTBPROM2.0, represents a substantial knowledge base update and extension of simulation capability. We incorporated a recent ChIP-seq based binding network of 2555 interactions linking to 104 transcription factors (TFs) (representing a 3.5-fold expansion of TF coverage). We integrated this expanded regulatory network with a refined genome-scale metabolic model that can correctly predict growth viability over 69 source metabolite conditions and predict metabolic gene essentiality more accurately than the original model. We used MTBPROM2.0 to simulate the metabolic consequences of knocking out and overexpressing each of the 104 TFs in the model. MTBPROM2.0 improves performance of knockout growth defect predictions compared to the original PROM MTB model, and it can successfully predict growth defects associated with TF overexpression. Moreover, condition-specific models of MTBPROM2.0 successfully predicted synergistic growth consequences of overexpressing the TF whiB4 in the presence of two standard anti-TB drugs. MTBPROM2.0 can screen in silico condition-specific transcription factor perturbations to generate putative targets of interest that can help prioritize future experiments for therapeutic development efforts.

Cell scale host-pathogen modeling: another branch in the evolution of constraint-based methods

Constraint-based models have become popular methods for systems biology as they enable the integration of complex, disparate datasets in a biologically cohesive framework that also supports the description of biological processes in terms of basic physicochemical constraints and relationships. The scope, scale, and application of genome scale models have grown from single cell bacteria to multi-cellular interaction modeling; host-pathogen modeling represents one of these examples at the current horizon of constraint-based methods. There are now a small number of examples of host-pathogen constraint-based models in the literature, however there has not yet been a definitive description of the methodology required for the functional integration of genome scale models in order to generate simulation capable host-pathogen models. Herein we outline a systematic procedure to produce functional host-pathogen models, highlighting steps which require debugging and iterative revisions in order to successfully build a functional model. The construction of such models will enable the exploration of host-pathogen interactions by leveraging the growing wealth of omic data in order to better understand mechanism of infection and identify novel therapeutic strategies.

A multi-scale approach to designing therapeutics for tuberculosis

Approximately one third of the world's population is infected with Mycobacterium tuberculosis. Limited information about how the immune system fights M. tuberculosis and what constitutes protection from the bacteria impact our ability to develop effective therapies for tuberculosis. We present an in vivo systems biology approach that integrates data from multiple model systems and over multiple length and time scales into a comprehensive multi-scale and multi-compartment view of the in vivo immune response to M. tuberculosis. We describe computational models that can be used to study (a) immunomodulation with the cytokines tumor necrosis factor and interleukin 10, (b) oral and inhaled antibiotics, and

A review on computational systems biology of pathogen-host interactions

Pathogens manipulate the cellular mechanisms of host organisms via pathogen-host interactions (PHIs) in order to take advantage of the capabilities of host cells, leading to infections. The crucial role of these interspecies molecular interactions in initiating and sustaining infections necessitates a thorough understanding of the corresponding mechanisms. Unlike the traditional approach of considering the host or pathogen separately, a systems-level approach, considering the PHI system as a whole is indispensable to elucidate the mechanisms of infection. Following the technological advances in the post-genomic era, PHI data have been produced in large-scale within the last decade. Systems biology-based methods for the inference and analysis of PHI regulatory, metabolic, and protein-protein networks to shed light on infection mechanisms are gaining increasing demand thanks to the availability of omics data. The knowledge derived from the PHIs may largely contribute to the identification of new and more efficient therapeutics to prevent or cure infections. There are recent efforts for the detailed documentation of these experimentally verified PHI data through Web-based databases. Despite these advances in data archiving, there are still large amounts of PHI data in the biomedical literature yet to be discovered, and novel text mining methods are in development to unearth such hidden data. Here, we review a collection of recent studies on computational systems biology of PHIs with a special focus on the methods for the inference and analysis of PHI networks, covering also the Web-based databases and text-mining efforts to unravel the data hidden in the literature.

PhIN: A Protein Pharmacology Interaction Network Database

Network pharmacology is a new and hot concept in drug discovery for its ability to investigate the complexity of polypharmacology, and becomes more and more important in drug development. Here we report a protein pharmacology interaction network database (PhIN), aiming to assist multitarget drug discovery by providing comprehensive and flexible network pharmacology analysis. Overall, PhIN contains 1,126,060 target–target interaction pairs in terms of shared compounds and 3,428,020 pairs in terms of shared scaffolds, which involve 12,419,700 activity data, 9,414 targets, 314 viral targets, 652 pathways, 1,359,400 compounds, and 309,556 scaffolds. Using PhIN, users can obtain interacting target networks within or across human pathways, between human and virus, by defining the number of shared compounds or scaffolds under an activity cutoff. We expect PhIN to be a useful tool for multitarget drug development. PhIN is freely available at http://cadd.pharmacy.nankai.edu.cn/phin/.

On correlated reaction sets and coupled reaction sets in metabolic networks

Two reactions are in the same "correlated reaction set" (or "Co-Set") if their fluxes are linearly correlated. On the other hand, two reactions are "coupled" if nonzero flux through one reaction implies nonzero flux through the other reaction. Flux correlation analysis has been previously used in the analysis of enzyme dysregulation and enzymopathy, while flux coupling analysis has been used to predict co-expression of genes and to model network evolution. The goal of this paper is to emphasize, through a few examples, that these two concepts are inherently different. In other words, except for the case of full coupling, which implies perfect correlation between two fluxes (R(2) = 1), there are no constraints on Pearson correlation coefficients (CC) in case of any other type of (un)coupling relations. In other words, Pearson CC can take any value between 0 and 1 in other cases. Furthermore, by analyzing genome-scale metabolic networks, we confirm that there are some examples in real networks of bacteria, yeast and human, which approve that flux coupling and flux correlation cannot be used interchangeably.

Using a genome-scale metabolic model of Enterococcus faecalis V583 to assess amino acid uptake and its impact on central metabolism

Increasing antibiotic resistance in pathogenic bacteria necessitates the development of new medication strategies. Interfering with the metabolic network of the pathogen can provide novel drug targets but simultaneously requires a deeper and more detailed organism-specific understanding of the metabolism, which is often surprisingly sparse. In light of this, we reconstructed a genome-scale metabolic model of the pathogen Enterococcus faecalis V583. The manually curated metabolic network comprises 642 metabolites and 706 reactions. We experimentally determined metabolic profiles of E. faecalis grown in chemically defined medium in an anaerobic chemostat setup at different dilution rates and calculated the net uptake and product fluxes to constrain the model. We computed growth-associated energy and maintenance parameters and studied flux distributions through the metabolic network. Amino acid auxotrophies were identified experimentally for model validation and revealed seven essential amino acids. In addition, the important metabolic hub of glutamine/glutamate was altered by constructing a glutamine synthetase knockout mutant. The metabolic profile showed a slight shift in the fermentation pattern toward ethanol production and increased uptake rates of multiple amino acids, especially l-glutamine and l-glutamate. The model was used to understand the altered flux distributions in the mutant and provided an explanation for the experimentally observed redirection of the metabolic flux. We further highlighted the importance of gene-regulatory effects on the redirection of the metabolic fluxes upon perturbation. The genome-scale metabolic model presented here includes gene-protein-reaction associations, allowing a further use for biotechnological applications, for studying essential genes, proteins, or reactions, and the search for novel drug targets.

Methods for integration of transcriptomic data in genome-scale metabolic models

Several computational methods have been developed that integrate transcriptomic data with genome-scale metabolic reconstructions to infer condition-specific system-wide intracellular metabolic flux distributions. In this mini-review, we describe each of these methods published to date with categorizing them based on four different grouping criteria (requirement for multiple gene expression datasets as input, requirement for a threshold to define a gene's high and low expression, requirement for a priori assumption of an appropriate objective function, and validation of predicted fluxes directly against measured intracellular fluxes). Then, we recommend which group of methods would be more suitable from a practical perspective.

Computational nanomedicine: modeling of nanoparticle-mediated hyperthermal cancer therapy

Nanoparticle-mediated hyperthermia for cancer therapy is a growing area of cancer nanomedicine because of the potential for localized and targeted destruction of cancer cells. Localized hyperthermal effects are dependent on many factors, including nanoparticle size and shape, excitation wavelength and power, and tissue properties. Computational modeling is an important tool for investigating and optimizing these parameters. In this review, we focus on computational modeling of magnetic and gold nanoparticle-mediated hyperthermia, followed by a discussion of new opportunities and challenges.