The roles of pteridine reductase 1 and dihydrofolate reductase-thymidylate synthase in pteridine metabolism in the protozoan parasite Leishmania major

Genomic insights into virulence mechanisms of Leishmania donovani: evidence from an atypical strain

Background

Leishmaniasis is a neglected tropical disease with diverse clinical phenotypes, determined by parasite, host and vector interactions. Despite the advances in molecular biology and the availability of more Leishmania genome references in recent years, the association between parasite species and distinct clinical phenotypes remains poorly understood. We present a genomic comparison of an atypical variant of Leishmania donovani from a South Asian focus, where it mostly causes cutaneous form of leishmaniasis.

Results

Clinical isolates from six cutaneous leishmaniasis patients (CL-SL); 2 of whom were poor responders to antimony (CL-PR), and two visceral leishmaniasis patients (VL-SL) were sequenced on an Illumina MiSeq platform. Chromosome aneuploidy was observed in both groups but was more frequent in CL-SL. 248 genes differed by 2 fold or more in copy number among the two groups. Genes involved in amino acid use (LdBPK_271940) and energy metabolism (LdBPK_271950), predominated the VL-SL group with the same distribution pattern reflected in gene tandem arrays. Genes encoding amastins were present in higher copy numbers in VL-SL and CL-PR as well as being among predicted pseudogenes in CL-SL. Both chromosome and SNP profiles showed CL-SL and VL-SL to form two distinct groups. While expected heterozygosity was much higher in VL-SL, SNP allele frequency patterns did not suggest potential recent recombination breakpoints. The SNP/indel profile obtained using the more recently generated PacBio sequence did not vary markedly from that based on the standard LdBPK282A1 reference. Several genes previously associated with resistance to antimonials were observed in higher copy numbers in the analysis of CL-PR. H-locus amplification was seen in one cutaneous isolate which however did not belong to the CL-PR group.

Conclusions

The data presented suggests that intra species variations at chromosome and gene level are more likely to influence differences in tropism as well as response to treatment, and contributes to greater understanding of parasite molecular mechanisms underpinning these differences. These findings should be substantiated with a larger sample number and expression/functional studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5271-z) contains supplementary material, which is available to authorized users.

Genomic Analysis of Colombian Leishmania panamensis strains with different level of virulence

The establishment of Leishmania infection in mammalian hosts and the subsequent manifestation of clinical symptoms require internalization into macrophages, immune evasion and parasite survival and replication. Although many of the genes involved in these processes have been described, the genetic and genomic variability associated to differences in virulence is largely unknown. Here we present the genomic variation of four Leishmania (Viannia) panamensis strains exhibiting different levels of virulence in BALB/c mice and its application to predict novel genes related to virulence. De novo DNA sequencing and assembly of the most virulent strain allowed comparative genomics analysis with sequenced L. (Viannia) panamensis and L. (Viannia) braziliensis strains, and showed important variations at intra and interspecific levels. Moreover, the mutation detection and a CNV search revealed both base and structural genomic variation within the species. Interestingly, we found differences in the copy number and protein diversity of some genes previously related to virulence. Several machine-learning approaches were applied to combine previous knowledge with features derived from genomic variation and predict a curated set of 66 novel genes related to virulence. These genes can be prioritized for validation experiments and could potentially become promising drug and immune targets for the development of novel prophylactic and therapeutic interventions.

Homology modelling, molecular docking, and molecular dynamics simulations reveal the inhibition of Leishmania donovani dihydrofolate reductase-thymidylate synthase enzyme by Withaferin-A

Objective

Present in silico study was carried out to explore the mode of inhibition of Leishmania donovani dihydrofolate reductase-thymidylate synthase (Ld DHFR-TS) enzyme by Withaferin-A, a withanolide isolated from Withania somnifera. Withaferin-A (WA) is known for its profound multifaceted properties, but its antileishmanial activity is not well understood. The parasite’s DHFR-TS enzyme is diverse from its mammalian host and could be a potential drug target in parasites.

Results

A 3D model of Ld DHFR-TS enzyme was built and verified using Ramachandran plot and SAVES tools. The protein was docked with WA-the ligand, methotrexate (MTX)-competitive inhibitor of DHFR, and dihydrofolic acid (DHFA)-substrate for DHFR-TS. Molecular docking studies reveal that WA competes for active sites of both Hu DHFR and TS enzymes whereas it binds to a site other than active site in Ld DHFR-TS. Moreover, Lys 173 residue of DHFR-TS forms a H-bond with WA and has higher binding affinity to Ld DHFR-TS than Hu DHFR and Hu TS. The MD simulations confirmed the H-bonding interactions were stable. The binding energies of WA with Ld DHFR-TS were calculated using MM-PBSA. Homology modelling, molecular docking and MD simulations of Ld DHFR-TS revealed that WA could be a potential anti-leishmanial drug.

Electronic supplementary material

The online version of this article (10.1186/s13104-018-3354-1) contains supplementary material, which is available to authorized users.

Protozoan Parasite Auxotrophies and Metabolic Dependencies

Diseases caused by protozoan parasites have a major impact on world health. These early branching eukaryotes cause significant morbidity and mortality in humans and livestock. During evolution, protozoan parasites have evolved toward complex life cycles in multiple host organisms with different nutritional resources. The conservation of functional metabolic pathways required for these successive environments is therefore a prerequisite for parasitic lifestyle. Nevertheless, parasitism drives genome evolution toward gene loss and metabolic dependencies (including strict auxotrophy), especially for obligatory intracellular parasites. In this chapter, we will compare and contrast how protozoan parasites have perfected this metabolic adaptation by focusing on specific auxotrophic pathways and scavenging strategies used by clinically relevant apicomplexan and trypanosomatid parasites to access host's nutritional resources. We will further see how these metabolic dependencies have in turn been exploited for therapeutic purposes against these human pathogens.

Discovery of novel, orally bioavailable, antileishmanial compounds using phenotypic screening

Leishmaniasis is a parasitic infection that afflicts approximately 12 million people worldwide. There are several limitations to the approved drug therapies for leishmaniasis, including moderate to severe toxicity, growing drug resistance, and the need for extended dosing. Moreover, miltefosine is currently the only orally available drug therapy for this infection. We addressed the pressing need for new therapies by pursuing a two-step phenotypic screen to discover novel, potent, and orally bioavailable antileishmanials. First, we conducted a high-throughput screen (HTS) of roughly 600,000 small molecules for growth inhibition against the promastigote form of the parasite life cycle using the nucleic acid binding dye SYBR Green I. This screen identified approximately 2,700 compounds that inhibited growth by over 65% at a single point concentration of 10 μM. We next used this 2700 compound focused library to identify compounds that were highly potent against the disease-causing intra-macrophage amastigote form and exhibited limited toxicity toward the host macrophages. This two-step screening strategy uncovered nine unique chemical scaffolds within our collection, including two previously described antileishmanials. We further profiled two of the novel compounds for in vitro absorption, distribution, metabolism, excretion, and in vivo pharmacokinetics. Both compounds proved orally bioavailable, affording plasma exposures above the half-maximal effective concentration (EC50) concentration for at least 12 hours. Both compounds were efficacious when administered orally in a murine model of cutaneous leishmaniasis. One of the two compounds exerted potent activity against trypanosomes, which are kinetoplastid parasites related to Leishmania species. Therefore, this compound could help control multiple parasitic diseases. The promising pharmacokinetic profile and significant in vivo efficacy observed from our HTS hits highlight the utility of our two-step phenotypic screening strategy and strongly suggest that medicinal chemistry optimization of these newly identified scaffolds will lead to promising candidates for an orally available anti-parasitic drug.

