A new approach for potential drug target discovery through in silico metabolic pathway analysis using Trypanosoma cruzi genome information

Extensive Translational Regulation through the Proliferative Transition of Trypanosoma cruzi Revealed by Multi-Omics

Trypanosoma cruzi is the etiological agent for Chagas disease, a neglected parasitic disease in Latin America. Gene transcription control governs the eukaryotic cell replication but is absent in trypanosomatids; thus, it must be replaced by posttranscriptional regulatory events. We investigated the entrance into the T. cruzi replicative cycle using ribosome profiling and proteomics on G₁/S epimastigote cultures synchronized with hydroxyurea. We identified 1,784 translationally regulated genes (change > 2, false-discovery rate [FDR] < 0.05) and 653 differentially expressed proteins (change > 1.5, FDR < 0.05), respectively. A major translational remodeling accompanied by an extensive proteome change is found, while the transcriptome remains largely unperturbed at the replicative entrance of the cell cycle. The differentially expressed genes comprise specific cell cycle processes, confirming previous findings while revealing candidate cell cycle regulators that undergo previously unnoticed translational regulation. Clusters of genes showing a coordinated regulation at translation and protein abundance share related biological functions such as cytoskeleton organization and mitochondrial metabolism; thus, they may represent posttranscriptional regulons. The translatome and proteome of the coregulated clusters change in both coupled and uncoupled directions, suggesting that complex cross talk between the two processes is required to achieve adequate protein levels of different regulons. This is the first simultaneous assessment of the transcriptome, translatome, and proteome of trypanosomatids, which represent a paradigm for the absence of transcriptional control. The findings suggest that gene expression chronology along the T. cruzi cell cycle is controlled mainly by translatome and proteome changes coordinated using different mechanisms for specific gene groups. IMPORTANCE Trypanosoma cruzi is an ancient eukaryotic unicellular parasite causing Chagas disease, a potentially life-threatening illness that affects 6 to 7 million people, mostly in Latin America. The antiparasitic treatments for the disease have incomplete efficacy and adverse reactions; thus, improved drugs are needed. We study the mechanisms governing the replication of the parasite, aiming to find differences with the human host, valuable for the development of parasite-specific antiproliferative drugs. Transcriptional regulation is essential for replication in most eukaryotes, but in trypanosomatids, it must be replaced by subsequent gene regulation steps since they lack transcription initiation control. We identified the genome-wide remodeling of mRNA translation and protein abundance during the entrance to the replicative phase of the cell cycle. We found that translation is strongly regulated, causing variation in protein levels of specific cell cycle processes, representing the first simultaneous study of the translatome and proteome in trypanosomatids.

A novel receptor for platelet-activating factor and lysophosphatidylcholine in Trypanosoma cruzi

The lipid mediators, platelet-activating factor (PAF) and lysophosphatidylcholine (LPC), play relevant pathophysiological roles in Trypanosoma cruzi infection. Several species of LPC, including C18:1 LPC, which mimics the effects of PAF, are synthesized by T. cruzi. The present study identified a receptor in T. cruzi, which was predicted to bind to PAF, and found it to be homologous to members of the progestin and adiponectin family of receptors (PAQRs). We constructed a three-dimensional model of the T. cruzi PAQR (TcPAQR) and performed molecular docking to predict the interactions of the TcPAQR model with C16:0 PAF and C18:1 LPC. We knocked out T. cruzi PAQR (TcPAQR) gene and confirmed the identity of the expressed protein through immunoblotting and immunofluorescence assays using an anti-human PAQR antibody. Wild-type and knockout (KO) parasites were also used to investigate the in vitro cell differentiation and interactions with peritoneal mouse macrophages; TcPAQR KO parasites were unable to react to C16:0 PAF or C18:1 LPC. Our data are highly suggestive that PAF and LPC act through TcPAQR in T. cruzi, triggering its cellular differentiation and ability to infect macrophages.

Systematic in silico Evaluation of Leishmania spp. Proteomes for Drug Discovery

