The enigmatic RNase MRP of kinetoplastids.

The ribonucleoprotein RNase MRP is responsible for the processing of ribosomal RNA precursors. It is found in virtually all eukaryotes that have been examined. In the Euglenozoa, including the genera Euglena, Diplonema and kinetoplastids, MRP RNA and protein subunits have so far escaped detection using bioinformatic methods. However, we now demonstrate that the RNA component is widespread among the Euglenozoa and that these RNAs have secondary structures that conform to the structure of all other phylogenetic groups. In Euglena, we identified the same set of P/MRP protein subunits as in many other protists. However, we failed to identify any of these proteins in the kinetoplastids. This finding poses interesting questions regarding the structure and function of RNase MRP in these species.

Phage Origin of Mitochondrion-Localized Family A DNA Polymerases in Kinetoplastids and Diplonemids.

Mitochondria retain their own genomes as other bacterial endosymbiont-derived organelles. Nevertheless, no protein for DNA replication and repair is encoded in any mitochondrial genomes (mtDNAs) assessed to date, suggesting that the nucleus primarily governs the maintenance of mtDNA. As the proteins of diverse evolutionary origins occupy a large proportion of the current mitochondrial proteomes, we anticipate finding the same evolutionary trend in the nucleus-encoded machinery for mtDNA maintenance. Indeed, none of the DNA polymerases (DNAPs) in the mitochondrial endosymbiont, a putative α-proteobacterium, seemingly had been inherited by their descendants (mitochondria), as none of the known types of mitochondrion-localized DNAP showed a specific affinity to the α-proteobacterial DNAPs. Nevertheless, we currently have no concrete idea of how and when the known types of mitochondrion-localized DNAPs emerged. We here explored the origins of mitochondrion-localized DNAPs after the improvement of the samplings of DNAPs from bacteria and phages/viruses. Past studies have revealed that a set of mitochondrion-localized DNAPs in kinetoplastids and diplonemids, namely PolIB, PolIC, PolID, PolI-Perk1/2, and PolI-dipl (henceforth designated collectively as "PolIBCD+") have emerged from a single DNAP. In this study, we recovered an intimate connection between PolIBCD+ and the DNAPs found in a particular group of phages. Thus, the common ancestor of kinetoplastids and diplonemids most likely converted a laterally acquired phage DNAP into a mitochondrion-localized DNAP that was ancestral to PolIBCD+. The phage origin of PolIBCD+ hints at a potentially large contribution of proteins acquired via nonvertical processes to the machinery for mtDNA maintenance in kinetoplastids and diplonemids.

Phosphoglycerate kinase: structural aspects and functions, with special emphasis on the enzyme from Kinetoplastea.

Phosphoglycerate kinase (PGK) is a glycolytic enzyme that is well conserved among the three domains of life. PGK is usually a monomeric enzyme of about 45 kDa that catalyses one of the two ATP-producing reactions in the glycolytic pathway, through the conversion of 1,3-bisphosphoglycerate (1,3BPGA) to 3-phosphoglycerate (3PGA). It also participates in gluconeogenesis, catalysing the opposite reaction to produce 1,3BPGA and ADP. Like most other glycolytic enzymes, PGK has also been catalogued as a moonlighting protein, due to its involvement in different functions not associated with energy metabolism, which include pathogenesis, interaction with nucleic acids, tumorigenesis progression, cell death and viral replication. In this review, we have highlighted the overall aspects of this enzyme, such as its structure, reaction kinetics, activity regulation and possible moonlighting functions in different protistan organisms, especially both free-living and parasitic Kinetoplastea. Our analysis of the genomes of different kinetoplastids revealed the presence of open-reading frames (ORFs) for multiple PGK isoforms in several species. Some of these ORFs code for unusually large PGKs. The products appear to contain additional structural domains fused to the PGK domain. A striking aspect is that some of these PGK isoforms are predicted to be catalytically inactive enzymes or 'dead' enzymes. The roles of PGKs in kinetoplastid parasites are analysed, and the apparent significance of the PGK gene duplication that gave rise to the different isoforms and their expression in Trypanosoma cruzi is discussed.

