An atypical orthologue of 6-pyruvoyltetrahydropterin synthase can provide the missing link in the folate biosynthesis pathway of malaria parasites

CRISPR-based functional profiling of the Toxoplasma gondii genome during acute murine infection

Examining host-pathogen interactions in animals can capture aspects of infection that are obscured in cell culture. Using CRISPR-based screens, we functionally profile the entire genome of the apicomplexan parasite Toxoplasma gondii during murine infection. Barcoded gRNAs enabled bottleneck detection and mapping of population structures within parasite lineages. Over 300 genes with previously unknown roles in infection were found to modulate parasite fitness in mice. Candidates span multiple axes of host-parasite interaction. Rhoptry Apical Surface Protein 1 was characterized as a mediator of host-cell tropism that facilitates repeated invasion attempts. GTP cyclohydrolase I was also required for fitness in mice and druggable through a repurposed compound, 2,4-diamino-6-hydroxypyrimidine. This compound synergized with pyrimethamine against T. gondii and malaria-causing Plasmodium falciparum parasites. This work represents a complete survey of an apicomplexan genome during infection of an animal host and points to novel interfaces of host-parasite interaction.

Computational Studies on 6-Pyruvoyl Tetrahydropterin Synthase (6-PTPS) of Plasmodium falciparum

6-Pyruvoyl tetrahydropterin synthase (6-PTPS) is a lyase involved in the synthesis of tetrahydrobiopterin. In Plasmodium species where dihydroneopterin aldolase (DHNA) is absent, it acts in the folate biosynthetic pathway necessary for the growth and survival of the parasite. The 6-pyruvoyl tetrahydropterin synthase of Plasmodium falciparum (PfPTPS) has been identified as a potential antimalarial drug target. This study identified potential inhibitors of PfPTPS using molecular docking techniques. Molecular docking and virtual screening of 62 compounds including the control to the deposited Protein Data Bank (PDB) structure was carried out using AutoDock Vina in PyRx. Five of the compounds, N,N-dimethyl-N'-[4-oxo-6-(2,2,5-trimethyl-1,3-dioxolan-4-yl)-3H-pteridin-2-yl]methanimidamide (140296439), 2-amino-6-[(1R)-3-cyclohexyl-1-hydroxypropyl]-3H-pteridin-4-one (140296495), 2-(2,3-dihydroxypropyl)-8,9-dihydro-6H-pyrimido[2,1-b]pteridine-7,11-dione (144380406), 2-(dimethylamino)-6-[(2,2-dimethyl-1,3-dioxolan-4-yl)-hydroxymethyl]-3H-pteridin-4-one (135573878), and [1-acetyloxy-1-(2-methyl-4-oxo-3H-pteridin-6-yl)propan-2-yl] acetate (136075207), showed better binding affinity than the control ligand, biopterin (135449517), and were selected and screened. Three conformers of 140296439 with the binding energy of -7.2, -7.1, and -7.0 kcal/mol along with 140296495 were better than the control at -5.7 kcal/mol. In silico absorption, distribution, metabolism, excretion, and toxicity (ADMET) studies predicted good pharmacokinetic properties of all the compounds while reporting a high risk of irritant toxicity in 140296439 and 144380406. The study highlights the five compounds, 140296439, 140296495, 144380406, 135573878 and 136075207, as potential inhibitors of PfPTPS and possible compounds for antimalarial drug development.

Functional profiling of the Toxoplasma genome during acute mouse infection

Within a host, pathogens encounter a diverse and changing landscape of cell types, nutrients, and immune responses. Examining host-pathogen interactions in animal models can therefore reveal aspects of infection absent from cell culture. We use CRISPR-based screens to functionally profile the entire genome of the model apicomplexan parasite Toxoplasma gondii during mouse infection. Barcoded gRNAs were used to track mutant parasite lineages, enabling detection of bottlenecks and mapping of population structures. We uncovered over 300 genes that modulate parasite fitness in mice with previously unknown roles in infection. These candidates span multiple axes of host-parasite interaction, including determinants of tropism, host organelle remodeling, and metabolic rewiring. We mechanistically characterized three novel candidates, including GTP cyclohydrolase I, against which a small-molecule inhibitor could be repurposed as an antiparasitic compound. This compound exhibited antiparasitic activity against T. gondii and Plasmodium falciparum, the most lethal agent of malaria. Taken together, we present the first complete survey of an apicomplexan genome during infection of an animal host, and point to novel interfaces of host-parasite interaction that may offer new avenues for treatment.

The genomic architecture of antimalarial drug resistance

Plasmodium falciparum and Plasmodium vivax, the two protozoan parasite species that cause the majority of cases of human malaria, have developed resistance to nearly all known antimalarials. The ability of malaria parasites to develop resistance is primarily due to the high numbers of parasites in the infected person's bloodstream during the asexual blood stage of infection in conjunction with the mutability of their genomes. Identifying the genetic mutations that mediate antimalarial resistance has deepened our understanding of how the parasites evade our treatments and reveals molecular markers that can be used to track the emergence of resistance in clinical samples. In this review, we examine known genetic mutations that lead to resistance to the major classes of antimalarial medications: the 4-aminoquinolines (chloroquine, amodiaquine and piperaquine), antifolate drugs, aryl amino-alcohols (quinine, lumefantrine and mefloquine), artemisinin compounds, antibiotics (clindamycin and doxycycline) and a napthoquinone (atovaquone). We discuss how the evolution of antimalarial resistance informs strategies to design the next generation of antimalarial therapies.

Folate biosynthesis pathway: mechanisms and insights into drug design for infectious diseases

