Reconstruction and flux-balance analysis of the Plasmodium falciparum metabolic network

Host cell CRISPR genomics and modelling reveal shared metabolic vulnerabilities in the intracellular development of Plasmodium falciparum and related hemoparasites

Parasitic diseases, particularly malaria (caused by Plasmodium falciparum) and theileriosis (caused by Theileria spp.), profoundly impact global health and the socioeconomic well-being of lower-income countries. Despite recent advances, identifying host metabolic proteins essential for these auxotrophic pathogens remains challenging. Here, we generate a novel metabolic model of human hepatocytes infected with P. falciparum and integrate it with a genome-wide CRISPR knockout screen targeting Theileria-infected cells to pinpoint shared vulnerabilities. We identify key host metabolic enzymes critical for the intracellular survival of both of these lethal hemoparasites. Remarkably, among the metabolic proteins identified by our synergistic approach, we find that host purine and heme biosynthetic enzymes are essential for the intracellular survival of P. falciparum and Theileria, while other host enzymes are only essential under certain metabolic conditions, highlighting P. falciparum's adaptability and ability to scavenge nutrients selectively. Unexpectedly, host porphyrins emerge as being essential for both parasites. The shared vulnerabilities open new avenues for developing more effective therapies against these debilitating diseases, with the potential for broader applicability in combating apicomplexan infections.

Comparison of extraction methods in vitro Plasmodium falciparum: A 1H NMR and LC-MS joined approach

Malaria is a parasitic disease that remains a global concern and the subject of many studies. Metabolomics has emerged as an approach to better comprehend complex pathogens and discover possible drug targets, thus giving new insights that can aid in the development of antimalarial therapies. However, there is no standardized method to extract metabolites from in vitro Plasmodium falciparum intraerythrocytic parasites, the stage that causes malaria. Additionally, most methods are developed with either LC-MS or NMR analysis in mind, and have rarely been evaluated with both tools. In this work, three extraction methods frequently found in the literature were reproduced and samples were analyzed through both LC-MS and 1H NMR, and evaluated in order to reveal which is the most repeatable and consistent through an array of different tools, including chemometrics, peak detection and annotation. The most reliable method in this study proved to be a double extraction with methanol and methanol/water (80:20, v/v). Metabolomic studies in the field should move towards standardization of methodologies and the use of both LC-MS and 1H NMR in order to make data more comparable between studies and facilitate the achievement of biologically interpretable information.

Inhibitors of NAD+ Production in Cancer Treatment: State of the Art and Perspectives

The addiction of tumors to elevated nicotinamide adenine dinucleotide (NAD+) levels is a hallmark of cancer metabolism. Obstructing NAD+ biosynthesis in tumors is a new and promising antineoplastic strategy. Inhibitors developed against nicotinamide phosphoribosyltransferase (NAMPT), the main enzyme in NAD+ production from nicotinamide, elicited robust anticancer activity in preclinical models but not in patients, implying that other NAD+-biosynthetic pathways are also active in tumors and provide sufficient NAD+ amounts despite NAMPT obstruction. Recent studies show that NAD+ biosynthesis through the so-called "Preiss-Handler (PH) pathway", which utilizes nicotinate as a precursor, actively operates in many tumors and accounts for tumor resistance to NAMPT inhibitors. The PH pathway consists of three sequential enzymatic steps that are catalyzed by nicotinate phosphoribosyltransferase (NAPRT), nicotinamide mononucleotide adenylyltransferases (NMNATs), and NAD+ synthetase (NADSYN1). Here, we focus on these enzymes as emerging targets in cancer drug discovery, summarizing their reported inhibitors and describing their current or potential exploitation as anticancer agents. Finally, we also focus on additional NAD+-producing enzymes acting in alternative NAD+-producing routes that could also be relevant in tumors and thus become viable targets for drug discovery.

Genome-scale metabolic modeling and in silico analysis of opportunistic skin pathogen Cutibacterium acnes

Cutibacterium acnes, one of the most abundant skin microbes found in the sebaceous gland, is known to contribute to the development of acne vulgaris when its strains become imbalanced. The current limitations of acne treatment using antibiotics have caused an urgent need to develop a systematic strategy for selectively targeting C. acnes, which can be achieved by characterizing their cellular behaviors under various skin environments. To this end, we developed a genome-scale metabolic model (GEM) of virulent C. acnes, iCA843, based on the genome information of a relevant strain from ribotype 5 to comprehensively understand the pathogenic traits of C. acnes in the skin environment. We validated the model qualitatively by demonstrating its accuracy prediction of propionate and acetate production patterns, which were consistent with experimental observations. Additionally, we identified unique biosynthetic pathways for short-chain fatty acids in C. acnes compared to other GEMs of acne-inducing skin pathogens. By conducting constraint-based flux analysis under endogenous carbon sources in human skin, we discovered that the Wood-Werkman cycle is highly activated under acnes-associated skin condition for the regeneration of NAD, resulting in enhanced propionate production. Finally, we proposed potential anti-C. acnes targets by using the model-guided systematic framework based on gene essentiality analysis and protein sequence similarity search with abundant skin microbiome taxa.

Recent metabolomic developments for antimalarial drug discovery

Malaria is a parasitic disease that remains a global health issue, responsible for a significant death and morbidity toll. Various factors have impacted the use and delayed the development of antimalarial therapies, such as the associated financial cost and parasitic resistance. In order to discover new drugs and validate parasitic targets, a powerful omics tool, metabolomics, emerged as a reliable approach. However, as a fairly recent method in malaria, new findings are timely and original practices emerge frequently. This review aims to discuss recent research towards the development of new metabolomic methods in the context of uncovering antiplasmodial mechanisms of action in vitro and to point out innovative metabolic pathways that can revitalize the antimalarial pipeline.