DNA Topoisomerases of Kinetoplastid Parasites: Brief Overview and Recent Perspectives.

Topoisomerases are a group of enzymes that resolve DNA topological problems and aid in different DNA transaction processes viz. replication, transcription, recombination, etc. inside cells. These proteins accomplish their feats by steps of DNA strand(s) scission, strand passage or rotation and subsequent rejoining activities. Topoisomerases of kinetoplastid parasites have been extensively studied because of their unusual features. The unique presence of heterodimeric Type IB topoisomerase and prokaryotic 'TopA homologue' Type IA topoisomerase in kinetoplastids still generates immense interest among scientists. Moreover, because of their structural dissimilarity with the host enzymes, topoisomerases of kinetoplastid parasites are attractive targets for chemotherapeutic interventions to kill these deadly parasites. In this review, we summarize historical perspectives and recent advances in kinetoplastid topoisomerase research and how these proteins are exploited for drug targeting.

Novel 8-nitroquinolin-2(1H)-ones as NTR-bioactivated antikinetoplastid molecules: Synthesis, electrochemical and SAR study.

To study the antiparasitic 8-nitroquinolin-2(1H)-one pharmacophore, a series of 31 derivatives was synthesized in 1-5 steps and evaluated in vitro against both Leishmania infantum and Trypanosoma brucei brucei. In parallel, the reduction potential of all molecules was measured by cyclic voltammetry. Structure-activity relationships first indicated that antileishmanial activity depends on an intramolecular hydrogen bond (described by X-ray diffraction) between the lactam function and the nitro group, which is responsible for an important shift of the redox potential (+0.3 V in comparison with 8-nitroquinoline). With the assistance of computational chemistry, a set of derivatives presenting a large range of redox potentials (from -1.1 to -0.45 V) was designed and provided a list of suitable molecules to be synthesized and tested. This approach highlighted that, in this series, only substrates with a redox potential above -0.6 V display activity toward L. infantum. Nevertheless, such relation between redox potentials and in vitro antiparasitic activities was not observed in T. b. brucei. Compound 22 is a new hit compound in the series, displaying both antileishmanial and antitrypanosomal activity along with a low cytotoxicity on the human HepG2 cell line. Compound 22 is selectively bioactivated by the type 1 nitroreductases (NTR1) of L. donovani and T. brucei brucei. Moreover, despite being mutagenic in the Ames test, as most of nitroaromatic derivatives, compound 22 was not genotoxic in the comet assay. Preliminary in vitro pharmacokinetic parameters were finally determined and pointed out a good in vitro microsomal stability (half-life > 40 min) and a 92% binding to human albumin.

Shedding light on lipid metabolism in Kinetoplastida: A phylogenetic analysis of phospholipase D protein homologs.

