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

Multitarget, Selective Compound Design Yields Potent Inhibitors of a Kinetoplastid Pteridine Reductase 1

The optimization of compounds with multiple targets is a difficult multidimensional problem in the drug discovery cycle. Here, we present a systematic, multidisciplinary approach to the development of selective antiparasitic compounds. Computational fragment-based design of novel pteridine derivatives along with iterations of crystallographic structure determination allowed for the derivation of a structure–activity relationship for multitarget inhibition. The approach yielded compounds showing apparent picomolar inhibition of T. brucei pteridine reductase 1 (PTR1), nanomolar inhibition of L. major PTR1, and selective submicromolar inhibition of parasite dihydrofolate reductase (DHFR) versus human DHFR. Moreover, by combining design for polypharmacology with a property-based on-parasite optimization, we found three compounds that exhibited micromolar EC₅₀ values against T. brucei brucei while retaining their target inhibition. Our results provide a basis for the further development of pteridine-based compounds, and we expect our multitarget approach to be generally applicable to the design and optimization of anti-infective agents.

Repurposing the Trypanosomatidic GSK Kinetobox for the Inhibition of Parasitic Pteridine and Dihydrofolate Reductases

Three open-source anti-kinetoplastid chemical boxes derived from a whole-cell phenotypic screening by GlaxoSmithKline (Tres Cantos Anti-Kinetoplastid Screening, TCAKS) were exploited for the discovery of a novel core structure inspiring new treatments of parasitic diseases targeting the trypansosmatidic pteridine reductase 1 (PTR1) and dihydrofolate reductase (DHFR) enzymes. In total, 592 compounds were tested through medium-throughput screening assays. A subset of 14 compounds successfully inhibited the enzyme activity in the low micromolar range of at least one of the enzymes from both Trypanosoma brucei and Lesihmania major parasites (pan-inhibitors), or from both PTR1 and DHFR-TS of the same parasite (dual inhibitors). Molecular docking studies of the protein–ligand interaction focused on new scaffolds not reproducing the well-known antifolate core clearly explaining the experimental data. TCMDC-143249, classified as a benzenesulfonamide derivative by the QikProp descriptor tool, showed selective inhibition of PTR1 and growth inhibition of the kinetoplastid parasites in the 5 μM range. In our work, we enlarged the biological profile of the GSK Kinetobox and identified new core structures inhibiting selectively PTR1, effective against the kinetoplastid infectious protozoans. In perspective, we foresee the development of selective PTR1 and DHFR inhibitors for studies of drug combinations.

Characterization of the Trypanosoma brucei Pteridine Reductase Active- Site using Computational Docking and Virtual Screening Techniques

Background:

Human African trypanosomiasis is a fatal disease prevalent in approximately 36 sub-Saharan countries. Emerging reports of drug resistance in Trypanosoma brucei are a serious cause of concern as only limited drugs are available for the treatment of the disease. Pteridine reductase is an enzyme of Trypanosoma brucei.

Methods:

It plays a critical role in the pterin metabolic pathway that is absolutely essential for its survival in the human host. The success of finding a potent inhibitor in structure-based drug design lies within the ability of computational tools to efficiently and accurately dock a ligand into the binding cavity of the target protein. Here we report the computational characterization of Trypanosoma brucei pteridine reductase (Tb-PR) active-site using twenty-four high-resolution co-crystal structures with various drugs. Structurally, the Tb-PR active site can be grouped in two clusters; one with high Root Mean Square Deviation (RMSD) of atomic positions and another with low RMSD of atomic positions. These clusters provide fresh insight for rational drug design against Tb-PR. Henceforth, the effect of several factors on docking accuracy, including ligand and protein flexibility were analyzed using Fred.

Results:

The online server was used to analyze the side chain flexibility and four proteins were selected on the basis of results. The proteins were subjected to small-scale virtual screening using 85 compounds, and statistics were calculated using Bedroc and roc curves. The enrichment factor was also calculated for the proteins and scoring functions. The best scoring function was used to understand the ligand protein interactions with top common compounds of four proteins. In addition, we made a 3D structural comparison between the active site of Tb-PR and Leishmania major pteridine reductase (Lm- PR). We described key structural differences between Tb-PR and Lm-PR that can be exploited for rational drug design against these two human parasites.

Conclusion:

