Evolutionary aspects of trypanosomes: analysis of genes

Recent Advances in the Development of Triose Phosphate Isomerase Inhibitors as Antiprotozoal Agents

BACKGROUND

Parasitic diseases caused by protozoa such as Chagas disease, leishmaniasis, malaria, African trypanosomiasis, amebiasis, trichomoniasis, and giardiasis are considered serious public health problems in developing countries. Drug-resistance among parasites justifies the search for new therapeutic drugs and the identification of new targets becomes a valuable approach. In this scenario, glycolysis pathway which consists of the conversion of glucose into pyruvate plays an important role in the protozoa energy supply and it is therefore considered as a promising target. In this pathway, triose phosphate isomerase (TIM) plays an essential role in efficient energy production. Furthermore, protozoa TIM show structural differences with human enzyme counterparts suggesting the possibility of obtaining selective inhibitors. Therefore, TIM is considered a valid approach to develop new antiprotozoal agents, inhibiting the glycolysis in the parasite.

OBJECTIVE

In this review, we discuss the drug design strategies, structure-activity relationship, and binding modes of outstanding TIM inhibitors against Trypanosoma cruzi, Trypanosoma brucei, Plasmodium falciparum, Giardia lamblia, Leishmania mexicana, Trichomonas vaginalis, and Entamoeba histolytica.

RESULTS

TIM inhibitors showed mainly aromatic systems and symmetrical structure, where the size and type of heteroatom are important for enzyme inhibition. This inhibition is mainly based on the interaction with i) the interfacial region of TIM inducing changes on the quaternary and tertiary structure or ii) with the TIM catalytic region were the main pathways that disabled the catalytic activity of the enzyme.

CONCLUSION

Benzothiazole, benzoxazole, benzimidazole, and sulfhydryl derivatives stand out as TIM inhibitors. In silico and in vitro studies demonstrate that the inhibitors bind mainly at the TIM dimer interface. In this review, the development of new TIM inhibitors as antiprotozoal drugs is demonstrated as an important pharmaceutical strategy that may lead to new therapies for these ancient parasitic diseases.

Transcriptome analysis of the bloodstream stage from the parasite Trypanosoma vivax

Background

Trypanosoma vivax is the earliest branching African trypanosome. This crucial phylogenetic position makes T. vivax a fascinating model to tackle fundamental questions concerning the origin and evolution of several features that characterize African trypanosomes, such as the Variant Surface Glycoproteins (VSGs) upon which antibody clearing and antigenic variation are based. Other features like gene content and trans-splicing patterns are worth analyzing in this species for comparative purposes.

Results

We present a RNA-seq analysis of the bloodstream stage of T. vivax from data obtained using two complementary sequencing technologies (454 Titanium and Illumina).

Assembly of 454 reads yielded 13385 contigs corresponding to proteins coding genes (7800 of which were identified). These sequences, their annotation and other features are available through an online database presented herein. Among these sequences, about 1000 were found to be species specific and 50 exclusive of the T. vivax strain analyzed here. Expression patterns and levels were determined for VSGs and the remaining genes. Interestingly, VSG expression level, although being high, is considerably lower than in Trypanosoma brucei. Indeed, the comparison of surface protein composition between both African trypanosomes (as inferred from RNA-seq data), shows that they are substantially different, being VSG absolutely predominant in T. brucei, while in T. vivax it represents only about 55%. This raises the question concerning the protective role of VSGs in T. vivax, hence their ancestral role in immune evasion.

It was also found that around 600 genes have their unique (or main) trans-splice site very close (sometimes immediately before) the start codon. Gene Ontology analysis shows that this group is enriched in proteins related to the translation machinery (e.g. ribosomal proteins, elongation factors).

Conclusions

This is the first RNA-seq data study in trypanosomes outside the model species T. brucei, hence it provides the possibility to conduct comparisons that allow drawing evolutionary and functional inferences. This analysis also provides several insights on the expression patterns and levels of protein coding sequences (such as VSG gene expression), trans-splicing, codon patterns and regulatory mechanisms. An online T. vivax RNA-seq database described herein could be a useful tool for parasitologists working with trypanosomes.