Folate pathway is a key target for the development of new drugs against infectious diseases since the discovery of sulfa drugs and trimethoprim. The knowledge about this pathway has increased in the last years and the catalytic mechanism and structures of all enzymes of the pathway are fairly understood. In addition, differences among enzymes from prokaryotes and eukaryotes could be used for the design of specific inhibitors. In this review, we show a panorama of progress that has been achieved within the folate pathway obtained in the last years. We explored the structure and mechanism of enzymes, several genetic features, strategies, and approaches used in the design of new inhibitors that have been used as targets in pathogen chemotherapy.

Wholly Rickettsia! Reconstructed Metabolic Profile of the Quintessential Bacterial Parasite of Eukaryotic Cells

Reductive genome evolution has purged many metabolic pathways from obligate intracellular Rickettsia (Alphaproteobacteria; Rickettsiaceae). While some aspects of host-dependent rickettsial metabolism have been characterized, the array of host-acquired metabolites and their cognate transporters remains unknown. This dearth of information has thwarted efforts to obtain an axenic Rickettsia culture, a major impediment to conventional genetic approaches. Using phylogenomics and computational pathway analysis, we reconstructed the Rickettsia metabolic and transport network, identifying 51 host-acquired metabolites (only 21 previously characterized) needed to compensate for degraded biosynthesis pathways. In the absence of glycolysis and the pentose phosphate pathway, cell envelope glycoconjugates are synthesized from three imported host sugars, with a range of additional host-acquired metabolites fueling the tricarboxylic acid cycle. Fatty acid and glycerophospholipid pathways also initiate from host precursors, and import of both isoprenes and terpenoids is required for the synthesis of ubiquinone and the lipid carrier of lipid I and O-antigen. Unlike metabolite-provisioning bacterial symbionts of arthropods, rickettsiae cannot synthesize B vitamins or most other cofactors, accentuating their parasitic nature. Six biosynthesis pathways contain holes (missing enzymes); similar patterns in taxonomically diverse bacteria suggest alternative enzymes that await discovery. A paucity of characterized and predicted transporters emphasizes the knowledge gap concerning how rickettsiae import host metabolites, some of which are large and not known to be transported by bacteria. Collectively, our reconstructed metabolic network offers clues to how rickettsiae hijack host metabolic pathways. This blueprint for growth determinants is an important step toward the design of axenic media to rescue rickettsiae from the eukaryotic cell.IMPORTANCE A hallmark of obligate intracellular bacteria is the tradeoff of metabolic genes for the ability to acquire host metabolites. For species of Rickettsia, arthropod-borne parasites with the potential to cause serious human disease, the range of pilfered host metabolites is unknown. This information is critical for dissociating rickettsiae from eukaryotic cells to facilitate rickettsial genetic manipulation. In this study, we reconstructed the Rickettsia metabolic network and identified 51 host metabolites required to compensate patchwork Rickettsia biosynthesis pathways. Remarkably, some metabolites are not known to be transported by any bacteria, and overall, few cognate transporters were identified. Several pathways contain missing enzymes, yet similar pathways in unrelated bacteria indicate convergence and possible novel enzymes awaiting characterization. Our work illuminates the parasitic nature by which rickettsiae hijack host metabolism to counterbalance numerous disintegrated biosynthesis pathways that have arisen through evolution within the eukaryotic cell. This metabolic blueprint reveals what a Rickettsia axenic medium might entail.

Tetrahydrobiopterin Supplementation Improves Phenylalanine Metabolism in a Murine Model of Severe Malaria

Tetrahydrobiopterin (BH4) is an essential cofactor for both phenylalanine hydroxylase and nitric oxide synthase. Patients with severe malaria have low urinary BH4, elevated plasma phenylalanine, and impaired endothelium-dependent vasodilation, suggesting that BH4 depletion may limit phenylalanine metabolism and nitric oxide synthesis. We infected C57BL/6 mice with Plasmodium berghei ANKA to characterize BH4 availability and to investigate the effects of BH4 supplementation. P. berghei ANKA infection lowered BH4 in plasma, erythrocytes, and brain tissue but raised it in aorta and liver tissue. The ratio of BH4 to 7,8-BH2 (the major product of BH4 oxidation) was decreased in plasma, erythrocytes, and brain tissue, suggesting that oxidation contributes to BH4 depletion. The continuous infusion of sepiapterin (a BH4 precursor) and citrulline (an arginine precursor) raised the concentrations of BH4 and arginine in both blood and tissue compartments. The restoration of systemic BH4 and arginine availability in infected mice produced only a minor improvement in whole blood nitrite concentrations, a biomarker of NO synthesis, and failed to prevent the onset of severe disease symptoms. However, sepiapterin and citrulline infusion reduced the ratio of phenylalanine to tyrosine in plasma, aortic tissue, and brain tissue. In summary, BH4 depletion in P. berghei infection may compromise both nitric oxide synthesis and phenylalanine metabolism; however, these findings require further investigation in human patients with severe malaria.

Exceptionally long-range haplotypes in Plasmodium falciparum chromosome 6 maintained in an endemic African population

BACKGROUND

Previous genome-wide analyses of single nucleotide variation in Plasmodium falciparum identified evidence of an extended haplotype region on chromosome 6 in West Africa, suggesting recent positive selection. Such a pattern is not seen in samples from East Africa or South East Asia, so it could be marking a selective process specific to West Africa. Analyses of the haplotype structure in samples taken at different times could give clues to possible causes of selection.

METHODS