Unicellular flagellates that make up the class Kinetoplastida include multiple parasites responsible for public health concerns, including Trypanosoma brucei and T. cruzi (agents of African sleeping sickness and Chagas disease, respectively), and various Leishmania species, which cause leishmaniasis. These diseases are generally difficult to eradicate, with treatments often having lethal side effects and/or being effective only during the acute phase of the diseases, when most patients are still asymptomatic. Phospholipid signaling and metabolism are important in the different life stages of Trypanosoma, including playing a role in transitions between stages and in immune system evasion, thus, making the responsible enzymes into potential therapeutic targets. However, relatively little is understood about how the pathways function in these pathogens. Thus, in this study we examined evolutionary history of proteins from one such signaling pathway, namely phospholipase D (PLD) homologs. PLD is an enzyme responsible for synthesizing phosphatidic acid (PA) from membrane phospholipids. PA is not only utilized for phospholipid synthesis, but is also involved in many other signaling pathways, including biotic and abiotic stress response. 37 different representative Kinetoplastida genomes were used for an exhaustive search to identify putative PLD homologs. The genome of Bodo saltans was the only one of surveyed Kinetoplastida genomes that encoded a protein that clustered with plant PLDs. The representatives from other Kinetoplastida species clustered together in two different clades, thought to be homologous to the PLD superfamily, but with shared sequence similarity with cardiolipin synthases (CLS), and phosphatidylserine synthases (PSS). The protein structure predictions showed that most Kinetoplastida sequences resemble CLS and PSS, with the exception of 5 sequences from Bodo saltans that shared significant structural similarities with the PLD sequences, suggesting the loss of PLD-like sequences during the evolution of parasitism in kinetoplastids. On the other hand, diacylglycerol kinase (DGK) homologs were identified for all species examined in this study, indicating that DGK could be the only pathway for the synthesis of PA involved in lipid signaling in these organisms due to genome streamlining during transition to parasitic lifestyle. Our findings offer insights for development of potential therapeutic and/or intervention approaches, particularly those focused on using PA, PLD and/or DGK related pathways, against trypanosomiasis, leishmaniasis, and Chagas disease.

Novel scaffolds for inhibition of Cruzipain identified from high-throughput screening of anti-kinetoplastid chemical boxes.

American Trypanosomiasis or Chagas disease is a prevalent, neglected and serious debilitating illness caused by the kinetoplastid protozoan parasite Trypanosoma cruzi. The current chemotherapy is limited only to nifurtimox and benznidazole, two drugs that have poor efficacy in the chronic phase and are rather toxic. In this scenario, more efficacious and safer drugs, preferentially acting through a different mechanism of action and directed against novel targets, are particularly welcome. Cruzipain, the main papain-like cysteine peptidase of T. cruzi, is an important virulence factor and a chemotherapeutic target with excellent pre-clinical validation evidence. Here, we present the identification of new Cruzipain inhibitory scaffolds within the GlaxoSmithKline HAT (Human African Trypanosomiasis) and Chagas chemical boxes, two collections grouping 404 non-cytotoxic compounds with high antiparasitic potency, drug-likeness, structural diversity and scientific novelty. We have adapted a continuous enzymatic assay to a medium-throughput format and carried out a primary screening of both collections, followed by construction and analysis of dose-response curves of the most promising hits. Using the identified compounds as a starting point a substructure directed search against CHEMBL Database revealed plausible common scaffolds while docking experiments predicted binding poses and specific interactions between Cruzipain and the novel inhibitors.

Proteasome inhibition for treatment of leishmaniasis, Chagas disease and sleeping sickness.

Chagas disease, leishmaniasis and sleeping sickness affect 20 million people worldwide and lead to more than 50,000 deaths annually. The diseases are caused by infection with the kinetoplastid parasites Trypanosoma cruzi, Leishmania spp. and Trypanosoma brucei spp., respectively. These parasites have similar biology and genomic sequence, suggesting that all three diseases could be cured with drugs that modulate the activity of a conserved parasite target. However, no such molecular targets or broad spectrum drugs have been identified to date. Here we describe a selective inhibitor of the kinetoplastid proteasome (GNF6702) with unprecedented in vivo efficacy, which cleared parasites from mice in all three models of infection. GNF6702 inhibits the kinetoplastid proteasome through a non-competitive mechanism, does not inhibit the mammalian proteasome or growth of mammalian cells, and is well-tolerated in mice. Our data provide genetic and chemical validation of the parasite proteasome as a promising therapeutic target for treatment of kinetoplastid infections, and underscore the possibility of developing a single class of drugs for these neglected diseases.

Comparative Metabolism of Free-living Bodo saltans and Parasitic Trypanosomatids.

