Network-based assessment of the selectivity of metabolic drug targets in Plasmodium falciparum with respect to human liver metabolism

Structure-based pharmacophore modeling, virtual screening, and molecular dynamics simulation studies for identification of Plasmodium falciparum 5-aminolevulinate synthase inhibitors

Plasmodium falciparum (Pf) 5-aminolevulinic acid synthase (5-ALAS) is an essential enzyme with high selectivity during liver stage development, signifying its potential as a prophylactic antimalarial drug target. The aim of this study was to identify important potential lead compounds which can serve as inhibitors of Pf 5-ALAS using pharmacophore modeling, virtual screening, qualitative structural assessment, in silico ADMET (Absorption, Distribution, Metabolism, Excretion and Toxicity) evaluation and molecular dynamics simulation. The best model of the tertiary structure of Pf 5-ALAS was obtained using MolProbity, while the following databases were explored for the pharmacophore-based virtual screening: CHEMBL, ChemDiv, ChemSpace, MCULE, MCULE-ULTIMATE, MolPort, NCI Open Chemical Repository, LabNetwork and ZINC databases. 2,621 compounds were screened against the modeled Pf 5-ALAS using AutoDock vina. The post-screening analysis was carried out using Discovery Studio while molecular dynamics simulation was performed on the best hits using NAMD-VMD and Galaxy Europe platform. Compound CSMS00081585868 was observed as the best hit with a binding affinity of -9.9 kcal/mol and predicted Ki of 52.10 nM, engaging in seven hydrogen bonds with the target's active site amino acid residues. The in silico ADMET prediction showed that all ten best hits possessed relatively good pharmacokinetic properties. The qualitative structural assessment of the best hit, CSMS00081585868, revealed that the presence of two pyridine scaffolds bearing hydroxy and fluorine groups linked by a pyrrolidine scaffold contributed significantly to its ability to have a strong binding affinity with the receptor. The best hit also showed stability in the active site of Pf 5-ALAS as confirmed from the RMSD obtained during the MD simulation.

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.

Genome reconstructions of metabolism of Plasmodium RBC and liver stages

Genome scale metabolic models (GEMs) offer a powerful means of integrating genome and biochemical information on an organism to make testable predictions of metabolic functions at different conditions and to systematically predict essential genes that may be targeted by drugs. This review describes how Plasmodium GEMs have become increasingly more accurate through the integration of omics and experimental genetic data. We also discuss how GEMs contribute to our increasing understanding of how Plasmodium metabolism is reprogrammed between life cycle stages.

Non-canonical metabolic pathways in the malaria parasite detected by isotope-tracing metabolomics

The malaria parasite, Plasmodium falciparum, proliferates rapidly in human erythrocytes by actively scavenging multiple carbon sources and essential nutrients from its host cell. However, a global overview of the metabolic capacity of intraerythrocytic stages is missing. Using multiplex 13 C-labelling coupled with untargeted mass spectrometry and unsupervised isotopologue grouping, we have generated a draft metabolome of P. falciparum and its host erythrocyte consisting of 911 and 577 metabolites, respectively, corresponding to 41% of metabolites and over 70% of the metabolic reaction predicted from the parasite genome. An additional 89 metabolites and 92 reactions were identified that were not predicted from genomic reconstructions, with the largest group being associated with metabolite damage-repair systems. Validation of the draft metabolome revealed four previously uncharacterised enzymes which impact isoprenoid biosynthesis, lipid homeostasis and mitochondrial metabolism and are necessary for parasite development and proliferation. This study defines the metabolic fate of multiple carbon sources in P. falciparum, and highlights the activity of metabolite repair pathways in these rapidly growing parasite stages, opening new avenues for drug discovery.

Genome-Scale Metabolic Modeling for Unraveling Molecular Mechanisms of High Threat Pathogens