Leishmaniasis is a group of neglected infectious diseases, with approximately 1. 3 million new cases each year, for which the available therapies have serious limitations. Therefore, it is extremely important to apply efficient and low-cost methods capable of selecting the best therapeutic targets to speed up the development of new therapies against those diseases. Thus, we propose the use of integrated computational methods capable of evaluating the druggability of the predicted proteomes of Leishmania braziliensis and Leishmania infantum, species responsible for the different clinical manifestations of leishmaniasis in Brazil. The protein members of those proteomes were assessed based on their structural, chemical, and functional contexts applying methods that integrate data on molecular function, biological processes, subcellular localization, drug binding sites, druggability, and gene expression. These data were compared to those extracted from already known drug targets (BindingDB targets), which made it possible to evaluate Leishmania proteomes for their biological relevance and treatability. Through this methodology, we identified more than 100 proteins of each Leishmania species with druggability characteristics, and potential interaction with available drugs. Among those, 31 and 37 proteins of L. braziliensis and L. infantum, respectively, have never been tested as drug targets, and they have shown evidence of gene expression in the evolutionary stage of pharmacological interest. Also, some of those Leishmania targets showed an alignment similarity of <50% when compared to the human proteome, making these proteins pharmacologically attractive, as they present a reduced risk of side effects. The methodology used in this study also allowed the evaluation of opportunities for the repurposing of compounds as anti-leishmaniasis drugs, inferring potential interaction between Leishmania proteins and ~1,000 compounds, of which only 15 have already been tested as a treatment for leishmaniasis. Besides, a list of potential Leishmania targets to be tested using drugs described at BindingDB, such as the potential interaction of the DEAD box RNA helicase, TRYR, and PEPCK proteins with the Staurosporine compound, was made available to the public.

NMR structure and dynamics of Q4DY78, a conserved kinetoplasid-specific protein from Trypanosoma cruzi

The 106-residue protein Q4DY78 (UniProt accession number) from Trypanosoma cruzi is highly conserved in the related kinetoplastid pathogens Trypanosoma brucei and Leishmania major. Given the essentiality of its orthologue in T. brucei, the high sequence conservation with other trypanosomatid proteins, and the low sequence similarity with mammalian proteins, Q4DY78 is an attractive protein for structural characterization. Here, we solved the structure of Q4DY78 by solution NMR and evaluated its backbone dynamics. Q4DY78 is composed of five α -helices and a small, two-stranded antiparallel β-sheet. The backbone RMSD is 0.22 ± 0.05 Å for the representative ensemble of the 20 lowest-energy structures. Q4DY78 is overall rigid, except for N-terminal residues (V⁸ to I¹⁰), residues at loop 4 (K⁵⁷ to G⁶⁵) and residues at the C-terminus (F⁸⁹ to F¹¹²). Q4DY78 has a short motif FPCAP that could potentially mediate interactions with the host cytoskeleton via interaction with EVH1 (Drosophila Enabled (Ena)/Vasodilator-stimulated phosphoprotein (VASP) homology 1) domains. Albeit Q4DY78 lacks calcium-binding motifs, its fold resembles that of eukaryotic calcium-binding proteins such as calcitracin, calmodulin, and polcacin Bet V4. We characterized this novel protein with a calcium binding fold without the capacity to bind calcium.

Of Drugs and Trypanosomatids: New Tools and Knowledge to Reduce Bottlenecks in Drug Discovery

Leishmaniasis (Leishmania species), sleeping sickness (Trypanosoma brucei), and Chagas disease (Trypanosoma cruzi) are devastating and globally spread diseases caused by trypanosomatid parasites. At present, drugs for treating trypanosomatid diseases are far from ideal due to host toxicity, elevated cost, limited access, and increasing rates of drug resistance. Technological advances in parasitology, chemistry, and genomics have unlocked new possibilities for novel drug concepts and compound screening technologies that were previously inaccessible. In this perspective, we discuss current models used in drug-discovery cascades targeting trypanosomatids (from in vitro to in vivo approaches), their use and limitations in a biological context, as well as different examples of recently discovered lead compounds.

Insights about resveratrol analogs against trypanothione reductase of Leishmania braziliensis: Molecular modeling, computational docking and in vitro antileishmanial studies