Triosephosphate isomerase: a highly evolved biocatalyst

Triosephosphate isomerase (TIM) is a perfectly evolved enzyme which very fast interconverts dihydroxyacetone phosphate and D: -glyceraldehyde-3-phosphate. Its catalytic site is at the dimer interface, but the four catalytic residues, Asn11, Lys13, His95 and Glu167, are from the same subunit. Glu167 is the catalytic base. An important feature of the TIM active site is the concerted closure of loop-6 and loop-7 on ligand binding, shielding the catalytic site from bulk solvent. The buried active site stabilises the enediolate intermediate. The catalytic residue Glu167 is at the beginning of loop-6. On closure of loop-6, the Glu167 carboxylate moiety moves approximately 2 Å to the substrate. The dynamic properties of the Glu167 side chain in the enzyme substrate complex are a key feature of the proton shuttling mechanism. Two proton shuttling mechanisms, the classical and the criss-cross mechanism, are responsible for the interconversion of the substrates of this enolising enzyme.

Codon usage suggests that translational selection has a major impact on protein expression in trypanosomatids

Background

Different proteins are required in widely different quantities to build a living cell. In most organisms, transcription control makes a major contribution to differential expression. This is not the case in trypanosomatids where most genes are transcribed at an equivalent rate within large polycistronic clusters. Thus, trypanosomatids must use post-transcriptional control mechanisms to balance gene expression requirements.

Results

Here, the evidence for translational selection, the enrichment of 'favoured' codons in more highly expressed genes, is explored. A set of highly expressed, tandem-repeated genes display codon bias in Trypanosoma cruzi, Trypanosoma brucei and Leishmania major. The tRNA complement reveals forty-five of the sixty-one possible anticodons indicating widespread use of 'wobble' tRNAs. Consistent with translational selection, cognate tRNA genes for favoured codons are over-represented. Importantly, codon usage (Codon Adaptation Index) correlates with predicted and observed expression level. In addition, relative codon bias is broadly conserved among syntenic genes from different trypanosomatids.

Conclusion

Synonymous codon bias is correlated with tRNA gene copy number and with protein expression level in trypanosomatids. Taken together, the results suggest that translational selection is the dominant mechanism underlying the control of differential protein expression in these organisms. The findings reveal how trypanosomatids may compensate for a paucity of canonical Pol II promoters and subsequent widespread constitutive RNA polymerase II transcription.

Characterisation of the gene encoding type II DNA topoisomerase from Leishmania donovani: a key molecular target in antileishmanial therapy

The gene encoding type II DNA topoisomerase from the kinetoplastid hemoflagellated protozoan parasite Leishmania donovani (LdTOP2) was isolated from a genomic DNA library of this parasite. DNA sequence analysis revealed an ORF of 3711 bp encoding a putative protein of 1236 amino acids with no introns. The deduced amino acid sequence of LdTOP2 showed strong homologies to TOP2 sequences from other kinetoplastids, namely Crithidia and Trypanosoma spp. with estimated identities of 86 and 68%, respectively. LdTOP2 shares a much lower identity of 32% with its human homologue. LdTOP2 is located as a single copy on a chromosome in the 0.7 Mb region in the L.donovani genome and is expressed as a 5 kb transcript. 5'-Mapping studies indicate that the LdTOP2 gene transcript is matured post-transcriptionally with the trans-splicing of the mini-exon occurring at -639 from the predicted initiation site. Antiserum raised in rabbit against glutathione S-transferase fusion protein containing the major catalytic portion of the recombinant L.donovani topoisomerase II protein could detect a band on western blots at approximately 132 kDa, the expected size of the entire protein. Use of the same antiserum for immunolocalisation analysis led to the identification of nuclear, as well as kinetoplast, antigens for L.donovani topoisomerase II. The in vitro biochemical properties of the full-length recombinant LdTOP2 when overexpressed in E.coli were similar to the Mg(II) and ATP-dependent activity found in cell extracts of L.donovani.