This study investigates chromosome 6 extended haplotypes in The Gambia by analysing alleles at multiple microsatellite loci using genome sequence data previously obtained from clinical isolates collected in 2008, followed by genotyping of 13 loci in 439 isolates from 1984, 1991, 2008 and 2014. Temporal changes in haplotype structure and frequencies were determined.

RESULTS

A region of high linkage disequilibrium spanning over 170 kilobases (kb) was identified with both NGS and laboratory determined microsatellite alleles. Multiple long haplotypes were found in all temporal populations from The Gambia. Two of the haplotypes were detected in samples from 1984 and 1991. The frequency of long-range haplotypes increased in 2008 and 2014 populations. There was higher Fst between older and more recent populations at loci in proximity to genes involved in drug metabolism pathways.

CONCLUSIONS

The occurrence of several long haplotypes at intermediate frequencies suggests an unusual mode of selection in chromosome 6, possibly combined with recombination suppression on specific haplotypes. Such selection apparently occurred before the emergence of known anti-malarial drug resistance alleles, and could be due to effects of other drugs or unknown processes that have long been operating in this endemic region.

Functional Annotations of Paralogs: A Blessing and a Curse

Gene duplication followed by mutation is a classic mechanism of neofunctionalization, producing gene families with functional diversity. In some cases, a single point mutation is sufficient to change the substrate specificity and/or the chemistry performed by an enzyme, making it difficult to accurately separate enzymes with identical functions from homologs with different functions. Because sequence similarity is often used as a basis for assigning functional annotations to genes, non-isofunctional gene families pose a great challenge for genome annotation pipelines. Here we describe how integrating evolutionary and functional information such as genome context, phylogeny, metabolic reconstruction and signature motifs may be required to correctly annotate multifunctional families. These integrative analyses can also lead to the discovery of novel gene functions, as hints from specific subgroups can guide the functional characterization of other members of the family. We demonstrate how careful manual curation processes using comparative genomics can disambiguate subgroups within large multifunctional families and discover their functions. We present the COG0720 protein family as a case study. We also discuss strategies to automate this process to improve the accuracy of genome functional annotation pipelines.

Possible selection of host folate pathway gene polymorphisms in patients with malaria from a malaria endemic region in North East India

BACKGROUND

Recent studies in experimental mice have shown that mild deficiency of methylenetetrahydrofolate reductase (MTHFR) enzyme confers protection against malaria, thus providing an important basis for the hypothesis that MTHFR polymorphism, i.e. C677T, might have been subjected to selection pressure against malaria. The present study was undertaken in a malaria endemic region in North East India to assess whether a similar selection advantage exists for other genes in folate metabolism pathway.

METHODS

A total of 401 subjects including 131 symptomatic malaria, 97 asymptomatic malaria and 173 normal healthy controls were analysed for nine polymorphisms (single-nucleotide polymorphisms [SNPs] in eight genes and insertion/deletion in one gene): MTHFR C677T, methionine synthase reductase (MTRR) A66G, glutamate carboxypeptidase II (GCPII) C1561T, cystathionine beta-synthase (CBS) 844ins68, reduced folate carrier-1 (RFC-1) G80A, serine hydroxymethyltransferase (SHMT) C1420T, methionine synthase (MTR) A2756G, MTHFR G1793A (D 919G), glycine N-methyltransferase (GNMT) 1289 by PCR-RFLP technique. Differences in frequencies of genotype distribution of each polymorphic marker between these groups were evaluated.

RESULTS

MTRR A2756G, SHMT C1420T, GCPII C1561T, MTRR A2756G and GNMT C1289T and RFC1 G80A polymorphisms showed significantly different prevalence between different groups analyzed. No significant differences were seen in the distribution of other polymorphisms.

CONCLUSIONS

The study gives a clue for the possible selection of specific polymorphisms in the genes involved in the folate metabolism pathway by malaria parasite.

MRP1 mediates folate transport and antifolate sensitivity in Plasmodium falciparum

Multidrug resistance-associated proteins (MRP) of Plasmodium falciparum have been associated with altered drug sensitivity. Knowledge on MRP substrate specificity is indispensible for the characterization of resistance mechanisms and identifying its physiological roles. An untargeted metabolomics approach detected decreased folate concentrations in red blood cells infected with schizont stage parasites lacking expression of MRP1. Furthermore, a tenfold decrease in sensitivity toward the folate analog methotrexate was detected for parasites lacking MRP1. PfMRP1 is involved in the export of folate from parasites into red blood cells and is therefore a relevant factor for efficient malaria treatment through the folate pathway.

A mathematical model of microbial folate biosynthesis and utilisation: implications for antifolate development

The metabolic biochemistry of folate biosynthesis and utilisation has evolved into a complex network of reactions. Although this complexity represents challenges to the field of folate research it has also provided a renewed source for antimetabolite targets. A range of improved folate chemotherapy continues to be developed and applied particularly to cancer and chronic inflammatory diseases. However, new or better antifolates against infectious diseases remain much more elusive. In this paper we describe the assembly of a generic deterministic mathematical model of microbial folate metabolism. Our aim is to explore how a mathematical model could be used to explore the dynamics of this inherently complex set of biochemical reactions. Using the model it was found that: (1) a particular small set of folate intermediates are overrepresented, (2) inhibitory profiles can be quantified by the level of key folate products, (3) using the model to scan for the most effective combinatorial inhibitions of folate enzymes we identified specific targets which could complement current antifolates, and (4) the model substantiates the case for a substrate cycle in the folinic acid biosynthesis reaction. Our model is coded in the systems biology markup language and has been deposited in the BioModels Database (MODEL1511020000), this makes it accessible to the community as a whole.