In this work, we combined molecular modeling, computational docking and in vitro analysis to explore the antileishmanial effect of some resveratrol analogs (ResAn), focusing on their pro-oxidant effect. The molecular target was the trypanothione reductase of Leishmania braziliensis (LbTryR), an essential component of the antioxidant defenses in trypanosomatid parasites. Three-dimensional structures of LbTryR were modeled and molecular docking studies of ResAn1-5 compounds showed the following affinity: ResAn1 > ResAn2 > ResAn4 > ResAn5 > ResAn3. Positive correlation was observed between these compounds' affinity to the LbTryR and the IC₅₀ values against Leishmania sp (ResAn1 < ResAn2 < ResAn4), which allows for TryR being considered an important target for them. As the compound ResAn1 showed the best antileishmanial activity, and docking studies showed its high affinity for NADP binding site (NS) of TryR, plus having been able to induce ROS production in L. braziliensis promastigotes treated, ResAn1 probably occupies NS interfering in the electron transfer processes responsible for the catalytic reaction. The in silico prediction of ADMET properties suggests that ResAn1 may be a promising drug candidate with properties to cross biological membranes and high gastrointestinal absorption, not violating Lipinski's rules. Ultimately, the antileishmanial effect of ResAn can be associated with a pro-oxidant effect which, in turn, can be exploited as an antimicrobial agent. Communicated by Ramaswamy H. Sarma.

Mining of potential drug targets through the identification of essential and analogous enzymes in the genomes of pathogens of Glycine max, Zea mays and Solanum lycopersicum

Pesticides are one of the most widely used pest and disease control measures in plant crops and their indiscriminate use poses a direct risk to the health of populations and environment around the world. As a result, there is a great need for the development of new, less toxic molecules to be employed against plant pathogens. In this work, we employed an in silico approach to study the genes coding for enzymes of the genomes of three commercially important plants, soybean (Glycine max), tomato (Solanum lycopersicum) and corn (Zea mays), as well as 15 plant pathogens (4 bacteria and 11 fungi), focusing on revealing a set of essential and non-homologous isofunctional enzymes (NISEs) that could be prioritized as drug targets. By combining sequence and structural data, we obtained an initial set of 568 cases of analogy, of which 97 were validated and further refined, revealing a subset of 29 essential enzymatic activities with a total of 119 different structural forms, most belonging to central metabolic routes, including the carbohydrate metabolism, the metabolism of amino acids, among others. Further, another subset of 26 enzymatic activities possess a tertiary structure specific for the pathogen, not present in plants, men and Apis mellifera, which may be of importance for the development of specific enzymatic inhibitors against plant diseases that are less harmful to humans and the environment.

In silico identification of inhibitors of ribose 5-phosphate isomerase from Trypanosoma cruzi using ligand and structure based approaches

Chagas disease, caused by the protozoan Trypanosoma cruzi, affects approximately seven million people, mainly in Latin America, and causes about 7000 deaths annually. The available treatments are unsatisfactory and search for more effective drugs against this pathogen is critical. In this context, the ribose 5-phosphate isomerase (Rpi) enzyme is a potential drug target mainly due to its function in the pentose phosphate pathway and its essentiality (previously shown in other trypanosomatids). In this study, we propose novel potential inhibitors for the Rpi of T. cruzi (TcRpi) based on a computer-aided approach, including structure-based and ligand-based pharmacophore modeling. Along with a substructural and similarity search, the selected pharmacophore hypotheses were used to screen the purchasable subset of the ZINC Database, yielding 20,183 candidate compounds. These compounds were submitted to molecular docking studies in the TcRpi and Human Rpi (HsRpi) active sites in order to identify potential selective inhibitors for the T. cruzi enzyme. After the molecular docking and ADME-T (absorption, distribution, metabolism, excretion and toxicity)/PAINS (pan-assay interference compounds) screenings, 211 molecules were selected as potential TcRpi inhibitors. Out of these, three compounds - ZINC36975961, ZINC63480117, and ZINC43763931 - were submitted to molecular dynamics simulations and two of them - ZINC36975961 and ZINC43763931- had good performance and made interactions with important active site residues over all the simulation time. These compounds could be considered potential TcRpi inhibitors candidates and also may be used as leads for developing new TcRpi inhibitors.

A Computational Methodology to Overcome the Challenges Associated With the Search for Specific Enzyme Targets to Develop Drugs Against Leishmania major

