PlasmoDB: exploring genomics and post-genomics data of the malaria parasite,Plasmodium falciparum

Single-genome sequencing reveals within-host evolution of human malaria parasites

Population genomics of bulk malaria infections is unable to examine intrahost evolution; therefore, most work has focused on the role of recombination in generating genetic variation. We used single-cell sequencing protocol for low-parasitaemia infections to generate 406 near-complete single Plasmodium vivax genomes from 11 patients sampled during sequential febrile episodes. Parasite genomes contain hundreds of de novo mutations, showing strong signatures of selection, which are enriched in the ApiAP2 family of transcription factors, known targets of adaptation. Comparing 315 P. falciparum single-cell genomes from 15 patients with our P. vivax data, we find broad complementary patterns of de novo mutation at the gene and pathway level, revealing the importance of within-host evolution during malaria infections.

Atomic Resolution Homology Models and Molecular Dynamics Simulations of Plasmodium falciparum Tubulins

Microtubules are tubulin polymers present in the eukaryotic cytoskeleton essential for structural stability and cell division that are also roadways for intracellular transport of vesicles and organelles. In the human malaria parasite Plasmodium falciparum, apart from providing structural stability and cell division, microtubules also facilitate important biological activities crucial for parasite survival in hosts, such as egression and motility. Hence, parasite structures and processes involving microtubules are among the most important drug targets for discovering much-needed novel Plasmodium inhibitors. The current study aims to construct reliable and high-quality 3D models of α-, β-, and γ-tubulins using various modeling techniques. We identified a common binding pocket specific to Plasmodium α-, β-, and γ-tubulins. Molecular dynamics simulations confirmed the stability of the Plasmodium tubulin 3D structures. The models generated in the present study may be used for protein-protein and protein-drug interaction investigations targeted toward designing malaria parasite tubulin-specific inhibitors.

Sequence analysis and characterization of pyruvate kinase from Clonorchis sinensis, a 53.1-kDa homopentamer, implicated immune protective efficacy against clonorchiasis

BACKGROUND

Clonorchis sinensis, the causative agent of clonorchiasis, is classified as one of the most neglected tropical diseases and affects more than 15 million people globally. This hepatobiliary disease is highly associated with cholangiocarcinoma. As key molecules in the infectivity and subsistence of trematodes, glycolytic enzymes have been targets for drug and vaccine development. Clonorchis sinensis pyruvate kinase (CsPK), a crucial glycolytic enzyme, was characterized in this research.

RESULTS

Differences were observed in the sequences and spatial structures of CsPK and PKs from humans, rats, mice and rabbits. CsPK possessed a characteristic active site signature (IKLIAKIENHEGV) and some unique sites but lacked the N-terminal domain. The predicted subunit molecular mass (Mr) of CsPK was 53.1 kDa. Recombinant CsPK (rCsPK) was a homopentamer with a Mr. of approximately 290 kDa by both native PAGE and gel filtration chromatography. Significant differences in the protein and mRNA levels of CsPK were observed among four life stages of C. sinensis (egg, adult worm, excysted metacercaria and metacercaria), suggesting that these developmental stages may be associated with diverse energy demands. CsPK was widely distributed in adult worms. Moreover, an intense Th1-biased immune response was persistently elicited in rats immunized with rCsPK. Also, rat anti-rCsPK sera suppressed C. sinensis adult subsistence both in vivo and in vitro.

CONCLUSIONS

The sequences and spatial structures, molecular mass, and expression profile of CsPK have been characterized. rCsPK was indicated to be a homopentamer. Rat anti-rCsPK sera suppressed C. sinensis adult subsistence both in vivo and in vitro. CsPK is worthy of further study as a promising target for drug and vaccine development.

Inhibition of the SR protein-phosphorylating CLK kinases of Plasmodium falciparum impairs blood stage replication and malaria transmission

Cyclin-dependent kinase-like kinases (CLKs) are dual specificity protein kinases that phosphorylate Serine/Arginine-rich (SR) proteins involved in pre-mRNA processing. Four CLKs, termed PfCLK-1-4, can be identified in the human malaria parasite Plasmodium falciparum, which show homology with the yeast SR protein kinase Sky1p. The four PfCLKs are present in the nucleus and cytoplasm of the asexual blood stages and of gametocytes, sexual precursor cells crucial for malaria parasite transmission from humans to mosquitoes. We identified three plasmodial SR proteins, PfSRSF12, PfSFRS4 and PfSF-1, which are predominantly present in the nucleus of blood stage trophozoites, PfSRSF12 and PfSF-1 are further detectable in the nucleus of gametocytes. We found that recombinantly expressed SR proteins comprising the Arginine/Serine (RS)-rich domains were phosphorylated by the four PfCLKs in in vitro kinase assays, while a recombinant PfSF-1 peptide lacking the RS-rich domain was not phosphorylated. Since it was hitherto not possible to knock-out the pfclk genes by conventional gene disruption, we aimed at chemical knock-outs for phenotype analysis. We identified five human CLK inhibitors, belonging to the oxo-β-carbolines and aminopyrimidines, as well as the antiseptic chlorhexidine as PfCLK-targeting compounds. The six inhibitors block P. falciparum blood stage replication in the low micromolar to nanomolar range by preventing the trophozoite-to-schizont transformation. In addition, the inhibitors impair gametocyte maturation and gametogenesis in in vitro assays. The combined data show that the four PfCLKs are involved in phosphorylation of SR proteins with essential functions for the blood and sexual stages of the malaria parasite, thus pointing to the kinases as promising targets for antimalarial and transmission blocking drugs.