Predicting novel substrates for enzymes with minimal experimental effort with active learning

Enzymatic substrate promiscuity is more ubiquitous than previously thought, with significant consequences for understanding metabolism and its application to biocatalysis. This realization has given rise to the need for efficient characterization of enzyme promiscuity. Enzyme promiscuity is currently characterized with a limited number of human-selected compounds that may not be representative of the enzyme's versatility. While testing large numbers of compounds may be impractical, computational approaches can exploit existing data to determine the most informative substrates to test next, thereby more thoroughly exploring an enzyme's versatility. To demonstrate this, we used existing studies and tested compounds for four different enzymes, developed support vector machine (SVM) models using these datasets, and selected additional compounds for experiments using an active learning approach. SVMs trained on a chemically diverse set of compounds were discovered to achieve maximum accuracies of ~80% using ~33% fewer compounds than datasets based on all compounds tested in existing studies. Active learning-selected compounds for testing resolved apparent conflicts in the existing training data, while adding diversity to the dataset. The application of these algorithms to wide arrays of metabolic enzymes would result in a library of SVMs that can predict high-probability promiscuous enzymatic reactions and could prove a valuable resource for the design of novel metabolic pathways.

Application of a simple quantum chemical approach to ligand fragment scoring for Trypanosoma brucei pteridine reductase 1 inhibition

There is a need for improved and generally applicable scoring functions for fragment-based approaches to ligand design. Here, we evaluate the performance of a computationally efficient model for inhibitory activity estimation, which is composed only of multipole electrostatic energy and dispersion energy terms that approximate long-range ab initio quantum mechanical interaction energies. We find that computed energies correlate well with inhibitory activity for a compound series with varying substituents targeting two subpockets of the binding site of Trypanosoma brucei pteridine reductase 1. For one subpocket, we find that the model is more predictive for inhibitory activity than the ab initio interaction energy calculated at the MP2 level. Furthermore, the model is found to outperform a commonly used empirical scoring method. Finally, we show that the results for the two subpockets can be combined, which suggests that this simple nonempirical scoring function could be applied in fragment–based drug design.

Essential multimeric enzymes in kinetoplastid parasites: A host of potentially druggable protein-protein interactions

Parasitic diseases caused by kinetoplastid parasites of the genera Trypanosoma and Leishmania are an urgent public health crisis in the developing world. These closely related species possess a number of multimeric enzymes in highly conserved pathways involved in vital functions, such as redox homeostasis and nucleotide synthesis. Computational alanine scanning of these protein-protein interfaces has revealed a host of potentially ligandable sites on several established and emerging anti-parasitic drug targets. Analysis of interfaces with multiple clustered hotspots has suggested several potentially inhibitable protein-protein interactions that may have been overlooked by previous large-scale analyses focusing solely on secondary structure. These protein-protein interactions provide a promising lead for the development of new peptide and macrocycle inhibitors of these enzymes.

Structure-guided discovery of thiazolidine-2,4-dione derivatives as a novel class of Leishmania major pteridine reductase 1 inhibitors

Leishmania major, as other protozoan parasites, plague human kind since pre-historic times but it remains a worldwide ailment for which the therapeutic arsenal remains scarce. Although L. major is pteridine- and purine-auxotroph, well-established folate biosynthesis inhibitors, such as methotrexate, have poor effect over the parasite survival. The lack of efficiency is related to an alternative biochemical pathway in which pteridine reductase 1 (PTR1) plays a major role. For this reason, this enzyme has been considered a promising target for anti-leishmanial drug development and several inhibitors that share the substrate scaffold have been reported. In order to design a novel class of PTR1 inhibitors, we employed the thiazolidinone ring as a bioisosteric replacement for pteridine/purine ring. Among seven novel thiazolidine-2,4-dione derivatives reported herein, 2d was identified as the most promising lead by thermal shift assays (ΔTm = 11 °C, p = 0,01). Kinetic assays reveal that 2d has IC50 = 44.67 ± 1.74 μM and shows a noncompetitive behavior. This information guided docking studies and molecular dynamics simulations (50 000 ps) that supports 2d putative binding profile (H-bonding to Ser-111 and Leu-66) and shall be useful to design more potent inhibitors.

Trypanosoma brucei DHFR-TS Revisited: Characterisation of a Bifunctional and Highly Unstable Recombinant Dihydrofolate Reductase-Thymidylate Synthase

Bifunctional dihydrofolate reductase-thymidylate synthase (DHFR-TS) is a chemically and genetically validated target in African trypanosomes, causative agents of sleeping sickness in humans and nagana in cattle. Here we report the kinetic properties and sensitivity of recombinant enzyme to a range of lipophilic and classical antifolate drugs. The purified recombinant enzyme, expressed as a fusion protein with elongation factor Ts (Tsf) in ThyA- Escherichia coli, retains DHFR activity, but lacks any TS activity. TS activity was found to be extremely unstable (half-life of 28 s) following desalting of clarified bacterial lysates to remove small molecules. Stability could be improved 700-fold by inclusion of dUMP, but not by other pyrimidine or purine (deoxy)-nucleosides or nucleotides. Inclusion of dUMP during purification proved insufficient to prevent inactivation during the purification procedure. Methotrexate and trimetrexate were the most potent inhibitors of DHFR (Ki 0.1 and 0.6 nM, respectively) and FdUMP and nolatrexed of TS (Ki 14 and 39 nM, respectively). All inhibitors showed a marked drop-off in potency of 100- to 1,000-fold against trypanosomes grown in low folate medium lacking thymidine. The most potent inhibitors possessed a terminal glutamate moiety suggesting that transport or subsequent retention by polyglutamylation was important for biological activity. Supplementation of culture medium with folate markedly antagonised the potency of these folate-like inhibitors, as did thymidine in the case of the TS inhibitors raltitrexed and pemetrexed.

Leishmania infantum Asparagine Synthetase A Is Dispensable for Parasites Survival and Infectivity