The Rickettsia Endosymbiont of Ixodes pacificus Contains All the Genes of De Novo Folate Biosynthesis

Ticks and other arthropods often are hosts to nutrient providing bacterial endosymbionts, which contribute to their host’s fitness by supplying nutrients such as vitamins and amino acids. It has been detected, in our lab, that Ixodes pacificus is host to Rickettsia species phylotype G021. This endosymbiont is predominantly present, and 100% maternally transmitted in I. pacificus. To study roles of phylotype G021 in I. pacificus, bioinformatic and molecular approaches were carried out. MUMmer genome alignments of whole genome sequence of I. scapularis, a close relative to I. pacificus, against completely sequenced genomes of R. bellii OSU85-389, R. conorii, and R. felis, identified 8,190 unique sequences that are homologous to Rickettsia sequences in the NCBI Trace Archive. MetaCyc metabolic reconstructions revealed that all folate gene orthologues (folA, folC, folE, folKP, ptpS) required for de novo folate biosynthesis are present in the genome of Rickettsia buchneri in I. scapularis. To examine the metabolic capability of phylotype G021 in I. pacificus, genes of the folate biosynthesis pathway of the bacterium were PCR amplified using degenerate primers. BLAST searches identified that nucleotide sequences of the folA, folC, folE, folKP, and ptpS genes possess 98.6%, 98.8%, 98.9%, 98.5% and 99.0% identity respectively to the corresponding genes of Rickettsia buchneri. Phylogenetic tree constructions show that the folate genes of phylotype G021 and homologous genes from various Rickettsia species are monophyletic. This study has shown that all folate genes exist in the genome of Rickettsia species phylotype G021 and that this bacterium has the genetic capability for de novo folate synthesis.

Structural basis of a novel activity of bacterial 6-pyruvoyltetrahydropterin synthase homologues distinct from mammalian 6-pyruvoyltetrahydropterin synthase activity

Escherichia coli 6-carboxytetrahydropterin synthase (eCTPS), a homologue of 6-pyruvoyltetrahydropterin synthase (PTPS), possesses a much stronger catalytic activity to cleave the side chain of sepiapterin in vitro compared with genuine PTPS activity and catalyzes the conversion of dihydroneopterin triphosphate to 6-carboxy-5,6,7,8-tetrahydropterin in vivo. Crystal structures of wild-type apo eCTPS and of a Cys27Ala mutant eCTPS complexed with sepiapterin have been determined to 2.3 and 2.5 Å resolution, respectively. The structures are highly conserved at the active site and the Zn(2+) binding site. However, comparison of the eCTPS structures with those of mammalian PTPS homologues revealed that two specific residues, Trp51 and Phe55, that are not found in mammalian PTPS keep the substrate bound by stacking it with their side chains. Replacement of these two residues by site-directed mutagenesis to the residues Met and Leu, which are only found in mammalian PTPS, converted eCTPS to the mammalian PTPS activity. These studies confirm that these two aromatic residues in eCTPS play an essential role in stabilizing the substrate and in the specific enzyme activity that differs from the original PTPS activity. These aromatic residues Trp51 and Phe55 are a key signature of bacterial PTPS enzymes that distinguish them from mammalian PTPS homologues.

Transition-state inhibitors of purine salvage and other prospective enzyme targets in malaria

Malaria is a leading cause of human death within the tropics. The gradual generation of drug resistance imposes an urgent need for the development of new and selective antimalarial agents. Kinetic isotope effects coupled to computational chemistry have provided the relevant details on geometry and charge of enzymatic transition states to facilitate the design of transition-state analogs. These features have been reproduced into chemically stable mimics through synthetic chemistry, generating inhibitors with dissociation constants in the pico- to femto-molar range. Transition-state analogs are expected to contribute to the control of malaria.

Biochemical fossils of the ancient transition from geoenergetics to bioenergetics in prokaryotic one carbon compound metabolism

The deep dichotomy of archaea and bacteria is evident in many basic traits including ribosomal protein composition, membrane lipid synthesis, cell wall constituents, and flagellar composition. Here we explore that deep dichotomy further by examining the distribution of genes for the synthesis of the central carriers of one carbon units, tetrahydrofolate (H4F) and tetrahydromethanopterin (H4MPT), in bacteria and archaea. The enzymes underlying those distinct biosynthetic routes are broadly unrelated across the bacterial-archaeal divide, indicating that the corresponding pathways arose independently. That deep divergence in one carbon metabolism is mirrored in the structurally unrelated enzymes and different organic cofactors that methanogens (archaea) and acetogens (bacteria) use to perform methyl synthesis in their H4F- and H4MPT-dependent versions, respectively, of the acetyl-CoA pathway. By contrast, acetyl synthesis in the acetyl-CoA pathway - from a methyl group, CO2 and reduced ferredoxin - is simpler, uniform and conserved across acetogens and methanogens, and involves only transition metals as catalysts. The data suggest that the acetyl-CoA pathway, while being the most ancient of known CO2 assimilation pathways, reflects two phases in early evolution: an ancient phase in a geochemically confined and non-free-living universal common ancestor, in which acetyl thioester synthesis proceeded spontaneously with the help of geochemically supplied methyl groups, and a later phase that reflects the primordial divergence of the bacterial and archaeal stem groups, which independently invented genetically-encoded means to synthesize methyl groups via enzymatic reactions. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference.

Utility of the Biosynthetic Folate Pathway for Targets in Antimicrobial Discovery