The results indicated that relying just on re-docking and cross-docking experiments for virtual screening of libraries isn’t enough and results might be misleading. Hence it has been suggested that small scale virtual screening should be performed prior to large scale screening.

Accelerating Drug Discovery Efforts for Trypanosomatidic Infections Using an Integrated Transnational Academic Drug Discovery Platform

According to the World Health Organization, more than 1 billion people are at risk of or are affected by neglected tropical diseases. Examples of such diseases include trypanosomiasis, which causes sleeping sickness; leishmaniasis; and Chagas disease, all of which are prevalent in Africa, South America, and India. Our aim within the New Medicines for Trypanosomatidic Infections project was to use (1) synthetic and natural product libraries, (2) screening, and (3) a preclinical absorption, distribution, metabolism, and excretion-toxicity (ADME-Tox) profiling platform to identify compounds that can enter the trypanosomatidic drug discovery value chain. The synthetic compound libraries originated from multiple scaffolds with known antiparasitic activity and natural products from the Hypha Discovery MycoDiverse natural products library. Our focus was first to employ target-based screening to identify inhibitors of the protozoan Trypanosoma brucei pteridine reductase 1 ( TbPTR1) and second to use a Trypanosoma brucei phenotypic assay that made use of the T. brucei brucei parasite to identify compounds that inhibited cell growth and caused death. Some of the compounds underwent structure-activity relationship expansion and, when appropriate, were evaluated in a preclinical ADME-Tox assay panel. This preclinical platform has led to the identification of lead-like compounds as well as validated hits in the trypanosomatidic drug discovery value chain.

Identification of a 2,4-diaminopyrimidine scaffold targeting Trypanosoma brucei pteridine reductase 1 from the LIBRA compound library screening campaign

The LIBRA compound library is a collection of 522 non-commercial molecules contributed by various Italian academic laboratories. These compounds have been designed and synthesized during different medicinal chemistry programs and are hosted by the Italian Institute of Technology. We report the screening of the LIBRA compound library against Trypanosoma brucei and Leishmania major pteridine reductase 1, TbPTR1 and LmPTR1. Nine compounds were active against parasitic PTR1 and were selected for cell-based parasite screening, as single agents and in combination with methotrexate (MTX). The most interesting TbPTR1 inhibitor identified was 4-(benzyloxy)pyrimidine-2,6-diamine (LIB_66). Subsequently, six new LIB_66 derivatives were synthesized to explore its Structure-Activity-Relationship (SAR) and absorption, distribution, metabolism, excretion and toxicity (ADMET) properties. The results indicate that PTR1 has a preference to bind inhibitors, which resemble its biopterin/folic acid substrates, such as the 2,4-diaminopyrimidine derivatives.

Structural Insights into the Development of Cycloguanil Derivatives as Trypanosoma brucei Pteridine-Reductase-1 Inhibitors

Cycloguanil is a known dihydrofolate-reductase (DHFR) inhibitor, but there is no evidence of its activity on pteridine reductase (PTR), the main metabolic bypass to DHFR inhibition in trypanosomatid parasites. Here, we provide experimental evidence of cycloguanil as an inhibitor of Trypanosoma brucei PTR1 (TbPTR1). A small library of cycloguanil derivatives was developed, resulting in 1 and 2a having IC₅₀ values of 692 and 186 nM, respectively, toward TbPTR1. Structural analysis revealed that the increased potency of 1 and 2a is due to the combined contributions of hydrophobic interactions, H-bonds, and halogen bonds. Moreover, in vitro cell-growth-inhibition tests indicated that 2a is also effective on T. brucei. The simultaneous inhibition of DHFR and PTR1 activity in T. brucei is a promising new strategy for the treatment of human African trypanosomiasis. For this purpose, 1,6-dihydrotriazines represent new molecular tools to develop potent dual PTR and DHFR inhibitors.

Mutualist-Provisioned Resources Impact Vector Competency

Parasites elicit several physiological changes in their host to enhance transmission. Little is known about the functional association between parasitism and microbiota-provisioned resources typically dedicated to animal hosts and how these goods may be rerouted to optimize parasite development. This study is the first to identify a specific symbiont-generated metabolite that impacts insect vector competence by facilitating parasite establishment and, thus, eventual transmission. Specifically, we demonstrate that the tsetse fly obligate mutualist Wigglesworthia provisions folate (vitamin B₉) that pathogenic African trypanosomes exploit in an effort to successfully establish an infection in the vector’s MG. This process is essential for the parasite to complete its life cycle and be transmitted to a new vertebrate host. Disrupting metabolic contributions provided by the microbiota of arthropod disease vectors may fuel future innovative control strategies while also offering minimal nontarget effects.

