Genetically Validated Drug Targets in Leishmania: Current Knowledge and Future Prospects

Py-CoMFA, docking, and molecular dynamics simulations of Leishmania (L.) amazonensis arginase inhibitors

Leishmaniasis is a disease caused by a protozoan of the genus Leishmania, affecting millions of people, mainly in tropical countries, due to poor social conditions and low economic development. First-line chemotherapeutic agents involve highly toxic pentavalent antimonials, while treatment failure is mainly due to the emergence of drug-resistant strains. Leishmania arginase (ARG) enzyme is vital in pathogenicity and contributes to a higher infection rate, thus representing a potential drug target. This study helps in designing ARG inhibitors for the treatment of leishmaniasis. Py-CoMFA (3D-QSAR) models were constructed using 34 inhibitors from different chemical classes against ARG from L. (L.) amazonensis (LaARG). The 3D-QSAR predictions showed an excellent correlation between experimental and calculated pIC50 values. The molecular docking study identified the favorable hydrophobicity contribution of phenyl and cyclohexyl groups as substituents in the enzyme allosteric site. Molecular dynamics simulations of selected protein-ligand complexes were conducted to understand derivatives' interaction modes and affinity in both active and allosteric sites. Two cinnamide compounds, 7g and 7k, were identified, with similar structures to the reference 4h allosteric site inhibitor. These compounds can guide the development of more effective arginase inhibitors as potential antileishmanial drugs.

Molecular tools to regulate gene expression in Trypanosoma cruzi

Developing molecular strategies to manipulate gene expression in trypanosomatids is challenging, particularly with respect to the unique gene expression mechanisms adopted by these unicellular parasites, such as polycistronic mRNA transcription and multi-gene families. In the case of Trypanosoma cruzi (T. cruzi), the causative agent of Chagas Disease, the lack of RNA interference machinery further complicated functional genetic studies important for understanding parasitic biology and developing biomarkers and potential therapeutic targets. Therefore, alternative methods of performing knockout and/or endogenous labelling experiments were developed to identify and understand the function of proteins for survival and interaction with the host. In this review, we present the main tools for the genetic manipulation of T. cruzi, focusing on the Clustered Regularly Interspaced Short Palindromic Repeats Cas9-associated system technique widely used in this organism. Moreover, we highlight the importance of using these tools to elucidate the function of uncharacterized and glycosylated proteins. Further developments of these technologies will allow the identification of new biomarkers, therapeutic targets and potential vaccines against Chagas disease with greater efficiency and speed.

TORC1 is an essential regulator of nutrient-controlled proliferation and differentiation in Leishmania

Leishmania parasites undergo differentiation between various proliferating and non-dividing forms to adapt to changing host environments. The mechanisms that link environmental cues with the parasite's developmental changes remain elusive. Here, we report that Leishmania TORC1 is a key environmental sensor for parasite proliferation and differentiation in the sand fly-stage promastigotes and for replication of mammalian-stage amastigotes. We show that Leishmania RPTOR1, interacts with TOR1 and LST8, and identify new parasite-specific proteins that interact in this complex. We investigate TORC1 function by conditional deletion of RPTOR1, where under nutrient-rich conditions RPTOR1 depletion results in decreased protein synthesis and growth, G1 cell cycle arrest and premature differentiation from proliferative promastigotes to non-dividing mammalian-infective metacyclic forms. These parasites are unable to respond to nutrients to differentiate into proliferative retroleptomonads, which are required for their blood-meal induced amplification in sand flies and enhanced mammalian infectivity. We additionally show that RPTOR1-/- metacyclic promastigotes develop into amastigotes but do not proliferate in the mammalian host to cause pathology. RPTOR1-dependent TORC1 functionality represents a critical mechanism for driving parasite growth and proliferation.

Fragment Merging, Growing, and Linking Identify New Trypanothione Reductase Inhibitors for Leishmaniasis

Trypanothione reductase (TR) is a suitable target for drug discovery approaches against leishmaniasis, although the identification of potent inhibitors is still challenging. Herein, we harnessed a fragment-based drug discovery (FBDD) strategy to develop new TR inhibitors. Previous crystallographic screening identified fragments 1-3, which provided ideal starting points for a medicinal chemistry campaign. In silico investigations revealed critical hotspots in the TR binding site, guiding our structure- and ligand-based structure-actvity relationship (SAR) exploration that yielded fragment-derived compounds 4-14. A trend of improvement in Leishmania infantum TR inhibition was detected along the optimization and confirmed by the crystal structures of 9, 10, and 14 in complex with Trypanosoma brucei TR. Compound 10 showed the best TR inhibitory profile (Ki = 0.2 μM), whereas 9 was the best one in terms of in vitro and ex vivo activity. Although further fine-tuning is needed to improve selectivity, we demonstrated the potentiality of FBDD on a classic but difficult target for leishmaniasis.

Repurposing of conformationally-restricted cyclopentane-based AKT-inhibitors leads to discovery of potential and more selective antileishmanial agents than miltefosine

Conformational restriction was addressed towards the development of more selective and effective antileishmanial agents than currently used drugs for treatment of Leishmania donovani; the causative parasite of the fatal visceral leishmaniasis. Five types of cyclopentane-based conformationally restricted miltefosine analogs that were previously explored in literature as anticancer AKT-inhibitors were reprepared and repurposed as antileishmanial agents. Amongst, positions-1 and 2 cis-conformationally-restricted compound 1a and positions-2 and 3 trans-conformationally-restricted compound 3b were highly potent eliciting sub-micromolar IC50 values for inhibition of infection and inhibition of parasite number compared with the currently used miltefosine drug that showed low micromolar IC50 values for inhibition of infection and inhibition of parasite number. Compounds 1a and 3b eradicated the parasite without triggering host cells cytotoxicity over more than one log concentration interval which is a superior performance compared to miltefosine. In silico studies suggested that conformational restriction conserved the conformer capable of binding LdAKT-like kinase while it might be possible that it excludes other conformers mediating undesirable effects and/or toxicity of miltefosine. Together, this study presents compounds 1a and 3b as antileishmanial agents with superior performance over the currently used miltefosine drug.

Bromodomain Factor 5 as a Target for Antileishmanial Drug Discovery

Leishmaniases are a collection of neglected tropical diseases caused by kinetoplastid parasites in the genus Leishmania. Current chemotherapies are severely limited, and the need for new antileishmanials is of pressing international importance. Bromodomains are epigenetic reader domains that have shown promising therapeutic potential for cancer therapy and may also present an attractive target to treat parasitic diseases. Here, we investigate Leishmania donovani bromodomain factor 5 (LdBDF5) as a target for antileishmanial drug discovery. LdBDF5 contains a pair of bromodomains (BD5.1 and BD5.2) in an N-terminal tandem repeat. We purified recombinant bromodomains of L. donovani BDF5 and determined the structure of BD5.2 by X-ray crystallography. Using a histone peptide microarray and fluorescence polarization assay, we identified binding interactions of LdBDF5 bromodomains with acetylated peptides derived from histones H2B and H4. In orthogonal biophysical assays including thermal shift assays, fluorescence polarization, and NMR, we showed that BDF5 bromodomains bind to human bromodomain inhibitors SGC-CBP30, bromosporine, and I-BRD9; moreover, SGC-CBP30 exhibited activity against Leishmania promastigotes in cell viability assays. These findings exemplify the potential BDF5 holds as a possible drug target in Leishmania and provide a foundation for the future development of optimized antileishmanial compounds targeting this epigenetic reader protein.

In Silico Exploration of the Trypanothione Reductase (TryR) of L. mexicana