Pathogens give rise to a wide range of diseases threatening global health and hence drawing public health agencies' attention to establish preventative and curative solutions. Genome-scale metabolic modeling is ever increasingly used tool for biomedical applications including the elucidation of antibiotic resistance, virulence, single pathogen mechanisms and pathogen-host interaction systems. With this approach, the sophisticated cellular system of metabolic reactions inside the pathogens as well as between pathogen and host cells are represented in conjunction with their corresponding genes and enzymes. Along with essential metabolic reactions, alternate pathways and fluxes are predicted by performing computational flux analyses for the growth of pathogens in a very short time. The genes or enzymes responsible for the essential metabolic reactions in pathogen growth are regarded as potential drug targets, as a priori guide to researchers in the pharmaceutical field. Pathogens alter the key metabolic processes in infected host, ultimately the objective of these integrative constraint-based context-specific metabolic models is to provide novel insights toward understanding the metabolic basis of the acute and chronic processes of infection, revealing cellular mechanisms of pathogenesis, identifying strain-specific biomarkers and developing new therapeutic approaches including the combination drugs. The reaction rates predicted during different time points of pathogen development enable us to predict active pathways and those that only occur during certain stages of infection, and thus point out the putative drug targets. Among others, fatty acid and lipid syntheses reactions are recent targets of new antimicrobial drugs. Genome-scale metabolic models provide an improved understanding of how intracellular pathogens utilize the existing microenvironment of the host. Here, we reviewed the current knowledge of genome-scale metabolic modeling in pathogen cells as well as pathogen host interaction systems and the promising applications in the extension of curative strategies against pathogens for global preventative healthcare.

New Therapeutic Candidates for the Treatment of Malassezia pachydermatis -Associated Infections

The opportunistic pathogen Malassezia pachydermatis causes bloodstream infections in preterm infants or individuals with immunodeficiency disorders and has been associated with a broad spectrum of diseases in animals such as seborrheic dermatitis, external otitis and fungemia. The current approaches to treat these infections are failing as a consequence of their adverse effects, changes in susceptibility and antifungal resistance. Thus, the identification of novel therapeutic targets against M. pachydermatis infections are highly relevant. Here, Gene Essentiality Analysis and Flux Variability Analysis was applied to a previously reported M. pachydermatis metabolic network to identify enzymes that, when absent, negatively affect biomass production. Three novel therapeutic targets (i.e., homoserine dehydrogenase (MpHSD), homocitrate synthase (MpHCS) and saccharopine dehydrogenase (MpSDH)) were identified that are absent in humans. Notably, L-lysine was shown to be an inhibitor of the enzymatic activity of MpHCS and MpSDH at concentrations of 1 mM and 75 mM, respectively, while L-threonine (1 mM) inhibited MpHSD. Interestingly, L- lysine was also shown to inhibit M. pachydermatis growth during in vitro assays with reference strains and canine isolates, while it had a negligible cytotoxic activity on HEKa cells. Together, our findings form the bases for the development of novel treatments against M. pachydermatis infections.

Novel Approaches for Systems Biology of Metabolism-Oriented Pathogen-Human Interactions: A Mini-Review

Pathogenic microorganisms exploit host metabolism for sustained survival by rewiring its metabolic interactions. Therefore, several metabolic changes are induced in both pathogen and host cells in the course of infection. A systems-based approach to elucidate those changes includes the integrative use of genome-scale metabolic networks and molecular omics data, with the overall goal of better characterizing infection mechanisms for novel treatment strategies. This review focuses on novel aspects of metabolism-oriented systems-based investigation of pathogen-human interactions. The reviewed approaches are the generation of dual-omics data for the characterization of metabolic signatures of pathogen-host interactions, the reconstruction of pathogen-host integrated genome-scale metabolic networks, which has a high potential to be applied to pathogen-gut microbiota interactions, and the structure-based analysis of enzymes playing role in those interactions. The integrative use of those approaches will pave the way for the identification of novel biomarkers and drug targets for the prediction and prevention of infectious diseases.

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.

In Silico Knockout Screening of Plasmodium falciparum Reactions and Prediction of Novel Essential Reactions by Analysing the Metabolic Network

Malaria is an infectious disease that affects close to half a million individuals every year and Plasmodium falciparum is a major cause of malaria. The treatment of this disease could be done effectively if the essential enzymes of this parasite are specifically targeted. Nevertheless, the development of the parasite in resisting existing drugs now makes discovering new drugs a core responsibility. In this study, a novel computational model that makes the prediction of new and validated antimalarial drug target cheaper, easier, and faster has been developed. We have identified new essential reactions as potential targets for drugs in the metabolic network of the parasite. Among the top seven (7) predicted essential reactions, four (4) have been previously identified in earlier studies with biological evidence and one (1) has been with computational evidence. The results from our study were compared with an extensive list of seventy-seven (77) essential reactions with biological evidence from a previous study. We present a list of thirty-one (31) potential candidates for drug targets in Plasmodium falciparum which includes twenty-four (24) new potential candidates for drug targets.