The need for new antimicrobials is great in face of a growing pool of resistant pathogenic organisms. This review will address the potential for antimicrobial therapy based on polypharmacological activities within the currently utilized bacterial biosynthetic folate pathway. The folate metabolic pathway leads to synthesis of required precursors for cellular function and contains a critical node, dihydrofolate reductase (DHFR), which is shared between prokaryotes and eukaryotes. The DHFR enzyme is currently targeted by methotrexate in anti-cancer therapies, by trimethoprim for antibacterial uses, and by pyrimethamine for anti-protozoal applications. An additional anti-folate target is dihyropteroate synthase (DHPS), which is unique to prokaryotes as they cannot acquire folate through dietary means. It has been demonstrated as a primary target for the longest standing antibiotic class, the sulfonamides, which act synergistically with DHFR inhibitors. Investigations have revealed most DHPS enzymes possess the ability to utilize sulfa drugs metabolically, producing alternate products that presumably inhibit downstream enzymes requiring the produced dihydropteroate. Recent work has established an off-target effect of sulfonamide antibiotics on a eukaryotic enzyme, sepiapterin reductase, causing alterations in neurotransmitter synthesis. Given that inhibitors of both DHFR and DHPS are designed to mimic their cognate substrate, which contain shared substructures, it is reasonable to expect such “off-target” effects. These inhibitors are also likely to interact with the enzymatic neighbors in the folate pathway that bind products of the DHFR or DHPS enzymes and/or substrates of similar substructure. Computational studies designed to assess polypharmacology reiterate these conclusions. This leads to hypotheses exploring the vast utility of multiple members of the folate pathway for modulating cellular metabolism, and includes an appealing capacity for prokaryotic-specific polypharmacology for antimicrobial applications.

Simple topological properties predict functional misannotations in a metabolic network

Motivation: Misannotation in sequence databases is an important obstacle for automated tools for gene function annotation, which rely extensively on comparison with sequences with known function. To improve current annotations and prevent future propagation of errors, sequence-independent tools are, therefore, needed to assist in the identification of misannotated gene products. In the case of enzymatic functions, each functional assignment implies the existence of a reaction within the organism’s metabolic network; a first approximation to a genome-scale metabolic model can be obtained directly from an automated genome annotation. Any obvious problems in the network, such as dead end or disconnected reactions, can, therefore, be strong indications of misannotation.

Results: We demonstrate that a machine-learning approach using only network topological features can successfully predict the validity of enzyme annotations. The predictions are tested at three different levels. A random forest using topological features of the metabolic network and trained on curated sets of correct and incorrect enzyme assignments was found to have an accuracy of up to 86% in 5-fold cross-validation experiments. Further cross-validation against unseen enzyme superfamilies indicates that this classifier can successfully extrapolate beyond the classes of enzyme present in the training data. The random forest model was applied to several automated genome annotations, achieving an accuracy of in most cases when validated against recent genome-scale metabolic models. We also observe that when applied to draft metabolic networks for multiple species, a clear negative correlation is observed between predicted annotation quality and phylogenetic distance to the major model organism for biochemistry (Escherichia coli for prokaryotes and Homo sapiens for eukaryotes).

Contact:j.pinney@imperial.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Targeting the vitamin biosynthesis pathways for the treatment of malaria

The most severe form of malaria is Malaria tropica, caused by Plasmodium falciparum. There are more than 1 billion people that are exposed to malaria parasites leading to more than 500,000 deaths annually. Vaccines are not available and the increasing drug resistance of the parasite prioritizes the need for novel drug targets and chemotherapeutics, which should be ideally designed to target selectively the parasite. In this sense, parasite-specific pathways, such as the vitamin biosyntheses, represent perfect drug-target characteristics because they are absent in humans. In the past, the vitamin B9 (folate) metabolism has been exploited by antifolates to treat infections caused by malaria parasites. Recently, two further vitamin biosynthesis pathways - for the vitamins B6 (pyridoxine) and B1 (thiamine) - have been identified in Plasmodium and analyzed for their suitability to discover new drugs. In this review, the current status of the druggability of plasmodial vitamin biosynthesis pathways is summarized.

Organic cofactors in the metabolism of Dehalococcoides mccartyi strains

Dehalococcoides mccartyi strains are strictly anaerobic organisms specialized to grow with halogenated compounds as electron acceptor via a respiratory process. Their genomes are among the smallest known for free-living organisms, and the embedded gene set reflects their strong specialization. Here, we briefly review main characteristics of published Dehalococcoides genomes and show how genome information together with cultivation and biochemical experiments have contributed to our understanding of Dehalococcoides physiology and biochemistry. We extend this approach by the detailed analysis of cofactor metabolism in Dehalococcoides strain CBDB1. Dehalococcoides genomes were screened for encoded proteins annotated to contain or interact with organic cofactors, and the expression of these proteins was analysed by shotgun proteomics to shed light on cofactor requirements. In parallel, cultivation experiments testing for vitamin requirements showed that cyanocobalamin (vitamin B12), thiamine and biotin were essential supplements and that cyanocobalamin could be substituted by dicyanocobinamide and dimethylbenzimidazole. Dehalococcoides genome analysis, detection of single enzymes by shotgun proteomics and inhibition studies confirmed the expression of the biosynthetic pathways for pyridoxal-5-phosphate, flavin nucleotides, folate, S-adenosylmethionine, pantothenate and nicotinic acids in strain CBDB1. Haem/cytochromes, quinones and lipoic acids were not necessary for cultivation or dechlorination activity and no biosynthetic pathways were identified in the genomes.

The folate metabolic network of Falciparum malaria