We present an approach for detecting enzymes that are specific of Leishmania major compared with Homo sapiens and provide targets that may assist research in drug development. This approach is based on traditional techniques of sequence homology comparison by similarity search and Markov modeling; it integrates the characterization of enzymatic functionality, secondary and tertiary protein structures, protein domain architecture, and metabolic environment. From 67 enzymes represented by 42 enzymatic activities classified by AnEnPi (Analogous Enzymes Pipeline) as specific for L major compared with H sapiens, only 40 (23 Enzyme Commission [EC] numbers) could actually be considered as strictly specific of L major and 27 enzymes (19 EC numbers) were disregarded for having ambiguous homologies or analogies with H sapiens. Among the 40 strictly specific enzymes, we identified sterol 24-C-methyltransferase, pyruvate phosphate dikinase, trypanothione synthetase, and RNA-editing ligase as 4 essential enzymes for L major that may serve as targets for drug development.

The Comparative Genomics and Phylogenomics of Leishmania amazonensis Parasite

Leishmaniasis is an infectious disease caused by Leishmania species. Leishmania amazonensis is a New World Leishmania species belonging to the Mexicana complex, which is able to cause all types of leishmaniasis infections. The L. amazonensis reference strain MHOM/BR/1973/M2269 was sequenced identifying 8,802 codifying sequences (CDS), most of them of hypothetical function. Comparative analysis using six Leishmania species showed a core set of 7,016 orthologs. L. amazonensis and Leishmania mexicana share the largest number of distinct orthologs, while Leishmania braziliensis presented the largest number of inparalogs. Additionally, phylogenomic analysis confirmed the taxonomic position for L. amazonensis within the “Mexicana complex”, reinforcing understanding of the split of New and Old World Leishmania. Potential non-homologous isofunctional enzymes (NISE) were identified between L. amazonensis and Homo sapiens that could provide new drug targets for development.

In vivo imaging of mice infected with bioluminescent Trypanosoma cruzi unveils novel sites of infection

Background

The development of techniques that allow the imaging of animals infected with parasites expressing luciferase opens up new possibilities for following the fate of parasites in infected mammals.

Methods

D-luciferin potassium salt stock solution was prepared in phosphate-buffered saline (PBS) at 15 mg/ml. To produce bioluminescence, infected and control mice received an intraperitoneal injection of luciferin stock solution (150 mg/kg). All mice were immediately anesthetized with 2% isofluorane, and after 10 minutes were imaged. Ex vivo evaluation of infected tissues and organs was evaluated in a 24-well plate in 150 μg/ml D-luciferin diluted in PBS. Images were captured using the IVIS Lumina image system (Xenogen). Dissected organs were also evaluated by microscopy of hematoxylin-eosin stained sections.

Results

Here we describe the results obtained using a genetically modified Dm28c strain of T. cruzi expressing the firefly luciferase to keep track of infection by bioluminescence imaging. Progression of infection was observed in vivo in BALB/c mice at various intervals after infection with transgenic Dm28c-luc. The bioluminescent signal was immediately observed at the site of T. cruzi inoculation, and one day post infection (dpi) it was disseminated in the peritoneal cavity. A similar pattern in the cavity was observed on 7 dpi, but the bioluminescence was more intense in the terminal region of the large intestine, rectum, and gonads. On 14 and 21 dpi, bioluminescent parasites were also observed in the heart, snout, paws, hind limbs, and forelimbs. From 28 dpi to 180 dpi in chronically infected mice, bioluminescence declined in regions of the body but was concentrated in the gonad region. Ex vivo evaluation of dissected organs and tissues by bioluminescent imaging confirmed the in vivo bioluminescent foci. Histopathological analysis of dissected organs demonstrated parasite nests at the rectum and snout, in muscle fibers of mice infected with Dm28c-WT and with Dm28c-luc, corroborating the bioluminescent imaging.

Conclusion

Bioluminescence imaging is accurate for tracking parasites in vivo, and this methodology is important to gain a better understanding of the infection, tissue inflammation, and parasite biology regarding host cell interaction, proliferation, and parasite clearance to subpatent levels.

Selection of molecular targets for drug development against trypanosomatids

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

Metabolic signatures of triatomine vectors of Trypanosoma cruzi unveiled by metabolomics