Comparison of the genomes of free-living Bodo saltans and those of parasitic trypanosomatids reveals that the transition from a free-living to a parasitic life style has resulted in the loss of approximately 50% of protein-coding genes. Despite this dramatic reduction in genome size, B. saltans and trypanosomatids still share a significant number of common metabolic traits: glycosomes; a unique set of the pyrimidine biosynthetic pathway genes; an ATP-PFK which is homologous to the bacterial PPi -PFKs rather than to the canonical eukaryotic ATP-PFKs; an alternative oxidase; three phosphoglycerate kinases and two GAPDH isoenzymes; a pyruvate kinase regulated by fructose-2,6-bisphosphate; trypanothione as a substitute for glutathione; synthesis of fatty acids via a unique set of elongase enzymes; and a mitochondrial acetate:succinate coenzyme A transferase. B. saltans has lost the capacity to synthesize ubiquinone. Among genes that are present in B. saltans and lost in all trypanosomatids are those involved in the degradation of mureine, tryptophan and lysine. Novel acquisitions of trypanosomatids are components of pentose sugar metabolism, pteridine reductase and bromodomain-factor proteins. In addition, only the subfamily Leishmaniinae has acquired a gene for catalase and the capacity to convert diaminopimelic acid to lysine.

Evolutionary comparison of prenylation pathway in kinetoplastid Leishmania and its sister Leptomonas.

BACKGROUND: Leptomonas is monogenetic kinetoplastid parasite of insects and is primitive in comparison to Leishmania. Comparative studies of these two kinetoplastid may share light on the evolutionary transition to dixenous parasitism in Leishmania. In order to adapt and survive within two hosts, Leishmania species must have acquired virulence factors in addition to mechanisms that mediate susceptibility/resistance to infection in the pathology associated with disease. Rab proteins are key mediators of vesicle transport and contribute greatly to the evolution of complexity of membrane transport system. In this study we used our whole genome sequence data of these two divergent kinetoplastids to analyze the orthologues/paralogues of Rab proteins. RESULTS: During change of lifestyle from monogenetic (Leptomonas) to digenetic (Leishmania), we found that the prenyl machinery remained unchanged. Geranylgeranyl transferase-I (GGTase-I) was absent in both Leishmania and its sister Leptomonas. Farnesyltransferase (FTase) and geranylgeranyl transferase-II (GGTase-II) were identified for protein prenylation. We predict that activity of the missing alpha-subunit (α-subunit) of GGTase-II in Leptomonas was probably contributed by the α-subunit of FTase, while beta-subunit (ß-subunit) of GGTase-II was conserved and indicated functional conservation in the evolution of these two kinetoplastids. Therefore the ß-subunit emerges as an excellent target for compounds inhibiting parasite activity in clinical cases of co-infections. We also confirmed that during the evolution to digenetic life style in Leishmania, the parasite acquired capabilities to evade drug action and maintain parasite virulence in the host with the incorporation of short-chain dehydrogenase/reductase (SDR/MDR) superfamily in Rab genes. CONCLUSION: Our study based on whole genome sequences is the first to build comparative evolutionary analysis and identification of prenylation proteins in Leishmania and its sister Leptomonas. The information presented in our present work has importance for drug design targeted to kill L. donovani in humans but not affect the human form of the prenylation enzymes.

Targeting Lysine Deacetylases (KDACs) in Parasites.