Computational Systems Biology of Metabolism in Infection

A systems approach to elucidate the effect of infection on cell metabolism provides several opportunities from a better understanding of molecular mechanisms to the identification of potential biomarkers and drug targets. This is obvious from the fact that we have witnessed the accelerated use of computational systems biology in the last five years to study metabolic changes in pathogen and/or host cells in response to infection. In this chapter, we aim to present a comprehensive review of the recent research by focusing on genome-scale metabolic network models of pathogen-host systems and genome-wide metabolomics and fluxomics analysis of infected cells.

Novel Plasmodium falciparum metabolic network reconstruction identifies shifts associated with clinical antimalarial resistance

Background

Malaria remains a major public health burden and resistance has emerged to every antimalarial on the market, including the frontline drug, artemisinin. Our limited understanding of Plasmodium biology hinders the elucidation of resistance mechanisms. In this regard, systems biology approaches can facilitate the integration of existing experimental knowledge and further understanding of these mechanisms.

Results

Here, we developed a novel genome-scale metabolic network reconstruction, iPfal17, of the asexual blood-stage P. falciparum parasite to expand our understanding of metabolic changes that support resistance. We identified 11 metabolic tasks to evaluate iPfal17 performance. Flux balance analysis and simulation of gene knockouts and enzyme inhibition predict candidate drug targets unique to resistant parasites. Moreover, integration of clinical parasite transcriptomes into the iPfal17 reconstruction reveals patterns associated with antimalarial resistance. These results predict that artemisinin sensitive and resistant parasites differentially utilize scavenging and biosynthetic pathways for multiple essential metabolites, including folate and polyamines. Our findings are consistent with experimental literature, while generating novel hypotheses about artemisinin resistance and parasite biology. We detect evidence that resistant parasites maintain greater metabolic flexibility, perhaps representing an incomplete transition to the metabolic state most appropriate for nutrient-rich blood.

Conclusion

Using this systems biology approach, we identify metabolic shifts that arise with or in support of the resistant phenotype. This perspective allows us to more productively analyze and interpret clinical expression data for the identification of candidate drug targets for the treatment of resistant parasites.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3905-1) contains supplementary material, which is available to authorized users.

Bioenergetics-based modeling of Plasmodium falciparum metabolism reveals its essential genes, nutritional requirements, and thermodynamic bottlenecks

Novel antimalarial therapies are urgently needed for the fight against drug-resistant parasites. The metabolism of malaria parasites in infected cells is an attractive source of drug targets but is rather complex. Computational methods can handle this complexity and allow integrative analyses of cell metabolism. In this study, we present a genome-scale metabolic model (iPfa) of the deadliest malaria parasite, Plasmodium falciparum, and its thermodynamics-based flux analysis (TFA). Using previous absolute concentration data of the intraerythrocytic parasite, we applied TFA to iPfa and predicted up to 63 essential genes and 26 essential pairs of genes. Of the 63 genes, 35 have been experimentally validated and reported in the literature, and 28 have not been experimentally tested and include previously hypothesized or novel predictions of essential metabolic capabilities. Without metabolomics data, four of the genes would have been incorrectly predicted to be non-essential. TFA also indicated that substrate channeling should exist in two metabolic pathways to ensure the thermodynamic feasibility of the flux. Finally, analysis of the metabolic capabilities of P. falciparum led to the identification of both the minimal nutritional requirements and the genes that can become indispensable upon substrate inaccessibility. This model provides novel insight into the metabolic needs and capabilities of the malaria parasite and highlights metabolites and pathways that should be measured and characterized to identify potential thermodynamic bottlenecks and substrate channeling. The hypotheses presented seek to guide experimental studies to facilitate a better understanding of the parasite metabolism and the identification of targets for more efficient intervention.

Almost half of the world population is at risk of infection by malaria parasites. The rise in drug-resistant parasites requires better understanding and targeting of their metabolism. In this study, we present a genome-scale metabolic reconstruction (iPfa) of the deadliest malaria parasite, Plasmodium falciparum, and its thermodynamics-based flux analysis (TFA). Our results support and extend the available experimental evidence on the essential genes and nutritional requirements of this organism. Besides, we identify metabolites that give rise to thermodynamic bottlenecks and suggest substrate channeling. Overall, these results provide novel insight into the metabolism of P. falciparum and may guide experimental studies to develop a better characterization of the parasite metabolism and the identification of antimalarial drug targets.