Two nucleus-localized CDK-like kinases with crucial roles for malaria parasite erythrocytic replication are involved in phosphorylation of splicing factor

The kinome of the human malaria parasite Plasmodium falciparum comprises representatives of most eukaryotic protein kinase groups, including kinases which regulate proliferation and differentiation processes. Despite extensive research on most plasmodial enzymes, little information is available regarding the four identified members of the cyclin-dependent kinase-like kinase (CLK) family. In other eukaryotes, CLKs regulate mRNA splicing through phosphorylation of Serine/Arginine-rich proteins. Here, we investigate two of the PfCLKs, the Lammer kinase homolog PfCLK-1, and PfCLK-2. Both PfCLKs show homology with the yeast Serine/Arginine protein kinase Sky1p and are transcribed throughout the asexual blood stages and in gametocytes. PfCLK-1/Lammer possesses two nuclear localization signal sites and PfCLK-2 possesses one of these signal sites upstream of the C-terminal catalytic domains. Indirect immunofluorescence, Western blot, and electron microscopy data confirm that the kinases are primarily localized in the parasite nucleus, and PfCLK-2 is further present in the cytoplasm. The two kinases are important for completion of the asexual replication cycle of P. falciparum, as demonstrated by reverse genetics approaches. In vitro kinase assays show substrate phosphorylation by the PfCLKs, including the Sky1p substrate, splicing factor Npl3p, and the plasmodial alternative splicing factor PfASF-1. Mass spectrometric analysis of co-immunoprecipitated proteins indicates assembly of the two PfCLKs with proteins with predicted nuclease, phosphatase, or helicase functions. Our data indicate a crucial role of PfCLKs for malaria blood stage parasites, presumably by participating in gene regulation through the post-transcriptional modification of mRNA.

Inconsistencies of genome annotations in apicomplexan parasites revealed by 5'-end-one-pass and full-length sequences of oligo-capped cDNAs

Background

Apicomplexan parasites are causative agents of various diseases including malaria and have been targets of extensive genomic sequencing. We generated 5'-EST collections for six apicomplexa parasites using our full-length oligo-capping cDNA library method. To improve upon the current genome annotations, as well as to validate the importance for physical cDNA clone resources, we generated a large-scale collection of full-length cDNAs for several apicomplexa parasites.

Results

In this study, we used a total of 61,056 5'-end-single-pass cDNA sequences from Plasmodium falciparum, P. vivax, P. yoelii, P. berghei, Cryptosporidium parvum, and Toxoplasma gondii. We compared these partially sequenced cDNA sequences with the currently annotated gene models and observed significant inconsistencies between the two datasets. In particular, we found that on average 14% of the exons in the current gene models were not supported by any cDNA evidence, and that 16% of the current gene models may contain at least one mis-annotation and should be re-evaluated. We also identified a large number of transcripts that had been previously unidentified. For 732 cDNAs in T. gondii, the entire sequences were determined in order to evaluate the annotated gene models at the complete full-length transcript level. We found that 41% of the T. gondii gene models contained at least one inconsistency. We also identified and confirmed by RT-PCR 140 previously unidentified transcripts found in the intergenic regions of the current gene annotations. We show that the majority of these discrepancies are due to questionable predictions of one or two extra exons in the upstream or downstream regions of the genes.

Conclusion

Our data indicates that the current gene models are likely to still be incomplete and have much room for improvement. Our unique full-length cDNA information is especially useful for further refinement of the annotations for the genomes of apicomplexa parasites.

A complex-centric view of protein network evolution

The recent availability of protein-protein interaction networks for several species makes it possible to study protein complexes in an evolutionary context. In this article, we present a novel network-based framework for reconstructing the evolutionary history of protein complexes. Our analysis is based on generalizing evolutionary measures for single proteins to the level of whole subnetworks, comprehensively considering a broad set of computationally derived complexes and accounting for both sequence and interaction changes. Specifically, we compute sets of orthologous complexes across species, and use these to derive evolutionary rate and age measures for protein complexes. We observe significant correlations between the evolutionary properties of a complex and those of its member proteins, suggesting that protein complexes form early in evolution and evolve as coherent units. Additionally, our approach enables us to directly quantify the extent to which gene duplication has played a role in the evolution of complexes. We find that about one quarter of the sets of orthologous complexes have originated from evolutionary cores of homodimers that underwent duplication and divergence, testifying to the important role of gene duplication in protein complex evolution.