Due to an increasing problem of drug resistance among almost all parasites species ranging from protists to worms, there is an urgent need to explore new drug targets and their inhibitors to provide new and effective parasitic therapeutics. In this regard, there is growing interest in exploring known drug leads of human epigenetic enzymes as potential starting points to develop novel treatments for parasitic diseases. This approach of repurposing (starting with validated targets and inhibitors) is quite attractive since it has the potential to reduce the expense of drug development and accelerate the process of developing novel drug candidates for parasite control. Lysine deacetylases (KDACs) are among the most studied epigenetic drug targets of humans, and a broad range of small-molecule inhibitors for these enzymes have been reported. In this work, we identify the KDAC protein families in representative species across important classes of parasites, screen a compound library of 23 hydroxamate- or benzamide-based small molecules KDAC inhibitors, and report their activities against a range of parasitic species, including the pathogen of malaria (Plasmodium falciparum), kinetoplastids (Trypanosoma brucei and Leishmania donovani), and nematodes (Brugia malayi, Dirofilaria immitis and Haemonchus contortus). Compound activity against parasites is compared to that observed against the mammalian cell line (L929 mouse fibroblast) in order to determine potential parasite-versus-host selectivity). The compounds showed nanomolar to sub-nanomolar potency against various parasites, and some selectivity was observed within the small panel of compounds tested. The possible binding modes of the active compounds at the different protein target sites within different species were explored by docking to homology models to help guide the discovery of more selective, parasite-specific inhibitors. This current work supports previous studies that explored the use of KDAC inhibitors in targeting Plasmodium to develop new anti-malarial treatments, and also pioneers experiments with these KDAC inhibitors as potential new anthelminthics. The selectivity observed begins to address the challenges of targeting specific parasitic diseases while limiting host toxicity.

Trypanosoma brucei harbours a divergent XPB helicase paralogue that is specialized in nucleotide excision repair and conserved among kinetoplastid organisms.

Conserved from yeast to humans, TFIIH is essential for RNA polymerase II transcription and nucleotide excision repair (NER). TFIIH consists of a core that includes the DNA helicase Xeroderma pigmentosum B (XPB) and a kinase subcomplex. Trypanosoma brucei TFIIH harbours all core complex components and is indispensable for RNA polymerase II transcription of spliced leader RNA genes (SLRNAs). Kinetoplastid organisms, however, possess two highly divergent XPB paralogues with only the larger being identified as a TFIIH subunit in T. brucei. Here we show that a knockout of the gene for the smaller paralogue, termed XPB-R (R for repair) resulted in viable cultured trypanosomes that grew slower than normal. XPB-R depletion did not affect transcription in vivo or in vitro and XPB-R was not found to occupy the SLRNA promoter which assembles a RNA polymerase II transcription pre-initiation complex including TFIIH. However, XPB-R(-/-) cells were much less tolerant than wild-type cells to UV light- and cisplatin-induced DNA damage, which require NER. Since XPB-R(-/-) cells were not impaired in DNA base excision repair, XPB-R appears to function specifically in NER. Interestingly, several other protists possess highly divergent XPB paralogues suggesting that XPBs specialized in transcription or NER exist beyond the Kinetoplastida.

Lipid synthesis in protozoan parasites: a comparison between kinetoplastids and apicomplexans.

Lipid metabolism is of crucial importance for pathogens. Lipids serve as cellular building blocks, signalling molecules, energy stores, posttranslational modifiers, and pathogenesis factors. Parasites rely on a complex system of uptake and synthesis mechanisms to satisfy their lipid needs. The parameters of this system change dramatically as the parasite transits through the various stages of its life cycle. Here we discuss the tremendous recent advances that have been made in the understanding of the synthesis and uptake pathways for fatty acids and phospholipids in apicomplexan and kinetoplastid parasites, including Plasmodium, Toxoplasma, Cryptosporidium, Trypanosoma and Leishmania. Lipid synthesis differs in significant ways between parasites from both phyla and the human host. Parasites have acquired novel pathways through endosymbiosis, as in the case of the apicoplast, have dramatically reshaped substrate and product profiles, and have evolved specialized lipids to interact with or manipulate the host. These differences potentially provide opportunities for drug development. We outline the lipid pathways for key species in detail as they progress through the developmental cycle and highlight those that are of particular importance to the biology of the pathogens and/or are the most promising targets for parasite-specific treatment.

RNA-seq-based metatranscriptomic and microscopic investigation reveals novel metalloproteases of Neobodo sp. as potential virulence factors for soft tunic syndrome in Halocynthia roretzi.