Complementation of a glucose transporter mutant of Schizosaccharomyces pombe by a novel Trypanosoma brucei gene

The African trypanosome Trypanosoma brucei has a digenetic life cycle that involves the insect vector and the mammalian host. This is underscored by biochemical switches in its nutritional requirements. In the insect vector, the parasite relies on amino acid catabolism, but in the mammalian host, it derives its energy exclusively from blood glucose. Glucose transport is facilitated, and constitutes the rate-limiting step in ATP synthesis. Here, we report the cloning of a novel glucose transporter-related gene by heterologous screening of a lambdaEMBL4 genomic library of T. brucei EATRO 164 using a rat liver glucose transporter cDNA clone. Genomic analysis shows that the gene is present as a single copy within the parasite genome. The gene encodes a protein with an estimated molecular mass of 55.9 kDa, which shares only segmental homology with members of the glucose transporter superfamily. Several potential post-translational modification sites including phosphorylation, N-glycosylation, and cotranslational myristoylation sites also punctuate the sequence. It is distinguished from classical transporter proteins by the absence of putative hydrophobic membrane-spanning domains. However, this protein was capable of complementing Schizosaccharomyces pombe glucose transporter mutants. The rescued phenotype conferred the ability of the cells to grow on a broad range of sugars, both monosaccharides and disaccharides. The kinetics of glucose uptake reflected those in T. brucei. In addition to complementation in yeast, we also showed that the gene enhanced glucose uptake in cultured mammalian cells.

Compositional similarities between the human immunodeficiency virus and surface antigens of pathogens

The genome of the human immunodeficiency virus (HIV) is rich in A but not U and deficient in C but not G. This asymmetric nucleotide bias is the major factor in determining the unusual composition of HIV proteins. In this report, we have identified the cellular genes in the GenBank database that are compositionally similar to HIV in order to further understand the significance of the nucleotide bias of the viral genome. A total of 101 genes in the bacterial and invertebrate subdivisions of the database were found to have a base composition that is similar to the composition of the HIV genome. The identified cellular sequences represent a discrete subset of the database since 81 of the 101 entries code for antigens from pathogens and nearly all of these organisms infect humans. The amino acid compositions of these surface antigens are also similar to the unusual composition of HIV proteins, which are deficient in proline and rich in lysine and other polar residues encoded by A-rich codons. The similarities between the HIV proteins and the immunodominant antigens from other pathogens may indicate a common pathogenic strategy for the promotion of immune dysregulation.

The nuclear genomes of African and American trypanosomes are strikingly different

We have investigated the compositional distributions of exons and their different codon positions, as well as the codon usage and amino-acid (aa) composition of the nuclear genomes of the African and American trypanosomes Trypanosoma brucei and T. cruzi. Very large differences between the two species were found in all the properties investigated. The most striking differences concern the compositional distributions of third codon positions and the extremely large nucleotide divergence of third codon position for homologous genes encoding proteins that are highly conserved in their aa sequences. Moreover, if coding sequences from each species are divided into two groups according to the GC levels in third codon positions, very different codon usages and aa compositions are found. This indicates a compositional compartmentalization in both genomes which had previously been detected in T. brucei (and T. equiperdum) by compositional fractionation.

Compositional compartmentalization of the nuclear genomes of Trypanosoma brucei and Trypanosoma equiperdum

High molecular weight DNA preparations from Trypanosoma brucei and Trypanosoma equiperdum were fractionated by preparative centrifugation in a Cs2SO4 density gradient in the presence of BAMD, bis(acetatomercurimethyl)dioxane, a sequence-specific DNA ligand. Analytical centrifugation in CsCl of the DNA fractions so obtained showed that both DNAs had a bimodal distribution with two major peaks banding at 1.702-1.703 and 1.708 g/cm3 and representing 1/3 and 2/3 of total DNA, respectively. Several minor components were also detected. These results indicate that a compositional compartmentalization is not only found in the genome of vertebrates and plants, as already described, but also in those of protozoa such as Trypanosomes.