Metabolic host responses to malarial infection during the intraerythrocytic developmental cycle

BACKGROUND

The malarial parasite Plasmodium falciparum undergoes a complex life cycle, including an intraerythrocytic developmental cycle, during which it is metabolically dependent on the infected human red blood cell (RBC). To describe whole cell metabolic activity within both P. falciparum and RBCs during the asexual reproduction phase of the intraerythrocytic developmental cycle, we developed an integrated host-parasite metabolic modeling framework driven by time-dependent gene expression data.

RESULTS

We validated the model by reproducing the experimentally determined 1) stage-specific production of biomass components and their precursors in the parasite and 2) metabolite concentration changes in the medium of P. falciparum-infected RBC cultures. The model allowed us to explore time- and strain-dependent P. falciparum metabolism and hypothesize how host cell metabolism alters in response to malarial infection. Specifically, the metabolic analysis showed that uninfected RBCs that coexist with infected cells in the same culture decrease their production of 2,3-bisphosphoglycerate, an oxygen-carrying regulator, reducing the ability of hemoglobin in these cells to release oxygen. Furthermore, in response to parasite-induced oxidative stress, infected RBCs downgraded their glycolytic flux by using the pentose phosphate pathway and secreting ribulose-5-phosphate. This mechanism links individually observed experimental phenomena, such as glycolytic inhibition and ribulose-5-phosphate secretion, to the oxidative stress response.

CONCLUSIONS

Although the metabolic model does not incorporate regulatory mechanisms per se, alterations in gene expression levels caused by regulatory mechanisms are manifested in the model as altered metabolic states. This provides the model the capability to capture complex multicellular host-pathogen metabolic interactions of the infected RBC culture. The system-level analysis revealed complex relationships such as how the parasite can reduce oxygen release in uninfected cells in the presence of infected RBCs as well as the role of different metabolic pathways involved in the oxidative stress response of infected RBCs.

Mathematical Modeling of Cellular Metabolism

Cellular metabolism basically consists of the conversion of chemical compounds taken up from the extracellular environment into energy (conserved in energy-rich bonds of organic phosphates) and a wide array of organic molecules serving as catalysts (enzymes), information carriers (nucleic acids), and building blocks for cellular structures such as membranes or ribosomes. Metabolic modeling aims at the construction of mathematical representations of the cellular metabolism that can be used to calculate the concentration of cellular molecules and the rates of their mutual chemical interconversion in response to varying external conditions as, for example, hormonal stimuli or supply of essential nutrients. Based on such calculations, it is possible to quantify complex cellular functions as cellular growth, detoxification of drugs and xenobiotic compounds or synthesis of exported molecules. Depending on the specific questions to metabolism addressed, the methodological expertise of the researcher, and available experimental information, different conceptual frameworks have been established, allowing the usage of computational methods to condense experimental information from various layers of organization into (self-) consistent models. Here, we briefly outline the main conceptual frameworks that are currently exploited in metabolism research.

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.

Drug target identification through systems biology

To rationalise drug target selection, we should understand the role of putative targets in biological pathways quantitatively. We review here how experimental and computational network-based approaches can aid more rational drug target selection and illustrate this with results obtained for microbes and for cancer. Comparison of the drug response of biochemical networks in target cells and (healthy) host cells can reveal network-selective targets.

Phylogenetic Analysis of Nucleus-Encoded Acetyl-CoA Carboxylases Targeted at the Cytosol and Plastid of Algae