Bodonids and trypanosomatids are derived from a common ancestor with the bodonids being a more primitive lineage. The Neobodonida, one of the three clades of bodonids, can be free-living, commensal or parasitic. Despite the ecological and evolutionary significance of these organisms, however, many of their biological and pathological features are currently unknown. Here, we employed metatranscriptomics using RNA-seq technology combined with field-emission microscopy to reveal the virulence factors of a recently described genus of Neobodonida that is considered to be responsible for ascidian soft tunic syndrome (AsSTS), but whose pathogenesis is unclear. Our microscopic observation of infected tunic tissues suggested putative virulence factors, enabling us to extract novel candidate transcripts; these included cysteine proteases of the families C1 and C2, serine proteases of S51 and S9 families, and metalloproteases grouped into families M1, M3, M8, M14, M16, M17, M24, M41, and M49. Protease activity/inhibition assays and the estimation of expression levels within gene clusters allowed us to identify metalloprotease-like enzymes as potential virulence attributes for AsSTS. Furthermore, a multimarker-based phylogenetic analysis using 1,184 concatenated amino acid sequences clarified the order Neobodo sp. In sum, we herein used metatranscriptomics to elucidate the in situ expression profiles of uncharacterized putative transcripts of Neobodo sp., combined these results with microscopic observation to select candidate genes relevant to pathogenesis, and used empirical screening to define important virulence factors.

Therapeutic potential of phosphodiesterase inhibitors in parasitic diseases.

Protozoan parasites of the order kinetoplastida are the causative agents of three of the world's most important neglected human diseases: African trypanosomiasis, American trypanosomiasis, and leishmaniasis. Current therapies are limited, with some treatments having serious and sometimes lethal side effects. The growing number of cases that are refractory to treatment is also of concern. With few new drugs in development, there is an unmet medical need for new, more effective, and safer medications. Recent studies employing genetic and pharmacological techniques have begun to shed light on the role of the cyclic nucleotide phosphodiesterases in the life cycle of these pathogens and suggest that these important regulators of cyclic nucleotide signaling may be promising new targets for the treatment of parasitic diseases.

Crystal structure of Leishmania major ADP-ribosylation factor-like 1 and a classification of related GTPase family members in this Kinetoplastid.

ADP-ribosylation factor-like (ARL) proteins are small GTPases that undergo conformational changes upon nucleotide binding, and which regulate the affinity of ARLs for binding other proteins, lipids or membranes. There is a paucity of structural data on this family of proteins in the Kinetoplastida, despite studies implicating them in key events related to vesicular transport and regulation of microtubule-dependent processes. The crystal structure of Leishmania major ARL1 in complex with GDP has been determined to 2.1 Å resolution and reveals a high degree of structural conservation with human ADP-ribosylation factor 1 (ARF1). Putative L. major and Trypanosoma brucei ARF/ARL family members have been classified based on structural considerations, amino acid sequence conservation combined with functional data on Kinetoplastid and human orthologues. This classification may guide future studies designed to elucidate the function of specific family members.

Kinetoplastid RNA editing ligases 1 and 2 exhibit different electrostatic properties.

Kinetoplastid RNA editing ligases 1 and 2 (KREL1 and KREL2) share a significant degree of sequence homology. However, biochemical experiments have reported that KREL1 and KREL2 differ in their functional roles during the RNA editing process. In this study, we hypothesize that dissimilar roles for KREL1 and KREL2 proteins arise from their different physicochemical characteristics. To test our hypothesis at sequence level, we plotted theoretical titration curves for KREL1, KREL2 and their binding partner proteins. The plots showed a lower isoelectric point for KREL1 compared to that for KREL2 as well as more relative alkalinity and acidity for binding partner proteins of KREL1 and KREL2 at net charge zero, respectively. At structure level, based on the available high resolution structure of KREL1 N-terminal domain and strong sequence similarity between KRELs and other ligases, we built the homology model of KREL2 N-terminal domain. Using Poisson-Boltzmann continuum approach, we calculated the electrostatic potential isosurfaces of KREL1 structure and KREL2 model. KREL1 and KREL2 coordinates differed in their electrostatic isopotential patterns. A wider negative patch on the surface of KREL1 suggests differential affinity for another protein compared to KREL2. In contrast, a larger positive patch on the KREL2 surface predicts its differential affinity and/or specificity for its RNA substrate. Subsequently, we employed in silico mutational scanning and identified the surface-exposed residues contributing to the long-range electrostatic energy of KRELs. We predict that two structurally conserved loops of KRELs, not previously reported in the literature, also recognize their RNA substrates. Our results provide important information about the physicochemical properties of RNA editing ligases that could contribute to the ligation step of RNA editing.