Many symbionts supplement their host’s diet with essential nutrients. However, whether these nutrients also enhance parasitism is unknown. In this study, we investigated whether folate (vitamin B₉) production by the tsetse fly (Glossina spp.) essential mutualist, Wigglesworthia, aids auxotrophic African trypanosomes in completing their life cycle within this obligate vector. We show that the expression of Wigglesworthia folate biosynthesis genes changes with the progression of trypanosome infection within tsetse. The disruption of Wigglesworthia folate production caused a reduction in the percentage of flies that housed midgut (MG) trypanosome infections. However, decreased folate did not prevent MG trypanosomes from migrating to and establishing an infection in the fly’s salivary glands, thus suggesting that nutrient requirements vary throughout the trypanosome life cycle. We further substantiated that trypanosomes rely on symbiont-generated folate by feeding this vitamin to Glossina brevipalpis, which exhibits low trypanosome vector competency and houses Wigglesworthia incapable of producing folate. Folate-supplemented G. brevipalpis flies were significantly more susceptible to trypanosome infection, further demonstrating that this vitamin facilitates parasite infection establishment. Our cumulative results provide evidence that Wigglesworthia provides a key metabolite (folate) that is “hijacked” by trypanosomes to enhance their infectivity, thus indirectly impacting tsetse species vector competency. Parasite dependence on symbiont-derived micronutrients, which likely also occurs in other arthropod vectors, represents a relationship that may be exploited to reduce disease transmission.

Development of Chemical Proteomics for the Folateome and Analysis of the Kinetoplastid Folateome

The folate pathway has been extensively studied in a number of organisms, with its essentiality exploited by a number of drugs. However, there has been little success in developing drugs that target folate metabolism in the kinetoplastids. Despite compounds being identified which show significant inhibition of the parasite enzymes, this activity does not translate well into cellular and animal models of disease. Understanding to which enzymes antifolates bind under physiological conditions and how this corresponds to the phenotypic response could provide insight on how to target the folate pathway in these organisms. To facilitate this, we have adopted a chemical proteomics approach to study binding of compounds to enzymes of folate metabolism. Clinical and literature antifolate compounds were immobilized onto resins to allow for “pull down” of the proteins in the “folateome”. Using competition studies, proteins, which bind the beads specifically and nonspecifically, were identified in parasite lysate (Trypanosoma brucei and Leishmania major) for each antifolate compound. Proteins were identified through tryptic digest, tandem mass tag (TMT) labeling of peptides followed by LC-MS/MS. This approach was further exploited by creating a combined folate resin (folate beads). The resin could pull down up to 9 proteins from the folateome. This information could be exploited in gaining a better understanding of folate metabolism in kinetoplastids and other organisms.

Exploiting the 2-Amino-1,3,4-thiadiazole Scaffold To Inhibit Trypanosoma brucei Pteridine Reductase in Support of Early-Stage Drug Discovery

Pteridine reductase-1 (PTR1) is a promising drug target for the treatment of trypanosomiasis. We investigated the potential of a previously identified class of thiadiazole inhibitors of Leishmania major PTR1 for activity against Trypanosoma brucei (Tb). We solved crystal structures of several TbPTR1-inhibitor complexes to guide the structure-based design of new thiadiazole derivatives. Subsequent synthesis and enzyme- and cell-based assays confirm new, mid-micromolar inhibitors of TbPTR1 with low toxicity. In particular, compound 4m, a biphenyl-thiadiazole-2,5-diamine with IC₅₀ = 16 μM, was able to potentiate the antitrypanosomal activity of the dihydrofolate reductase inhibitor methotrexate (MTX) with a 4.1-fold decrease of the EC₅₀ value. In addition, the antiparasitic activity of the combination of 4m and MTX was reversed by addition of folic acid. By adopting an efficient hit discovery platform, we demonstrate, using the 2-amino-1,3,4-thiadiazole scaffold, how a promising tool for the development of anti-T. brucei agents can be obtained.

The Potential of Secondary Metabolites from Plants as Drugs or Leads against Protozoan Neglected Diseases-Part III: In-Silico Molecular Docking Investigations

