Tunahan Çakır

PHISTO: pathogen–host interaction search tool

SUMMARY Knowledge of pathogen-host protein interactions is required to better understand infection mechanisms. The pathogen-host interaction search tool (PHISTO) is a web-accessible platform that provides relevant information about pathogen-host interactions (PHIs). It enables access to the most up-to-date PHI data for all pathogen types for which experimentally verified protein interactions with human are available. The platform also offers integrated tools for visualization of PHI networks, graph-theoretical analysis of targeted human proteins, BLAST search and text mining for detecting missing experimental methods. PHISTO will facilitate PHI studies that provide potential therapeutic targets for infectious diseases. AVAILABILITY http://www.phisto.org. CONTACT saliha.durmus@boun.edu.tr SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

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.

Infection Strategies of Bacterial and Viral Pathogens through Pathogen–Human Protein–Protein Interactions

Since ancient times, even in today’s modern world, infectious diseases cause lots of people to die. Infectious organisms, pathogens, cause diseases by physical interactions with human proteins. A thorough analysis of these interspecies interactions is required to provide insights about infection strategies of pathogens. Here we analyzed the most comprehensive available pathogen–human protein interaction data including 23,435 interactions, targeting 5,210 human proteins. The data were obtained from the newly developed pathogen–host interaction search tool, PHISTO. This is the first comprehensive attempt to get a comparison between bacterial and viral infections. We investigated human proteins that are targeted by bacteria and viruses to provide an overview of common and special infection strategies used by these pathogen types. We observed that in the human protein interaction network the proteins targeted by pathogens have higher connectivity and betweenness centrality values than those proteins not interacting with pathogens. The preference of interacting with hub and bottleneck proteins is found to be a common infection strategy of all types of pathogens to manipulate essential mechanisms in human. Compared to bacteria, viruses tend to interact with human proteins of much higher connectivity and centrality values in the human network. Gene Ontology enrichment analysis of the human proteins targeted by pathogens indicates crucial clues about the infection mechanisms of bacteria and viruses. As the main infection strategy, bacteria interact with human proteins that function in immune response to disrupt human defense mechanisms. Indispensable viral strategy, on the other hand, is the manipulation of human cellular processes in order to use that transcriptional machinery for their own genetic material transcription. A novel observation about pathogen–human systems is that the human proteins targeted by both pathogens are enriched in the regulation of metabolic processes.

Effect of carbon source perturbations on transcriptional regulation of metabolic fluxes in Saccharomyces cerevisiae

BackgroundControl effective flux (CEF) of a reaction is the weighted sum of all fluxes through that reaction, derived from elementary flux modes (EFM) of a metabolic network. Change in CEFs under different environmental conditions has earlier been proven to be correlated with the corresponding changes in the transcriptome. Here we use this to investigate the degree of transcriptional regulation of fluxes in the metabolism of Saccharomyces cerevisiae. We do this by quantifying correlations between changes in CEFs and changes in transcript levels for shifts in carbon source, i.e. between the fermentative carbon source glucose and nonfermentative carbon sources like ethanol, acetate, and lactate. The CEF analysis is based on a simple stoichiometric model that includes reactions of the central carbon metabolism and the amino acid metabolism.ResultsThe effect of the carbon shift on the metabolic fluxes was investigated for both batch and chemostat cultures. For growth on glucose in batch (respiro-fermentative) cultures, EFMs with no by-product formation were removed from the analysis of the CEFs, whereas those including any by-products (ethanol, glycerol, acetate, succinate) were omitted in the analysis of growth on glucose in chemostat (respiratory) cultures. This resulted in improved correlations between CEF changes and transcript levels. A regression correlation coefficient of 0.60 was obtained between CEF changes and gene expression changes in the central carbon metabolism for the analysis of 5 different perturbations. Out of 45 data points there were no more than 6 data points deviating from the correlation. Additionally, up- or down-regulation of at least 75% of the genes were in qualitative agreement with the CEF changes for all perturbations studied.ConclusionThe analysis indicates that changes in carbon source are associated with a high degree of hierarchical regulation of metabolic fluxes in the central carbon metabolism as the change in fluxes are correlating directly with the change in transcript levels of genes encoding their corresponding enzymes. For amino acid biosynthesis there was, however, not found to exist a similar correlation, and this may point to either post-transcriptional and/or metabolic regulation, or be due to the absence of a direct perturbation on the amino acid pathways in these experiments.