Human leishmaniasis is a neglected tropical disease which affects nearly 1.5 million people every year, with Mexico being an important endemic region. One of the major defense mechanisms of these parasites is based in the polyamine metabolic pathway, as it provides the necessary compounds for its survival. Among the enzymes in this route, trypanothione reductase (TryR), an oxidoreductase enzyme, is crucial for the Leishmania genus' survival against oxidative stress. Thus, it poses as an attractive drug target, yet due to the size and features of its catalytic pocket, modeling techniques such as molecular docking focusing on that region is not convenient. Herein, we present a computational study using several structure-based approaches to assess the druggability of TryR from L. mexicana, the predominant Leishmania species in Mexico, beyond its catalytic site. Using this consensus methodology, three relevant pockets were found, of which the one we call σ-site promises to be the most favorable one. These findings may help the design of new drugs of trypanothione-related diseases.

Gene editing of putative cAMP and Ca2+ -regulated proteins using an efficient cloning-free CRISPR/Cas9 system in Trypanosoma cruzi

Trypanosoma cruzi, the agent of Chagas disease, must adapt to a diversity of environmental conditions that it faces during its life cycle. The adaptation to these changes is mediated by signaling pathways that coordinate the cellular responses to the new environmental settings. Cyclic AMP (cAMP) and Calcium (Ca2+ ) signaling pathways regulate critical cellular processes in this parasite, such as differentiation, osmoregulation, host cell invasion and cell bioenergetics. Although the use of CRISPR/Cas9 technology prompted reverse genetics approaches for functional analysis in T. cruzi, it is still necessary to expand the toolbox for genome editing in this parasite, as for example to perform multigene analysis. Here we used an efficient T7RNAP/Cas9 strategy to tag and delete three genes predicted to be involved in cAMP and Ca2+ signaling pathways: a putative Ca2+ /calmodulin-dependent protein kinase (CAMK), Flagellar Member 6 (FLAM6) and Cyclic nucleotide-binding domain/C2 domain-containing protein (CC2CP). We endogenously tagged these three genes and determined the subcellular localization of the tagged proteins. Furthermore, the strategy used to knockout these genes allows us to presume that TcCC2CP is an essential gene in T. cruzi epimastigotes. Our results will open new venues for future research on the role of these proteins in T. cruzi.

Gene editing of putative cAMP and Ca 2+ -regulated proteins using an efficient cloning-free CRISPR/Cas9 system in Trypanosoma cruzi

Trypanosoma cruzi , the agent of Chagas disease, must adapt to a diversity of environmental conditions that it faces during its life cycle. The adaptation to these changes is mediated by signaling pathways that coordinate the cellular responses to the new environmental settings. Cyclic AMP (cAMP) and Calcium (Ca 2+ ) signaling pathways regulate critical cellular processes in this parasite, such as differentiation, osmoregulation, host cell invasion and cell bioenergetics. Although the use of CRISPR/Cas9 technology prompted reverse genetics approaches for functional analysis in T. cruzi , it is still necessary to expand the toolbox for genome editing in this parasite, as for example to perform multigene analysis. Here we used an efficient T7RNAP/Cas9 strategy to tag and delete three genes predicted to be involved in cAMP and Ca 2+ signaling pathways: a putative Ca 2+ /calmodulin-dependent protein kinase ( CAMK ), Flagellar Member 6 ( FLAM6 ) and Cyclic nucleotide-binding domain/C2 domain-containing protein ( CC2CP ). We endogenously tagged these three genes and determined the subcellular localization of the tagged proteins. Furthermore, the strategy used to knockout these genes allow us to presume that TcCC2CP is an essential gene in T. cruzi epimastigotes. Our results will open new venues for future research on the role of these proteins in T. cruzi .

CYP5122A1 encodes an essential sterol C4-methyl oxidase in Leishmania donovani and determines the antileishmanial activity of antifungal azoles

Visceral leishmaniasis, caused by Leishmania donovani, is a life-threatening parasitic disease, but current antileishmanial drugs are limited and have severe drawbacks. There have been efforts to repurpose antifungal azole drugs for the treatment of Leishmania infection. Antifungal azoles are known to potently inhibit the activity of cytochrome P450 (CYP) 51 enzymes which are responsible for removing the C14α-methyl group of lanosterol, a key step in ergosterol biosynthesis in Leishmania. However, they exhibit varying degrees of antileishmanial activities in culture, suggesting the existence of unrecognized molecular targets for these compounds. Our previous study reveals that, in Leishmania, lanosterol undergoes parallel C4- and C14-demethylation reactions to form 4α,14α-dimethylzymosterol and T-MAS, respectively. In the current study, CYP5122A1 is identified as a sterol C4-methyl oxidase that catalyzes the sequential oxidation of lanosterol to form C4-oxidation metabolites. CYP5122A1 is essential for both L. donovani promastigotes in culture and intracellular amastigotes in infected mice. Overexpression of CYP5122A1 results in growth delay, differentiation defects, increased tolerance to stress, and altered expression of lipophosphoglycan and proteophosphoglycan. CYP5122A1 also helps to determine the antileishmanial effect of antifungal azoles in vitro. Dual inhibitors of CYP51 and CYP5122A1, e.g., clotrimazole and posaconazole, possess superior antileishmanial activity against L. donovani promastigotes whereas CYP51-selective inhibitors, e.g., fluconazole and voriconazole, have little effect on promastigote growth. Our findings uncover the critical biochemical and biological role of CYP5122A1 in L. donovani and provide an important foundation for developing new antileishmanial drugs by targeting both CYP enzymes.

Next-Generation Leishmanization: Revisiting Molecular Targets for Selecting Genetically Engineered Live-Attenuated Leishmania

Despite decades of research devoted to finding a vaccine against leishmaniasis, we are still lacking a safe and effective vaccine for humans. Given this scenario, the search for a new prophylaxis alternative for controlling leishmaniasis should be a global priority. Inspired by leishmanization-a first generation vaccine strategy where live L. major parasites are inoculated in the skin to protect against reinfection-live-attenuated Leishmania vaccine candidates are promising alternatives due to their robust elicited protective immune response. In addition, they do not cause disease and could provide long-term protection upon challenge with a virulent strain. The discovery of a precise and easy way to perform CRISPR/Cas-based gene editing allowed the selection of safer null mutant live-attenuated Leishmania parasites obtained by gene disruption. Here, we revisited molecular targets associated with the selection of live-attenuated vaccinal strains, discussing their function, their limiting factors and the ideal candidate for the next generation of genetically engineered live-attenuated Leishmania vaccines to control leishmaniasis.

Inhibiting Leishmania donovani Sterol Methyltransferase to Identify Lead Compounds Using Molecular Modelling

The recent outlook of leishmaniasis as a global public health concern coupled with the reportage of resistance and lack of efficacy of most antileishmanial drugs calls for a concerted effort to find new leads. The study combined In silico and in vitro approaches to identify novel potential synthetic small-molecule inhibitors targeting the Leishmania donovani sterol methyltransferase (LdSMT). The LdSMT enzyme in the ergosterol biosynthetic pathway is required for the parasite’s membrane fluidity, distribution of membrane proteins, and control of the cell cycle. The lack of LdSMT homologue in the human host and its conserved nature among all Leishmania parasites makes it a viable target for future antileishmanial drugs. Initially, six known inhibitors of LdSMT with IC50 < 10 μM were used to generate a pharmacophore model with a score of 0.9144 using LigandScout. The validated model was used to screen a synthetic library of 95,630 compounds obtained from InterBioScreen limited. Twenty compounds with pharmacophore fit scores above 50 were docked against the modelled three-dimensional structure of LdSMT using AutoDock Vina. Consequently, nine compounds with binding energies ranging from −7.5 to −8.7 kcal/mol were identified as potential hit molecules. Three compounds comprising STOCK6S-06707, STOCK6S-84928, and STOCK6S-65920 with respective binding energies of −8.7, −8.2, and −8.0 kcal/mol, lower than 22,26-azasterol (−7.6 kcal/mol), a known LdSMT inhibitor, were selected as plausible lead molecules. Molecular dynamics simulation studies and molecular mechanics Poisson–Boltzmann surface area calculations showed that the residues Asp25 and Trp208 were critical for ligand binding. The compounds were also predicted to have antileishmanial activity with reasonable pharmacological and toxicity profiles. When the antileishmanial activity of the three hits was evaluated in vitro against the promastigotes of L. donovani, mean half-maximal inhibitory concentrations (IC50) of 21.9 ± 1.5 μM (STOCK6S-06707), 23.5 ± 1.1 μM (STOCK6S-84928), and 118.3 ± 5.8 μM (STOCK6S-65920) were obtained. Furthermore, STOCK6S-84928 and STOCK6S-65920 inhibited the growth of Trypanosoma brucei, with IC50 of 14.3 ± 2.0 μM and 18.1 ± 1.4 μM, respectively. The identified compounds could be optimised to develop potent antileishmanial therapeutic agents.