Chagas disease is a trypanosomiasis whose causative agent is the protozoan parasite Trypanosoma cruzi, which is transmitted to humans by hematophagous insects known as triatomines and affects a large proportion of South America. The digestive tract of the insect vectors in which T. cruzi develops constitutes a dynamic environment that affects the development of the parasite. Thus, we set out to investigate the chemical composition of the triatomine intestinal tract through a metabolomics approach. We performed Direct Infusion Fourier Transform Ion Cyclotron Resonance Mass Spectrometry on fecal samples of three triatomine species (Rhodnius prolixus, Triatoma infestans, Panstrongylus megistus) fed with rabbit blood. We then identified groups of metabolites whose frequencies were either uniform in all species or enriched in each of them. By querying the Human Metabolome Database, we obtained putative identities of the metabolites of interest. We found that a core group of metabolites with uniform frequencies in all species represented approximately 80% of the molecules detected, whereas the other 20% varied among triatomine species. The uniform core was composed of metabolites of various categories, including fatty acids, steroids, glycerolipids, nucleotides, sugars, and others. Nevertheless, the metabolic fingerprint of triatomine feces differs depending on the species considered. The variable core was mainly composed of prenol lipids, amino acids, glycerolipids, steroids, phenols, fatty acids and derivatives, benzoic acid and derivatives, flavonoids, glycerophospholipids, benzopyrans, and quinolines. Triatomine feces constitute a rich and varied chemical medium whose constituents are likely to affect T. cruzi development and infectivity. The complexity of the fecal metabolome of triatomines suggests that it may affect triatomine vector competence for specific T. cruzi strains. Knowledge of the chemical environment of T. cruzi in its invertebrate host is likely to generate new ways to understand the factors influencing parasite proliferation as well as methods to control Chagas disease.

Specific and nonhomologous isofunctional enzymes of the genetic information processing pathways as potential therapeutical targets for tritryps

Leishmania major, Trypanosoma brucei, and Trypanosoma cruzi (Tritryps) are unicellular protozoa that cause leishmaniasis, sleeping sickness and Chagas' disease, respectively. Most drugs against them were discovered through the screening of large numbers of compounds against whole parasites. Nonhomologous isofunctional enzymes (NISEs) may present good opportunities for the identification of new putative drug targets because, though sharing the same enzymatic activity, they possess different three-dimensional structures thus allowing the development of molecules against one or other isoform. From public data of the Tritryps' genomes, we reconstructed the Genetic Information Processing Pathways (GIPPs). We then used AnEnPi to look for the presence of these enzymes between Homo sapiens and Tritryps, as well as specific enzymes of the parasites. We identified three candidates (ECs 3.1.11.2 and 6.1.1.-) in these pathways that may be further studied as new therapeutic targets for drug development against these parasites.

Structural modelling and comparative analysis of homologous, analogous and specific proteins from Trypanosoma cruzi versus Homo sapiens: putative drug targets for chagas' disease treatment

Background

Trypanosoma cruzi is the etiological agent of Chagas' disease, an endemic infection that causes thousands of deaths every year in Latin America. Therapeutic options remain inefficient, demanding the search for new drugs and/or new molecular targets. Such efforts can focus on proteins that are specific to the parasite, but analogous enzymes and enzymes with a three-dimensional (3D) structure sufficiently different from the corresponding host proteins may represent equally interesting targets. In order to find these targets we used the workflows MHOLline and AnEnΠ obtaining 3D models from homologous, analogous and specific proteins of Trypanosoma cruzi versus Homo sapiens.

Results

We applied genome wide comparative modelling techniques to obtain 3D models for 3,286 predicted proteins of T. cruzi. In combination with comparative genome analysis to Homo sapiens, we were able to identify a subset of 397 enzyme sequences, of which 356 are homologous, 3 analogous and 38 specific to the parasite.

Conclusions

In this work, we present a set of 397 enzyme models of T. cruzi that can constitute potential structure-based drug targets to be investigated for the development of new strategies to fight Chagas' disease. The strategies presented here support the concept of structural analysis in conjunction with protein functional analysis as an interesting computational methodology to detect potential targets for structure-based rational drug design. For example, 2,4-dienoyl-CoA reductase (EC 1.3.1.34) and triacylglycerol lipase (EC 3.1.1.3), classified as analogous proteins in relation to H. sapiens enzymes, were identified as new potential molecular targets.