A growing interest in asparagine (Asn) metabolism has currently been observed in cancer and infection fields. Asparagine synthetase (AS) is responsible for the conversion of aspartate into Asn in an ATP-dependent manner, using ammonia or glutamine as a nitrogen source. There are two structurally distinct AS: the strictly ammonia dependent, type A, and the type B, which preferably uses glutamine. Absent in humans and present in trypanosomatids, AS-A was worthy of exploring as a potential drug target candidate. Appealingly, it was reported that AS-A was essential in Leishmania donovani, making it a promising drug target. In the work herein we demonstrate that Leishmania infantum AS-A, similarly to Trypanosoma spp. and L. donovani, is able to use both ammonia and glutamine as nitrogen donors. Moreover, we have successfully generated LiASA null mutants by targeted gene replacement in L. infantum, and these parasites do not display any significant growth or infectivity defect. Indeed, a severe impairment of in vitro growth was only observed when null mutants were cultured in asparagine limiting conditions. Altogether our results demonstrate that despite being important under asparagine limitation, LiAS-A is not essential for parasite survival, growth or infectivity in normal in vitro and in vivo conditions. Therefore we exclude AS-A as a suitable drug target against L. infantum parasites.

Molecular characterization and functional analysis of pteridine reductase in wild-type and antimony-resistant Leishmania lines

Pteridine reductase (PTR1) is an NADPH-dependent reductase that participates in the salvage of pteridines, which are essential to maintain growth of Leishmania. In this study, we performed the molecular characterization of ptr1 gene in wild-type (WTS) and SbIII-resistant (SbR) lines from Leishmania guyanensis (Lg), Leishmania amazonensis (La), Leishmania braziliensis (Lb) and Leishmania infantum (Li), evaluating the chromosomal location, mRNA levels of the ptr1 gene and PTR1 protein expression. PFGE results showed that the ptr1 gene is located in a 797 kb chromosomal band in all Leishmania lines analyzed. Interestingly, an additional chromosomal band of 1070 kb was observed only in LbSbR line. Northern blot results showed that the levels of ptr1 mRNA are increased in the LgSbR, LaSbR and LbSbR lines. Western blot assays using the polyclonal anti-LmPTR1 antibody demonstrated that PTR1 protein is more expressed in the LgSbR, LaSbR and LbSbR lines compared to their respective WTS counterparts. Nevertheless, no difference in the level of mRNA and protein was observed between the LiWTS and LiSbR lines. Functional analysis of PTR1 enzyme was performed to determine whether the overexpression of ptr1 gene in the WTS L. braziliensis and L. infantum lines would change the SbIII-resistance phenotype of transfected parasites. Western blot results showed that the expression level of PTR1 protein was increased in the transfected parasites compared to the non-transfected ones. IC50 analysis revealed that the overexpression of ptr1 gene in the WTS L. braziliensis line increased 2-fold the SbIII-resistance phenotype compared to the non-transfected counterpart. Furthermore, the overexpression of ptr1 gene in the WTS L. infantum line did not change the SbIII-resistance phenotype. These results suggest that the PTR1 enzyme may be implicated in the SbIII-resistance phenotype in L. braziliensis line.

Folate Biosynthesis, Reduction, and Polyglutamylation and the Interconversion of Folate Derivatives

Many microorganisms and plants possess the ability to synthesize folic acid derivatives de novo, initially forming dihydrofolate. All the folic acid derivatives that serve as recipients and donors of one-carbon units are derivatives of tetrahydrofolate, which is formed from dihydrofolate by an NADPH-dependent reduction catalyzed by dihydrofolate reductase (FolA). This review discusses the biosynthesis of dihydrofolate monoglutamate, its reduction to tetrahydrofolate monoglutamate, and the addition of glutamyl residues to form folylpolyglutamates. Escherichia coli and Salmonella, like many microorganisms that can synthesize folate de novo, appear to lack the ability to transport folate into the cell and are thus highly susceptible to inhibitors of folate biosynthesis. The review includes a brief discussion of the inhibition of folate biosynthesis by sulfa drugs. The folate biosynthetic pathway can be divided into two sections. First, the aromatic precursor chorismate is converted to paminobenzoic acid (PABA) by the action of three proteins. Second, the pteridine portion of folate is made from GTP and coupled to PABA to generate dihydropteroate, and the bifunctional protein specified by folC, dihydrofolate synthetase, or folylpolyglutamate synthetase, adds the initial glutamate molecule to form dihydrofolate (H2PteGlu1, or dihydropteroylmonoglutamate). Bacteriophage T4 infection of E. coli has been shown to cause alterations in the metabolism of folate derivatives. Infection is associated with an increase in the chain lengths in folylpolyglutamates and particularly the accumulation of hexaglutamate derivatives.

Design, synthesis and biological evaluation of quinazoline derivatives as anti-trypanosomatid and anti-plasmodial agents

In this paper, the design, synthesis and biological evaluation of a set of quinazoline-2,4,6-triamine derivatives (1-9) as trypanocidal, antileishmanial and antiplasmodial agents are explained. The compounds were rationalized basing on docking studies of the dihydrofolate reductase (DHFR from Trypanosoma cruzi, Leishmania major and Plasmodium vivax) and pteridin reductase (PTR from T. cruzi and L. major) structures. All compounds were in vitro screened against both bloodstream trypomastigotes of T. cruzi (NINOA and INC-5 strains) and promatigotes of Leishmania mexicana (MHOM/BZ/61/M379 strain), and also for cytotoxicity using Vero cell line. Against T. cruzi, three compounds (5, 6 and 8) were the most effective showing a better activity profile than nifurtimox and benznidazole (reference drugs). Against L. mexicana, four compounds (5, 6, 8, and 9) exhibited the highest activity, even than glucantime (reference drug). In the cytotoxicity assay, protozoa were more susceptible than Vero cells. In vivo Plasmodium berghei assay (ANKA strain), the compounds 1, 5, 6 and 8 showed a more comparable activity than chloroquine and pyrimethamine (reference drugs) when they were administrated by the oral route. The antiprotozoal activity of these substances, endowed with redox properties, represented a good starting point for a medicinal chemistry program aiming for chemotherapy of Chagas' disease, leishmaniosis and malaria.

Structure, Activity and Stereoselectivity of NADPH-Dependent Oxidoreductases Catalysing the S-Selective Reduction of the Imine Substrate 2-Methylpyrroline

Oxidoreductases from Streptomyces sp. GF3546 [3546-IRED], Bacillus cereus BAG3X2 (BcIRED) and Nocardiopsis halophila (NhIRED) each reduce prochiral 2-methylpyrroline (2MPN) to (S)-2-methylpyrrolidine with >95 % ee and also a number of other imine substrates with good selectivity. Structures of BcIRED and NhIRED have helped to identify conserved active site residues within this subgroup of imine reductases that have S selectivity towards 2MPN, including a tyrosine residue that has a possible role in catalysis and superimposes with an aspartate in related enzymes that display R selectivity towards the same substrate. Mutation of this tyrosine residue-Tyr169-in 3546-IRED to Phe resulted in a mutant of negligible activity. The data together provide structural evidence for the location and significance of the Tyr residue in this group of imine reductases, and permit a comparison of the active sites of enzymes that reduce 2MPN with either R or S selectivity.