Anti-trypanosomatid drug discovery: progress and challenges

Leishmaniasis (visceral and cutaneous), Chagas disease and human African trypanosomiasis cause substantial death and morbidity, particularly in low- and middle-income countries. Although the situation has improved for human African trypanosomiasis, there remains an urgent need for new medicines to treat leishmaniasis and Chagas disease; the clinical development pipeline is particularly sparse for Chagas disease. In this Review, we describe recent advances in our understanding of the biology of the causative pathogens, particularly from the drug discovery perspective, and we explore the progress that has been made in the development of new drug candidates and the identification of promising molecular targets. We also explore the challenges in developing new clinical candidates and discuss potential solutions to overcome such hurdles.

Molecular Modeling, Virtual Screening, and Molecular Dynamics for Leishmania infantum Methionyl-tRNA Synthetase

Visceral leishmaniasis is a neglected tropical disease (NTD) caused by Leishmania infantum and L. donovani that is lethal in cases of nontreatment. The treatments are limited by serious drawbacks involving safety, resistance, stability, and high costs. In this work, we aimed to identify inhibitors of Leishmania infantum methionyl-tRNA synthetase (LiMetRS), a validated molecular target for leishmaniasis drug discovery, using a combination of strategies. A virtual database of compounds was organized by filtering compounds from the ZINC15 database. Homology modeling was used to obtain the structure of LiMetRS based on the crystal coordinates of the enzyme from Trypanosoma brucei (TbMetRS). A virtual screening using molecular docking identified 10 candidate compounds from among more than 5 million that were included in the initial database. The selected hits were further evaluated using a script created in this work to select only the ligands that interacted with specific amino acids in the catalytic site of the enzyme. Furthermore, suitable pharmacokinetic profiles were predicted for the selected compounds, especially a good balance between aqueous solubility and lipophilic character, no ability to cross the blood-brain barrier, good oral absorption, and no liability toward P-gp efflux for most compounds. Six compounds were then subjected to all-atom molecular dynamics. Two compounds showed good stability when bound to the leishmanial enzyme, which provided a deeper understanding of the structural differences between TbMetRS and LiMetRS that can guide further drug discovery efforts for visceral leishmaniasis.

Expanded Proteomic Survey of the Human Parasite Leishmania major Focusing on Changes in Null Mutants of the Golgi GDP-Mannose/Fucose/Arabinopyranose Transporter LPG2 and of the Mitochondrial Fucosyltransferase FUT1

The trypanosomatid protozoan parasite Leishmania has a significant impact on human health globally. Understanding the pathways associated with virulence within this significant pathogen is critical for identifying novel vaccination and chemotherapy targets. Within this study we leverage an ultradeep proteomic approach to improve our understanding of two virulence-associated genes in Leishmania, encoding the Golgi mannose/arabinopyranose/fucose nucleotide-sugar transporter (LPG2) and the mitochondrial fucosyltransferase (FUT1). Using deep peptide fractionation followed by complementary fragmentation approaches with higher-energy collisional dissociation (HCD) and electron transfer dissociation (ETD) allowed the identification of over 6,500 proteins, nearly doubling the experimentally known Leishmania major proteome. This deep proteomic analysis revealed significant quantitative differences in both Δlpg2- and Δfut1s mutants with FUT1-dependent changes linked to marked alterations within mitochondrion-associated proteins, while LPG2-dependent changes impacted many pathways, including the secretory pathway. While the FUT1 enzyme has been shown to fucosylate peptides in vitro, no evidence for protein fucosylation was identified within our ultradeep analysis, nor did we observe fucosylated glycans within Leishmania glycopeptides isolated using hydrophilic interaction liquid chromatography (HILIC) enrichment. This work provides a critical resource for the community on the observable Leishmania proteome as well as highlighting phenotypic changes associated with LPG2 or FUT1, ablation of which may guide the development of future therapeutics. IMPORTANCE Leishmania is a widespread trypanosomatid protozoan parasite of humans, with ~12 million cases currently, ranging from mild to fatal, and hundreds of millions asymptomatically infected. This work advances knowledge of the experimental proteome by nearly 2-fold, to more than 6,500 proteins and thus provides a great resource to investigators seeking to decode how this parasite is transmitted and causes disease and to identify new targets for therapeutic intervention. The ultradeep proteomics approach identified potential proteins underlying the "persistence-without-pathology" phenotype of mutants with deletion of the Golgi nucleotide transporter LPG2, showing many alterations and several candidates. Studies of a rare mutant with deletion of the mitochondrial fucosyltransferase FUT1 revealed changes underlying its strong mitochondrial dysfunction but did not reveal examples of fucosylation of either peptides or N-glycans. This suggests that this vital protein's elusive target(s) may be more complex than the methods used could detect or that this target may not be a protein but perhaps another glycoconjugate or glycolipid.

Genome deletions to overcome the directed loss of gene function in Leishmania

With the global reach of the Neglected Tropical Disease leishmaniasis increasing, coupled with a tiny armory of therapeutics which all have problems with resistance, cost, toxicity and/or administration, the validation of new drug targets in the causative insect vector borne protozoa Leishmania spp is more important than ever. Before the introduction of CRISPR Cas9 technology in 2015 genetic validation of new targets was carried out largely by targeted gene knockout through homologous recombination, with the majority of genes targeted (~70%) deemed non-essential. In this study we exploit the ready availability of whole genome sequencing technology to reanalyze one of these historic cell lines, a L. major knockout in the catalytic subunit of serine palmitoyltransferase (LCB2), which causes a complete loss of sphingolipid biosynthesis but remains viable and infective. This revealed a number of Single Nucleotide Polymorphisms, but also the complete loss of several coding regions including a gene encoding a putative ABC3A orthologue, a putative sterol transporter. Hypothesizing that the loss of such a transporter may have facilitated the directed knockout of the catalytic subunit of LCB2 and the complete loss of de novo sphingolipid biosynthesis, we re-examined LCB2 in a L. mexicana line engineered for straightforward CRISPR Cas9 directed manipulation. Strikingly, LCB2 could not be knocked out indicating essentiality. However, simultaneous deletion of LCB2 and the putative ABC3A was possible. This indicated that the loss of the putative ABC3A facilitated the loss of sphingolipid biosynthesis in Leishmania, and suggested that we should re-examine the many other Leishmania knockout lines where genes were deemed non-essential.

Prediction of Antileishmanial Compounds: General Model, Preparation, and Evaluation of 2-Acylpyrrole Derivatives

In this work, the SOFT.PTML tool has been used to pre-process a ChEMBL dataset of pre-clinical assays of antileishmanial compound candidates. A comparative study of different ML algorithms, such as logistic regression (LOGR), support vector machine (SVM), and random forests (RF), has shown that the IFPTML-LOGR model presents excellent values of specificity and sensitivity (81-98%) in training and validation series. The use of this software has been illustrated with a practical case study focused on a series of 28 derivatives of 2-acylpyrroles 5a,b, obtained through a Pd(II)-catalyzed C-H radical acylation of pyrroles. Their in vitro leishmanicidal activity against visceral (L. donovani) and cutaneous (L. amazonensis) leishmaniasis was evaluated finding that compounds 5bc (IC₅₀ = 30.87 μM, SI > 10.17) and 5bd (IC₅₀ = 16.87 μM, SI > 10.67) were approximately 6-fold more selective than the drug of reference (miltefosine) in in vitro assays against L. amazonensis promastigotes. In addition, most of the compounds showed low cytotoxicity, CC₅₀ > 100 μg/mL in J774 cells. Interestingly, the IFPMTL-LOGR model predicts correctly the relative biological activity of these series of acylpyrroles. A computational high-throughput screening (cHTS) study of 2-acylpyrroles 5a,b has been performed calculating >20,700 activity scores vs a large space of 647 assays involving multiple Leishmania species, cell lines, and potential target proteins. Overall, the study demonstrates that the SOFT.PTML all-in-one strategy is useful to obtain IFPTML models in a friendly interface making the work easier and faster than before. The present work also points to 2-acylpyrroles as new lead compounds worthy of further optimization as antileishmanial hits.