The understanding of algal phylogeny is being impeded by an unknown number of events of horizontal gene transfer (HGT), and primary and secondary/tertiary endosymbiosis. Through these events, previously heterotrophic eukaryotes developed photosynthesis and acquired new biochemical pathways. Acetyl-CoA carboxylase (ACCase) is a key enzyme in the fatty acid synthesis and elongation pathways in algae, where ACCase exists in two locations (cytosol and plastid) and in two forms (homomeric and heteromeric). All algae contain nucleus-encoded homomeric ACCase in the cytosol, independent of the origin of the plastid. Nucleus-encoded homomeric ACCase is also found in plastids of algae that arose from a secondary/tertiary endosymbiotic event. In contrast, plastids of algae that arose from a primary endosymbiotic event contain heteromeric ACCase, which consists of three nucleus-encoded and one plastid-encoded subunits. These properties of ACCase provide the potential to inform on the phylogenetic relationships of hosts and their plastids, allowing different hypothesis of endosymbiotic events to be tested. Alveolata (Dinoflagellata and Apicomplexa) and Chromista (Stramenopiles, Haptophyta and Cryptophyta) have traditionally been grouped together as Chromalveolata, forming the red lineage. However, recent genetic evidence groups the Stramenopiles, Alveolata and green plastid containing Rhizaria as SAR, excluding Haptophyta and Cryptophyta. Sequences coding for plastid and cytosol targeted homomeric ACCases were isolated from Isochrysis aff. galbana (TISO), Chromera velia and Nannochloropsis oculata, representing three taxonomic groups for which sequences were lacking. Phylogenetic analyses show that cytosolic ACCase strongly supports the SAR grouping. Conversely, plastidial ACCase groups the SAR with the Haptophyta, Cryptophyta and Prasinophyceae (Chlorophyta). These two ACCase based, phylogenetic relationships suggest that the plastidial homomeric ACCase was acquired by the Haptophyta, Cryptophyta and SAR, before the photosynthetic Rhizaria acquired their green plastid. Additionally, plastidial ACCase was derived by HGT from an ancestor or relative of the Prasinophyceae and not by duplication of cytosolic ACCase.

From next-generation sequencing to systematic modeling of the gut microbiome

Changes in the human gut microbiome are associated with altered human metabolism and health, yet the mechanisms of interactions between microbial species and human metabolism have not been clearly elucidated. Next-generation sequencing has revolutionized the human gut microbiome research, but most current applications concentrate on studying the microbial diversity of communities and have at best provided associations between specific gut bacteria and human health. However, little is known about the inner metabolic mechanisms in the gut ecosystem. Here we review recent progress in modeling the metabolic interactions of gut microbiome, with special focus on the utilization of metabolic modeling to infer host–microbe interactions and microbial species interactions. The systematic modeling of metabolic interactions could provide a predictive understanding of gut microbiome, and pave the way to synthetic microbiota design and personalized-microbiome medicine and healthcare. Finally, we discuss the integration of metabolic modeling and gut microbiome engineering, which offer a new way to explore metabolic interactions across members of the gut microbiota.

Modeling metabolism and stage-specific growth of Plasmodium falciparum HB3 during the intraerythrocytic developmental cycle

The human malaria parasite Plasmodium falciparum goes through a complex life cycle, including a roughly 48-hour-long intraerythrocytic developmental cycle (IDC) in human red blood cells. A better understanding of the metabolic processes required during the asexual blood-stage reproduction will enhance our basic knowledge of P. falciparum and help identify critical metabolic reactions and pathways associated with blood-stage malaria. We developed a metabolic network model that mechanistically links time-dependent gene expression, metabolism, and stage-specific growth, allowing us to predict the metabolic fluxes, the biomass production rates, and the timing of production of the different biomass components during the IDC. We predicted time- and stage-specific production of precursors and macromolecules for P. falciparum (strain HB3), allowing us to link specific metabolites to specific physiological functions. For example, we hypothesized that coenzyme A might be involved in late-IDC DNA replication and cell division. Moreover, the predicted ATP metabolism indicated that energy was mainly produced from glycolysis and utilized for non-metabolic processes. Finally, we used the model to classify the entire tricarboxylic acid cycle into segments, each with a distinct function, such as superoxide detoxification, glutamate/glutamine processing, and metabolism of fumarate as a byproduct of purine biosynthesis. By capturing the normal metabolic and growth progression in P. falciparum during the IDC, our model provides a starting point for further elucidation of strain-specific metabolic activity, host-parasite interactions, stress-induced metabolic responses, and metabolic responses to antimalarial drugs and drug candidates.

Metabolic Needs and Capabilities of Toxoplasma gondii through Combined Computational and Experimental Analysis