Flux balance analysis of a genome-scale yeast model constrained by exometabolomic data allows metabolic system identification of genetically different strains

A systems approach to biology requires a principled approach to pathway identification. In this study, the two nuclear petite yeast mutants K1Δpet191a and K1Δpet191ab and their parental industrial strain K1 were cultured in glucose‐containing microaerobic chemostats. Exometabolomic profiles were used to infer the differences in the fermentation characteristics and respiration capacity of the strains. The ability of the metabolite measurement information to describe genetically different strains was investigated using a genome‐scale yeast model. Flux balance analysis (FBA) of the model reveals that the objective function of minimal oxygen consumption enables the identification of the effect of genotypic differences when combined with the knowledge of the extracellular state of metabolism. The predicted decrease in oxygen consumption flux of K1Δpet191a and K1Δpet191ab strains with respect to the parental strain is about 80% and 100%, respectively, which coincides with the respiratory deficiencies of the strains. The expected increase in ethanol production rates in response to the decrease in the respiratory capacity was also predicted to be very close to the experimental values. This study shows the predictive power of the integrated analysis of genome‐scale models with exometabolomic profiles, since accurate predictions could be made without any information about the respiration capacity of the strains. The FBA approach thereby enables identification of responsive pathways and so permits the elucidation of the genetic characteristics of strains in terms of expressed metabolite profiles.

Systematic analysis of transcription-level effects of neurodegenerative diseases on human brain metabolism by a newly reconstructed brain-specific metabolic network

Network‐oriented analysis is essential to identify those parts of a cell affected by a given perturbation. The effect of neurodegenerative perturbations in the form of diseases of brain metabolism was investigated by using a newly reconstructed brain‐specific metabolic network. The developed stoichiometric model correctly represents healthy brain metabolism, and includes 630 metabolic reactions in and between astrocytes and neurons, which are controlled by 570 genes. The integration of transcriptome data of six neurodegenerative diseases (Alzheimers disease, Parkinsons disease, Huntingtons disease, amyotrophic lateral sclerosis, multiple sclerosis, schizophrenia) with the model was performed to identify reporter features specific and common for these diseases, which revealed metabolites and pathways around which the most significant changes occur. The identified metabolites are potential biomarkers for the pathology of the related diseases. Our model indicated perturbations in oxidative stress, energy metabolism including TCA cycle and lipid metabolism as well as several amino acid related pathways, in agreement with the role of these pathways in the studied diseases. The computational prediction of transcription factors that commonly regulate the reporter metabolites was achieved through binding‐site analysis. Literature support for the identified transcription factors such as USF1, SP1 and those from FOX families are known from the literature to have regulatory roles in the identified reporter metabolic pathways as well as in the neurodegenerative diseases. In essence, the reconstructed brain model enables the elucidation of effects of a perturbation on brain metabolism and the illumination of possible machineries in which a specific metabolite or pathway acts as a regulatory spot for cellular reorganization.

Metabolic network discovery by top-down and bottom-up approaches and paths for reconciliation.

The primary focus in the network-centric analysis of cellular metabolism by systems biology approaches is to identify the active metabolic network for the condition of interest. Two major approaches are available for the discovery of the condition-specific metabolic networks. One approach starts from genome-scale metabolic networks, which cover all possible reactions known to occur in the related organism in a condition-independent manner, and applies methods such as the optimization-based Flux-Balance Analysis to elucidate the active network. The other approach starts from the condition-specific metabolome data, and processes the data with statistical or optimization-based methods to extract information content of the data such that the active network is inferred. These approaches, termed bottom-up and top-down, respectively, are currently employed independently. However, considering that both approaches have the same goal, they can both benefit from each other paving the way for the novel integrative analysis methods of metabolome data- and flux-analysis approaches in the post-genomic era. This study reviews the strengths of constraint-based analysis and network inference methods reported in the metabolic systems biology field; then elaborates on the potential paths to reconcile the two approaches to shed better light on how the metabolism functions.