Genome evolution in trypanosomatid parasites

A decade of genome sequencing has transformed our understanding of how trypanosomatid parasites have evolved and provided fresh impetus to explaining the origins of parasitism in the Kinetoplastida. In this review, I will consider the many ways in which genome sequences have influenced our view of genomic reduction in trypanosomatids; how species-specific genes, and the genomic domains they occupy, have illuminated the innovations in trypanosomatid genomes; and how comparative genomics has exposed the molecular mechanisms responsible for innovation and adaptation to a parasitic lifestyle.

Investigation of osmolyte effects on FolM: comparison with other dihydrofolate reductases

A weak association between osmolytes and dihydrofolate (DHF) decreases the affinity of the substrate for the Escherichia coli chromosomal and R67 plasmid dihydrofolate reductase (DHFR) enzymes. To test whether the osmolyte-DHF association also interferes with binding of DHF to FolM, an E. coli enzyme that possesses weak DHFR activity, ligand binding was monitored in the presence of osmolytes. The affinity of FolM for DHF, measured by kcat/Km(DHF), was decreased by the addition of an osmolyte. Additionally, binding of the antifolate drug, methotrexate, to FolM was weakened by the addition of an osmolyte. The changes in ligand binding with water activity were unique for each osmolyte, indicating preferential interaction between the osmolyte and folate and its derivatives; however, additional evidence provided support for further interactions between FolM and osmolytes. Binding of the reduced nicotinamide adenine dinucleotide phosphate (NADPH) cofactor to FolM was monitored by isothermal titration calorimetry as a control for protein-osmolyte association. In the presence of betaine (proposed to be the osmolyte most excluded from protein surfaces), the NADPH Kd decreased, consistent with dehydration effects. However, other osmolytes did not tighten binding to the cofactor. Rather, dimethyl sulfoxide (DMSO) had no effect on the NADPH Kd, while ethylene glycol and polyethylene glycol 400 weakened cofactor binding. Differential scanning calorimetry of FolM in the presence of osmolytes showed that both DMSO and ethylene glycol decreased the stability of FolM, while betaine increased the stability of the protein. These results suggest that some osmolytes can destabilize FolM by preferentially interacting with the protein. Further, these weak attractions can impede ligand binding. These various contributions have to be considered when interpreting osmotic pressure results.

Selection of molecular targets for drug development against trypanosomatids

Trypanosomatid parasites are a group of flagellated protozoa that includes the genera Leishmania and Trypanosoma, which are the causative agents of diseases (leishmaniases, sleeping sickness and Chagas disease) that cause considerable morbidity and mortality, affecting more than 27 million people worldwide. Today no effective vaccines for the prevention of these diseases exist, whereas current chemotherapy is ineffective, mainly due to toxic side effects of current drugs and to the emergence of drug resistance and lack of cost effectiveness. For these reasons, rational drug design and the search of good candidate drug targets is of prime importance. The search for drug targets requires a multidisciplinary approach. To this end, the completion of the genome project of many trypanosomatid species gives a vast amount of new information that can be exploited for the identification of good drug candidates with a prediction of "druggability" and divergence from mammalian host proteins. In addition, an important aspect in the search for good drug targets is the "target identification" and evaluation in a biological pathway, as well as the essentiality of the gene in the mammalian stage of the parasite, which is provided by basic research and genetic and proteomic approaches. In this chapter we will discuss how these bioinformatic tools and experimental evaluations can be integrated for the selection of candidate drug targets, and give examples of metabolic and signaling pathways in the parasitic protozoa that can be exploited for rational drug design.

Antimicrobial Dihydrofolate Reductase Inhibitors - Achievements and Future Options: Review

Despite all progress made in the fight against infections caused by bacteria, fungi, protozoa or viruses, there is a need for more and new active agents. Intensive efforts are currently directed against many new and attractive targets, and are hoped to result in new useful agents. The opportunities offered by some known and validated targets are, however, by far not exhausted. Dihydrofolate reductase (DHFR, EC 1.5.1.3) attracted much attention over several decades, which yielded several useful agents. There are excellent chances for new drugs in this field, and they are thought to increase by limiting the spectrum of activity. Whereas trimethoprim seems to present the optimum which can be achieved for a broad spectrum antibacterial agent, specific agents could probably be designed for well defined groups or specific organisms, such as staphylococci among the bacteria, or for a number of parasites, such as Plasmodium falciparum, the fungus Pneumocystis carinii, and several protozoa, such as Trypanosoma, Toxoplasma, and others. This would even extend to herbicides or specific plant pathogens. Achievements and current efforts directed against new DHFR-inhibitors are reviewed, considering only the most recent literature.

Biochemical effects of riluzole on Leishmania parasites

We have previously shown that riluzole (6-(trifluoromethoxy)benzothiazol-2-amine), an agent used to treat CNS disorders, possesses inhibitory activity against pteridine reductase (PTR1) in pathogenic protists at low micromolar concentrations. Therefore, the potential use of this drug in anti-parasitic chemotherapy deserves evaluation. In this study, we report the effect of this compound on cell cultures of Leishmania mexicana and L. major. The anti-parasitic activity of riluzole was confirmed, with the largest effect observed when the drug was administered to cells during their exponential growth phase. Moreover, a remarkable decrease in PTR1 activity was observed in the lysates of cells pretreated with the compound, which is due to impairment of the enzyme's preferential reaction with biopterin as a cofactor. In addition, the treatment increased the parasites' susceptibility to oxidative stress, affecting the ability of Leishmania to survive under severe oxidative conditions. These results suggest that the inhibitory effect of riluzole on PTR1 is not the only mechanism through which it induces the death of Leishmania parasites.

Leishmania donovani pteridine reductase 1: comparative protein modeling and protein-ligand interaction studies of the leishmanicidal constituents isolated from the fruits of Piper longum