Plasmodium falciparum Nicotinamidase as A Novel Antimalarial Target

Inhibition of Plasmodium falciparum nicotinamidase could represent a potential antimalarial since parasites require nicotinic acid to successfully recycle nicotinamide to NAD⁺, and importantly, humans lack this biosynthetic enzyme. Recently, mechanism-based inhibitors of nicotinamidase have been discovered. The most potent compound inhibits both recombinant P. falciparum nicotinamidase and parasites replication in infected human red blood cells (RBCs). These studies provide evidence for the importance of nicotinamide salvage through nicotinamidase as a central master player of NAD⁺ homeostasis in P. falciparum.

ELIMINATOR: essentiality analysis using multisystem networks and integer programming

A gene is considered as essential when it is indispensable for cells to grow and replicate in a certain environment. However, gene essentiality is not a structural property but rather a contextual one, which depends on the specific biological conditions affecting the cell. This circumstantial essentiality of genes is what brings the attention of scientist since we can identify genes essential for cancer cells but not essential for healthy cells. This same contextuality makes their identification extremely challenging. Huge experimental efforts such as Project Achilles where the essentiality of thousands of genes is measured together with a plethora of molecular data (transcriptomics, copy number, mutations, etc.) in over one thousand cell lines can shed light on the causality behind the essentiality of a gene in a given environment. Here, we present an in-silico method for the identification of patient-specific essential genes using constraint-based modelling (CBM). Our method expands the ideas behind traditional CBM to accommodate multisystem networks. In essence, it first calculates the minimum number of lowly expressed genes required to be activated by the cell to sustain life as defined by a set of requirements; and second, it performs an exhaustive in-silico gene knockout to find those that lead to the need of activating additional lowly expressed genes. We validated the proposed methodology using a set of 452 cancer cell lines derived from the Cancer Cell Line Encyclopedia where an exhaustive experimental large-scale gene knockout study using CRISPR (Achilles Project) evaluates the impact of each removal. We also show that the integration of different essentiality predictions per gene, what we called Essentiality Congruity Score, reduces the number of false positives. Finally, we explored our method in a breast cancer patient dataset, and our results showed high concordance with previous publications. These findings suggest that identifying genes whose activity is fundamental to sustain cellular life in a patient-specific manner is feasible using in-silico methods. The patient-level gene essentiality predictions can pave the way for precision medicine by identifying potential drug targets whose deletion can induce death in tumour cells.

Dual transcriptome based reconstruction of Salmonella-human integrated metabolic network to screen potential drug targets

Salmonella enterica serovar Typhimurium (S. Typhimurium) is a highly adaptive pathogenic bacteria with a serious public health concern due to its increasing resistance to antibiotics. Therefore, identification of novel drug targets for S. Typhimurium is crucial. Here, we first created a pathogen-host integrated genome-scale metabolic network by combining the metabolic models of human and S. Typhimurium, which we further tailored to the pathogenic state by the integration of dual transcriptome data. The integrated metabolic model enabled simultaneous investigation of metabolic alterations in human cells and S. Typhimurium during infection. Then, we used the tailored pathogen-host integrated genome-scale metabolic network to predict essential genes in the pathogen, which are candidate novel drug targets to inhibit infection. Drug target prioritization procedure was applied to these targets, and pabB was chosen as a putative drug target. It has an essential role in 4-aminobenzoic acid (PABA) synthesis, which is an essential biomolecule for many pathogens. A structure based virtual screening was applied through docking simulations to predict candidate compounds that eliminate S. Typhimurium infection by inhibiting pabB. To our knowledge, this is the first comprehensive study for predicting drug targets and drug like molecules by using pathogen-host integrated genome-scale models, dual RNA-seq data and structure-based virtual screening protocols. This framework will be useful in proposing novel drug targets and drugs for antibiotic-resistant pathogens.

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.

rMTA: robust metabolic transformation analysis

MOTIVATION

The development of computational tools exploiting -omics data and high-quality genome-scale metabolic networks for the identification of novel drug targets is a relevant topic in Systems Medicine. Metabolic Transformation Algorithm (MTA) is one of these tools, which aims to identify targets that transform a disease metabolic state back into a healthy state, with potential application in any disease where a clear metabolic alteration is observed.

RESULTS

