Skip to main content


Data mining of the transcriptome of Plasmodium falciparum: the pentose phosphate pathway and ancillary processes


The general paradigm that emerges from the analysis of the transcriptome of the malaria parasite Plasmodium falciparum is that the expression clusters of genes that code for enzymes engaged in the same cellular function is coordinated. Here the consistency of this perception is examined by analysing specific pathways that metabolically-linked. The pentose phosphate pathway (PPP) is a fundamental element of cell biochemistry since it is the major pathway for the recycling of NADP+ to NADPH and for the production of ribose-5-phosphate that is needed for the synthesis of nucleotides. The function of PPP depends on the synthesis of NADP+ and thiamine pyrophosphate, a co-enzyme of the PPP enzyme transketolase. In this essay, the transcription of gene coding for enzymes involved in the PPP, thiamine and NAD(P)+ syntheses are analysed. The genes coding for two essential enzymes in these pathways, transaldolase and NAD+ kinase could not be found in the genome of P. falciparum. It is found that the transcription of the genes of each pathway is not always coordinated and there is usually a gene whose transcription sets the latest time for the full deployment of the pathway's activity. The activity of PPP seems to involve only the oxidative arm of PPP that is geared for maximal NADP+ reduction and ribose-5-phosphate production during the early stages of parasite development. The synthesis of thiamine diphosphate is predicted to occur much later than the expression of transketolase. Later in the parasite cycle, the non-oxidative arm of PPP that can use fructose-6-phosphate and glyceraldehyde-3-phosphate supplied by glycolysis, becomes fully deployed allowing to maximize the production of ribose-5-phosphate. These discrepancies require direct biochemical investigations to test the activities of the various enzymes in the developing parasite. Notably, several transcripts of PPP enzyme-coding genes display biphasic pattern of transcription unlike most transcripts that peak only once during the parasite cycle. The physiological meaning of this pattern requires further investigation.


The analysis of the transcriptome of Plasmodium falciparum has revealed that during the intraerythrocytic development of the parasite, genes coding for enzymes and proteins that are involved in complex cellular functions such as transcription, replication or energy metabolism, each involving many gene products, are transcribed in a coordinated fashion, supporting the notion that all components must be present at the right time to allow for optimal function [13]. While this may be true in general, it has already been found that scrutinizing the details of specific metabolic functions reveal some significant departures from this paradigm [4]. Such scrutiny has also provided some intriguing peculiarities that through detailed biochemical studies may reveal some parasite-specific functions that may enlighten our understanding of parasitism or even provide for novel targets for chemotherapeutic intervention. The present analysis explores of the transcriptome to fathom additional metabolic pathways.

During the erythrocytic stage the malaria parasite is engaged in intensive synthesis of nucleotides and is subjected to endogenously produced oxidative radicals that must be detoxified. Like other cells, in order to perform its anabolism, the parasite needs not only energy (ATP): it also needs reducing power, under the form of NADPH. Enzymes that function primarily in the reductive direction utilize NADP+/NADPH pair as co-factors as opposed to oxidative enzymes that utilize the NAD+/NADH cofactor pair. The conversion of ribonucleotides to deoxyribonucleotides (through the action of ribonucleotide reductase) requires NADPH as the electron source. Thus, any cell that proliferates rapidly requires large quantities of NADPH. NADPH can be produced during glucose-6-phosphate oxidation through the pentose-phosphate pathway (PPP; Figure 1). This pathway also produces ribose-5-phosphate (R5P), the sugar component of nucleic acids.

Figure 1

Pentose phosphate pathway. The no-entry symbol indicates that the gene coding for this enzyme could not be found in the genome of P. falciparum, but biochemical evidence suggests that the enzyme activity must be present. P represents phosphate.