Inhibition of the glycolytic enzymes in the trypanosome: an approach in the development of new leads in the therapy of parasitic diseases

Glycolysis in the trypanosome represents an important target for the development of new therapeutic agents due to the fact that this metabolism is essential for the parasite, glucose being its sole source of energy. In addition, different features of this metabolism and those associated with glycolytic enzymes offer opportunities for the development of efficient and selective compounds. Examples are given in this work of inhibitors directed to the enzymes aldolase and glyceraldehyde-phosphate-dehydrogenase and also of molecules acting specifically on the clusters of basic amino-acids present at the surfaces of the glycolytic enzymes in the parasite.

Isolation and characterization of a gene from Trypanosoma cruzi encoding a 46-kilodalton protein with homology to human and rat tyrosine aminotransferase

The complete sequence of a gene encoding a 46-kDa protein of Trypanosoma cruzi is presented. The first ATG complies with the consensus sequence for initiation of translation. A single band of 2 kb was highlighted by hybridizing a probe from the 46-kDa protein gene to a Northern filter containing total T. cruzi RNA. The gene is present in 50-80 copies per cell and most of them are contained in 2 tandem arrays on large T. cruzi chromosomes (> 2000 kb). A strong homology with rat and human tyrosine aminotransferase was detected. Homology with a Trypanosoma brucei retrotransposon was found in the nonsense strand of the intergenic region.

Molecular analysis of the cytosolic and glycosomal glyceraldehyde-3-phosphate dehydrogenase in Leishmania mexicana

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) activity was detected in two cell compartments of Leishmania mexicana promastigotes. These activities could be attributed to two different isoenzymes, one residing in glycosomes, the other in the cytosol. We have cloned and sequenced the genes for both isoenzymes. The glycosomal enzyme is encoded by two tandemly linked genes of identical sequence and contains features frequently found in glycosomal enzymes: the presence of peptide insertions, a small carboxy-terminal extension with a potential glycosomal targeting signal (-SKM) and an excess of positively charged residues (net charge +7). Only one open reading frame was detected for the cytosolic enzyme. The amino acid sequences of the two proteins are only 55% identical. We discuss some evolutionary aspects of the observed organization of the GAPDH genes in the Trypanosomatidae and the role of the two isoenzymes in the metabolism of these organisms. The possibility to develop GAPDH-specific inhibitors that will be effective against the enzyme of various parasitic members of this family is explored.

Characterization of trypanosome protein phosphatase 1 and 2A catalytic subunits

Oligonucleotides corresponding to highly conserved regions of mammalian protein phosphatase catalytic subunits were used in the polymerase chain reaction (PCR) to generate an amplification product from genomic DNA of Trypanosoma brucei rhodesiense. The PCR product was used to screen a T. b. rhodesiense cDNA library for cDNA clones encoding putative protein phosphatase catalytic subunits. Two cDNA clones, (TPP1A and TPP1B) representing two distinct type 1 catalytic subunit isotypes, encode 39-kDa proteins of 346 amino acids that show 66% and 40% identity, respectively, to mammalian protein phosphatase 1 and 2A catalytic subunits. Both cDNAs are derived from 2.3-kb mRNAs, and Northern blot analysis has provided indirect evidence that these mRNAs are part of the same transcription unit as mRNAs for RNA polymerase II largest subunit. Another cDNA, TPP2, represents the type 2A class of phosphatases and codes for a 34.5-kDa protein of 303 amino acids. The deduced amino acid sequence has 39% and 55% identity, respectively, to the catalytic subunits of mammalian protein phosphatase 1 and 2A. Southern and Northern blot analyses are consistent with TPP2 being encoded by a single copy gene from which is derived a mRNA of 2.5 kb. This finding constitutes the first example in eukaryotes in which a single gene encodes the type 2A class of protein phosphatases. Sera from mice immunized with TPP1A fusion protein reacted with the catalytic subunits of mammalian types 1, 2A and 2B protein phosphatases. However, antisera to TPP2 fusion protein was specific for the type 2A catalytic subunit and recognized a polypeptide of 35 kDa in a Western blot of crude trypanosomal lysate.