Here, we present a robust extension to MTA (rMTA), which additionally incorporates a worst-case scenario analysis and minimization of metabolic adjustment to evaluate the beneficial effect of gene knockouts. We show that rMTA complements MTA in the different datasets analyzed (gene knockout perturbations in different organisms, Alzheimer's disease and prostate cancer), bringing a more accurate tool for predicting therapeutic targets.

AVAILABILITY AND IMPLEMENTATION

rMTA is freely available on The Cobra Toolbox: https://opencobra.github.io/cobratoolbox/latest/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Network-Based Metabolism-Centered Screening of Potential Drug Targets in Klebsiella pneumoniae at Genome Scale

Klebsiella pneumoniae is an opportunistic bacterial pathogen leading to life-threatening nosocomial infections. Emergence of highly resistant strains poses a major challenge in the management of the infections by healthcare-associated K. pneumoniae isolates. Thus, despite intensive efforts, the current treatment strategies remain insufficient to eradicate such infections. Failure of the conventional infection-prevention and treatment efforts explicitly indicates the requirement of new therapeutic approaches. This prompted us to systematically analyze the K. pneumoniae metabolism to investigate drug targets. Genome-scale metabolic networks (GMNs) facilitating the systematic analysis of the metabolism are promising platforms. Thus, we used a GMN of K. pneumoniae MGH 78578 to determine putative targets through gene- and metabolite-centric approaches. To develop more realistic infection models, we performed the bacterial growth simulations within different host-mimicking media, using an improved biomass formation reaction. We selected more suitable targets based on several property-based prioritization procedures. KdsA was identified as the high-ranked putative target satisfying most of the target prioritization criteria specified under the gene-centric approach. Through a structure-based virtual screening protocol, we identified potential KdsA inhibitors. In addition, the metabolite-centric approach extended the drug target list based on synthetic lethality. This revealed the importance of combined metabolic analyses for a better understanding of the metabolism. To our knowledge, this is the first comprehensive effort on the investigation of the K. pneumoniae metabolism for drug target prediction through the constraint-based analysis of its GMN in conjunction with several bioinformatic approaches. This study can guide the researchers for the future drug designs by providing initial findings regarding crucial components of the Klebsiella metabolism.

Rapid Discrimination of Malaria- and Dengue-Infected Patients Sera Using Raman Spectroscopy

Malaria and dengue have overlapping clinical symptoms and are prevalent in the same geographic region (tropical and subtropical), hence precise diagnosis is challenging. The high mortality rate associated with both malaria and dengue could be attributed to "false", "delayed", or "missed" diagnosis. The present study thus aims to stratify malaria and dengue using Raman spectroscopy (RS). In total, 130 human sera were analyzed for model development and double-blinded testing. Principal components linear discriminant analysis (PC-LDA) of acquired RS-spectra could classify malaria and dengue with a minor overlap of 16.7%. Receiver operating characteristic (ROC) analysis of test samples showed sensitivity/specificity of 0.9529 for malaria vs healthy controls (HC) and 0.9584 for dengue vs HC. The Raman findings were complemented by mass spectroscopy (MS)-based metabolite analysis of 8 individuals, each from malaria, dengue, and HC. Several of the metabolites, including amino acids, cell-free DNA, creatinine, and bilirubin, assigned for the predominant RS-bands were also identified by MS and showed similar trends. Our data clearly indicates that RS-based serum analysis using a microprobe has immense potential for early, accurate, and automated detection and discrimination of malaria and dengue, and in the future, it could be extrapolated in field-settings combined with hand-held RS. Further, this approach might be extended to diagnose other closely related infections with similar clinical manifestations.

Efflux pumps and antimicrobial resistance: Paradoxical components in systems genomics

Efflux pumps play a major role in the increasing antimicrobial resistance rendering a large number of drugs of no use. Large numbers of pathogens are becoming multidrug resistant due to inadequate dosage and use of the existing antimicrobials. This leads to the need for identifying new efflux pump inhibitors. Design of novel targeted therapies using inherent complexity involved in the biological network modeling has gained increasing importance in recent times. The predictive approaches should be used to determine antimicrobial activities with high pathogen specificity and microbicidal potency. Antimicrobial peptides, which are part of our innate immune system, have the ability to respond to infections and have gained much attention in making resistant strain sensitive to existing drugs. In this review paper, we outline evidences linking host-directed therapy with the efflux pump activity to infectious disease.

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.

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.

Functional interrogation of Plasmodium genus metabolism identifies species- and stage-specific differences in nutrient essentiality and drug targeting

Several antimalarial drugs exist, but differences between life cycle stages among malaria species pose challenges for developing more effective therapies. To understand the diversity among stages and species, we reconstructed genome-scale metabolic models (GeMMs) of metabolism for five life cycle stages and five species of Plasmodium spanning the blood, transmission, and mosquito stages. The stage-specific models of Plasmodium falciparum uncovered stage-dependent changes in central carbon metabolism and predicted potential targets that could affect several life cycle stages. The species-specific models further highlight differences between experimental animal models and the human-infecting species. Comparisons between human- and rodent-infecting species revealed differences in thiamine (vitamin B1), choline, and pantothenate (vitamin B5) metabolism. Thus, we show that genome-scale analysis of multiple stages and species of Plasmodium can prioritize potential drug targets that could be both anti-malarials and transmission blocking agents, in addition to guiding translation from non-human experimental disease models.

Malaria kills nearly one-half million people a year and over 1 billion people are at risk of becoming infected by the parasite. Plasmodial infections are difficult to treat for a myriad of reasons, but the ability of the organism to remain latent in hosts and the complex life cycles greatly contributed to the difficulty in treat malaria. Genome-scale metabolic models (GeMMs) enable hierarchical integration of disparate data types into a framework amenable to computational simulations enabling deeper mechanistic insights from high-throughput data measurements. In this study, GeMMs of multiple Plasmodium species are used to study metabolic similarities and differences across the Plasmodium genus. In silico gene-knock out simulations across species and stages uncovered functional metabolic differences between human- and rodent-infecting species as well as across the parasite’s life-cycle stages. These findings may help identify drug regimens that are more effective in targeting human-infecting species across multiple stages of the organism.

In silico engineering of Pseudomonas metabolism reveals new biomarkers for increased biosurfactant production

BACKGROUND

Rhamnolipids, biosurfactants with a wide range of biomedical applications, are amphiphilic molecules produced on the surfaces of or excreted extracellularly by bacteria including Pseudomonas aeruginosa. However, Pseudomonas putida is a non-pathogenic model organism with greater metabolic versatility and potential for industrial applications.

METHODS

We investigate in silico the metabolic capabilities of P. putida for rhamnolipids biosynthesis using statistical, metabolic and synthetic engineering approaches after introducing key genes (RhlA and RhlB) from P. aeruginosa into a genome-scale model of P. putida. This pipeline combines machine learning methods with multi-omic modelling, and drives the engineered P. putida model toward an optimal production and export of rhamnolipids out of the membrane.

RESULTS

We identify a substantial increase in synthesis of rhamnolipids by the engineered model compared to the control model. We apply statistical and machine learning techniques on the metabolic reaction rates to identify distinct features on the structure of the variables and individual components driving the variation of growth and rhamnolipids production. We finally provide a computational framework for integrating multi-omics data and identifying latent pathways and genes for the production of rhamnolipids in P. putida.

CONCLUSIONS

We anticipate that our results will provide a versatile methodology for integrating multi-omics data for topological and functional analysis of P. putida toward maximization of biosurfactant production.

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.

In-silico gene essentiality analysis of polyamine biosynthesis reveals APRT as a potential target in cancer

Constraint-based modeling for genome-scale metabolic networks has emerged in the last years as a promising approach to elucidate drug targets in cancer. Beyond the canonical biosynthetic routes to produce biomass, it is of key importance to focus on metabolic routes that sustain the proliferative capacity through the regulation of other biological means in order to improve in-silico gene essentiality analyses. Polyamines are polycations with central roles in cancer cell proliferation, through the regulation of transcription and translation among other things, but are typically neglected in in silico cancer metabolic models. In this study, we analysed essential genes for the biosynthesis of polyamines. Our analysis corroborates the importance of previously known regulators of the pathway, such as Adenosylmethionine Decarboxylase 1 (AMD1) and uncovers novel enzymes predicted to be relevant for polyamine homeostasis. We focused on Adenine Phosphoribosyltransferase (APRT) and demonstrated the detrimental consequence of APRT gene silencing on different leukaemia cell lines. Our results highlight the importance of revisiting the metabolic models used for in-silico gene essentiality analyses in order to maximize the potential for drug target identification in cancer.

Revealing the mystery of metabolic adaptations using a genome scale model of Leishmania infantum

Human macrophage phagolysosome and sandfly midgut provide antagonistic ecological niches for Leishmania parasites to survive and proliferate. Parasites optimize their metabolism to utilize the available inadequate resources by adapting to those environments. Lately, a number of metabolomics studies have revived the interest to understand metabolic strategies utilized by the Leishmania parasite for optimal survival within its hosts. For the first time, we propose a reconstructed genome-scale metabolic model for Leishmania infantum JPCM5, the analyses of which not only captures observations reported by metabolomics studies in other Leishmania species but also divulges novel features of the L. infantum metabolome. Our results indicate that Leishmania metabolism is organized in such a way that the parasite can select appropriate alternatives to compensate for limited external substrates. A dynamic non-essential amino acid motif exists within the network that promotes a restricted redistribution of resources to yield required essential metabolites. Further, subcellular compartments regulate this metabolic re-routing by reinforcing the physiological coupling of specific reactions. This unique metabolic organization is robust against accidental errors and provides a wide array of choices for the parasite to achieve optimal survival.

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.

Genomics of apicomplexan parasites

The increasing prevalence of infections involving intracellular apicomplexan parasites such as Plasmodium, Toxoplasma, and Cryptosporidium (the causative agents of malaria, toxoplasmosis, and cryptosporidiosis, respectively) represent a significant global healthcare burden. Despite their significance, few treatments are available; a situation that is likely to deteriorate with the emergence of new resistant strains of parasites. To lay the foundation for programs of drug discovery and vaccine development, genome sequences for many of these organisms have been generated, together with large-scale expression and proteomic datasets. Comparative analyses of these datasets are beginning to identify the molecular innovations supporting both conserved processes mediating fundamental roles in parasite survival and persistence, as well as lineage-specific adaptations associated with divergent life-cycle strategies. The challenge is how best to exploit these data to derive insights into parasite virulence and identify those genes representing the most amenable targets. In this review, we outline genomic datasets currently available for apicomplexans and discuss biological insights that have emerged as a consequence of their analysis. Of particular interest are systems-based resources, focusing on areas of metabolism and host invasion that are opening up opportunities for discovering new therapeutic targets.

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.

Biological characterization of chemically diverse compounds targeting the Plasmodium falciparum coenzyme A synthesis pathway

Background

In the fight against malaria, the discovery of chemical compounds with a novel mode of action and/or chemistry distinct from currently used drugs is vital to counteract the parasite’s known ability to develop drug resistance. Another desirable aspect is efficacy against gametocytes, the sexual developmental stage of the parasite which enables the transmission through Anopheles vectors. Using a chemical rescue approach, we previously identified compounds targeting Plasmodium falciparum coenzyme A (CoA) synthesis or utilization, a promising target that has not yet been exploited in anti-malarial drug development.

Results

We report on the outcomes of a series of biological tests that help to define the species- and stage-specificity, as well as the potential targets of these chemically diverse compounds. Compound activity against P. falciparum gametocytes was determined to assess stage-specificity and transmission-reducing potential. Against early stage gametocytes IC₅₀ values ranging between 60 nM and 7.5 μM were obtained. With the exception of two compounds with sub-micromolar potencies across all intra-erythrocytic stages, activity against late stage gametocytes was lower. None of the compounds were specific pantothenate kinase inhibitors. Chemical rescue profiling with CoA pathway intermediates demonstrated that most compounds acted on either of the two final P. falciparum CoA synthesis enzymes, phosphopantetheine adenylyltransferase (PPAT) or dephospho CoA kinase (DPCK). The most active compound targeted either phosphopantothenoylcysteine synthetase (PPCS) or phosphopantothenoylcysteine decarboxylase (PPCDC). Species-specificity was evaluated against Trypanosoma cruzi and Trypanosoma brucei brucei. No specific activity against T. cruzi amastigotes was observed; however three compounds inhibited the viability of trypomastigotes with sub-micromolar potencies and were confirmed to act on T. b. brucei CoA synthesis.

Conclusions

Utilizing the compounds we previously identified as effective against asexual P. falciparum, we demonstrate for the first time that gametocytes, like the asexual stages, depend on CoA, with two compounds exhibiting sub-micromolar potencies across asexual forms and all gametocytes stages tested. Furthermore, three compounds inhibited the viability of T. cruzi and T. b. brucei trypomastigotes with sub-micromolar potencies and were confirmed to act on T. b. brucei CoA synthesis, indicating that the CoA synthesis pathway might represent a valuable new drug target in these parasite species.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1860-3) contains supplementary material, which is available to authorized users.