Malaria, leishmaniasis, Chagas disease, and human African trypanosomiasis continue to cause considerable suffering and death in developing countries. Current treatment options for these parasitic protozoal diseases generally have severe side effects, may be ineffective or unavailable, and resistance is emerging. There is a constant need to discover new chemotherapeutic agents for these parasitic infections, and natural products continue to serve as a potential source. This review presents molecular docking studies of potential phytochemicals that target key protein targets in Leishmania spp., Trypanosoma spp., and Plasmodium spp.

The Role of Folate Transport in Antifolate Drug Action in Trypanosoma brucei

The aim of this study was to identify and characterize mechanisms of resistance to antifolate drugs in African trypanosomes. Genome-wide RNAi library screens were undertaken in bloodstream form Trypanosoma brucei exposed to the antifolates methotrexate and raltitrexed. In conjunction with drug susceptibility and folate transport studies, RNAi knockdown was used to validate the functions of the putative folate transporters. The transport kinetics of folate and methotrexate were further characterized in whole cells. RNA interference target sequencing experiments identified a tandem array of genes encoding a folate transporter family, TbFT1-3, as major contributors to antifolate drug uptake. RNAi knockdown of TbFT1-3 substantially reduced folate transport into trypanosomes and reduced the parasite's susceptibly to the classical antifolates methotrexate and raltitrexed. In contrast, knockdown of TbFT1-3 increased susceptibly to the non-classical antifolates pyrimethamine and nolatrexed. Both folate and methotrexate transport were inhibited by classical antifolates but not by non-classical antifolates or biopterin. Thus, TbFT1-3 mediates the uptake of folate and classical antifolates in trypanosomes, and TbFT1-3 loss-of-function is a mechanism of antifolate drug resistance.

Profiling of Flavonol Derivatives for the Development of Antitrypanosomatidic Drugs

Flavonoids represent a potential source of new antitrypanosomatidic leads. Starting from a library of natural products, we combined target-based screening on pteridine reductase 1 with phenotypic screening on Trypanosoma brucei for hit identification. Flavonols were identified as hits, and a library of 16 derivatives was synthesized. Twelve compounds showed EC50 values against T. brucei below 10 μM. Four X-ray crystal structures and docking studies explained the observed structure-activity relationships. Compound 2 (3,6-dihydroxy-2-(3-hydroxyphenyl)-4H-chromen-4-one) was selected for pharmacokinetic studies. Encapsulation of compound 2 in PLGA nanoparticles or cyclodextrins resulted in lower in vitro toxicity when compared to the free compound. Combination studies with methotrexate revealed that compound 13 (3-hydroxy-6-methoxy-2-(4-methoxyphenyl)-4H-chromen-4-one) has the highest synergistic effect at concentration of 1.3 μM, 11.7-fold dose reduction index and no toxicity toward host cells. Our results provide the basis for further chemical modifications aimed at identifying novel antitrypanosomatidic agents showing higher potency toward PTR1 and increased metabolic stability.

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.

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.

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.

Trypanosoma brucei pteridine reductase 1 is essential for survival in vitro and for virulence in mice

Gene knockout and knockdown methods were used to examine essentiality of pteridine reductase (PTR1) in pterin metabolism in the African trypanosome. Attempts to generate PTR1 null mutants in bloodstream form Trypanosoma brucei proved unsuccessful; despite integration of drug selectable markers at the target locus, the gene for PTR1 was either retained at the same locus or elsewhere in the genome. However, RNA interference (RNAi) resulted in complete knockdown of endogenous protein after 48 h, followed by cell death after 4 days. This lethal phenotype was reversed by expression of enzymatically active Leishmania major PTR1 in RNAi lines ((oe)RNAi) or by addition of tetrahydrobiopterin to cultures. Loss of PTR1 was associated with gross morphological changes due to a defect in cytokinesis, resulting in cells with multiple nuclei and kinetoplasts, as well as multiple detached flagella. Electron microscopy also revealed increased numbers of glycosomes, while immunofluorescence microscopy showed increased and more diffuse staining for glycosomal matrix enzymes, indicative of mis-localisation to the cytosol. Mis-localisation was confirmed by digitonin fractionation experiments. RNAi cell lines were markedly less virulent than wild-type parasites in mice and virulence was restored in the (oe)RNAi line. Thus, PTR1 may be a drug target for human African trypanosomiasis.