Discovery of Leishmania Druggable Serine Proteases by Activity-Based Protein Profiling

Leishmaniasis are a group of diseases caused by parasitic protozoa of the genus Leishmania. Current treatments are limited by difficult administration, high cost, poor efficacy, toxicity, and growing resistance. New agents, with new mechanisms of action, are urgently needed to treat the disease. Although extensively studied in other organisms, serine proteases (SPs) have not been widely explored as antileishmanial drug targets. Herein, we report for the first time an activity-based protein profiling (ABPP) strategy to investigate new therapeutic targets within the SPs of the Leishmania parasites. Active-site directed fluorophosphonate probes (rhodamine and biotin-conjugated) were used for the detection and identification of active Leishmania serine hydrolases (SHs). Significant differences were observed in the SHs expression levels throughout the Leishmania life cycle and between different Leishmania species. Using iTRAQ-labelling-based quantitative proteomic mass spectrometry, we identified two targetable SPs in Leishmania mexicana: carboxypeptidase LmxM.18.0450 and prolyl oligopeptidase LmxM.36.6750. Druggability was ascertained by selective inhibition using the commercial serine protease inhibitors chymostatin, lactacystin and ZPP, which represent templates for future anti-leishmanial drug discovery programs. Collectively, the use of ABPP method complements existing genetic methods for target identification and validation in Leishmania.

Bromodomain factor 5 is an essential regulator of transcription in Leishmania

Leishmania are unicellular parasites that cause human and animal diseases. Like other kinetoplastids, they possess large transcriptional start regions (TSRs) which are defined by histone variants and histone lysine acetylation. Cellular interpretation of these chromatin marks is not well understood. Eight bromodomain factors, the reader modules for acetyl-lysine, are found across Leishmania genomes. Using L. mexicana, Cas9-driven gene deletions indicate that BDF1–5 are essential for promastigotes. Dimerisable, split Cre recombinase (DiCre)-inducible gene deletion of BDF5 show it is essential for both promastigotes and murine infection. ChIP-seq identifies BDF5 as enriched at TSRs. XL-BioID proximity proteomics shows the BDF5 landscape is enriched for BDFs, HAT2, proteins involved in transcriptional activity, and RNA processing; revealing a Conserved Regulators of Kinetoplastid Transcription (CRKT) Complex. Inducible deletion of BDF5 causes global reduction in RNA polymerase II transcription. Our results indicate the requirement of Leishmania to interpret histone acetylation marks through the bromodomain-enriched CRKT complex for normal gene expression and cellular viability.

Leishmania use large (5–10 kb) transcriptional start regions, where the chromatin is highly enriched for acetylated histones, to drive the expression of polycistronic gene arrays. Here the authors show bromodomain-containing protein BDF5 is enriched at transcriptional start sites and its depletion leads to cell death in vitro and in murine infections, and they identify its interactors.

In silico analysis of Leishmania proteomes and protein-protein interaction network: Prioritizing therapeutic targets and drugs for repurposing to treat leishmaniasis