Model-based transcriptome engineering promotes a fermentative transcriptional state in yeast

The ability to rationally manipulate the transcriptional states of cells would be of great use in medicine and bioengineering. We have developed an algorithm, NetSurgeon, which uses genome-wide gene-regulatory networks to identify interventions that force a cell toward a desired expression state. We first validated NetSurgeon extensively on existing datasets. Next, we used NetSurgeon to select transcription factor deletions aimed at improving ethanol production in Saccharomyces cerevisiae cultures that are catabolizing xylose. We reasoned that interventions that move the transcriptional state of cells using xylose toward that of cells producing large amounts of ethanol from glucose might improve xylose fermentation. Some of the interventions selected by NetSurgeon successfully promoted a fermentative transcriptional state in the absence of glucose, resulting in strains with a 2.7-fold increase in xylose import rates, a 4-fold improvement in xylose integration into central carbon metabolism, or a 1.3-fold increase in ethanol production rate. We conclude by presenting an integrated model of transcriptional regulation and metabolic flux that will enable future efforts aimed at improving xylose fermentation to prioritize functional regulators of central carbon metabolism.

Metabolomic Profiling of the Malaria Box Reveals Antimalarial Target Pathways

The threat of widespread drug resistance to frontline antimalarials has renewed the urgency for identifying inexpensive chemotherapeutic compounds that are effective against Plasmodium falciparum, the parasite species responsible for the greatest number of malaria-related deaths worldwide. To aid in the fight against malaria, a recent extensive screening campaign has generated thousands of lead compounds with low micromolar activity against blood stage parasites. A subset of these leads has been compiled by the Medicines for Malaria Venture (MMV) into a collection of structurally diverse compounds known as the MMV Malaria Box. Currently, little is known regarding the activity of these Malaria Box compounds on parasite metabolism during intraerythrocytic development, and a majority of the targets for these drugs have yet to be defined. Here we interrogated the in vitro metabolic effects of 189 drugs (including 169 of the drug-like compounds from the Malaria Box) using ultra-high-performance liquid chromatography-mass spectrometry (UHPLC-MS). The resulting metabolic fingerprints provide information on the parasite biochemical pathways affected by pharmacologic intervention and offer a critical blueprint for selecting and advancing lead compounds as next-generation antimalarial drugs. Our results reveal several major classes of metabolic disruption, which allow us to predict the mode of action (MoA) for many of the Malaria Box compounds. We anticipate that future combination therapies will be greatly informed by these results, allowing for the selection of appropriate drug combinations that simultaneously target multiple metabolic pathways, with the aim of eliminating malaria and forestalling the expansion of drug-resistant parasites in the field.