The cytosolic and glycosomal isoenzymes of glyceraldehyde-3-phosphate dehydrogenase in Trypanosoma brucei have a distant evolutionary relationship

Trypanosoma brucei contains two isoenzymes for glyceraldehyde-3-phosphate dehydrogenase: one enzyme resides in a microbody-like organelle, the glycosome; the other is found in the cytosol. Previously we have reported the characterization of the gene for the glycosomal enzyme [Michels, P. A. M., Poliszczak, A., Osinga, K. A., Misset, O., Van Beeumen, J., Wierenga, R. K., Borst, P. & Opperdoes, F. R. (1986) EMBO J. 5, 1049-1056]. Here we describe the cloning and analysis of the gene that codes for the cytosolic isoenzyme. The gene encodes a polypeptide of 330 amino acids, with a calculated molecular mass of 35440 Da. The two isoenzymes are only 55% identical. The cytosolic glyceraldehyde-3-phosphate dehydrogenase differs from the glycosomal enzyme in the following respects: (a) its subunit molecular mass is 3.4 kDa smaller due to the absence of insertions and a small C-terminal extension which are unique to the glycosomal protein; (b) the cytosolic enzyme has a lower pI (7.9, as compared to 9.3 for the glycosomal isoenzyme), which is due to a reduction in the excess of positively charged amino acids (the calculated net charges of the polypeptides are +2 and +11, respectively). We have compared the amino acid sequences of the two T. brucei glyceraldehyde-3-phosphate dehydrogenases, with 24 available sequences of the corresponding enzyme of other organisms from various phylogenetic groups. On the basis of this comparison an evolutionary tree was constructed. This analysis strongly supports the theory that T. brucei diverged early in evolution from the main eukaryotic branch of the phylogenetic tree. Further, two separate branches for the lineages leading to Trypanosoma are inferred from the amino acid sequences, suggesting that the genes for the two glyceraldehyde-3-phosphate dehydrogenases of the trypanosome are distantly related and must have been acquired independently by the trypanosomal ancestor. The branching determined with the glycosomal enzyme precedes that found with the cytosolic enzyme. The available data do not allow us to decide which of the two genes originally belonged to the trypanosome lineage and which entered the cell later by horizontal gene transfer.

Genomic organisation of nuclear tRNAGly and tRNALeu genes in Trypanosoma brucei

We have isolated a 0.3-kb HaeIII restriction fragment from Trypanosoma brucei which contains two tRNA genes. Secondary structure models predict that the two genes identified encode tRNA molecules which specify glycine (anticodon UCC) and leucine (anticodon CAG). The two genes are separated by 86 nucleotides, transcribed in the same direction and contain features of conventional RNA polymerase III transcription units. Southern blot analysis indicates the presence of multicopy tRNa gene families in T. brucei.

Evidence for gene conversion between the phosphoglycerate kinase genes of Trypanosoma brucei

Trypanosoma brucei contains a tandem array of three genes for phosphoglycerate kinase (PGKase), genes A, B and C, each coding for a different protein. We have compared allelic variants of this gene array and find evidence for gene conversion between the three genes. Near the 3' end, the different alleles and gene B contain a variable sequence that is similar to the corresponding sequence in either gene A or gene C. This sequence is flanked by glycine triplets that are conserved in all PGKases from bacteria to mammals. The triplets are encoded by (GGT)n, resulting in sequences that resemble the recombination-promoting chi-sites of Escherichia coli. Upstream of the variable sequence, there is an area of 800 base-pairs in which genes A, B and C are highly homologous; in all three genes this region ends with a sharp boundary at which gene B again shows segmental homology with both genes A and C. These results suggest that repeated gene conversion events partially erase the differences between genes A, B and C that arise in evolution and suggest that chi-like sequences may act as recombinational hotspots in protozoa such as T. brucei.