Visceral leishmaniasis or kala-azar is caused by the dimorphic parasite Leishmania donovani in the Indian subcontinent. Treatment options for kala-azar are currently inadequate due to various limitations. Currently, drug discovery for leishmaniases is oriented towards rational drug design; the aim is to identify specific inhibitors that target particular metabolic activities as a possible means of controlling the parasites without affecting the host. Leishmania salvages pteridin from its host and reduces it using pteridine reductase 1 (PTR1, EC 1.5.1.33), which makes this reductase an excellent drug target. Recently, we identified six alkamides and one benzenoid compound from the n-hexane fraction of the fruit of Piper longum that possess potent leishmanicidal activity against promastigotes as well as axenic amastigotes. Based on a homology model derived for recombinant pteridine reductase isolated from a clinical isolate of L. donovani, we carried out molecular modeling and docking studies with these compounds to evaluate their binding affinity. A fairly good agreement between experimental data and the results of molecular modeling investigation of the bioactive and inactive compounds was observed. The amide group in the conjugated alkamides and the 3,4-methylenedioxystyrene moiety in the benzenoid compound acts as heads and the long aliphatic chain acts as a tail, thus playing important roles in the binding of the inhibitor to the appropriate position at the active site. The remarkably high activity of a component containing piperine and piperine isomers (3.36:1) as observed by our group prompted us to study the activities of all four isomers of piperine-piperine (2E,4E), isopiperine (2Z,4E), isochavicine (2E,4Z), and chavicine (2Z,4Z)-against LdPTR1. The maximum inhibitory effect was demonstrated by isochavicine. The identification of these predicted inhibitors of LdPTR1 allowed us to build up a stereoview of the structure of the binding site in relation to activity, affording significant information that should prove useful during the structure-based design of leishmanicidal drugs.

Folate metabolic pathways in Leishmania

Trypanosomatid parasitic protozoans of the genus Leishmania are autotrophic for both folate and unconjugated pteridines. Leishmania salvage these metabolites from their mammalian hosts and insect vectors through multiple transporters. Within the parasite, folates are reduced by a bifunctional DHFR (dihydrofolate reductase)-TS (thymidylate synthase) and by a novel PTR1 (pteridine reductase 1), which reduces both folates and unconjugated pteridines. PTR1 can act as a metabolic bypass of DHFR inhibition, reducing the effectiveness of existing antifolate drugs. Leishmania possess a reduced set of folate-dependent metabolic reactions and can salvage many of the key products of folate metabolism from their hosts. For example, they lack purine synthesis, which normally requires 10-formyltetrahydrofolate, and instead rely on a network of purine salvage enzymes. Leishmania elaborate at least three pathways for the synthesis of the key metabolite 5,10-methylene-tetrahydrofolate, required for the synthesis of thymidylate, and for 10-formyltetrahydrofolate, whose presumptive function is for methionyl-tRNAMet formylation required for mitochondrial protein synthesis. Genetic studies have shown that the synthesis of methionine using 5-methyltetrahydrofolate is dispensable, as is the activity of the glycine cleavage complex, probably due to redundancy with serine hydroxymethyltransferase. Although not always essential, the loss of several folate metabolic enzymes results in attenuation or loss of virulence in animal models, and a null DHFR-TS mutant has been used to induce protective immunity. The folate metabolic pathway provides numerous opportunities for targeted chemotherapy, with strong potential for 'repurposing' of compounds developed originally for treatment of human cancers or other infectious agents.

Discovery of safe and orally effective 4-aminoquinaldine analogues as apoptotic inducers with activity against experimental visceral leishmaniasis

Novel antileishmanials are urgently required to overcome emergence of drug resistance, cytotoxic effects, and difficulties in oral delivery. Toward this, we investigated a series of novel 4-aminoquinaldine derivatives, a new class of molecules, as potential antileishmanials. 4-Aminoquinaldine derivatives presented inhibitory effects on L. donovani promastigotes and amastigotes (50% inhibitory concentration range, 0.94 to 127 μM). Of these, PP-9 and PP-10 were the most effective in vitro and demonstrated strong efficacies in vivo through the intraperitoneal route. They were also found to be effective against both sodium antimony gluconate-sensitive and -resistant Leishmania donovani strains in BALB/c mice when treated orally, resulting in more than 95% protection. Investigation of their mode of action revealed that killing by PP-10 involved moderate inhibition of dihydrofolate reductase and elicitation of the apoptotic cascade. Our studies implicate that PP-10 augments reactive oxygen species generation, evidenced from decreased glutathione levels and increased lipid peroxidation. Subsequent disruption of Leishmania promastigote mitochondrial membrane potential and activation of cytosolic proteases initiated the apoptotic pathway, resulting in DNA fragmentation and parasite death. Our results demonstrate that PP-9 and PP-10 are promising lead compounds with the potential for treating visceral leishmaniasis (VL) through the oral route.

Trypanosomal dihydrofolate reductase reveals natural antifolate resistance

Dihydrofolate reductase (DHFR) is a potential drug target for Trypanosoma brucei, a human parasite, which is the causative agent for African sleeping sickness. No drug is available against this target, since none of the classical antifolates such as pyrimethamine (PYR), cycloguanil, or trimethoprim are effective as selective inhibitors of T. brucei DHFR (TbDHFR). In order to design effective drugs that target TbDHFR, co-crystal structures with bound antifolates were studied. On comparison with malarial Plasmodium falciparum DHFR (PfDHFR), the co-crystal structures of wild-type TbDHFR reveal greater structural similarities to a mutant PfDHFR causing antifolate resistance than the wild-type enzyme. TbDHFR imposes steric hindrance for rigid inhibitors like PYR around Thr86, which is equivalent to Ser108Asn of the malarial enzymes. In addition, a missing residue on TbDHFR active-site loop together with the presence of Ile51 widens its active site even further than the structural effect of Asn51Ile, which is observed in PfDHFR structures. The structural similarities are paralleled by the similarly poor affinities of the trypanosomal enzyme for rigid inhibitors. Mutations of TbDHFR at Thr86 resulted in 10-fold enhancement or 7-fold reduction in the rigid inhibitors affinities for Thr86Ser or Thr86Asn, respectively. The co-crystal structure of TbDHFR with a flexible antifolate WR99210 suggests that its greater affinity result from its ability to avoid such Thr86 clash and occupy the widened binding space similarly to what is observed in the PfDHFR structures. Natural resistance to antifolates of TbDHFR can therefore be explained, and potential antifolate chemotherapy of trypanosomiasis should be possible taking this into account.

3D-QSAR based pharmacophore modeling and virtual screening for identification of novel pteridine reductase inhibitors

Pteridine reductase is a promising target for development of novel therapeutic agents against Trypanosomatid parasites. A 3D-QSAR pharmacophore hypothesis has been generated for a series of L. major pteridine reductase inhibitors using Catalyst/HypoGen algorithm for identification of the chemical features that are responsible for the inhibitory activity. Four pharmacophore features, namely: two H-bond donors (D), one Hydrophobic aromatic (H) and one Ring aromatic (R) have been identified as key features involved in inhibitor-PTR1 interaction. These features are able to predict the activity of external test set of pteridine reductase inhibitors with a correlation coefficient (r) of 0.80. Based on the analysis of the best hypotheses, some potent Pteridine reductase inhibitors were screened out and predicted with anti-PTR1 activity. It turned out that the newly identified inhibitory molecules are at least 300 fold more potent than the current crop of existing inhibitors. Overall the current SAR study is an effort for elucidating quantitative structure-activity relationship for the PTR1 inhibitors. The results from the combined 3D-QSAR modeling and molecular docking approach have led to the prediction of new potent inhibitory scaffolds.