The targeting of key enzymes in the folate pathway continues to be an effective chemotherapeutic approach that has earned antifolate drugs a valuable position in the medical pharmacopoeia. The successful therapeutic use of antifolates as antimalarials has been a catalyst for ongoing research into the biochemistry of folate and pterin biosynthesis in malaria parasites. However, our understanding of the parasites folate metabolism remains partial and patchy, especially in relation to the shikimate pathway, the folate cycle, and folate salvage. A sizeable number of potential folate targets remain to be characterised. Recent reports on the parasite specific transport of folate precursors that would normally be present in the human host awaken previous hypotheses on the salvage of folate precursors or by-products. As the parasite progresses through its life-cycle it encounters very contrasting host cell environments that present radically different metabolic milieus and biochemical challenges. It would seem probable that as the parasite encounters differing environments it would need to modify its biochemistry. This would be reflected in the folate homeostasis in Plasmodium. Recent drug screening efforts and insights into folate membrane transport substantiate the argument that folate metabolism may still offer unexplored opportunities for therapeutic attack.

Comparative genomics guided discovery of two missing archaeal enzyme families involved in the biosynthesis of the pterin moiety of tetrahydromethanopterin and tetrahydrofolate

C-1 carriers are essential cofactors in all domains of life, and in Archaea, these can be derivatives of tetrahydromethanopterin (H₄-MPT) or tetrahydrofolate (H₄-folate). Their synthesis requires 6-hydroxymethyl-7,8-dihydropterin diphosphate (6-HMDP) as the precursor, but the nature of pathways that lead to its formation were unknown until the recent discovery of the GTP cyclohydrolase IB/MptA family that catalyzes the first step, the conversion of GTP to dihydroneopterin 2′,3′-cyclic phosphate or 7,8-dihydroneopterin triphosphate [El Yacoubi, B.; et al. (2006) J. Biol. Chem., 281, 37586–37593 and Grochowski, L. L.; et al. (2007) Biochemistry46, 6658–6667]. Using a combination of comparative genomics analyses, heterologous complementation tests, and in vitro assays, we show that the archaeal protein families COG2098 and COG1634 specify two of the missing 6-HMDP synthesis enzymes. Members of the COG2098 family catalyze the formation of 6-hydroxymethyl-7,8-dihydropterin from 7,8-dihydroneopterin, while members of the COG1634 family catalyze the formation of 6-HMDP from 6-hydroxymethyl-7,8-dihydropterin. The discovery of these missing genes solves a long-standing mystery and provides novel examples of convergent evolutions where proteins of dissimilar architectures perform the same biochemical function.

Functional promiscuity of the COG0720 family

The biosynthesis of GTP derived metabolites such as tetrahydrofolate (THF), biopterin (BH(4)), and the modified tRNA nucleosides queuosine (Q) and archaeosine (G(+)) relies on several enzymes of the Tunnel-fold superfamily. A subset of these proteins includes the 6-pyruvoyltetrahydropterin (PTPS-II), PTPS-III, and PTPS-I homologues, all members of the COG0720 family that have been previously shown to transform 7,8-dihydroneopterin triphosphate (H(2)NTP) into different products. PTPS-II catalyzes the formation of 6-pyruvoyltetrahydropterin in the BH(4) pathway, PTPS-III catalyzes the formation of 6-hydroxylmethyl-7,8-dihydropterin in the THF pathway, and PTPS-I catalyzes the formation of 6-carboxy-5,6,7,8-tetrahydropterin in the Q pathway. Genes of these three enzyme families are often misannotated as they are difficult to differentiate by sequence similarity alone. Using a combination of physical clustering, signature motif, phylogenetic codistribution analyses, in vivo complementation studies, and in vitro enzymatic assays, a complete reannotation of the COG0720 family was performed in prokaryotes. Notably, this work identified and experimentally validated dual function PTPS-I/III enzymes involved in both THF and Q biosynthesis. Both in vivo and in vitro analyses showed that the PTPS-I family could tolerate a translation of the active site cysteine and was inherently promiscuous, catalyzing different reactions on the same substrate or the same reaction on different substrates. Finally, the analysis and experimental validation of several archaeal COG0720 members confirmed the role of PTPS-I in archaeosine biosynthesis and resulted in the identification of PTPS-III enzymes with variant signature sequences in Sulfolobus species. This study reveals an expanded versatility of the COG0720 family members and illustrates that for certain protein families extensive comparative genomic analysis beyond homology is required to correctly predict function.

Quantitative mass spectrometry catalogues Salmonella pathogenicity island-2 effectors and identifies their cognate host binding partners

Gram-negative bacterial pathogens have developed specialized secretion systems to transfer bacterial proteins directly into host cells. These bacterial effectors are central to virulence and reprogram host cell processes to favor bacterial survival, colonization, and proliferation. Knowing the complete set of effectors encoded by a particular pathogen is the key to understanding bacterial disease. In addition, the identification of the molecular assemblies that these effectors engage once inside the host cell is critical to determining the mechanism of action of each effector. In this work we used stable isotope labeling of amino acids in cell culture (SILAC), a powerful quantitative proteomics technique, to identify the proteins secreted by the Salmonella pathogenicity island-2 type three secretion system (SPI-2 T3SS) and to characterize the host interaction partners of SPI-2 effectors. We confirmed many of the known SPI-2 effectors and were able to identify several novel substrate candidates of this secretion system. We verified previously published host protein-effector binding pairs and obtained 11 novel interactions, three of which were investigated further and confirmed by reciprocal co-immunoprecipitation. The host cell interaction partners identified here suggest that Salmonella SPI-2 effectors target, in a concerted fashion, cellular processes such as cell attachment and cell cycle control that are underappreciated in the context of infection. The technology outlined in this study is specific and sensitive and serves as a robust tool for the identification of effectors and their host targets that is readily amenable to the study of other bacterial pathogens.