The reactions of the PPP operate exclusively in the cytoplasm. PPP has both an oxidative and a non-oxidative arm. The oxidation steps, utilizing glucose-6-phosphate (G6P) as the substrate, occur at the beginning of the pathway and are the reactions that generate NADPH. Thus, the first carbon of glucose-6-phosphate is first oxidized to a lactone (catalyzed by glucose-6-phosphate dehydrogenase) concomitantly releasing two electrons that reduce one molecule of NADP+ to NADPH. The ensuing decarboxylation of 6-phospho-D-gluconate (catalyzed by 6-phosphogluconate dehydrogenase) releases two additional electrons, which reduce a second molecule of NADP+. A five-carbon sugar, D-ribulose-5-phosphate, is produced in the reaction. By isomerization, D-ribulose-5-phosphate is transformed into D-ribose-5-phosphate (R5P). To be used in nucleic acid synthesis, R5P is transformed into 5-Phospho-?-D-ribose 1-pyrophosphoric acid (PRPP) by ribose-phosphate diphosphokinase (EC:

The non-oxidative reactions of the PPP are primarily designed to generate R5P. Equally important reactions of the PPP are to convert dietary 5 carbon sugars or D-ribose-1-phosphate generated in the salvage of purines (that can be slowly converted to R5P by phosphoglucomutase; EC: into both 6 (fructose-6-phosphate) and 3 (glyceraldehyde-3-phosphate) carbon sugars which can then be utilized by the pathways of glycolysis. In the first reaction, R5P will accept two carbon atoms from xylulose-5-phosphate (obtained by epimerization of ribulose-5-P), yielding sedoheptulose-7-phosphate and glyceraldehyde-3-phosphate (catalyzed by transketolase). Sedoheptulose-7-phosphate transfers three carbons to glyceraldehyde-3-phosphate (catalyzed by transaldolase), yielding fructose-6-phosphate (F6P) and erythrose-4-phosphate. Erythrose-4-phosphate then accepts two carbon atoms from a second molecule of xylulose-5-phosphate (catalyzed again by transketolase), yielding a second molecule of F6P and a glyceraldehyde-3-P (GAP) molecule. F6P and a GAP can then enter glycolysis and eventually produce ATP. The intermediate erythrose-4-phosphate is a substrate for the shikimate pathway

The non-oxidative part of PPP can also work in the reverse direction utilizing fructose-6-phosphate glyceraldehyde-3-phosphate generated by glycolysis to produce ribose-5-phosphate. In essence, the PPP can operate in different modes (Figure 2). When both R5P and NADPH are needed (Figure 2 mode 1) only the oxidative path will operate: all the ribulose 5-phosphate is isomerized to R5P and the pathway is completed. When more ribose-5-phosphate than NADPH is needed (Figure 2 mode 2), R5P is synthesized by a reversion of the non-oxidative path from F6P and GAP generated by glycolysis. When the cell needs both NADPH and ATP but not R5P (Figure 2 mode 3): ribulose-5-phosphate is converted to F6P and GAP for glycolysis to make ATP and NADH.

Figure 2

Different modes of PPP action. The PPP can function in different modes depending on the needs of the cell. Mode 1: Both ribose-5-phosphate and NADPH needed – predominating reactions are shown. – All the ribulose 5-phosphate is isomerized to ribose 5-phosphate, which is used for the synthesis of PRPP. Mode 2: More ribose-5-phosphate is needed than NADPH Ribose 5-phosphate is synthesized by the non-oxidative arm using fructose-6-phosphate and glyceraldehydes-3-phosphate supplied by glycolysis. Mode 3: The cell needs NADPH and ATP but not ribose-5-phosphate. Ribulose-5-phosphate is converted to fructose-6-phosphate and glyceraldehydes-3-phosphate which are channeled into glycolysis

The genes that code for enzymes participating in ancillary processes that produce NADP+ and thiamine pyrophosphate, that serves as a co-factor for transketolase activity, should be expressed in coordination with the enzymes of the PPP. The details of the mentioned pathways can be also grasped at (;, respectively) and the time-dependent transcription of the genes coding for the different enzymes will be discussed below.

PPP activity in P. falciparum- infected erythrocytes has been measured [57]. The reverse activity of the non-oxidative arm of PPP has also been demonstrated by measuring the incorporation of radiolabel from [1-14C]glucose into nucleotides [6, 8]. The former investigators have shown that 4/5 of the glucose incorporated into parasite nucleic acids comes from the condensation of F6P and GAP in the reverse action of this arm. Atamna et al, reported that infected cells have large increase of PPP activity where 82 % is contributed by the parasite while the host cell's PPP activity is activated some 24-fold as a result of the oxidative stress that the parasite generates and impinges on the host cell [9].

The gene coding for G6PD has been cloned [10] and the biochemical properties of the isolated enzyme have been characterized [11, 12]. Molecular investigations have revealed that G6PD is coded by a hybrid gene that contains also the sequence of 6-phosphogluconolactonase (the second enzyme of PPP) [13, 14]. Activity of 6-phosphogluconate dehydrogenase in Plasmodium- infected erythrocytes has been detected indirectly and alluded to a parasite enzyme, but not characterized [1517]. Ribose-phosphate diphosphokinase of P. falciparum activity has been characterized and levels of its product 5'-phosphoribosyl-pyrophospate (PRPP) were measured in P. falciparum- infectederythrocytes [8, 18]. The levels of PRPP were found to be increased 56-fold in infected cells at the trohozoite stage compared to uninfected erythrocytes.

The parasite contains substantial levels of the pyridine nucleotides NAD+ and NADP+ and their reduced forms [19] and jointly with the presence of genes coding for the enzymes necessary for their synthesis in the genome of P. falciparum (but see below the discussion concerning the absence of NAD+ kinase-coding gene), indicate the ability of the parasite to produce these co-factors. NADP+/NADPH are used by several biochemical reactions in the parasite (Table 1), but the major role of NADPH is probably in the antioxidant defense of the parasite [20, 21] (see for details). Interestingly enough, the biochemical (such as the intermediate products, enzyme activities for each step, timing of production?) details of NAD(P)+ synthesis in the parasite such as the intermediate products, enzyme activities for each step, timing of production, were not investigated.

Table 1 NADP utilizing enzymes Enzymes are arranged by their sequential functional order. They are given by their name, their EC numbers, the gene identification (PfID) in the Plasmodium genome database (PlasmoDB), the time (in hours post invasion (HPI)) in the parasite's developmental cycle when they are maximally transcribed obtained from the IDC database and the metabolic function of the enzyme.

In this in silico analysis, the stage-dependent transcription of genes that code for enzymes that are involved in the PPP activity of the parasite will be analyzed in a functional context.

Materials and methods

The expression data used in this study was obtained from the transcriptome database of the P. falciparum intraerythrocytic developmental cycle as described [2]. This database contains the relative mRNA abundance for every hour of the intraerythrocytic cycle of parasite development based on the 70-mer oligonucleotide microarray [22]. The expression profile of each transcript is represented by an array of ratios between the mRNA level in the time point sample versus a fixed mRNA level in the control RNA pool. Loess smoothed expression profiles used in this study were calculated as described (see [2]). The peak hour induction was calculated from these profiles. This represents more accurate estimate of the timing of maximal expression since the loess smoothing provides and average mRNA abundance values over a time interval, correcting possible minor fluctuations in the raw value profiles. Each expression profile was subsequently normalized by its peak value. It is assumed that for each metabolic pathway there is a stoichiometric relationship between its individual enzymes and therefore the normalized values are more meaningful for functional evaluation. All enzymes discussed in this essay with their EC numbers and gene ID's are shown in Table 2. Amplification of transcription was calculated by dividing the maximum value by the minimum value of the profile of relative mRNA abundance. A fair agreement between the time-dependent transcription data of the DeRisi transcriptome and Winzeler's data has been observed (data not shown).

Table 2 Enzymes of the pentose phosphate pathway and pyridine nucleotide metabolism. Enzymes are grouped into pathways. They are given by their name, their EC numbers, the gene identification (PfID) in the Plasmodium genome database (PlasmoDB)), the time (in hours post invasion (HPI)) in the parasite's developmental cycle when they are maximally transcribed obtained from the IDC database and the metabolic function of the enzyme.

Results and Discussion

At the onset of the present analysis, it should be underscored that time-dependent transcription does not always overlaps translation. Hence, transcript levels cannot be directly extrapolated to levels of their translated product. Moreover, transient transcription does not divulge on the stability of the translated products. Indeed, a recent analysis has shown for several genes that while the transcript peaks at the trophozoite stage and declines thereafter, the translated protein continues to accumulate [23]. Significant discrepancies between mRNA and protein abundance in P. falciparum were also reported by LeRoche et al. [24]. This was shown to be due to a delay between the maximum detection of an mRNA transcript and that of its cognate protein. Surely, a protein cannot be produced if its cognate gene is not transcribed. Therefore, the present analysis can at best set a time for possible translation but not for its actual occurrence. Thus, while transcript levels are informative for the concerted action of metabolically related enzymes, post-transcriptional mechanisms for controlling protein levels and protein stability must also be considered.

In this analysis, the transcription of genes coding for enzymes that constitute the PPP will be discussed first followed by those that are involved in the synthesis of NAD(P)+ and of thiamine pyrophosphate.

The pentose phosphate pathway

The pathway is shown in Figure 1 and the time-dependent transcription of the different genes is depicted in Figure 3. Whereas G6PD (EC: and hence its conjoint 6-phosphogluconolactonase (EC: are transcribed at high levels during the early stages of parasite development, the transcription of 6-phosphogluconate dehydrogenase (EC: seems to be biphasic (first peak at 9 hours post invasion (HPI) while the major peak is observed at 26 HPI). The two other enzymes that need to be present for full synthesis of PRPP are ribose phosphate isomerase (EC: and ribose-phosphate diphosphokinase (EC: The first has a small peak at 12 HPI and a major one at 21 HPI. The parasite contains 2 genes that code for ribose-phosphate diphosphokinase (PF13_0143 and PF13_0157). Both peak at 21 HPI. Thus, the oxidative arm seems to be fully activated to meet the requirement of ribonucleotide synthesis according to mode 1 of Figure 2. The genes that code for this process are coordinately transcribed starting immediately after invasion and peak at 12 HPI [2]. Biochemical evidence indicates that ribonucleotide synthesis starts at the mature ring stage [25]. Functioning according to mode 1 (Figure 2) would also supply reducing power to neutralize the toxic reactive oxidative species that are maximally produced during the hemoglobin-rich cytosol of the host cell [9]. It seems that the synthesis of PRPP and its utilization in the synthesis purine and pyrimidine nucleotides serves as a metabolic sink that could accelerate the pace of the oxidative arm. Contrary to this judicious coordination, the transcription of genes coding for enzymes needed for the synthesis of NADP+ that is essential for the operation of the oxidative arm, lags considerably behind (see below). This is a crucial paradox that waits to be resolved.

Figure 3

Time dependent transcription of genes coding for enzymes involved in the pentose phosphate pathway. Glucose-6-phosphate dehydrogenase/6-phosphogluconolactoase – black filled circles; 6-phosphogluconate dehydrogenase – red filled circles; Ribose phosphate isomerase – blue filled circles; Ribulose 5-phosphate epimerase – green filled circles. Transketolase – purple filled circles; Ribose phosphate diphosphokinase (PF13_0157) – gold filled circles and (PF13_0143) – dark green filled circles. Glucose-6-phosphate dehydrogenase and 6-phosphogluconolactoase were grouped together because they constitute a chimerical protein.

The situation of the non-oxidative arm genes is different. The gene that codes for ribose phosphate isomerase (EC: peaks at 21 HPI and this seems to restrict the expression of the non-oxidative arm, since transketolase (EC: is transcribed maximally immediately following invasion. Most importantly, the gene that codes for transaldolase (EC: could not be found in the genome of P. falciparum (or in any other Plasmodium species sequenced so far, or in any other Apicomplexans for that matter) thus precluding the completion of the analysis of the non-oxidative arm transcription. As mentioned above, biochemical evidence indicates that this arm is active in the parasite. If it is indeed fully activated when ribose phosphate isomerase is transcribed (and supported by the transcription of the genes that code for ribose-phosphate diphosphokinase that peak at 21 HPI), it would match the transcription of the deoxyribonucleotide synthesis-related genes that starts at 18 HPI and peaks at 30 HPI [2] as well as those related to NADP+ synthesis. When the synthesis of ribonucleotides and deoxyribonucleotides ebbs off towards the end of the intraerythrocytic life cycle of the parasite, the PPP probably functions according to mode 3 (Figure 2) providing both NADPH and ATP. A feature that emerges from this analysis is that the transcription of the genes that code for enzymes acting in the non-oxidative arm are coordinated as shown for other clusters of genes whose products are functionally related [13]. However, it remains to be seen if the cognate enzymes of these transcripts are similarly coordinated since significant discrepancies between mRNA and protein abundance has been recently reported [24].

NAD(P) biosynthesis

The pathway is shown in Figure 4 and the time-dependent transcription of the different genes is depicted in Figure 5. Nicotine adenine dinucleotide (NAD+) and NAD phosphate (NADP+) play key roles in glycolysis and PPP as well as many other enzymatic reactions (see Tables 3 and 1, respectively). NAD synthesis starts by converting nicotinic acid to nicotinate D-ribonucleotide by nicotinate phosphoribosyltransferase (NAPRT; EC: Nicotinate D-ribonucleotide is then adenylated to deamino-NAD+ in a reaction catalyzed by nicotinate-nucleotide adenylyltransferase (EC: then this intermediate is amidated to NAD+ by NAD synthetase (EC: Nicotinamide can be used for NAD synthesis either by being deamidated to nicotininc acid by nicotinamidase (EC: or by being converted to nicotinamide mononucleotide by nicotinamide phosphoribosyltransferase (NPRT; EC: and then to NAD by NAD pyrophosphorylase (EC: Genes coding for the two latter enzymes could not be found in the genome of P. falciparum, although a measurable activity of NPRT has been detected [26]. These enzymes are present in normal erythrocytes and at physiological conditions the production of NAD+ from nicotinamide seems to be more important than that of NADP+ [27]. Zerez et al. [26] demonstrated 15-fold increase in the levels of NAD+ (NAD+ + NADH) in P. falciparum-infected erythrocytes, as well as 3-fold increase in NAPRT activity. These results indicate that the parasite is capable of NAD synthesis although neither the activity of nicotinate-nucleotide adenylyltransferase or of NAD synthetase were increased in infected cells compared to uninfected erythrocytes. This observation suggests that these enzymes must be active in the parasite since the digestion of host cell cytosol during parasite development should have reduced their activities. Abundant nicotinamidase activity was also detected in infected cells implying that the parasite can synthesize NAD+ from both nicotinic acid and nicotinamide (both present in RPMI-1640 medium used for parasite cultivation). The concentrations of NADP+ in infected cells are 10-fold lower that those of NAD+. A gene coding for NAD+ kinase could not be found in the genome of P. falciparum or in any other Plasmodium species or in other Apicomplexans, but could be detected in Giardia lamblia. No information about the action of NAD kinase in the parasite could be found.

Figure 4

Nicotine nucleotide metabolism. The no-entry symbol indicates that the gene coding for this enzyme could not be found in the genome of P. falciparum, but biochemical evidence suggests that the enzyme activity must be present.

Figure 5

Time dependent transcription of genes coding for enzymes involved in the biosynthesis of NAD(P)+. Nicotinamidase – purple filled circles; Nicotinate phosphoribosyl transferase – green filled circles; Nicotinate-nucleotide adenylyltransferase – red filled circles; NAD+ synthetase – black filled circles; NAD(P) transhydrogenase – blue filled circles.

Table 3 NAD+-utilizing enzymes. Enzymes are arranged by their sequential functional order. They are given by their name, their EC numbers, the gene identification (PfID) in the Plasmodium genome database (PlasmoDB), the time (in hours post invasion (HPI)) in the parasite's developmental cycle when they are maximally transcribed obtained from the IDC database, the metabolic function of the enzyme.

Inspection of the time-dependent transcription of genes coding enzymes that are involved in NAD(P)+ biosynthesis reveal an unusual pattern. Transcriptional profiles of at least three transcripts show two peaks, namely those of nicotinamidase, of nicotinate-nucleotide adenylyltransferase and of NAD(P)+ transhydrogenase, as compared to most other transcripts in the parasite transcriptome that are monophasic. The reason for this pattern is unclear as biochemical data are uanavilable for any physiological intepretation, but it may suggest that the the enzymes are not stable. Assuming that the first smaller peak of nicotinate-nucleotide adenylyltransferase (20 HPI) is sufficient for adequate expression of enzymatic activity, the biosynthesis of NAD(P)+ should peak when NAD+ synthetase is at its peak, i.e., at 28 HPI. Thus, by 28 HPI the transcripts of all enzymes necessary for NAD+ synthesis are fully deployed and somewhat later, synthesis itself is probably fully deployed. This time pattern means that important activities, such as glycolysis (the transcription of genes coding for glycolytic enzymes starts to peak at 9 HPI and starts to decrease at 24 HPI [2] At earlier stages glycolysis and PPP activity probably depend on the pool of NAD(P)+ that were present in the invading merozoite. For the full activation of glycolysis at the trophozoite stage [28], the pool of NAD(P)+ has to be significantly amplified. As can be seen in Table 3, several genes coding for enzymes requiring NAD(H) are transcribed earlier than the complete transcription of genes related to NAD+ synthesis. Such an outstanding exception is that of inosine-5'-monophosphate dehydrogenase, an essential enzyme in the purine metabolic pathway: it is transcribed too early to allow its translated product to be functionally useful. The functional meaning of the second peaks of nicotinamidase and of nicotinate-nucleotide adenylyltransferase is enigmatic since the transcription of all other enzyme-coding genes declines to a minimum when they peak. For being meaningful physiologically, the stability of all other enzymes must be sustained with time. Indeed, it has been recently shown that the levels of some proteins such as methionine adenosyltransferase, ornithine aminotransferase, lactate dehydrogenase, glyceraldehyde 3-phosphate dehydrogenase and enolase, are not only increased following transcription, but further increase after their transcript levels decline [23]. Finally, NAD+ kinase is essential for the production of NADP+. Its absence from the genome is indeed very perplexing because biochemical evidence indicate that the parasite is able to synthesize NADP+ [19] and the presence of many enzymes that need it as a cofactor, implicate that its synthesis is essential for parasite growth. As long as the gene (or some surrogate mecahnism) is not identified, nothing can be said about the coordination of NADP+ biosynthesis and the expression of enzymes that utilize it (Table 1).

Transhydrogenase operates at an important interface between NAD(H) and NADP(H) and between the mitochondrial proton electrochemical gradient ??; [29]. Under regular physiological conditions, the enzyme is a consumer of ??:

NADH + NADP+ + H+out_ NAD+ + NADPH + H+in.

The energy of the gradient can drive the [NADPH][NAD+]/[NADP+][NADH] ratio to values >400. Transhydrogenation can also function in the reverse direction from NADPH to NAD+. This is accompanied by outward proton translocation and formation of ??. In this mode, the enzyme utilizes substrate binding energy for proton pumping. Therefore, in terms of energy transduction, the transhydrogenase works in principle like the ATP synthase complex of mitochondria, the proton ATPase. Given the fact that the parasite genome does not have the full complement of genes coding for the mitochondrial proton ATPase [30], it is tempting to suggest that the transhydrogenase could fulfill such role. The transcription of the gene coding for transhydrogenase shows two very distinct peaks at 20 and 30 HPI, with no detectable transcription between them. It may well be that the second peak is adjusted to the time of elongation and division of the single mitochondrion [31].

Thiamine biosynthesis

The essential product of this pathway, thiamine diphosphate, is a cofactor of many enzymes (Table 4). The pathway is shown in Figure 6 and the time-dependent transcription of the different genes is depicted in Figure 7. Thiamine (vitamin B1) is the precursor of the coenzyme thiamine pyrophosphate that is involved in the action of many decarboxylating enzymes and relevant to the present analysis, in the function of transketolase. P. berghei-infected erythrocytes contain higher levels of thiamine than their uninfected counterparts [32] and thiamine deficiency retards the propagation of this parasite in vivo [33]. The incidence of thiamine deficiency in adults admitted to hospital with malaria in Thailand has been examined [34] and it was observed that in hospitalized thiamine deficiency commonly complicates acute falciparum malaria, particularly in severe infections, and could contribute to dysfunction of the central nervous system. Thiamine could be obtained by the parasite from the extracellular space (like in animal cells) and pyrophosphorylated by thiamine diphosphokinase (EC: However the presence of genes coding for enzymes that participate in the eukaryote biosynthetic pathway of thiamine phosphate, suggest that the parasite is able to synthesize its own precursor. Nevertheless, the absence of a gene coding for thiamine phosphate kinase (EC: must be resolved and direct biochemical demonstration of activity must be performed before a firm statement can be made on this metabolic pathway. It is not unlikely that thiamine phosphate kinase is not needed, as is the situation in yeast [35]. There it is suggested that thiamine phosphate is hydrolyzed to thiamine by a non-specific phosphatase and thiamine is then converted to thiamine pyrophosphate by thiamine diphosphokinase.

Table 4 Thiamine diphosphate utilizing enzymes. Enzymes are arranged by their sequential functional order. They are given by their name, their EC numbers, the gene identification (PfID) in the Plasmodium genome database (PlasmoDB), the time (in hours post invasion (HPI)) in the parasite's developmental cycle when they are maximally transcribed obtained from the IDC database, the metabolic function of the enzyme.
Figure 6

Biosynthesis of thiamine diphosphate. The no entry symbol indicates that the gene coding for this enzyme could not be found in the genome of P. falciparum. No biochemical evidence exists to confirm or refute enzyme activity. It has been added to the scheme because genes coding for enzymes preceding it in the pathway are present in the genome.

Figure 7

Time dependent transcription of genes coding for enzymes involved in the biosynthesis of thiamine diphosphate. Hydroxyethylthiazole kinase – red filled circles; Hydroxymethylpyrimidine kinase/phosphomethylpyrimidine kinase – black filled circles; Thiamine-phosphate diphosphorylase – blue filled circles; Thiamine diphosphokinase – green filled circles. Hydroxymethylpyrimidine kinase and phosphomethylpyrimidine kinase were grouped together because they probably constitute a chimerical protein.

The transcription of all genes coding for enzymes involved in thiamine pyrophosphate seems to be coordinated, single phased and peaking between 20 and 30 HPI. If thiamine is obtained from the host and thiamine pyrophosphate is synthesized by the single step mediated by thiamine diphosphokinase, the peak transcription of this gene at 28 HPI lags by several hours after that of transketolase (17 HPI), but since transketoalse is not the expression time-setter of PPP, this lag does not seem to limit the full activity of PPP. However, if most of the thiamine is generated endogenously by the parasite, the supply of this precursor will peak only at 33 HPI, thus limiting full PPP activity. It is not unlikely that both processes occur in tandem or that at advanced stage of parasite development, the need for thiamine cannot be met anymore by exogenous supply and the endogenous synthesis joins the game.


The analysis of time-dependence transcription of parasite genes concluded that the parasite has evolved a highly specialized mode of transcriptional regulation that produces a continuous cascade of gene expression, beginning with genes corresponding to general cellular processes, such as protein synthesis, and ending with Plasmodium-specific functionalities, such as genes involved in erythrocyte invasion [2]. However, a meticulous analysis shows marked and important deviations from this prototype that reveal a lack of coordinated transcription of genes coding for enzymes of the same metabolic pathway and between pathways. There are three most straightforward explanations for these apparent discrepancies. First, there are additional enzymes facilitating the "missing" activities and their identity was not revealed by the present annotations due to their diverse amino acid sequence. Second, the misaligned transcriptional regulation reflect an intricate interplay of the biosynthetic pathways where delayed production of metabolites in one pathway functions as a rate limiting factor for other pathway, which is otherwise fully deployed. Last but not least, post-transcriptional regulation may also play a role. All theories create an intriguing possibility for further studies. Clearly this effort will be enhanced by substantial progress in proteomics and most importantly, direct biochemical demonstrations of activities of individual enzymes and entire pathways.


  1. 1.

    Ben Mamoun C, Gluzman IY, Hott C, MacMillan SK, Amarakone AS, Anderson DL: Co-ordinated programme of gene expression during asexual intraerythrocytic development of the human malaria parasite Plasmodium falciparum revealed by microarray analysis. Mol Microbiol. 2001, 39: 26-36. 10.1046/j.1365-2958.2001.02222.x.

  2. 2.

    Bozdech Z, Llinas M, Pulliam BL, Wong ED, Zhu J, DeRisi JL: The transcriptome of the intraerythrocytic developmental cycle of Plasmodium falciparum. PLoS Biol. 2003, 1: E5-10.1371/journal.pbio.0000005.

  3. 3.

    Le Roch KG, Zhou Y, Blair PL, Grainger M, Moch JK, Haynes JD, De La Vega P, Holder AA, Batalov S, Carucci DJ, Winzeler EA: Discovery of gene function by expression profiling of the malaria parasite life cycle. Science. 2003, 301: 1503-1508. 10.1126/science.1087025.

  4. 4.

    Bozdech Z, Ginsburg H: Antioxidant defense in Plasmodium falciparum: data mining of the transcriptome. Malar J. 2004, 3: 23-10.1186/1475-2875-3-23.

  5. 5.

    Roth EF, Raventos-Suarez C, Perkins M, Nagel RL: Glutathione stability and oxidative stress in P. falciparum infection in vitro: responses of normal and G6PD deficient cells. Biochem Biophys Res Commun. 1982, 109: 355-362. 10.1016/0006-291X(82)91728-4.

  6. 6.

    Atamna H, Pescarmona G, Ginsburg H: Hexose-monophosphate shunt activity in intact Plasmodium falciparum infected erythrocytes and in free parasites. Mol Biochem Parasitol. 1994, 67: 79-89. 10.1016/0166-6851(94)90098-1.

  7. 7.

    Cappadoro M, Giribaldi G, Obrien E, Turrini F, Mannu F, Ulliers D, Simula G, Luzzatto L, Arese P: Early phagocytosis of glucose-6-phosphate dehydrogenase (G6PD)- deficient erythrocytes parasitized by Plasmodium falciparum may explain malaria protection in G6PD deficiency. Blood. 1998, 92: 2527-2534.

  8. 8.

    Roth E, Ruprecht R, Schulman S, Vanderberg J, Olson J: Ribose metabolism and nucleic acid synthesis in normal and glucose-6-phosphate dehydrogenase-deficient human erythrocytes infected with Plasmodium falciparum. J Clin Invest. 1986, 77: 1129-1135.

  9. 9.

    Atamna H, Ginsburg H: Origin of reactive oxygen species in erythrocytes infected with Plasmodium falciparum. Mol Biochem Parasitol. 1993, 61: 231-241. 10.1016/0166-6851(93)90069-A.

  10. 10.

    O'Brien E, Kurdi-Haidar B, Wanachiwanawin W, Carvajal JL, Vulliamy TJ, Cappadoro M, Mason PJ, Luzzatto L: Cloning of the glucose 6-phosphate dehydrogenase gene from Plasmodium falciparum. Mol Biochem Parasitol. 1994, 64: 313-326. 10.1016/0166-6851(94)00028-X.

  11. 11.

    Yoshida A, Roth E: Glucose-6-phosphate dehydrogenase of malaria parasite Plasmodium falciparum. Blood. 1987, 69: 1528-1530.

  12. 12.

    Kurdi-Haidar B, Luzzatto L: Expression and characterization of glucose-6-phosphate dehydrogenase of Plasmodium falciparum. Mol Biochem Parasitol. 1990, 41: 83-92. 10.1016/0166-6851(90)90099-8.

  13. 13.

    Sodeinde O, Clarke JL, Vulliamy TJ, Luzzatto L, Mason PJ: Expression of Plasmodium falciparum G6PD-6PGL in laboratory parasites and in patient isolates in G6PD- deficient and normal Nigerian children. Br J Haematol. 2003, 122: 662-668. 10.1046/j.1365-2141.2003.04397.x.

  14. 14.

    Clarke JL, Scopes DA, Sodeinde O, Mason PJ: Glucose-6-phosphate dehydrogenase-6-phosphogluconolactonase – A novel bifunctional enzyme in malaria parasites. Eur J Biochem. 2001, 268: 2013-2019. 10.1046/j.1432-1327.2001.02078.x.

  15. 15.

    Theakston R, Ali S, Moore G: Electron microscope autoradiographic studies on the effect of chloroquine on the uptake of tritiated nucleosides and methionine by Plasmodium berghei. Ann Trop Med Parasitol. 1972, 66: 295-302.

  16. 16.

    Theakston R, Fletcher K: An electron cytochemical study of 6-phosphogluconate dehydrogenase activity in infected erythrocytes during malaria. Life Sci. 1973, 13: 405-410. 10.1016/0024-3205(73)90233-6.

  17. 17.

    Fletcher KA, Canning MV, Theakston RD: Electrophoresis of glucose-6-phosphate and 6-phosphogluconate dehydrogenases in erythrocytes from malaria-infected animals. Ann Trop Med Parasitol. 1977, 71: 125-130.

  18. 18.

    Reyes P, Rathod PK, Sanchez DJ, Mrema JE, Rieckmann KH, Heidrich HG: Enzymes of purine and pyrimidine metabolism from the human malaria parasite, Plasmodium falciparum. Mol Biochem Parasitol. 1982, 5: 275-290. 10.1016/0166-6851(82)90035-4.

  19. 19.

    Sherman IW: Biochemistry of plasmodium (malarial parasites). Microbiol Rev. 1979, 43: 453-495.

  20. 20.

    Kanzok SM, Schirmer RH, Turbachova I, Iozef R, Becker K: The thioredoxin system of the malaria parasite Plasmodium falciparum – Glutathione reduction revisited. J Biol Chem. 2000, 275: 40180-40186. 10.1074/jbc.M007633200.

  21. 21.

    Becker K, Rahlfs S, Nickel C, Schirmer RH: Glutathione – functions and metabolism in the malarial parasite Plasmodium falciparum. Biol Chem. 2003, 384: 551-566. 10.1515/BC.2003.063.

  22. 22.

    Bozdech Z, Zhu J, Joachimiak MP, Cohen FE, Pulliam B, DeRisi JL: Expression profiling of the schizont and trophozoite stages of Plasmodium falciparum with a long-oligonucleotide microarray. Genome Biol. 2003, 4: R9-10.1186/gb-2003-4-2-r9.

  23. 23.

    Nirmalan N, Sims PF, Hyde JE: Quantitative proteomics of the human malaria parasite Plasmodium falciparum and its application to studies of development and inhibition. Mol Microbiol. 2004, 52: 1187-1199. 10.1111/j.1365-2958.2004.04049.x.

  24. 24.

    Le Roch KG, Johnson JR, Florens L, Zhou Y, Santrosyan A, Grainger M, Yan SF, Williamson KC, Holder AA, Carucci DJ, Yates JR, Winzeler EA: Global analysis of transcript and protein levels across the Plasmodium falciparum life cycle. Genome Res. 2004, 14: 2308-2318. 10.1101/gr.2523904.

  25. 25.

    Scott H, Gero AM, O'Sullivan WJ: In vitro inhibition of Plasmodium falciparum by pyrazofurin, an inhibitor of pyrimidine biosynthesis de novo. Mol Biochem Parasitol. 1986, 18: 3-15. 10.1016/0166-6851(86)90045-9.

  26. 26.

    Zerez C, Roth E, Schulman S, Tanaka K: Increased nicotinamide adenine dinucleotide content and synthesis in Plasmodium falciparum- infected human erythrocytes. Blood. 1990, 75: 1705-1710.

  27. 27.

    Micheli V, Simmonds HA, Sestini S, Ricci C: Importance of nicotinamide as an NAD precursor in the human erythrocyte. Arch Biochem Biophys. 1990, 283: 40-45. 10.1016/0003-9861(90)90609-3.

  28. 28.

    Pfaller M, Krogstad D, Parquette A, Nguyen-Dinh P: Plasmodium falciparum: Stage-specific lactate production in synchronized cultures. Exp Parasitol. 1982, 54: 391-396. 10.1016/0014-4894(82)90048-0.

  29. 29.

    Hatefi Y, Yamaguchi M: Nicotinamide nucleotide transhydrogenase: a model for utilization of substrate binding energy for proton translocation. FASEB J. 1996, 10: 444-452.

  30. 30.

    Krugkrai J: The multiple roles of the mitochondrion of the malarial parasite. Parasitology. 2004, 129: 511-524. 10.1017/S0031182004005888.

  31. 31.

    Divo AA, Geary TG, Jensen JB, Ginsburg H: The mitochondrion of Plasmodium falciparum visualized by rhodamine 123 fluorescence. J Protozool. 1985, 32: 442-446.

  32. 32.

    Singer I: Tissue thiamine changes in rats with experimental trypanosomiasis or malaria. Exp Parasitol. 1961, 11: 391-401. 10.1016/0014-4894(61)90043-1.

  33. 33.

    Singer I: Coenzyme A changes in liver, spleen and kidneys of rats with infections of Plasmodium berghei. Proc Soc Exp Biol Med. 1956, 91: 315-318.

  34. 34.

    Krishna S, Taylor AM, Supanaranond W, Pukrittayakamee S, ter Kuile F, Tawfiq KM, Holloway PA, White NJ: Thiamine deficiency and malaria in adults from southeast Asia. Lancet. 1999, 353: 546-549. 10.1016/S0140-6736(98)06316-8.

  35. 35.

    Hohmann S, Meacock PA: Thiamin metabolism and thiamin diphosphate-dependent enzymes in the yeast Saccharomyces cerevisiae: genetic regulation. Biochim Biophys Acta. 1998, 1385: 201-219.

Download references


We thank Professor I. Ohad for critical reading of the manuscript.

Author information

Correspondence to Hagai Ginsburg.

Additional information

Authors' contributions

Each author contributed equally to this investigation.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Authors’ original file for figure 4

Authors’ original file for figure 5

Authors’ original file for figure 6

Authors’ original file for figure 7

Rights and permissions

Reprints and Permissions

About this article


  • Thiamine
  • Pentose Phosphate Pathway
  • Thiamine Deficiency
  • Parasite Development
  • Thiamine Pyrophosphate


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.