Leishmania-macrophage interactions: insights into the redox biology

Leishmaniasis is a neglected tropical disease that affects about 350 million individuals worldwide. The protozoan parasite has a relatively simple life cycle with two principal stages: the flagellated mobile promastigote living in the gut of the sandfly vector and the intracellular amastigote within phagolysosomal vesicles of the vertebrate host macrophage. This review presents a state-of-the-art overview of the redox biology at the parasite-macrophage interface. Although Leishmania species are susceptible in vitro to exogenous superoxide radical, hydrogen peroxide, nitric oxide, and peroxynitrite, they manage to survive the endogenous oxidative burst during phagocytosis and the subsequent elevated nitric oxide production in the macrophage. The parasite adopts various defense mechanisms to cope with oxidative stress: the lipophosphoglycan membrane decreases superoxide radical production by inhibiting NADPH oxidase assembly and the parasite also protects itself by expressing antioxidant enzymes and proteins. Some of these enzymes could be considered potential drug targets because they are not expressed in mammals. In respect to antileishmanial therapy, the effects of current drugs on parasite-macrophage redox biology and its involvement in the development of drug resistance and treatment failure are presented.

Contribution of proteomics of Leishmania spp. to the understanding of differentiation, drug resistance mechanisms, vaccine and drug development

Leishmania spp., protozoan parasites with a digenetic life cycle, cause a spectrum of diseases in humans. Recently several Leishmania spp. have been sequenced which significantly boosted the number and quality of proteomic studies conducted. Here a historic review will summarize work of the pre-genomic era and then focus on studies after genome information became available. Firstly works comparing the different life cycle stages, in order to identify stage specific proteins, will be discussed. Identifying post-translational modifications by proteomics especially phosphorylation events will be discussed. Further the contribution of proteomics to the understanding of the molecular mechanism of drug resistance and the investigation of immunogenic proteins for the identification of vaccine candidates will be summarized. Approaches of how potentially secreted proteins were identified are discussed. So far 30-35% of the total predicted proteome of Leishmania spp. have been identified. This comprises mainly the abundant proteins, therefore the last section will look into technological approaches on how this coverage may be increased and what the gel-free and gel-based proteomics have to offer will be compared.

Design, synthesis and biological evaluation of novel inhibitors of Trypanosoma brucei pteridine reductase 1

Genetic studies indicate that the enzyme pteridine reductase 1 (PTR1) is essential for the survival of the protozoan parasite Trypanosoma brucei. Herein, we describe the development and optimisation of a novel series of PTR1 inhibitors, based on benzo[d]imidazol-2-amine derivatives. Data are reported on 33 compounds. This series was initially discovered by a virtual screening campaign (J. Med. Chem., 2009, 52, 4454). The inhibitors adopted an alternative binding mode to those of the natural ligands, biopterin and dihydrobiopterin, and classical inhibitors, such as methotrexate. Using both rational medicinal chemistry and structure-based approaches, we were able to derive compounds with potent activity against T. brucei PTR1 (K(i)(app)=7 nM), which had high selectivity over both human and T. brucei dihydrofolate reductase. Unfortunately, these compounds displayed weak activity against the parasites. Kinetic studies and analysis indicate that the main reason for the lack of cell potency is due to the compounds having insufficient potency against the enzyme, which can be seen from the low K(m) to K(i) ratio (K(m)=25 nM and K(i)=2.3 nM, respectively).

Dissecting the metabolic roles of pteridine reductase 1 in Trypanosoma brucei and Leishmania major

Leishmania parasites are pteridine auxotrophs that use an NADPH-dependent pteridine reductase 1 (PTR1) and NADH-dependent quinonoid dihydropteridine reductase (QDPR) to salvage and maintain intracellular pools of tetrahydrobiopterin (H₄B). However, the African trypanosome lacks a credible candidate QDPR in its genome despite maintaining apparent QDPR activity. Here we provide evidence that the NADH-dependent activity previously reported by others is an assay artifact. Using an HPLC-based enzyme assay, we demonstrate that there is an NADPH-dependent QDPR activity associated with both TbPTR1 and LmPTR1. The kinetic properties of recombinant PTR1s are reported at physiological pH and ionic strength and compared with LmQDPR. Specificity constants (k_cat/K_m) for LmPTR1 are similar with dihydrobiopterin (H₂B) and quinonoid dihydrobiopterin (qH₂B) as substrates and about 20-fold lower than LmQDPR with qH₂B. In contrast, TbPTR1 shows a 10-fold higher k_cat/K_m for H₂B over qH₂B. Analysis of Trypanosoma brucei isolated from infected rats revealed that H₄B (430 nm, 98% of total biopterin) was the predominant intracellular pterin, consistent with a dual role in the salvage and regeneration of H₄B. Gene knock-out experiments confirmed this: PTR1-nulls could only be obtained from lines overexpressing LmQDPR with H₄B as a medium supplement. These cells grew normally with H₄B, which spontaneously oxidizes to qH₂B, but were unable to survive in the absence of pterin or with either biopterin or H₂B in the medium. These findings establish that PTR1 has an essential and dual role in pterin metabolism in African trypanosomes and underline its potential as a drug target.

Phenylalanine hydroxylase (PAH) from the lower eukaryote Leishmania major

Aromatic amino acid hydroxylases (AAAH) typically use tetrahydrobiopterin (H(4)B) as the cofactor. The protozoan parasite Leishmania major requires biopterin for growth and expresses strong salvage and regeneration systems to maintain H(4)B levels. Here we explored the consequences of genetic manipulation of the sole L. major phenylalanine hydroxylase (PAH) to explore whether it could account for the Leishmania H(4)B requirement. L. major PAH resembles AAAHs of other organisms, bearing eukaryotic-type domain organization, and conservation of key catalytic residues including those implicated in pteridine binding. A pah(-) null mutant and an episomal complemented overexpressing derivative (pah-/+PAH) were readily obtained, and metabolic labeling studies established that PAH was required to hydroxylate Phe to Tyr. Neither WT nor overexpressing lines were able to hydroxylate radiolabeled tyrosine or tryptophan, nor to synthesize catecholamines. WT but not pah(-) parasites showed reactivity with an antibody to melanin when grown with l-3,4-dihydroxyphenylalanine (L-DOPA), although the reactive product is unlikely to be melanin sensu strictu. WT was auxotrophic for Phe, Trp and Tyr, suggesting that PAH activity was insufficient to meet normal Tyr requirements. However, pah(-) showed an increased sensitivity to Tyr deprivation, while the pah(-)/+PAH overexpressor showed increased survival and could be adapted to grow well without added Tyr. pah(-) showed no alterations in H(4)B-dependent differentiation, as established by in vitro metacyclogenesis, or survival in mouse or macrophage infections. Thus Leishmania PAH may mitigate but not alleviate Tyr auxotrophy, but plays no essential role in the steps of the parasite infectious cycle. These findings suggest PAH is unlikely to explain the Leishmania requirement for biopterin.