The complete genome sequence of 'Candidatus Liberibacter solanacearum', the bacterium associated with potato zebra chip disease

Zebra Chip (ZC) is an emerging plant disease that causes aboveground decline of potato shoots and generally results in unusable tubers. This disease has led to multi-million dollar losses for growers in the central and western United States over the past decade and impacts the livelihood of potato farmers in Mexico and New Zealand. ZC is associated with 'Candidatus Liberibacter solanacearum', a fastidious alpha-proteobacterium that is transmitted by a phloem-feeding psyllid vector, Bactericera cockerelli Sulc. Research on this disease has been hampered by a lack of robust culture methods and paucity of genome sequence information for 'Ca. L. solanacearum'. Here we present the sequence of the 1.26 Mbp metagenome of 'Ca. L. solanacearum', based on DNA isolated from potato psyllids. The coding inventory of the 'Ca. L. solanacearum' genome was analyzed and compared to related Rhizobiaceae to better understand 'Ca. L. solanacearum' physiology and identify potential targets to develop improved treatment strategies. This analysis revealed a number of unique transporters and pathways, all potentially contributing to ZC pathogenesis. Some of these factors may have been acquired through horizontal gene transfer. Taxonomically, 'Ca. L. solanacearum' is related to 'Ca. L. asiaticus', a suspected causative agent of citrus huanglongbing, yet many genome rearrangements and several gene gains/losses are evident when comparing these two Liberibacter. species. Relative to 'Ca. L. asiaticus', 'Ca. L. solanacearum' probably has reduced capacity for nucleic acid modification, increased amino acid and vitamin biosynthesis functionalities, and gained a high-affinity iron transport system characteristic of several pathogenic microbes.

Use of thermal melt curves to assess the quality of enzyme preparations

This study sought to determine whether the quality of enzyme preparations can be determined from their melting curves, which may easily be obtained using a fluorescent probe and a standard reverse transcription-polymerase chain reaction (RT-PCR) machine. Thermal melt data on 31 recombinant enzymes from Plasmodium parasites were acquired by incrementally heating them to 90 degrees C and measuring unfolding with a fluorescent dye. Activity assays specific to each enzyme were also performed. Four of the enzymes were denatured to varying degrees with heat and sodium dodecyl sulfate (SDS) prior to the thermal melt and activity assays. In general, melting curve quality was correlated with enzyme activity; enzymes with high-quality curves were found almost uniformly to be active, whereas those with lower quality curves were more varied in their catalytic performance. Inspection of melting curves of bovine xanthine oxidase and Entamoeba histolytica cysteine protease 1 allowed active stocks to be distinguished from inactive stocks, implying that a relationship between melting curve quality and activity persists over a wide range of experimental conditions and species. Our data suggest that melting curves can help to distinguish properly folded proteins from denatured ones and, therefore, may be useful in selecting stocks for further study and in optimizing purification procedures for specific proteins.

Vitamin B metabolism in Plasmodium falciparum as a source of drug targets

The malaria parasite Plasmodium falciparum depends primarily on nutrient sources from its human host. Most compounds, such as glucose, purines, amino acids, as well as cofactors and vitamins, are abundantly available in the host cell, and can be readily salvaged by the parasite. However, in some cases the parasite can also synthesize cofactors de novo in reactions that appear to be essential. Importantly, the three biosynthetic pathways that produce vitamins B(1), B(6) and B(9) are absent from the host, but are well established in P. falciparum. This review summarizes and updates the current knowledge of vitamin B de novo synthesis and salvage in P. falciparum and focuses on their potential as targets for drug intervention.

6-pyruvoyltetrahydropterin synthase paralogs replace the folate synthesis enzyme dihydroneopterin aldolase in diverse bacteria

Dihydroneopterin aldolase (FolB) catalyzes conversion of dihydroneopterin to 6-hydroxymethyldihydropterin (HMDHP) in the classical folate biosynthesis pathway. However, folB genes are missing from the genomes of certain bacteria from the phyla Chloroflexi, Acidobacteria, Firmicutes, Planctomycetes, and Spirochaetes. Almost all of these folB-deficient genomes contain an unusual paralog of the tetrahydrobiopterin synthesis enzyme 6-pyruvoyltetrahydropterin synthase (PTPS) in which a glutamate residue replaces or accompanies the catalytic cysteine. A similar PTPS paralog from the malaria parasite Plasmodium falciparum is known to form HMDHP from dihydroneopterin triphosphate in vitro and has been proposed to provide a bypass to the FolB step in vivo. Bacterial genes encoding PTPS-like proteins with active-site glutamate, cysteine, or both residues were accordingly tested together with the P. falciparum gene for complementation of the Escherichia coli folB mutation. The P. falciparum sequence and bacterial sequences with glutamate or glutamate plus cysteine were active; those with cysteine alone were not. These results demonstrate that PTPS paralogs with an active-site glutamate (designated PTPS-III proteins) can functionally replace FolB in vivo. Recombinant bacterial PTPS-III proteins, like the P. falciparum enzyme, mediated conversion of dihydroneopterin triphosphate to HMDHP, but other PTPS proteins did not. Neither PTPS-III nor other PTPS proteins exhibited significant dihydroneopterin aldolase activity. Phylogenetic analysis indicated that PTPS-III proteins may have arisen independently in various PTPS lineages. Consistent with this possibility, merely introducing a glutamate residue into the active site of a PTPS protein conferred incipient activity in the growth complementation assay, and replacing glutamate with alanine in a PTPS-III protein abolished complementation.