Efficient Reconstruction of Predictive Consensus Metabolic Network Models

Understanding cellular function requires accurate, comprehensive representations of metabolism. Genome-scale, constraint-based metabolic models (GSMs) provide such representations, but their usability is often hampered by inconsistencies at various levels, in particular for concurrent models. COMMGEN, our tool for COnsensus Metabolic Model GENeration, automatically identifies inconsistencies between concurrent models and semi-automatically resolves them, thereby contributing to consolidate knowledge of metabolic function. Tests of COMMGEN for four organisms showed that automatically generated consensus models were predictive and that they substantially increased coherence of knowledge representation. COMMGEN ought to be particularly useful for complex scenarios in which manual curation does not scale, such as for eukaryotic organisms, microbial communities, and host-pathogen interactions.

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.

PSAMM: A Portable System for the Analysis of Metabolic Models

The genome-scale models of metabolic networks have been broadly applied in phenotype prediction, evolutionary reconstruction, community functional analysis, and metabolic engineering. Despite the development of tools that support individual steps along the modeling procedure, it is still difficult to associate mathematical simulation results with the annotation and biological interpretation of metabolic models. In order to solve this problem, here we developed a Portable System for the Analysis of Metabolic Models (PSAMM), a new open-source software package that supports the integration of heterogeneous metadata in model annotations and provides a user-friendly interface for the analysis of metabolic models. PSAMM is independent of paid software environments like MATLAB, and all its dependencies are freely available for academic users. Compared to existing tools, PSAMM significantly reduced the running time of constraint-based analysis and enabled flexible settings of simulation parameters using simple one-line commands. The integration of heterogeneous, model-specific annotation information in PSAMM is achieved with a novel format of YAML-based model representation, which has several advantages, such as providing a modular organization of model components and simulation settings, enabling model version tracking, and permitting the integration of multiple simulation problems. PSAMM also includes a number of quality checking procedures to examine stoichiometric balance and to identify blocked reactions. Applying PSAMM to 57 models collected from current literature, we demonstrated how the software can be used for managing and simulating metabolic models. We identified a number of common inconsistencies in existing models and constructed an updated model repository to document the resolution of these inconsistencies.