Structure of recombinant Leishmania donovani pteridine reductase reveals a disordered active site

Pteridine reductase (PTR1) is a potential target for drug development against parasitic Trypanosoma and Leishmania species, protozoa that are responsible for a range of serious diseases found in tropical and subtropical parts of the world. As part of a structure-based approach to inhibitor development, specifically targeting Leishmania species, well ordered crystals of L. donovani PTR1 were sought to support the characterization of complexes formed with inhibitors. An efficient system for recombinant protein production was prepared and the enzyme was purified and crystallized in an orthorhombic form with ammonium sulfate as the precipitant. Diffraction data were measured to 2.5 Å resolution and the structure was solved by molecular replacement. However, a sulfate occupies a phosphate-binding site used by NADPH and occludes cofactor binding. The nicotinamide moiety is a critical component of the active site and without it this part of the structure is disordered. The crystal form obtained under these conditions is therefore unsuitable for the characterization of inhibitor complexes.

High-resolution structures of Trypanosoma brucei pteridine reductase ligand complexes inform on the placement of new molecular entities in the active site of a potential drug target

Pteridine reductase (PTR1) is a potential target for drug development against parasitic Trypanosoma and Leishmania species. These protozoa cause serious diseases for which current therapies are inadequate. High-resolution structures have been determined, using data between 1.6 and 1.1 Å resolution, of T. brucei PTR1 in complex with pemetrexed, trimetrexate, cyromazine and a 2,4-diaminopyrimidine derivative. The structures provide insight into the interactions formed by new molecular entities in the enzyme active site with ligands that represent lead compounds for structure-based inhibitor development and to support early-stage drug discovery.

In silico screening, structure-activity relationship, and biologic evaluation of selective pteridine reductase inhibitors targeting visceral leishmaniasis

In this study we utilized the concept of rational drug design to identify novel compounds with optimal selectivity, efficacy and safety, which would bind to the target enzyme pteridine reductase 1 (PTR1) in Leishmania parasites. Twelve compounds afforded from Baylis-Hillman chemistry were docked by using the QUANTUM program into the active site of Leishmania donovani PTR1 homology model. The biological activity for these compounds was estimated in green fluorescent protein-transfected L. donovani promastigotes, and the most potential analogue was further investigated in intracellular amastigotes. Structure-activity relationship based on homology model drawn on our recombinant enzyme was substantiated by recombinant enzyme inhibition assay and growth of the cell culture. Flow cytometry results indicated that 7-(4-chlorobenzyl)-3-methyl-4-(4-trifluoromethyl-phenyl)-3,4,6,7,8,9-hexahydro-pyrimido[1,2-a]pyrimidin-2-one (compound 7) was 10 times more active on L. donovani amastigotes (50% inhibitory concentration [IC(50)] = 3 μM) than on promastigotes (IC(50) = 29 μM). Compound 7 exhibited a K(i) value of 0.72 μM in a recombinant enzyme inhibition assay. We discovered that novel pyrimido[1,2-a]pyrimidin-2-one systems generated from the allyl amines afforded from the Baylis-Hillman acetates could have potential as a valuable pharmacological tool against the neglected disease visceral leishmaniasis.

Frequency of drug resistance gene amplification in clinical leishmania strains

Experimental studies about Leishmania resistance to metal and antifolates have pointed out that gene amplification is one of the main mechanisms of drug detoxification. Amplified genes code for adenosine triphosphate-dependent transporters (multidrug resistance and P-glycoproteins P), enzymes involved in trypanothione pathway, particularly gamma glutamyl cysteine synthase, and others involved in folates metabolism, such as dihydrofolate reductase and pterine reductase. The aim of this study was to detect and quantify the amplification of these genes in clinical strains of visceral leishmaniasis agents: Leishmania infantum, L. donovani, and L. archibaldi. Relative quantification experiments by means of real-time polymerase chain reaction showed that multidrug resistance gene amplification is the more frequent event. For P-glycoproteins P and dihydrofolate reductase genes, level of amplification was comparable to the level observed after in vitro selection of resistant clones. Gene amplification is therefore a common phenomenon in wild strains concurring to Leishmania genomic plasticity. This finding, which corroborates results of experimental studies, supports a better understanding of metal resistance selection and spreading in endemic areas.

Structure-based design of pteridine reductase inhibitors targeting African sleeping sickness and the leishmaniases

Pteridine reductase (PTR1) is a target for drug development against Trypanosoma and Leishmania species, parasites that cause serious tropical diseases and for which therapies are inadequate. We adopted a structure-based approach to the design of novel PTR1 inhibitors based on three molecular scaffolds. A series of compounds, most newly synthesized, were identified as inhibitors with PTR1-species specific properties explained by structural differences between the T. brucei and L. major enzymes. The most potent inhibitors target T. brucei PTR1, and two compounds displayed antiparasite activity against the bloodstream form of the parasite. PTR1 contributes to antifolate drug resistance by providing a molecular bypass of dihydrofolate reductase (DHFR) inhibition. Therefore, combining PTR1 and DHFR inhibitors might improve therapeutic efficacy. We tested two new compounds with known DHFR inhibitors. A synergistic effect was observed for one particular combination highlighting the potential of such an approach for treatment of African sleeping sickness.

Development and validation of a cytochrome c-coupled assay for pteridine reductase 1 and dihydrofolate reductase

Activity of the pterin- and folate-salvaging enzymes pteridine reductase 1 (PTR1) and dihydrofolate reductase-thymidylate synthetase (DHFR-TS) is commonly measured as a decrease in absorbance at 340 nm, corresponding to oxidation of nicotinamide adenine dinucleotide phosphate (NADPH). Although this assay has been adequate to study the biology of these enzymes, it is not amenable to support any degree of routine inhibitor assessment because its restricted linearity is incompatible with enhanced throughput microtiter plate screening. In this article, we report the development and validation of a nonenzymatically coupled screening assay in which the product of the enzymatic reaction reduces cytochrome c, causing an increase in absorbance at 550 nm. We demonstrate this assay to be robust and accurate, and we describe its utility in supporting a structure-based design, small-molecule inhibitor campaign against Trypanosoma brucei PTR1 and DHFR-TS.

Removal of substrate inhibition and increase in maximal velocity in the short chain dehydrogenase/reductase salutaridine reductase involved in morphine biosynthesis