Toxoplasma gondii is a human pathogen prevalent worldwide that poses a challenging and unmet need for novel treatment of toxoplasmosis. Using a semi-automated reconstruction algorithm, we reconstructed a genome-scale metabolic model, ToxoNet1. The reconstruction process and flux-balance analysis of the model offer a systematic overview of the metabolic capabilities of this parasite. Using ToxoNet1 we have identified significant gaps in the current knowledge of Toxoplasma metabolic pathways and have clarified its minimal nutritional requirements for replication. By probing the model via metabolic tasks, we have further defined sets of alternative precursors necessary for parasite growth. Within a human host cell environment, ToxoNet1 predicts a minimal set of 53 enzyme-coding genes and 76 reactions to be essential for parasite replication. Double-gene-essentiality analysis identified 20 pairs of genes for which simultaneous deletion is deleterious. To validate several predictions of ToxoNet1 we have performed experimental analyses of cytosolic acetyl-CoA biosynthesis. ATP-citrate lyase and acetyl-CoA synthase were localised and their corresponding genes disrupted, establishing that each of these enzymes is dispensable for the growth of T. gondii, however together they make a synthetic lethal pair.

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.

Promise and reality in the expanding field of network interaction analysis: metabolic networks

In the last few decades, metabolic networks revealed their capabilities as powerful tools to analyze the cellular metabolism. Many research fields (eg, metabolic engineering, diagnostic medicine, pharmacology, biochemistry, biology and physiology) improved the understanding of the cell combining experimental assays and metabolic network-based computations. This process led to the rise of the “systems biology” approach, where the theory meets experiments and where two complementary perspectives cooperate in the study of biological phenomena. Here, the reconstruction of metabolic networks is presented, along with established and new algorithms to improve the description of cellular metabolism. Then, advantages and limitations of modeling algorithms and network reconstruction are discussed.

Partial inhibition and bilevel optimization in flux balance analysis

MOTIVATION

Within Flux Balance Analysis, the investigation of complex subtasks, such as finding the optimal perturbation of the network or finding an optimal combination of drugs, often requires to set up a bilevel optimization problem. In order to keep the linearity and convexity of these nested optimization problems, an ON/OFF description of the effect of the perturbation (i.e. Boolean variable) is normally used. This restriction may not be realistic when one wants, for instance, to describe the partial inhibition of a reaction induced by a drug.

RESULTS

In this paper we present a formulation of the bilevel optimization which overcomes the oversimplified ON/OFF modeling while preserving the linear nature of the problem. A case study is considered: the search of the best multi-drug treatment which modulates an objective reaction and has the minimal perturbation on the whole network. The drug inhibition is described and modulated through a convex combination of a fixed number of Boolean variables. The results obtained from the application of the algorithm to the core metabolism of E.coli highlight the possibility of finding a broader spectrum of drug combinations compared to a simple ON/OFF modeling.

CONCLUSIONS

The method we have presented is capable of treating partial inhibition inside a bilevel optimization, without loosing the linearity property, and with reasonable computational performances also on large metabolic networks. The more fine-graded representation of the perturbation allows to enlarge the repertoire of synergistic combination of drugs for tasks such as selective perturbation of cellular metabolism. This may encourage the use of the approach also for other cases in which a more realistic modeling is required.

Functional genomics of Plasmodium falciparum using metabolic modelling and analysis

Plasmodium falciparum is an obligate intracellular parasite and the leading cause of severe malaria responsible for tremendous morbidity and mortality particularly in sub-Saharan Africa. Successful completion of the P. falciparum genome sequencing project in 2002 provided a comprehensive foundation for functional genomic studies on this pathogen in the following decade. Over this period, a large spectrum of experimental approaches has been deployed to improve and expand the scope of functionally annotated genes. Meanwhile, rapidly evolving methods of systems biology have also begun to contribute to a more global understanding of various aspects of the biology and pathogenesis of malaria. Herein we provide an overview on metabolic modelling, which has the capability to integrate information from functional genomics studies in P. falciparum and guide future malaria research efforts towards the identification of novel candidate drug targets.

Docking studies and network analyses reveal capacity of compounds from Kandelia rheedii to strengthen cellular immunity by interacting with host proteins during tuberculosis infection

Kandelia rheedii (locally known as Guria or Rasunia), widely found and used in Indian subcontinent, is a well-known herbal cure to tuberculosis. However, neither the mechanism nor the active components of the plant extract responsible for mediating this action has yet been confirmed. Here in this study, molecular interactions of three compounds (emodin, fusaric acid and skyrin) from the plant extract with the host protein targets (casein kinase (CSNK), estrogen receptor (ERBB), dopamine β-hydroxylase (DBH) and glucagon receptor (Gcgr)) has been found. These protein targets are known to be responsible for strengthening cellular immunity against Mycobacteria tuberculosis. The specific interactions of these three compounds with the respective protein targets have been discussed here. The insights from study should further help us designing molecular medicines against tuberculosis.