Bodo sp., a free-living flagellate, expresses divergent proteolytic activities from the closely related parasitic trypanosomatids.

We report the characterization of cell-associated and extracellular peptidases of Bodo sp., a free-living flagellate of the Bodonidae family, order Kinetoplastida, which is considered ancestral to the trypanosomatids. This bodonid isolate is phylogenetically related to Bodo caudatus and Bodo curvifilus. The proteolytic activity profiles of Bodo sp. were determined by sodium dodecyl sulfate-polyacrylamide gel electrophoresis containing co-polymerized gelatin, casein, hemoglobin, or bovine serum albumin as substrates. The enzymatic complex degraded gelatin better in acidic pH, and under these conditions four proteolytic bands (120, 100, 90, and 75 kDa) were detected in the cellular or extracellular extracts. Two peptidases (250 and 200 kDa) were exclusively detected with the substrate casein. All these enzymes belong to the serine peptidase class, based on inhibition by aprotinin and phenylmethylsulfonyl fluoride. This is the first biochemical characterization of peptidases in a free-living Bodo sp., potentially providing insight into the physiology of these protozoa and the evolutionary importance of peptidases to the order Kinetoplastida as some of these enzymes are important virulence factors in pathogenic trypanosomatids.

Kinetoplastid papain-like cysteine peptidases.

Cysteine peptidases are important for growth and survival of kinetoplastid parasites. The best characterised are those homologous to mammalian cathepsins B and L. To address a somewhat confusing terminology, we introduce a unifying nomenclature for kinetoplastid CATB and CATL peptidases. We review their evolutionary relatedness, genomic organisation, developmental expression, subcellular location and physiological functions. In addition, the applications of kinetoplastid CATB and CATL enzymes as vaccine candidates, diagnostic markers and drug targets are discussed.

Unassigned MURF1 of kinetoplastids codes for NADH dehydrogenase subunit 2.

BACKGROUND: In a previous study, we conducted a large-scale similarity-free function prediction of mitochondrion-encoded hypothetical proteins, by which the hypothetical gene murf1 (maxicircle unidentified reading frame 1) was assigned as nad2, encoding subunit 2 of NADH dehydrogenase (Complex I of the respiratory chain). This hypothetical gene occurs in the mitochondrial genome of kinetoplastids, a group of unicellular eukaryotes including the causative agents of African sleeping sickness and leishmaniasis. In the present study, we test this assignment by using bioinformatics methods that are highly sensitive in identifying remote homologs and confront the prediction with available biological knowledge. RESULTS: Comparison of MURF1 profile Hidden Markov Model (HMM) against function-known profile HMMs in Pfam, Panther and TIGR shows that MURF1 is a Complex I protein, but without specifying the exact subunit. Therefore, we constructed profile HMMs for each individual subunit, using all available sequences clustered at various identity thresholds. HMM-HMM comparison of these individual NADH subunits against MURF1 clearly identifies this hypothetical protein as NAD2. Further, we collected the relevant experimental information about kinetoplastids, which provides additional evidence in support of this prediction. CONCLUSION: Our in silico analyses provide convincing evidence for MURF1 being a highly divergent member of NAD2.