Leishmaniasis is a serious world health problem and its current therapies have several limitations demanding to develop novel therapeutics for this disease. The present study aims to prioritize novel broad-spectrum targets using proteomics and protein-protein interaction network (PPIN) data for 11 Leishmania species. Proteome comparison and host non-homology analysis resulted in 3605 pathogen-specific conserved core proteins. Gene ontology analysis indicated their involvement in major molecular functions like DNA binding, transportation, dioxygenase, and catalytic activity. PPIN analysis of these core proteins identified eight hub proteins (viz., vesicle-trafficking protein (LBRM2903_190011800), ribosomal proteins S17 (LBRM2903_34004790) and L2 (LBRM2903_080008100), eukaryotic translation initiation factor 3 (LBRM2903_350086700), replication factor A (LBRM2903_150008000), U3 small nucleolar RNA-associated protein (LBRM2903_340025600), exonuclease (LBRM2903_200021800), and mitochondrial RNA ligase (LBRM2903_200074100)). Among the hub proteins, six were classified as drug targets and two as vaccine candidates. Further, druggability analysis indicated three hub proteins, namely eukaryotic translation initiation factor 3, ribosomal proteins S17 and L2 as druggable. Their three-dimensional structures were modelled and docked with the identified ligands (2-methylthio-N6-isopentenyl-adenosine-5'-monophosphate, artenimol and omacetaxine mepesuccinate). These ligands could be experimentally validated (in vitro and in vivo) and repurposed for the development of novel antileishmanial agents.

Applied Machine Learning Toward Drug Discovery Enhancement: Leishmaniases as a Case Study

Drug discovery (DD) research is a complex field with a high attrition rate. Machine learning (ML) approaches combined to chemoinformatics are of valuable input to this field. We, herein, focused on implementing multiple ML algorithms that shall learn from different molecular fingerprints (FPs) of 65 057 molecules that have been identified as active or inactive against Leishmania major promastigotes. We sought to build a classifier able to predict whether a given molecule has the potential of being anti-leishmanial or not. Using the RDkit library, we calculated 5 molecular FPs of the molecules. Then, we implemented 4 ML algorithms that we trained and tested for their ability to classify the molecules into active/inactive classes based on their chemical structure, encoded by the molecular FPs. Best performers were random forest (RF) and support vector machine (SVM), while atom-pair and topology torsion FPs were the best embedding functions. Both models were further assessed on different stratification levels of the dataset and showed stable performances. At last, we used them to predict the potential of molecules within the Food and Drug Administration (FDA)-approved drugs collection to present anti-Leishmania effects. We ranked these drugs according to their anti-Leishmanial probability and obtained in total seven anti-Leishmania agents, previously described in the literature, within the top 10 of each model. This validates the robustness of the approach, the algorithms, and FPs choices as well as the importance of the dataset size and content. We further engaged these molecules into reverse docking experiments on 3D crystal structures of seven well-studied Leishmania drug targets and could predict the molecular targets for 4 drugs. The results bring novel insights into anti-Leishmania compounds.

Absence of DEATH kinesin is fatal for Leishmania mexicana amastigotes

Kinesins are motor proteins present in organisms from protists to mammals playing important roles in cell division, intracellular organisation and flagellum formation and maintenance. Leishmania mexicana is a protozoan parasite of the order Kinetoplastida causing human cutaneous leishmaniasis. Kinetoplastida genome sequence analyses revealed a large number of kinesins showing sequence and structure homology to eukaryotic kinesins. Here, we investigate the L. mexicana kinesin LmxKIN29 (LmxM.29.0350), also called DEATH kinesin. The activated MAP kinase LmxMPK3, a kinase affecting flagellum length in Leishmania, is able to phosphorylate recombinant full length LmxKIN29 at serine 554. Insect promastigote LmxKIN29 Leishmania null mutants showed no obvious phenotype. However, in mouse infection experiments, the null mutants were unable to cause the disease, whereas LmxKIN29 add-backs and single allele knockouts caused footpad lesions. Localisation using promastigotes expressing GFP-tagged LmxKIN29 revealed that the kinesin is predominantly found in between the nucleus and the flagellar pocket, while in dividing cells the GFP-fusion protein was found at the anterior and posterior ends of the cells indicating a role in cytokinesis. The inability to cause lesions in infected animals and the amino acid sequence divergence from mammalian kinesins suggests that LmxKIN29 is a potential drug target against leishmaniasis.

Challenges and Tools for In Vitro Leishmania Exploratory Screening in the Drug Development Process: An Updated Review

Leishmaniases are a group of vector-borne diseases caused by infection with the protozoan parasites Leishmania spp. Some of them, such as Mediterranean visceral leishmaniasis, are zoonotic diseases transmitted from vertebrate to vertebrate by a hematophagous insect, the sand fly. As there is an endemic in more than 90 countries worldwide, this complex and major health problem has different clinical forms depending on the parasite species involved, with the visceral form being the most worrying since it is fatal when left untreated. Nevertheless, currently available antileishmanial therapies are significantly limited (low efficacy, toxicity, adverse side effects, drug-resistance, length of treatment, and cost), so there is an urgent need to discover new compounds with antileishmanial activity, which are ideally inexpensive and orally administrable with few side effects and a novel mechanism of action. Therefore, various powerful approaches were recently applied in many interesting antileishmanial drug development programs. The objective of this review is to focus on the very first step in developing a potential drug and to identify the exploratory methods currently used to screen in vitro hit compounds and the challenges involved, particularly in terms of harmonizing the results of work carried out by different research teams. This review also aims to identify innovative screening tools and methods for more extensive use in the drug development process.

Effective Genome Editing in Leishmania (Viannia) braziliensis Stably Expressing Cas9 and T7 RNA Polymerase

Until 2015, loss-of-function studies to elucidate protein function in Leishmania relied on gene disruption through homologous recombination. Then, the CRISPR/Cas9 revolution reached these protozoan parasites allowing efficient genome editing with one round of transfection. In addition, the development of LeishGEdit, a PCR-based toolkit for generating knockouts and tagged lines using CRISPR/Cas9, allowed a more straightforward and effective genome editing. In this system, the plasmid pTB007 is delivered to Leishmania for episomal expression or integration in the β-tubulin locus and for the stable expression of T7 RNA polymerase and Cas9. In South America, and especially in Brazil, Leishmania (Viannia) braziliensis is the most frequent etiological agent of tegumentary leishmaniasis. The L. braziliensis β-tubulin locus presents significant sequence divergence in comparison with Leishmania major, which precludes the efficient integration of pTB007 and the stable expression of Cas9. To overcome this limitation, the L. major β-tubulin sequences, present in the pTB007, were replaced by a Leishmania (Viannia) β-tubulin conserved sequence generating the pTB007_Viannia plasmid. This modification allowed the successful integration of the pTB007_Viannia cassette in the L. braziliensis M2903 genome, and in silico predictions suggest that this can also be achieved in other Viannia species. The activity of Cas9 was evaluated by knocking out the flagellar protein PF16, which caused a phenotype of immobility in these transfectants. Endogenous PF16 was also successfully tagged with mNeonGreen, and an in-locus complementation strategy was employed to return a C-terminally tagged copy of the PF16 gene to the original locus, which resulted in the recovery of swimming capacity. The modified plasmid pTB007_Viannia allowed the integration and stable expression of both T7 RNA polymerase and Cas9 in L. braziliensis and provided an important tool for the study of the biology of this parasite.

The long and winding road of reverse genetics in Trypanosoma cruzi

Trypanosomes are early divergent protists with distinctive features among eukaryotic cells. Together with Trypanosoma brucei and Leishmania spp., Trypanosoma cruzi has been one of the most studied members of the group. This protozoan parasite is the causative agent of Chagas disease, a leading cause of heart disease in the Americas, for which there is no vaccine or satisfactory treatment available. Understanding T. cruzi biology is crucial to identify alternative targets for antiparasitic interventions. Genetic manipulation of T. cruzi has been historically challenging. However, the emergence of CRISPR/Cas9 technology has significantly improved the ability to generate genetically modified T. cruzi cell lines. Still, the system alone is not sufficient to answer all biologically relevant questions. In general, current genetic methods have limitations that should be overcome to advance in the study of this peculiar parasite. In this brief historic overview, we highlight the strengths and weaknesses of the molecular strategies that have been developed to genetically modify T. cruzi, emphasizing the future directions of the field.

The thermal proteome stability profile of Trypanosoma cruzi in epimastigote and trypomastigote life stages

Trypanosoma cruzi is a flagellate protozoa being the etiological agent of Chagas disease, a neglected tropical disease, which still poses a public health problem worldwide. The intricate molecular changes during T. cruzi-host interaction have been explored using different largescale omics techniques. However, protein stability is largely unknown. Thermal proteome profiling (TPP) methodology has the potential to characterize proteome-wide stability highlighting key proteins during T. cruzi infection and life stage transition from the invertebrate to the mammalian host. In the present work, T. cruzi epimastigotes and trypomastigotes cell lysates were subjected to TPP workflow and analyzed by quantitative large-scale mass spectrometry-based proteomics to fit a melting profile for each protein. A total of 2884 proteins were identified and associated to 1741 melting curves being 1370 in trypomastigotes (Tm_AVG 53.53 °C) and 1279 in epimastigotes (Tm_AVG 50.89 °C). A total of 453 proteins were identified with statistically different melting profiles between the two life stages. Proteins associated to pathogenesis and intracellular transport had regulated melting temperatures. Membrane and glycosylated proteins had a higher average Tm in trypomastigotes compared to epimastigotes. This study represents the first large-scale comparison of parasite protein stability between life stages. SIGNIFICANCE: Trypanosoma cruzi, a unicellular flagellate parasite, is the etiological agent of Chagas disease, endemic in South America and affecting more that 7 million people worldwide. There is an intense research to identify novel chemotherapeutic and diagnostic targets of Chagas disease. Proteomic approaches have helped in elucidating the quantitative proteome and PTMs changes of T. cruzi during life cycle transition and upon different biotic and abiotic stimuli. However, a comprehensive knowledge of the protein-protein interaction and protein conformation is still missing. In order to fill this gap, this manuscript elucidates the T. cruzi Y strain proteome-wide thermal stability map in the epimastigote and trypomastigote life stages. Comparison between life stages showed a higher average melting temperature stability for trypomastigotes than epimastigotes indicating a host temperature adaptation. Both presented a selective thermal stability shift for cellular compartments, molecular functions and biological processes based on the T. cruzi life stage. Membrane and glycosylated proteins presented a higher thermal stability in trypomastigotes when compared to the epimastigotes.

Palladium-mediated synthesis and biological evaluation of C-10b substituted Dihydropyrrolo[1,2-b]isoquinolines as antileishmanial agents

The development of new molecules for the treatment of leishmaniasis is, a neglected parasitic disease, is urgent as current anti-leishmanial therapeutics are hampered by drug toxicity and resistance. The pyrrolo[1,2-b]isoquinoline core was selected as starting point, and palladium-catalyzed Heck-initiated cascade reactions were developed for the synthesis of a series of C-10 substituted derivatives. Their in vitro leishmanicidal activity against visceral (L. donovani) and cutaneous (L. amazonensis) leishmaniasis was evaluated. The best activity was found, in general, for the 10-arylmethyl substituted pyrroloisoquinolines. In particular, 2ad (IC₅₀ = 3.30 μM, SI > 77.01) and 2bb (IC₅₀ = 3.93 μM, SI > 58.77) were approximately 10-fold more potent and selective than the drug of reference (miltefosine), against L. amazonensis on in vitro promastigote assays, while 2ae was the more active compound in the in vitro amastigote assays (IC₅₀ = 33.59 μM, SI > 8.93). Notably, almost all compounds showed low cytotoxicity, CC₅₀ > 100 μg/mL in J774 cells, highest tested dose. In addition, we have developed the first Perturbation Theory Machine Learning (PTML) algorithm able to predict simultaneously multiple biological activity parameters (IC₅₀, K_i, etc.) vs. any Leishmania species and target protein, with high values of specificity (>98%) and sensitivity (>90%) in both training and validation series. Therefore, this model may be useful to reduce time and assay costs (material and human resources) in the drug discovery process.

Systematic functional analysis of Leishmania protein kinases identifies regulators of differentiation or survival

Differentiation between distinct stages is fundamental for the life cycle of intracellular protozoan parasites and for transmission between hosts, requiring stringent spatial and temporal regulation. Here, we apply kinome-wide gene deletion and gene tagging in Leishmania mexicana promastigotes to define protein kinases with life cycle transition roles. Whilst 162 are dispensable, 44 protein kinase genes are refractory to deletion in promastigotes and are likely core genes required for parasite replication. Phenotyping of pooled gene deletion mutants using bar-seq and projection pursuit clustering reveal functional phenotypic groups of protein kinases involved in differentiation from metacyclic promastigote to amastigote, growth and survival in macrophages and mice, colonisation of the sand fly and motility. This unbiased interrogation of protein kinase function in Leishmania allows targeted investigation of organelle-associated signalling pathways required for successful intracellular parasitism.

The ultimate fate determinants of drug induced cell-death mechanisms in Trypanosomatids

Chemotherapy constitutes a major part of modern-day therapy for infectious and chronic diseases. A drug is said to be effective if it can inhibit its target, induce stress, and thereby trigger an array of cell death pathways in the form of programmed cell death, autophagy, necrosis, etc. Chemotherapy is the only treatment choice against trypanosomatid diseases like Leishmaniasis, Chagas disease, and sleeping sickness. Anti-trypanosomatid drugs can induce various cell death phenotypes depending upon the drug dose and growth stage of the parasites. The mechanisms and pathways triggering cell death in Trypanosomatids serve to help identify potential targets for the development of effective anti-trypanosomatids. Studies show that the key proteins involved in cell death of trypanosomatids are metacaspases, Endonuclease G, Apoptosis-Inducing Factor, cysteine proteases, serine proteases, antioxidant systems, etc. Unlike higher eukaryotes, these organisms either lack the complete set of effectors involved in cell death pathways, or are yet to be deciphered. A detailed summary of the existing knowledge of different drug-induced cell death pathways would help identify the lacuna in each of these pathways and therefore open new avenues for research and thereby new therapeutic targets to explore. The cell death pathway associated complexities in metazoans are absent in trypanosomatids; hence this summary can also help understand the trigger points as well as cross-talk between these pathways. Here we provide an in-depth overview of the existing knowledge of these drug-induced trypanosomatid cell death pathways, describe their associated physiological changes, and suggest potential interconnections amongst them.

Trypanothione Metabolism as Drug Target for Trypanosomatids

Chagas Disease, African sleeping sickness, and leishmaniasis are neglected diseases caused by pathogenic trypanosomatid parasites, which have a considerable impact on morbidity and mortality in poor countries. The available drugs used as treatment have high toxicity, limited access, and can cause parasite drug resistance. Long-term treatments, added to their high toxicity, result in patients that give up therapy. Trypanosomatids presents a unique trypanothione based redox system, which is responsible for maintaining the redox balance. Therefore, inhibition of these essential and exclusive parasite's metabolic pathways, absent from the mammalian host, could lead to the development of more efficient and safe drugs. The system contains different redox cascades, where trypanothione and tryparedoxins play together a central role in transferring reduced power to different enzymes, such as 2-Cys peroxiredoxins, non-selenium glutathione peroxidases, ascorbate peroxidases, glutaredoxins and methionine sulfoxide reductases, through NADPH as a source of electrons. There is sufficient evidence that this complex system is essential for parasite survival and infection. In this review, we explore what is known in terms of essentiality, kinetic and structural data, and the development of inhibitors of enzymes from this trypanothione-based redox system. The recent advances and limitations in the development of lead inhibitory compounds targeting these enzymes have been discussed. The combination of molecular biology, bioinformatics, genomics, and structural biology is fundamental since the knowledge of unique features of the trypanothione-dependent system will provide tools for rational drug design in order to develop better treatments for these diseases.

Effect of DODAB Nano-Sized Cationic Bilayer Fragments against Leishmania amazonensis

The dioctadecyldimethylammonium bromide (DODAB) is a double-chained cationic lipid with potent bactericide and fungistatic activities; however, its toxicity on protozoan parasites is still unknown. Here, we show the antileishmanial activity of DODAB nano-sized cationic bilayer fragments on stationary-phase promastigotes and amastigotes of Leishmania amazonensis, the causative agent of cutaneous leishmaniasis. Upon treatment with DODAB, we analyzed the parasite surface zeta-potential, parasite viability, cellular structural modifications, and intracellular proliferation. The DODAB cytotoxic effect was dose-dependent, with a median effective concentration (EC₅₀) of 25 µM for both life-cycle stages, comparable to the reported data for bacteria and fungi. The treatment with DODAB changed the membrane zeta-potential from negative to positive, compromised the parasite's morphology, affected the cell size regulation, caused a loss of intracellular organelles, and probably dysregulated the plasma membrane permeability without membrane disruption. Moreover, the parasites that survived after treatment induced small parasitophorous vacuoles and failed to proliferate inside macrophages. In conclusion, DODAB displayed antileishmanial activity, and it remains to be elucidated how DODAB acts on the protozoan membrane. Understanding this mechanism can provide insights into the development of new parasite-control strategies.

Essential gene prediction using limited gene essentiality information-An integrative semi-supervised machine learning strategy

Essential gene prediction helps to find minimal genes indispensable for the survival of any organism. Machine learning (ML) algorithms have been useful for the prediction of gene essentiality. However, currently available ML pipelines perform poorly for organisms with limited experimental data. The objective is the development of a new ML pipeline to help in the annotation of essential genes of less explored disease-causing organisms for which minimal experimental data is available. The proposed strategy combines unsupervised feature selection technique, dimension reduction using the Kamada-Kawai algorithm, and semi-supervised ML algorithm employing Laplacian Support Vector Machine (LapSVM) for prediction of essential and non-essential genes from genome-scale metabolic networks using very limited labeled dataset. A novel scoring technique, Semi-Supervised Model Selection Score, equivalent to area under the ROC curve (auROC), has been proposed for the selection of the best model when supervised performance metrics calculation is difficult due to lack of data. The unsupervised feature selection followed by dimension reduction helped to observe a distinct circular pattern in the clustering of essential and non-essential genes. LapSVM then created a curve that dissected this circle for the classification and prediction of essential genes with high accuracy (auROC > 0.85) even with 1% labeled data for model training. After successful validation of this ML pipeline on both Eukaryotes and Prokaryotes that show high accuracy even when the labeled dataset is very limited, this strategy is used for the prediction of essential genes of organisms with inadequate experimentally known data, such as Leishmania sp. Using a graph-based semi-supervised machine learning scheme, a novel integrative approach has been proposed for essential gene prediction that shows universality in application to both Prokaryotes and Eukaryotes with limited labeled data. The essential genes predicted using the pipeline provide an important lead for the prediction of gene essentiality and identification of novel therapeutic targets for antibiotic and vaccine development against disease-causing parasites.

Emerging Functions of Actins and Actin Binding Proteins in Trypanosomatids

Actin is the major protein constituent of the cytoskeleton that performs wide range of cellular functions. It exists in monomeric and filamentous forms, dynamics of which is regulated by a large repertoire of actin binding proteins. However, not much was known about existence of these proteins in trypanosomatids, till the genome sequence data of three important organisms of this class, viz. Trypanosoma brucei, Trypanosoma cruzi and Leishmania major, became available. Here, we have reviewed most of the findings reported to date on the intracellular distribution, structure and functions of these proteins and based on them, we have hypothesized some of their functions. The major findings are as follows: (1) All the three organisms encode at least a set of ten actin binding proteins (profilin, twinfilin, ADF/cofilin, CAP/srv2, CAPz, coronin, two myosins, two formins) and one isoform of actin, except that T. cruzi encodes for three formins and several myosins along with four actins. (2) Actin 1 and a few actin binding proteins (ADF/cofilin, profilin, twinfilin, coronin and myosin13 in L. donovani; ADF/cofilin, profilin and myosin1 in T. brucei; profilin and myosin-F in T.cruzi) have been identified and characterized. (3) In all the three organisms, actin cytoskeleton has been shown to regulate endocytosis and intracellular trafficking. (4) Leishmania actin1 has been the most characterized protein among trypanosomatid actins. (5) This protein is localized to the cytoplasm as well as in the flagellum, nucleus and kinetoplast, and in vitro, it binds to DNA and displays scDNA relaxing and kDNA nicking activities. (6) The pure protein prefers to form bundles instead of thin filaments, and does not bind DNase1 or phalloidin. (7) Myosin13, myosin1 and myosin-F regulate endocytosis and intracellular trafficking, respectively, in Leishmania, T. brucei and T. cruzi. (8) Actin-dependent myosin13 motor is involved in dynamics and assembly of Leishmania flagellum. (9) Leishmania twinfilin localizes mostly to the nucleolus and coordinates karyokinesis by effecting splindle elongation and DNA synthesis. (10) Leishmania coronin binds and promotes actin filament formation and exists in tetrameric form rather than trimeric form, like other coronins. (11) Trypanosomatid profilins are essential for survival of all the three parasites.

Alkyl and Aryl Derivatives Based on p-Coumaric Acid Modification and Inhibitory Action against Leishmania braziliensis and Plasmodium falciparum

In low-income populations, neglected diseases are the principal cause of mortality. Of these, leishmaniasis and malaria, being parasitic, protozoan infections, affect millions of people worldwide and are creating a public health problem. The present work evaluates the leishmanicidal and antiplasmodial action of a series of twelve p-coumaric acid derivatives. Of the tested derivatives, eight presented antiparasitic activities 1–3, 8–12. The hexyl p-coumarate derivative (9) (4.14 ± 0.55 μg/mL; selectivity index (SI) = 2.72) showed the highest leishmanicidal potency against the Leishmania braziliensis amastigote form. The results of the molecular docking study suggest that this compound inhibits aldehyde dehydrogenase (ALDH), mitogen-activated kinase protein (MPK4), and DNA topoisomerase 2 (TOP2), all of which are key enzymes in the development of Leishmania braziliensis. The data indicate that these enzymes interact via Van der Waals bonds, hydrophobic interactions, and hydrogen bonds with phenolic and aliphatic parts of this same compound. Of the other compounds analyzed, methyl p-coumarate (64.59 ± 2.89 μg/mL; IS = 0.1) demonstrated bioactivity against Plasmodium falciparum. The study reveals that esters presenting a p-coumarate substructure are promising for use in synthesis of derivatives with good antiparasitic profiles.

Bar-seq strategies for the LeishGEdit toolbox

The number of fully sequenced genomes increases steadily but the function of many genes remains unstudied. To accelerate dissection of gene function in Leishmania spp. and other kinetoplastids we previously developed a streamlined pipeline for CRISPR-Cas9 gene editing, which we termed LeishGEdit. To facilitate high-throughput mutant screens we have adapted this pipeline by barcoding mutants with unique 17-nucleotide barcodes, allowing loss-of-function screens in mixed populations. Here we present primer design and analysis tools that facilitate these bar-seq strategies. We have developed a standalone easy-to-use pipeline to design CRISPR primers suitable for the LeishGEdit toolbox for any given genome and have generated a list of 14,995 barcodes. Barcodes and oligo sequences are now accessible through our website www.leishgedit.net allowing researchers to pursue bar-seq experiments in all currently available TriTrypDB genomes (release 41). This will streamline CRISPR bar-seq assays in kinetoplastids, enabling pooled mutant screens across the community.

Defeating the trypanosomatid trio: proteomics of the protozoan parasites causing neglected tropical diseases

Mass spectrometry-based proteomics enables accurate measurement of the modulations of proteins on a large scale upon perturbation and facilitates the understanding of the functional roles of proteins in biological systems. It is a particularly relevant methodology for studying Leishmania spp., Trypanosoma cruzi and Trypanosoma brucei, as the gene expression in these parasites is primarily regulated by posttranscriptional mechanisms. Large-scale proteomics studies have revealed a plethora of information regarding modulated proteins and their molecular interactions during various life processes of the protozoans, including stress adaptation, life cycle changes and interactions with the host. Important molecular processes within the parasite that regulate the activity and subcellular localisation of its proteins, including several co- and post-translational modifications, are also accurately captured by modern proteomics mass spectrometry techniques. Finally, in combination with synthetic chemistry, proteomic techniques facilitate unbiased profiling of targets and off-targets of pharmacologically active compounds in the parasites. This provides important data sets for their mechanism of action studies, thereby aiding drug development programmes.

Essential roles for deubiquitination in Leishmania life cycle progression

The parasitic protozoan Leishmania requires proteasomal, autophagic and lysosomal proteolytic pathways to enact the extensive cellular remodelling that occurs during its life cycle. The proteasome is essential for parasite proliferation, yet little is known about the requirement for ubiquitination/deubiquitination processes in growth and differentiation. Activity-based protein profiling of L. mexicana C12, C19 and C65 deubiquitinating cysteine peptidases (DUBs) revealed DUB activity remains relatively constant during differentiation of procyclic promastigote to amastigote. However, when life cycle phenotyping (bar-seq) was performed on a pool including 15 barcoded DUB null mutants created in promastigotes using CRISPR-Cas9, significant loss of fitness was observed during differentiation and intracellular infection. DUBs 4, 7, and 13 are required for successful transformation from metacyclic promastigote to amastigote and DUBs 3, 5, 6, 8, 10, 11 and 14 are required for normal amastigote proliferation in mice. DUBs 1, 2, 12 and 16 are essential for promastigote viability and the essential role of DUB2 in establishing infection was demonstrated using DiCre inducible gene deletion in vitro and in vivo. DUB2 is found in the nucleus and interacts with nuclear proteins associated with transcription/chromatin dynamics, mRNA splicing and mRNA capping. DUB2 has broad linkage specificity, cleaving all the di-ubiquitin chains except for Lys27 and Met1. Our study demonstrates the crucial role that DUBs play in differentiation and intracellular survival of Leishmania and that amastigotes are exquisitely sensitive to disruption of ubiquitination homeostasis.

The complexity and diversity of the actin cytoskeleton of trypanosomatids

Trypanosomatids are a monophyletic group of parasitic flagellated protists belonging to the order Kinetoplastida. Their cytoskeleton is primarily made up of microtubules in which no actin microfilaments have been detected. Although all these parasites contain actin, it is widely thought that their actin cytoskeleton is reduced when compared to most eukaryotic organisms. However, there is increasing evidence that it is more complex than previously thought. As in other eukaryotic organisms, trypanosomatids encode for a conventional actin that is expected to form microfilament-like structures, and for members of three conserved actin-related proteins probably involved in microfilament nucleation (ARP2, ARP3) and in gene expression regulation (ARP6). In addition to these canonical proteins, also encode for an expanded set of actins and actin-like proteins that seem to be restricted to kinetoplastids. Analysis of their amino acid sequences demonstrated that, although very diverse in primary sequence when compared to actins of model organisms, modelling of their tertiary structure predicted the presence of the actin fold in all of them. Experimental characterization has been done for only a few of the trypanosomatid actins and actin-binding proteins. The most studied is the conventional actin of Leishmania donovani (LdAct), which unusually requires both ATP and Mg²⁺ for polymerization, unlike other conventional actins that do not require ATP. Additionally, polymerized LdAct tends to assemble in bundles rather than in single filaments. Regulation of actin polymerization depends on their interaction with actin-binding proteins. In trypanosomatids, there is a reduced but sufficient core of actin-binding proteins to promote microfilament nucleation, turnover and stabilization. There are also genes encoding for members of two families of myosin motor proteins, including one lineage-specific. Homologues to all identified actin-family proteins and actin-binding proteins of trypanosomatids are also present in Paratrypanosoma confusum (an early branching trypanosomatid) and in Bodo saltans (a closely related free-living organism belonging to the trypanosomatid sister order of Bodonida) suggesting they were all present in their common ancestor. Secondary losses of these genes may have occurred during speciation within the trypanosomatids, with salivarian trypanosomes having lost many of them and stercorarian trypanosomes retaining most.

The current drug discovery landscape for trypanosomiasis and leishmaniasis: Challenges and strategies to identify drug targets

Human trypanosomiasis and leishmaniasis are vector-borne neglected tropical diseases caused by infection with the protozoan parasites Trypanosoma spp. and Leishmania spp., respectively. Once restricted to endemic areas, these diseases are now distributed worldwide due to human migration, climate change, and anthropogenic disturbance, causing significant health and economic burden globally. The current chemotherapy used to treat these diseases has limited efficacy, and drug resistance is spreading. Hence, new drugs are urgently needed. Phenotypic compound screenings have prevailed as the leading method to discover new drug candidates against these diseases. However, the publication of the complete genome sequences of multiple strains, advances in the application of CRISPR/Cas9 technology, and in vivo bioluminescence-based imaging have set the stage for advancing target-based drug discovery. This review analyses the limitations of the narrow pool of available drugs presently used for treating these diseases. It describes the current drug-based clinical trials highlighting the most promising leads. Furthermore, the review presents a focused discussion on the most important biological and pharmacological challenges that target-based drug discovery programs must overcome to advance drug candidates. Finally, it examines the advantages and limitations of modern research tools designed to identify and validate essential genes as drug targets, including genomic editing applications and in vivo imaging.

Computational approaches for drug discovery against trypanosomatid-caused diseases

During three decades, only about 20 new drugs have been developed for malaria, tuberculosis and all neglected tropical diseases (NTDs). This critical situation was reached because NTDs represent only 10% of health research investments; however, they comprise about 90% of the global disease burden. Computational simulations applied in virtual screening (VS) strategies are very efficient tools to identify pharmacologically active compounds or new indications for drugs already administered for other diseases. One of the advantages of this approach is the low time-consuming and low-budget first stage, which filters for testing experimentally a group of candidate compounds with high chances of binding to the target and present trypanocidal activity. In this work, we review the most common VS strategies that have been used for the identification of new drugs with special emphasis on those applied to trypanosomiasis and leishmaniasis. Computational simulations based on the selected protein targets or their ligands are explained, including the method selection criteria, examples of successful VS campaigns applied to NTDs, a list of validated molecular targets for drug development and repositioned drugs for trypanosomatid-caused diseases. Thereby, here we present the state-of-the-art of VS and drug repurposing to conclude pointing out the future perspectives in the field.

Outwitting an Old Neglected Nemesis: A Review on Leveraging Integrated Data-Driven Approaches to Aid in Unraveling of Leishmanicides of Therapeutic Potential

The global prevalence of leishmaniasis has increased with skyrocketed mortality in the past decade. The causative agent of leishmaniasis is Leishmania species, which infects populations in almost all the continents. Prevailing treatment regimens are consistently inefficient with reported side effects, toxicity and drug resistance. This review complements existing ones by discussing the current state of treatment options, therapeutic bottlenecks including chemoresistance and toxicity, as well as drug targets. It further highlights innovative applications of nanotherapeutics-based formulations, inhibitory potential of leishmanicides, anti-microbial peptides and organometallic compounds on leishmanial species. Moreover, it provides essential insights into recent machine learning-based models that have been used to predict novel leishmanicides and also discusses other new models that could be adopted to develop fast, efficient, robust and novel algorithms to aid in unraveling the next generation of anti-leishmanial drugs. A plethora of enriched functional genomic, proteomic, structural biology, high throughput bioassay and drug-related datasets are currently warehoused in both general and leishmania-specific databases. The warehoused datasets are essential inputs for training and testing algorithms to augment the prediction of biotherapeutic entities. In addition, we demonstrate how pharmacoinformatics techniques including ligand-, structure- and pharmacophore-based virtual screening approaches have been utilized to screen ligand libraries against both modeled and experimentally solved 3D structures of essential drug targets. In the era of data-driven decision-making, we believe that highlighting intricately linked topical issues relevant to leishmanial drug discovery offers a one-stop-shop opportunity to decipher critical literature with the potential to unlock implicit breakthroughs.

Comparative high-throughput analysis of the Trypanosoma cruzi response to organometallic compounds

An in-depth, comparative look at the effects of two structurally related organometallic Pd and Pt compounds on the global gene expression pattern of T. cruzi epimastigotes. This parasite is the causative agent of Chagas disease.

Metabolic Control Analysis for Drug Target Prioritization in Trypanosomatids

To validate therapeutic targets in metabolic pathways of trypanosomatids, the criterion of enzyme essentiality determined by gene knockout or knockdown is usually being applied. Since, it is often found that most of the enzymes/proteins analyzed are essential, additional criteria have to be implemented for drug target prioritization. Metabolic control analysis (MCA), often in conjunction with kinetic pathway modeling, offers such possibility for prioritization. MCA is a theoretical and experimental approach to analyze how metabolic pathways are controlled. It involves strategies to perform quantitative analyses to determine the degree in which an enzyme controls a pathway flux, a value called flux control coefficient ([Formula: see text]). By determining the [Formula: see text] of individual steps in a metabolic pathway, the distribution of control of the pathway is established, that is, the identification of the main flux-controlling steps. Therefore, MCA can help in ranking pathway enzymes as drug targets from a metabolic perspective. In this chapter, three approaches to determine [Formula: see text] are reviewed: (1) In vitro pathway reconstitution, (2) manipulation of enzyme activities within parasites, and (3) in silico kinetic modeling of the metabolic pathway. To perform these methods, accurate experimental data of enzyme activities, metabolite concentrations and pathway fluxes are necessary. The methodology is illustrated with the example of trypanothione metabolism of Trypanosoma cruzi and protocols to determine such experimental data for this metabolic process are also described. However, the MCA strategy can be applied to any metabolic pathway in the parasite and general directions to perform it are provided in this chapter.

State-of-the-art CRISPR/Cas9 Technology for Genome Editing in Trypanosomatids

CRISPR/Cas9 technology has revolutionized biology. This prokaryotic defense system against foreign DNA has been repurposed for genome editing in a broad range of cell tissues and organisms. Trypanosomatids are flagellated protozoa belonging to the order Kinetoplastida. Some of its most representative members cause important human diseases affecting millions of people worldwide, such as Chagas disease, sleeping sickness and different forms of leishmaniases. Trypanosomatid infections represent an enormous burden for public health and there are no effective treatments for most of the diseases they cause. Since the emergence of the CRISPR/Cas9 technology, the genetic manipulation of these parasites has notably improved. As a consequence, genome editing is now playing a key role in the functional study of proteins, in the characterization of metabolic pathways, in the validation of alternative targets for antiparasitic interventions, and in the study of parasite biology and pathogenesis. In this work we review the different strategies that have been used to adapt the CRISPR/Cas9 system to Trypanosoma cruzi, Trypanosoma brucei, and Leishmania spp., as well as the research progress achieved using these approaches. Thereby, we will present the state-of-the-art molecular tools available for genome editing in trypanosomatids to finally point out the future perspectives in the field.

Genome and transcriptome analyses of Leishmania spp.: opening Pandora's box

In the last 30 years, significant advances in genetic manipulation tools along with complete genome and transcriptome sequencing have advanced our understanding of the biology of Leishmania parasites and their interplay with the sand fly and mammalian hosts. High-throughput sequencing in association with CRISPR/Cas9 have prepared the ground for significant advances. Given the richness of the progress made over the last decade, in this article, we focused on the most recent contributions of genome-wide and transcriptome analyses of Leishmania spp., which permit the comparison of life cycle stages, the evaluation of different strains and species in their natural niches and in the field and the simultaneously comparison of the gene expression profiles of parasites and hosts.