MADIBA: a web server toolkit for biological interpretation of Plasmodium and plant gene clusters

Background

Microarray technology makes it possible to identify changes in gene expression of an organism, under various conditions. Data mining is thus essential for deducing significant biological information such as the identification of new biological mechanisms or putative drug targets. While many algorithms and software have been developed for analysing gene expression, the extraction of relevant information from experimental data is still a substantial challenge, requiring significant time and skill.

Description

MADIBA (MicroArray Data Interface for Biological Annotation) facilitates the assignment of biological meaning to gene expression clusters by automating the post-processing stage. A relational database has been designed to store the data from gene to pathway for Plasmodium, rice and Arabidopsis. Tools within the web interface allow rapid analyses for the identification of the Gene Ontology terms relevant to each cluster; visualising the metabolic pathways where the genes are implicated, their genomic localisations, putative common transcriptional regulatory elements in the upstream sequences, and an analysis specific to the organism being studied.

Conclusion

MADIBA is an integrated, online tool that will assist researchers in interpreting their results and understand the meaning of the co-expression of a cluster of genes. Functionality of MADIBA was validated by analysing a number of gene clusters from several published experiments – expression profiling of the Plasmodium life cycle, and salt stress treatments of Arabidopsis and rice. In most of the cases, the same conclusions found by the authors were quickly and easily obtained after analysing the gene clusters with MADIBA.

Characterization of amino acid variation at strategic positions in parasite and human proteases for selective inhibition of falcipains in Plasmodium falciparum

Falcipains (FP) of Plasmodium falciparum are important virulence factors marked as potential targets for antimalarial drug discovery. In this study, the previously uncharacterized fp2B (PF11_0161) was shown to be highly expressed as an active enzyme during the erythrocytic stage. With three related proteases in the FP family and the existence of human homologues, it is prudent to identify clusters of residues unique to the parasite proteases that can be targeted selectively for drug design. Using bioinformatic tools, we have carefully mapped out a highly conserved and unique region constituted by I85, S149, and A151 in the plasmodial proteases that can influence the development of compounds capable of inhibiting the entire FP family. Taking drug interactions with the human homologues into consideration, these residues in FP2B were replaced with the cognate residues found in human cathepsin L (catL) for evaluation. Despite the high sequence similarity between the FP2 isozymes (97.5%), FP2B is found to be more tolerant to amino acid substitution at position 149 than FP2A. This structural disparity implied that residues mediating peptide substrate interactions are not fully conserved across the FP family and warrant attention in the design and evaluation of protease inhibitors focused on the FPs. The simultaneous substitution of the neighboring residues (I85 or A151) rendered the double mutants (S149A/I85M and S149A/A151D) completely inactive. Significantly, the mutations did not result in 'catL-like' specificity, suggesting that substrate-based inhibitors could be rationally designed against these important parasite-specific structural determinants.

Functional analysis, overexpression, and kinetic characterization of pyruvate kinase from Plasmodium falciparum

The important role of pyruvate kinase during malarial infection has prompted the cloning of a cDNA encoding Plasmodium falciparum pyruvate kinase (pfPyrK), using mRNA from intraerythrocytic-stage malaria parasites. The full-length cDNA encodes a protein with a computed molecular weight of 55.6 kDa and an isoelectric point of 7.5. The purified recombinant pfPyrK is enzymatically active and exists as a homotetramer in its active form. The enzyme exhibits hyperbolic kinetics with respect to phosphoenolpyruvate and ADP, with K(m) of 0.19 and 0.12 mM, respectively. pfPyrK is not affected by fructose-1,6-bisphosphate, a general activating factor of pyruvate kinase for most species. Glucose-6-phosphate, an activator of the Toxoplasma gondii enzyme, does not affect pfPyrK activity. Similar to rabbit pyruvate kinase, pfPyrK is susceptible to inactivation by 1mM pyridoxal-5'-phosphate, but to a lesser extent. A screen for inhibitors to pfPyrK revealed that it is markedly inhibited by ATP and citrate. Detailed kinetic analysis revealed a transition from hyperbolic to sigmoidal kinetics for PEP in the presence of citrate, as well as competitive inhibitory behavior for ATP with respect to PEP. Citrate exhibits non-competitive inhibition with respect to ADP with a K(i) of 0.8mM. In conclusion, P. falciparum expresses an active pyruvate kinase during the intraerythrocytic-stage of its developmental cycle that may play important metabolic roles during infection.