Sparsity as Cellular Objective to Infer Directed Metabolic Networks from Steady-State Metabolome Data: A Theoretical Analysis

Since metabolome data are derived from the underlying metabolic network, reverse engineering of such data to recover the network topology is of wide interest. Lyapunov equation puts a constraint to the link between data and network by coupling the covariance of data with the strength of interactions (Jacobian matrix). This equation, when expressed as a linear set of equations at steady state, constitutes a basis to infer the network structure given the covariance matrix of data. The sparse structure of metabolic networks points to reactions which are active based on minimal enzyme production, hinting at sparsity as a cellular objective. Therefore, for a given covariance matrix, we solved Lyapunov equation to calculate Jacobian matrix by a simultaneous use of minimization of Euclidean norm of residuals and maximization of sparsity (the number of zeros in Jacobian matrix) as objective functions to infer directed small-scale networks from three kingdoms of life (bacteria, fungi, mammalian). The inference performance of the approach was found to be promising, with zero False Positive Rate, and almost one True positive Rate. The effect of missing data on results was additionally analyzed, revealing superiority over similarity-based approaches which infer undirected networks. Our findings suggest that the covariance of metabolome data implies an underlying network with sparsest pattern. The theoretical analysis forms a framework for further investigation of sparsity-based inference of metabolic networks from real metabolome data.

Reconstructed Metabolic Network Models Predict Flux-Level Metabolic Reprogramming in Glioblastoma.

Developments in genome scale metabolic modeling techniques and omics technologies have enabled the reconstruction of context-specific metabolic models. In this study, glioblastoma multiforme (GBM), one of the most common and aggressive malignant brain tumors, is investigated by mapping GBM gene expression data on the growth-implemented brain specific genome-scale metabolic network, and GBM-specific models are generated. The models are used to calculate metabolic flux distributions in the tumor cells. Metabolic phenotypes predicted by the GBM-specific metabolic models reconstructed in this work reflect the general metabolic reprogramming of GBM, reported both in in-vitro and in-vivo experiments. The computed flux profiles quantitatively predict that major sources of the acetyl-CoA and oxaloacetic acid pool used in TCA cycle are pyruvate dehydrogenase from glycolysis and anaplerotic flux from glutaminolysis, respectively. Also, our results, in accordance with recent studies, predict a contribution of oxidative phosphorylation to ATP pool via a slightly active TCA cycle in addition to the major contributor aerobic glycolysis. We verified our results by using different computational methods that incorporate transcriptome data with genome-scale models and by using different transcriptome datasets. Correct predictions of flux distributions in glycolysis, glutaminolysis, TCA cycle and lipid precursor metabolism validate the reconstructed models for further use in future to simulate more specific metabolic patterns for GBM.

Transfer function approach in structured modeling of recombinant yeast utilizing starch

Abstract Compartment models are based on the division of biomass into a few pools (compartment). In this work, three biochemically structured (three-compartment, four-compartment, and eight-compartment) models are constructed to analyze anaerobic starch fermentation and growth characteristics of recombinant Saccharomyces cerevisiae . The common compartments in these models are the active compartment ( X A ), the plasmid product fusion protein ( X P ) and the structural part of the microorganism ( X G ). Each compartment is described by a kinetic expression depending on experimental observations and studies in the literature. Parameters of the resultant kinetic model are estimated by an optimization routine. The model equations are then solved simultaneously to produce time profiles of the constructed compartments. In addition, the transfer function approach is applied to the eight-compartment model to analyze its degree of complexity in terms of the characteristic time concept. The characteristic time constants of the compartments lead to a conclusion that the eight-compartment model can be reduced to a four-compartment model. Best simulation results were obtained with the four-compartment model, confirming the results of the transfer function approach.