The broad application of genome-scale metabolic modeling has made it a useful technique for tackling fundamental questions in biological research and engineering. Today over 100 models have been constructed for organisms that carry out a diverse array of metabolic activities spanning all three kingdoms of life. These models, however, have been curated independently following different conventions. The maintenance of model consistency has been challenging due to the lack of consensus in model representation and the absence of integrated modeling software for associating mathematical simulations with the annotation and biological interpretation of metabolic models. To solve this problem, we developed a new software package, PSAMM, and a new model format that incorporates heterogeneous, model-specific annotation information into modular representations of model definitions and simulation settings. PSAMM provides significant advances in standardizing the workflow of model annotation and consistency checking. Compared to existing tools, PSAMM supports more flexible configurations and is more efficient in running constraint-based simulations. All functions of PSAMM are freely available for academic users and can be downloaded from a public Git repository (https://zhanglab.github.io/psamm/) under the GNU General Public License.

Simulating Serial-Target Antibacterial Drug Synergies Using Flux Balance Analysis

Flux balance analysis (FBA) is an increasingly useful approach for modeling the behavior of metabolic systems. However, standard FBA modeling of genetic knockouts cannot predict drug combination synergies observed between serial metabolic targets, even though such synergies give rise to some of the most widely used antibiotic treatments. Here we extend FBA modeling to simulate responses to chemical inhibitors at varying concentrations, by diverting enzymatic flux to a waste reaction. This flux diversion yields very similar qualitative predictions to prior methods for single target activity. However, we find very different predictions for combinations, where flux diversion, which mimics the kinetics of competitive metabolic inhibitors, can explain serial target synergies between metabolic enzyme inhibitors that we confirmed in Escherichia coli cultures. FBA flux diversion opens the possibility for more accurate genome-scale predictions of drug synergies, which can be used to suggest treatments for infections and other diseases.

In silico multiple-targets identification for heme detoxification in the human malaria parasite Plasmodium falciparum

Detoxification of hemoglobin byproducts or free heme is an essential step and considered potential targets for anti-malaria drug development. However, most of anti-malaria drugs are no longer effective due to the emergence and spread of the drug resistant malaria parasites. Therefore, it is an urgent need to identify potential new targets and even for target combinations for effective malaria drug design. In this work, we reconstructed the metabolic networks of Plasmodium falciparum and human red blood cells for the simulation of steady mass and flux flows of the parasite's metabolites under the blood environment by flux balance analysis (FBA). The integrated model, namely iPF-RBC-713, was then adjusted into two stage-specific metabolic models, which first was for the pathological stage metabolic model of the parasite when invaded the red blood cell without any treatment and second was for the treatment stage of the parasite when a drug acted by inhibiting the hemozoin formation and caused high production rate of heme toxicity. The process of identifying target combinations consisted of two main steps. Firstly, the optimal fluxes of reactions in both the pathological and treatment stages were computed and compared to determine the change of fluxes. Corresponding enzymes of the reactions with zero fluxes in the treatment stage but non-zero fluxes in the pathological stage were predicted as a preliminary list of potential targets in inhibiting heme detoxification. Secondly, the combinations of all possible targets listed in the first step were examined to search for the best promising target combinations resulting in more effective inhibition of the detoxification to kill the malaria parasites. Finally, twenty-three enzymes were identified as a preliminary list of candidate targets which mostly were in pyruvate metabolism and citrate cycle. The optimal set of multiple targets for blocking the detoxification was a set of heme ligase, adenosine transporter, myo-inositol 1-phosphate synthase, ferrodoxim reductase-like protein and guanine transporter. In conclusion, the method has shown an effective and efficient way to identify target combinations which are obviously useful in the development of novel antimalarial drug combinations.

Quantitative analysis of drug effects at the whole-body level: a case study for glucose metabolism in malaria patients

We propose a hierarchical modelling approach to construct models for disease states at the whole-body level. Such models can simulate effects of drug-induced inhibition of reaction steps on the whole-body physiology. We illustrate the approach for glucose metabolism in malaria patients, by merging two detailed kinetic models for glucose metabolism in the parasite Plasmodium falciparum and the human red blood cell with a coarse-grained model for whole-body glucose metabolism. In addition we use a genome-scale metabolic model for the parasite to predict amino acid production profiles by the malaria parasite that can be used as a complex biomarker.

Metabolic modeling predicts metabolite changes in Mycobacterium tuberculosis

BACKGROUND

Mycobacterium tuberculosis (MTB) is the causal agent of the disease tuberculosis (TB). Metabolic adaptations are thought to be critical to the survival of MTB during pathogenesis. Computational tools that can be used to study MTB metabolism in silico and prioritize resource-intensive experimental work could significantly accelerate research.

RESULTS

We have developed E-Flux-MFC, an enhancement of our original E-Flux method that enables the prediction of changes in the production of external and internal metabolites corresponding to gene expression measurements. We have used this method to simulate the changes in the metabolic state of Mycobacterium tuberculosis (MTB). We have validated the accuracy of E-Flux-MFC for predicting changes in lipids and metabolites during a hypoxia time course using previously published metabolomics and transcriptomics data. We have further validated the accuracy of the method for predicting changes in MTB lipids following the deletion and induction of two well-studied transcription factors (TFs). We have applied the method to predict the metabolic impact of the induction of each of the approximately 180 MTB TFs using a previously generated and publically available expression data set.

CONCLUSIONS

E-flux-MFC can be used to study global changes in MTB metabolites from gene expression data associated with environmental and genetic perturbations. The application of this method to a data set of MTB TF perturbations provides a resource for studying the large number of TFs whose functions remain unknown. Most TFs impact metabolites indirectly through the propagation of gene expression changes through the regulatory network rather than through their direct regulons. E-Flux-MFC is also applicable to any organism for which accurate metabolic models are available.

Systems analysis of host-parasite interactions

Parasitic diseases caused by protozoan pathogens lead to hundreds of thousands of deaths per year in addition to substantial suffering and socioeconomic decline for millions of people worldwide. The lack of effective vaccines coupled with the widespread emergence of drug‐resistant parasites necessitates that the research community take an active role in understanding host–parasite infection biology in order to develop improved therapeutics. Recent advances in next‐generation sequencing and the rapid development of publicly accessible genomic databases for many human pathogens have facilitated the application of systems biology to the study of host–parasite interactions. Over the past decade, these technologies have led to the discovery of many important biological processes governing parasitic disease. The integration and interpretation of high‐throughput ‐omic data will undoubtedly generate extraordinary insight into host–parasite interaction networks essential to navigate the intricacies of these complex systems. As systems analysis continues to build the foundation for our understanding of host–parasite biology, this will provide the framework necessary to drive drug discovery research forward and accelerate the development of new antiparasitic therapies. WIREs Syst Biol Med 2015, 7:381–400. doi: 10.1002/wsbm.1311

For further resources related to this article, please visit the WIREs website.

Screening of potential targets in Plasmodium falciparum using stage-specific metabolic network analysis

The Apicomplexa parasite Plasmodium is a major cause of death in developing countries which are less equipped to bring new medicines to the market. Currently available drugs used for treatment of malaria are limited either by inadequate efficacy, toxicity and/or increased resistance. Availability of the genome sequence, microarray data and metabolic profile of Plasmodium parasite offers an opportunity for the identification of stage-specific genes important to the organism's lifecycle. In this study, microarray data were analysed for differential expression and overlapped onto metabolic pathways to identify differentially regulated pathways essential for transition to successive erythrocytic stages. The results obtained indicate that S-adenosylmethionine decarboxylase/ornithine decarboxylase, a bifunctional enzyme required for polyamine synthesis, is important for the Plasmodium cell growth in the absence of exogenous polyamines. S-adenosylmethionine decarboxylase/ornithine decarboxylase is a valuable target for designing therapeutically useful inhibitors. One such inhibitor, [Formula: see text]-difluoromethyl ornithine, is currently in use for the treatment of African sleeping sickness caused by Trypanosoma brucei. Structural studies of ornithine decarboxylase along with known inhibitors and their analogues were carried out to screen drug databases for more effective and less toxic compounds.

Fast-SL: an efficient algorithm to identify synthetic lethal sets in metabolic networks

MOTIVATION

Synthetic lethal sets are sets of reactions/genes where only the simultaneous removal of all reactions/genes in the set abolishes growth of an organism. Previous approaches to identify synthetic lethal genes in genome-scale metabolic networks have built on the framework of flux balance analysis (FBA), extending it either to exhaustively analyze all possible combinations of genes or formulate the problem as a bi-level mixed integer linear programming (MILP) problem. We here propose an algorithm, Fast-SL, which surmounts the computational complexity of previous approaches by iteratively reducing the search space for synthetic lethals, resulting in a substantial reduction in running time, even for higher order synthetic lethals.

RESULTS

We performed synthetic reaction and gene lethality analysis, using Fast-SL, for genome-scale metabolic networks of Escherichia coli, Salmonella enterica Typhimurium and Mycobacterium tuberculosis. Fast-SL also rigorously identifies synthetic lethal gene deletions, uncovering synthetic lethal triplets that were not reported previously. We confirm that the triple lethal gene sets obtained for the three organisms have a precise match with the results obtained through exhaustive enumeration of lethals performed on a computer cluster. We also parallelized our algorithm, enabling the identification of synthetic lethal gene quadruplets for all three organisms in under 6 h. Overall, Fast-SL enables an efficient enumeration of higher order synthetic lethals in metabolic networks, which may help uncover previously unknown genetic interactions and combinatorial drug targets.

AVAILABILITY AND IMPLEMENTATION

The MATLAB implementation of the algorithm, compatible with COBRA toolbox v2.0, is available at https://github.com/RamanLab/FastSL CONTACT: kraman@iitm.ac.in

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Predicting and exploring network components involved in pathogenesis in the malaria parasite via novel subnetwork alignments

BACKGROUND

Malaria is a major health threat, affecting over 40% of the world's population. The latest report released by the World Health Organization estimated about 207 million cases of malaria infection, and about 627,000 deaths in 2012 alone. During the past decade, new therapeutic targets have been identified and are at various stages of characterization, thanks to the emerging omics-based technologies. However, the mechanism of malaria pathogenesis remains largely unknown. In this paper, we employ a novel neighborhood subnetwork alignment approach to identify network components that are potentially involved in pathogenesis.

RESULTS

Our module-based subnetwork alignment approach identified 24 functional homologs of pathogenesis-related proteins in the malaria parasite P. falciparum, using the protein-protein interaction networks in Escherichia coli as references. Eighteen out of these 24 proteins are associated with 418 other proteins that are related to DNA replication, transcriptional regulation, translation, signaling, metabolism, cell cycle regulation, as well as cytoadherence and entry to the host.

CONCLUSIONS

The subnetwork alignments and subsequent protein-protein association network mining predicted a group of malarial proteins that may be involved in parasite development and parasite-host interaction, opening a new systems-level view of parasite pathogenesis and virulence.

Gene transfers shaped the evolution of de novo NAD+ biosynthesis in eukaryotes

NAD⁺ is an essential molecule for life, present in each living cell. It can function as an electron carrier or cofactor in redox biochemistry and energetics, and serves as substrate to generate the secondary messenger cyclic ADP ribose and nicotinic acid adenine dinucleotide phosphate. Although de novo NAD⁺ biosynthesis is essential, different metabolic pathways exist in different eukaryotic clades. The kynurenine pathway starting with tryptophan was most likely present in the last common ancestor of all eukaryotes, and is active in fungi and animals. The aspartate pathway, detected in most photosynthetic eukaryotes, was probably acquired from the cyanobacterial endosymbiont that gave rise to chloroplasts. An evolutionary analysis of enzymes catalyzing de novo NAD⁺ biosynthesis resulted in evolutionary trees incongruent with established organismal phylogeny, indicating numerous gene transfers. Endosymbiotic gene transfers probably introduced the aspartate pathway into eukaryotes and may have distributed it among different photosynthetic clades. In addition, several horizontal gene transfers substituted eukaryotic genes with bacterial orthologs. Although horizontal gene transfer is accepted as a key mechanism in prokaryotic evolution, it is supposed to be rare in eukaryotic evolution. The essential metabolic pathway of de novo NAD⁺ biosynthesis in eukaryotes was shaped by numerous gene transfers.

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.

Evolutionary Conservation of Bacterial Essential Metabolic Genes across All Bacterial Culture Media

One of the basic postulates of molecular evolution is that functionally important genes should evolve slower than genes of lesser significance. Essential genes, whose knockout leads to a lethal phenotype are considered of high functional importance, yet whether they are truly more conserved than nonessential genes has been the topic of much debate, fuelled by a host of contradictory findings. Here we conduct the first large-scale study utilizing genome-scale metabolic modeling and spanning many bacterial species, which aims to answer this question. Using the novel Media Variation Analysis, we examine the range of conservation of essential vs. nonessential metabolic genes in a given species across all possible media. We are thus able to obtain for the first time, exact upper and lower bounds on the levels of differential conservation of essential genes for each of the species studied. The results show that bacteria do exhibit an overall tendency for differential conservation of their essential genes vs. their non-essential ones, yet this tendency is highly variable across species. We show that the model bacterium E. coli K12 may or may not exhibit differential conservation of essential genes depending on its growth medium, shedding light on previous experimental studies showing opposite trends.

Predicting functional and regulatory divergence of a drug resistance transporter gene in the human malaria parasite

Background

The paradigm of resistance evolution to chemotherapeutic agents is that a key coding mutation in a specific gene drives resistance to a particular drug. In the case of resistance to the anti-malarial drug chloroquine (CQ), a specific mutation in the transporter pfcrt is associated with resistance. Here, we apply a series of analytical steps to gene expression data from our lab and leverage 3 independent datasets to identify pfcrt-interacting genes. Resulting networks provide insights into pfcrt’s biological functions and regulation, as well as the divergent phenotypic effects of its allelic variants in different genetic backgrounds.

Results

To identify pfcrt-interacting genes, we analyze pfcrt co-expression networks in 2 phenotypic states - CQ-resistant (CQR) and CQ-sensitive (CQS) recombinant progeny clones - using a computational approach that prioritizes gene interactions into functional and regulatory relationships. For both phenotypic states, pfcrt co-expressed gene sets are associated with hemoglobin metabolism, consistent with CQ’s expected mode of action. To predict the drivers of co-expression divergence, we integrate topological relationships in the co-expression networks with available high confidence protein-protein interaction data. This analysis identifies 3 transcriptional regulators from the ApiAP2 family and histone acetylation as potential mediators of these divergences. We validate the predicted divergences in DNA mismatch repair and histone acetylation by measuring the effects of small molecule inhibitors in recombinant progeny clones combined with quantitative trait locus (QTL) mapping.

Conclusions

This work demonstrates the utility of differential co-expression viewed in a network framework to uncover functional and regulatory divergence in phenotypically distinct parasites. pfcrt-associated co-expression in the CQ resistant progeny highlights CQR-specific gene relationships and possible targeted intervention strategies. The approaches outlined here can be readily generalized to other parasite populations and drug resistances.

Electronic supplementary material

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

Overcoming drug resistance through in silico prediction

Prediction tools are commonly used in pre-clinical research to assist target selection, to optimize drug potency or to predict the pharmacological profile of drug candidates. In silico prediction and overcoming drug resistance is a new opportunity that creates a high interest in pharmaceutical research. This review presents two main in silico strategies to meet this challenge: a structure-based approach to study the influence of mutations on the drug-target interaction and a system-biology approach to identify resistance pathways for a given drug. In silico screening of synergies between therapeutic and resistant pathways through biological network analysis is an example of technique to escape drug resistance. Structure-based drug design and in silico system biology are complementary approaches to reach few objectives at once: increase efficiency, reduce toxicity and overcoming drug resistance.

Computational prediction of essential metabolic genes using constraint-based approaches

In this chapter, we describe the application of constraint-based modeling to predict the impact of gene deletions on a metabolic phenotype. The metabolic reactions taking place inside cells form large networks, which have been reconstructed at a genome-scale for several organisms at increasing levels of detail. By integrating mathematical modeling techniques with biochemical principles, constraint-based approaches enable predictions of metabolite fluxes and growth under specific environmental conditions or for genetically modified microorganisms. Similar to the experimental knockout of a gene, predicting the essentiality of a metabolic gene for a phenotype further allows to generate hypotheses on its biological function and design of genetic engineering strategies for biotechnological applications. Here, we summarize the principles of constraint-based approaches and provide a detailed description of the procedure to predict the essentiality of metabolic genes with respect to a specific metabolic function. We exemplify the approach by predicting the essentiality of reactions in the citric acid cycle for the production of glucose from fatty acids.