A unique dual activity amino acid hydroxylase in Toxoplasma gondii

The genome of the protozoan parasite Toxoplasma gondii was found to contain two genes encoding tyrosine hydroxylase; that produces l-DOPA. The encoded enzymes metabolize phenylalanine as well as tyrosine with substrate preference for tyrosine. Thus the enzymes catabolize phenylalanine to tyrosine and tyrosine to l-DOPA. The catalytic domain descriptive of this class of enzymes is conserved with the parasite enzyme and exhibits similar kinetic properties to metazoan tyrosine hydroxylases, but contains a unique N-terminal extension with a signal sequence motif. One of the genes, TgAaaH1, is constitutively expressed while the other gene, TgAaaH2, is induced during formation of the bradyzoites of the cyst stages of the life cycle. This is the first description of an aromatic amino acid hydroxylase in an apicomplexan parasite. Extensive searching of apicomplexan genome sequences revealed an ortholog in Neospora caninum but not in Eimeria, Cryptosporidium, Theileria, or Plasmodium. Possible role(s) of these bi-functional enzymes during host infection are discussed.

Caveat emptor: limitations of the automated reconstruction of metabolic pathways in Plasmodium

The functional reconstruction of metabolic pathways from an annotated genome is a tedious and demanding enterprise. Automation of this endeavor using bioinformatics algorithms could cope with the ever-increasing number of sequenced genomes and accelerate the process. Here, the manual reconstruction of metabolic pathways in the functional genomic database of Plasmodium falciparum--Malaria Parasite Metabolic Pathways--is described and compared with pathways generated automatically as they appear in PlasmoCyc, metaSHARK and the Kyoto Encyclopedia for Genes and Genomes. A critical evaluation of this comparison discloses that the automatic reconstruction of pathways generates manifold paths that need an expert manual verification to accept some and reject most others based on manually curated gene annotation.

Adaptive copy number evolution in malaria parasites

Copy number polymorphism (CNP) is ubiquitous in eukaryotic genomes, but the degree to which this reflects the action of positive selection is poorly understood. The first gene in the Plasmodium folate biosynthesis pathway, GTP-cyclohydrolase I (gch1), shows extensive CNP. We provide compelling evidence that gch1 CNP is an adaptive consequence of selection by antifolate drugs, which target enzymes downstream in this pathway. (1) We compared gch1 CNP in parasites from Thailand (strong historical antifolate selection) with those from neighboring Laos (weak antifolate selection). Two percent of chromosomes had amplified copy number in Laos, while 72% carried multiple (2–11) copies in Thailand, and differentiation exceeded that observed at 73 synonymous SNPs. (2) We found five amplicon types containing one to greater than six genes and spanning 1 to >11 kb, consistent with parallel evolution and strong selection for this gene amplification. gch1 was the only gene occurring in all amplicons suggesting that this locus is the target of selection. (3) We observed reduced microsatellite variation and increased linkage disequilibrium (LD) in a 900-kb region flanking gch1 in parasites from Thailand, consistent with rapid recent spread of chromosomes carrying multiple copies of gch1. (4) We found that parasites bearing dhfr-164L, which causes high-level resistance to antifolate drugs, carry significantly (p = 0.00003) higher copy numbers of gch1 than parasites bearing 164I, indicating functional association between genes located on different chromosomes but linked in the same biochemical pathway. These results demonstrate that CNP at gch1 is adaptive and the associations with dhfr-164L strongly suggest a compensatory function. More generally, these data demonstrate how selection affects multiple enzymes in a single biochemical pathway, and suggest that investigation of structural variation may provide a fast-track to locating genes underlying adaptation.

Recent comparative genomic hybridization studies have revealed extensive copy number variation in eukaryotic genomes. The first gene in the Plasmodium folate biosynthesis pathway, GTP-cyclohydrolase I (gch1), shows extensive copy number polymorphism (CNP). We provide compelling evidence that gch1 CNP is adaptive and most likely results from selection by antifolate drugs, which target enzymes downstream in this pathway. Gch1 CNP shows extreme geographical differentiation; hitchhiking reduces diversity and increases LD in flanking sequence, indicating recent rapid spread within Thailand, while amplicon structure reveals multiple origins and parallel evolution. Furthermore, strong association between elevated copy number and a critical mutation dhfr-I164L that underlies high-level antifolate resistance indicates functional linkage and fitness epistasis between genes on different chromosomes. These data reveal hidden complexity in the evolutionary response to antifolate treatment and demonstrate that analysis of structural variation can provide a fast-track to locating genes that underlie adaptation.

Plasmodium falciparum: a paradigm for alternative folate biosynthesis in diverse microorganisms?

Folates have a key role in metabolism, and the folate-dependent generation of DNA precursors in the form of deoxythymidine 5'-phosphate is particularly important for the replication of malaria parasites. Although Plasmodium falciparum can synthesize folate derivatives de novo, a long-standing mystery has been the apparent absence of a key enzyme, dihydroneopterin aldolase, in the classical folate biosynthetic pathway of this organism. The discovery that a different enzyme, pyruvoyltetrahydropterin synthase, can produce the necessary substrate for the subsequent step in folate synthesis raises the question of whether this solution is unique to P. falciparum. Bioinformatic analyses suggest otherwise and indicate that an alternative route to folate could be widespread among diverse microorganisms and could be a target for novel drugs.