Salutaridine reductase (SalR, EC 1.1.1.248) catalyzes the stereospecific reduction of salutaridine to 7(S)-salutaridinol in the biosynthesis of morphine. It belongs to a new, plant-specific class of short-chain dehydrogenases, which are characterized by their monomeric nature and increased length compared with related enzymes. Homology modeling and substrate docking suggested that additional amino acids form a novel alpha-helical element, which is involved in substrate binding. Site-directed mutagenesis and subsequent studies on enzyme kinetics revealed the importance of three residues in this element for substrate binding. Further replacement of eight additional residues led to the characterization of the entire substrate binding pocket. In addition, a specific role in salutaridine binding by either hydrogen bond formation or hydrophobic interactions was assigned to each amino acid. Substrate docking also revealed an alternative mode for salutaridine binding, which could explain the strong substrate inhibition of SalR. An alternate arrangement of salutaridine in the enzyme was corroborated by the effect of various amino acid substitutions on substrate inhibition. In most cases, the complete removal of substrate inhibition was accompanied by a substantial loss in enzyme activity. However, some mutations greatly reduced substrate inhibition while maintaining or even increasing the maximal velocity. Based on these results, a double mutant of SalR was created that exhibited the complete absence of substrate inhibition and higher activity compared with wild-type SalR.

One scaffold, three binding modes: novel and selective pteridine reductase 1 inhibitors derived from fragment hits discovered by virtual screening

The enzyme pteridine reductase 1 (PTR1) is a potential target for new compounds to treat human African trypanosomiasis. A virtual screening campaign for fragments inhibiting PTR1 was carried out. Two novel chemical series were identified containing aminobenzothiazole and aminobenzimidazole scaffolds, respectively. One of the hits (2-amino-6-chloro-benzimidazole) was subjected to crystal structure analysis and a high resolution crystal structure in complex with PTR1 was obtained, confirming the predicted binding mode. However, the crystal structures of two analogues (2-amino-benzimidazole and 1-(3,4-dichloro-benzyl)-2-amino-benzimidazole) in complex with PTR1 revealed two alternative binding modes. In these complexes, previously unobserved protein movements and water-mediated protein-ligand contacts occurred, which prohibited a correct prediction of the binding modes. On the basis of the alternative binding mode of 1-(3,4-dichloro-benzyl)-2-amino-benzimidazole, derivatives were designed and selective PTR1 inhibitors with low nanomolar potency and favorable physicochemical properties were obtained.

LeishCyc: a biochemical pathways database for Leishmania major

Background

Leishmania spp. are sandfly transmitted protozoan parasites that cause a spectrum of diseases in more than 12 million people worldwide. Much research is now focusing on how these parasites adapt to the distinct nutrient environments they encounter in the digestive tract of the sandfly vector and the phagolysosome compartment of mammalian macrophages. While data mining and annotation of the genomes of three Leishmania species has provided an initial inventory of predicted metabolic components and associated pathways, resources for integrating this information into metabolic networks and incorporating data from transcript, protein, and metabolite profiling studies is currently lacking. The development of a reliable, expertly curated, and widely available model of Leishmania metabolic networks is required to facilitate systems analysis, as well as discovery and prioritization of new drug targets for this important human pathogen.

Description

The LeishCyc database was initially built from the genome sequence of Leishmania major (v5.2), based on the annotation published by the Wellcome Trust Sanger Institute. LeishCyc was manually curated to remove errors, correct automated predictions, and add information from the literature. The ongoing curation is based on public sources, literature searches, and our own experimental and bioinformatics studies. In a number of instances we have improved on the original genome annotation, and, in some ambiguous cases, collected relevant information from the literature in order to help clarify gene or protein annotation in the future. All genes in LeishCyc are linked to the corresponding entry in GeneDB (Wellcome Trust Sanger Institute).

Conclusion

The LeishCyc database describes Leishmania major genes, gene products, metabolites, their relationships and biochemical organization into metabolic pathways. LeishCyc provides a systematic approach to organizing the evolving information about Leishmania biochemical networks and is a tool for analysis, interpretation, and visualization of Leishmania Omics data (transcriptomics, proteomics, metabolomics) in the context of metabolic pathways. LeishCyc is the first such database for the Trypanosomatidae family, which includes a number of other important human parasites. Flexible query/visualization capabilities are provided by the Pathway Tools software and its Web interface. The LeishCyc database is made freely available over the Internet .

PTR1-dependent synthesis of tetrahydrobiopterin contributes to oxidant susceptibility in the trypanosomatid protozoan parasite Leishmania major

Leishmania must survive oxidative stress, but lack many classical antioxidant enzymes and rely heavily on trypanothione-dependent pathways. We used forward genetic screens to recover loci mediating oxidant resistance via overexpression in Leishmania major, which identified pteridine reductase 1 (PTR1). Comparisons of isogenic lines showed ptr1 (-) null mutants were 18-fold more sensitive to H(2)O(2) than PTR1-overproducing lines, and significant three- to fivefold differences were seen with a broad panel of oxidant-inducing agents. The toxicities of simple nitric oxide generators and other drug classes (except antifolates) were unaffected by PTR1 levels. H(2)O(2) susceptibility could be modulated by exogenous biopterin but not folate, in a PTR1- but not dihydrofolate reductase-dependent manner, implicating H(4)B metabolism specifically. Neither H(2)O(2) consumption nor the level of intracellular oxidative stress was affected by PTR1 levels. Coupled with the fact that reduced pteridines are at least 100-fold less abundant than cellular thiols, these data argue strongly that reduced pteridines act through a mechanism other than scavenging. The ability of unconjugated pteridines to counter oxidative stress has implications to infectivity and response to chemotherapy. Since the intracellular pteridine levels of Leishmania can be readily manipulated, these organisms offer a powerful setting for the dissection of pteridine-dependent oxidant susceptibility in higher eukaryotes.

The role of reduced pterins in resistance to reactive oxygen and nitrogen intermediates in the protozoan parasite Leishmania

During its life cycle, the protozoan parasite Leishmania experiences oxidative stress when interacting with macrophages. Reduced pterins are known scavengers of reactive oxygen and nitrogen intermediates. Leishmania has a pteridine reductase, PTR1, whose main function is to provide reduced pterins. We investigated the role of PTR1 in resistance to oxidative and nitrosative stress in Leishmania tarentolae, Leishmania infantum, and Leishmania major PTR1(-/-) mutants. The PTR1(-/-) cells of the three species were more sensitive to H2O2- and NO-induced stress. Using a fluorescent probe allowing ROI quantification, we demonstrated an increase in intracellular oxidant molecules in the PTR1(-/-) mutants. The disruption of PTR1 increased metacyclogenesis in L. infantum and L. major. We purified metacyclic parasites from PTR1(-/-) mutants and control cells and tested their intracellular survival in the J774 mouse cell line and in human monocyte-derived macrophages. Our results showed that PTR1(-/-) null mutants survived less in both macrophage models compared to control cells and this decrease was more pronounced in macrophages activated for oxidant production. This study demonstrates that one physiological role of reduced pterins in Leishmania is to deal with oxidative and nitrosative species, and a decreased ability to provide reduced pterins leads to decreased intracellular survival.