malaria , which is caused by a protozoan parasite of the plasmodium genus , represents a serious global health concern given that nearly half of the world population is at risk of infection ( who , 2014 ) .
although several anti - malarial drugs for treating the symptomatic stage ( or blood stage ) of malarial infection are commercially available ( antony and parija , 2016 ) , their efficacy has declined appreciably in the last few decades owing to widespread drug resistance developed by the parasite ( breman et  al . ,
chloroquine was the first - line malaria treatment for many decades until drug - resistant p.  falciparum strains became common .
the drug causes a dose - dependent decrease in hemozoin formation ( chou and fitch , 1992 , slater and cerami , 1992 ) and an associated increase in toxic free heme in the food vacuole of the parasite ( combrinck et  al .
, 2013 , loria et  al . , 1999 ) . over the past few decades
, researchers have proposed many different mechanisms for chloroquine action , including 1 ) dna intercalation ( meshnick , 1990 ) , 2 ) alteration of digestive food vacuole ph ( yayon et  al . , 1985 ) ,
3 ) inhibition of heme polymerase ( yayon et  al . , 1985 ) , and 4 ) formation of a toxic chloroquine - ferriprotoporphyrin ix complex ( sugioka et  al . ,
yet , there is no consensus on the exact mechanisms by which chloroquine kills the parasite and under what circumstances they operate .
shedding light on the possible actions of chloroquine is critical for developing more potent drugs or combination drug treatments against the parasite and for understanding how the parasite can develop resistance against them .
here we used systems biology model and transcriptional profiles of p.  falciparum to obtain information about metabolic and biological precursors that determine the physiological or pathophysiological state of the parasite .
specifically , we used the p.  falciparum transcriptomic data obtained by hu and colleagues during the intraerythrocytic development cycle ( idc ) at different concentrations of chloroquine ( hu et  al . , 2010 ) .
our primary goal was to use these data to investigate the effect of chloroquine on p.  falciparum dna replication .
we chose to study this effect because it has been proposed as a mechanism of chloroquine action in different plasmodium species ( gutteridge et  al .
, 1972 , meshnick , 1990 , polet and barr , 1968a , polet and barr , 1968b ) , in bacteria ( ciak and hahn , 1966 ) , and in bioassays examining the effect of chloroquine on rna and dna polymerases ( cohen and yielding , 1965 , o'brien et  al .
we first developed a method that integrates the information embedded in the transcriptome data with the whole - genome metabolic network model for p.  falciparum ( fang et  al . , 2014 ) to predict the metabolic phenotype corresponding to those transcription data .
we validated the developed method by using a data set from an independent study that jointly capture transcriptomic and metabolomic data from p.  falciparum , to correctly predict the metabolic phenotype tied to genetic perturbation of two krebs cycle enzymes ( ke et  al . , 2015 ) .
we then focused on the inhibition of dna replication in p.  falciparum ( as a consequence of chloroquine treatment ) a critical metabolic phenotype ( ciak and hahn , 1966 , cohen and yielding , 1965 , polet and barr , 1968a ) , and used our method to simulate the effect of chloroquine on the metabolic fluxes of p.  falciparum during the idc .
we identified genes that were substantially altered in response to chloroquine treatment and linked to the inhibition of p.  falciparum dna replication .
the cohort of identified genes suggests that dna replication is inhibited by the downstream effect of heme accumulation .
specifically , our analysis suggests that continuous accumulation of heme inhibits redox metabolism , carbon fixation , and pyrimidine metabolism , which leads to inhibition of dna replication and facilitates death of the parasite .
our results provide a mechanistic explanation for why parasites with an efficient redox metabolism may have a lower sensitivity to chloroquine ( kasozi et  al . , 2013 , lehane et  al . , 2012 ,
metabolic flux profiles of p.  falciparum strains 3d7 ( pf3d7 ) and dd2 ( pfdd2 ) during the idc were estimated by using transcriptome data previously reported by llinas and colleagues ( llinas et  al . , 2006 ) , who extracted cdna from parasites cultured in erythrocytes and hybridized them to dna microarrays to yield hourly levels of gene expression for the pf3d7 ( or pfdd2 ) strain during the idc .
the temporal resolution of these hourly gene expression profiles allowed us to make high - resolution time - dependent metabolic flux predictions during the idc .
we used previously reported gene expression data ( hu et  al . , 2010 ) to simulate the effect of chloroquine on the time - dependent metabolic fluxes of pf3d7 and pfdd2 during the idc .
briefly , hu and colleagues ( hu et  al . , 2010 ) administered chloroquine 18  h after erythrocyte invasion at concentrations of 41 , 72 , and 144  nm for pf3d7 and at 43  nm for pfdd2 , and obtained gene expression profiles every hour for 8  h during the idc following the chloroquine treatment .
the model was identical to a previously published model ( fang et  al . , 2014 ) , with one minor modification : we added reactions explicitly carried out by pdx2 ( pf11_0169 ) and pdx1 ( mal6p1.215 ) .
this was done by 1 ) modifying an old reaction that incorrectly assumed that pyridoxal 5-phosphate ( plp ) was synthesized from ribulose 5-phosphate and 2 ) adding a new ( and separate ) reaction for pdx1 because experiments suggest that it can synthesize plp independent of the ammonia generated by pdx2 ( hanes et  al . ,
2008 ) . in the updated model , pdx2 hydrolyzes l - glutamine to yield glutamate and ammonia while pdx1 incorporates the ammonia generated by pdx2 to yield plp from the substrates ribose 5-phosphate and glyceraldehyde 3-phosphate ( gengenbacher et  al . ,
2006 ) , with the isomerization of ribose 5-phosphate to ribulose 5-phosphate as a side reaction .
we used gene - to - reaction mappings available in the metabolic network model to obtain reaction expressions rexp , which depend on the expression of a gene ( or genes ) catalyzing a particular metabolic reaction . for a reaction catalyzed by a single enzyme , rexp was equal to the expression of the gene catalyzing that reaction
however , if a reaction was associated with more than one gene , then their gene - to - reaction mappings were defined by using the boolean operators  and  and  or .
 we implemented these operations by taking the minimal ( for the and operator ) and maximal ( for the or operator ) values of the associated gene expression data ( song et  al . , 2014 ) .
, 2014 ) that integrates hourly gene expression data with a whole - genome metabolic network model to estimate the hourly metabolic flux profiles of p.  falciparum during the 48-h long idc period . to take advantage of additional studies that collect gene expression data as a function of genetic or chemical perturbations
, we modified the method because such studies usually do not provide the same granularity and density of gene expression data to fully cover the entire idc .
therefore , we aimed to develop a method that could utilize a reduced transcriptomic data set to predict alterations in metabolic flux profiles during the idc .
briefly , we used the following steps to estimate metabolic fluxes under no - stress conditions : step 1 : estimate nutrient fluxes sufficient for the nominal growth rate ( 40% of the maximum growth rate).step 2 : estimate reaction fluxes vin for metabolite i corresponding to the nominal growth rate constrained by the nutrient fluxes obtained in step 1.step 3 : estimate reaction fluxes vit for reaction i at any time point t during the idc by incorporating time - dependent transcriptomic profiles under stoichiometric and nutrient uptake constraints .
step 1 : estimate nutrient fluxes sufficient for the nominal growth rate ( 40% of the maximum growth rate ) .
step 2 : estimate reaction fluxes vin for metabolite i corresponding to the nominal growth rate constrained by the nutrient fluxes obtained in step 1 .
step 3 : estimate reaction fluxes vit for reaction i at any time point t during the idc by incorporating time - dependent transcriptomic profiles under stoichiometric and nutrient uptake constraints .
the last step was achieved by globally minimizing for all reactions and time points the value of t that minimized the weighted absolute difference , ig|vitritvin| , where g represents the set of unidirectional reactions with known gene - to - reaction mapping .
the factor rit = rexp , itmin(rexp , it)mean(rexp , it)min(rexp , it ) ; where , rexp , it represents reaction expression of ith reaction at time t obtained from the gene - to - reaction mapping of the p.  falciparum metabolic network , according to the approach outlined in section 2.2 .
this factor , rit , ensured that the mean of the ith reaction during the idc was equal to vin .
, we estimated the organism - level impact of a drug on metabolism by using transcriptomic data obtained after drug treatment .
the original formulation outlined above ( steps 13 ) can not be directly translated to stress conditions unless the complete transcriptomic profile for each stress condition during the entire idc is available ( fang et  al . , 2014 ) .
therefore , we introduced a parameter i to account for up- or down - regulation of individual reactions in response to stress to simulate the effect of exogenous stress on metabolic flux profiles during the idc .
we modified the factor rit to iritmin(irit)mean(rit)min(irit ) , where i denotes the median relative change in the reaction expression of the ith reaction following chloroquine treatment , as observed in previous experiments ( hu et  al . , 2010 ) .
note that the mean of each individual modifying factor rit , i.e. , 1nj=1tnrij , obtained by summing the values across all time points in the idc , tracks with i ; that is , when the mean is less than one , greater than one , or equal to one , i<1 , i>1 , or i=1 , respectively . in other words , i increases or decreases the nominal flux of the ith reaction during the entire idc depending on the relative increase or decrease in the ith reaction expression . in principle ,
a drug may differentially affect gene transcription at different time points during the idc . here
, we assumed that the median change in gene transcription after a stress treatment was representative of the transcriptional response throughout the idc .
we used a previously published data set ( ke et  al . , 2015 ) to validate our methodology to simulate the effect of external stress on metabolic fluxes during the idc .
ke and colleagues ( ke et  al . , 2015 ) constructed a p.  falciparum mutant with two tricarboxylic acid ( tca ) cycle enzymes deleted ( idh/kdh ) , cultured wild - type and mutant parasites in erythrocytes , and obtained matched transcriptomic and metabolomic data .
we used this data set because it describes changes in gene expression ( transcriptomics ) and phenotype ( metabolomics ) under the same experimental stress conditions . if the reaction fluxes estimated from previously reported transcriptome data ( ke et  al . ,
2015 ) adequately represent the reaction fluxes in  vitro , then the predicted changes in reaction fluxes should be sufficient to predict the changes observed in the corresponding metabolomic data .
we used a computer model of mitochondrial oxidative phosphorylation and the tca cycle ( wu et  al . , 2007 ) to predict changes in tca cycle metabolites in response to changes in these estimated reaction fluxes ( ke et  al .
the kinetic parameters of this computer model were determined and validated by using extensive experimental data .
we simulated wild - type reactions by running this model , with its default parameters , for 2000  s ( time required to reach steady - state condition ) and stored the concentrations of tca cycle metabolites .
we performed simulations for the mutant parasite in the same way , albeit after scaling the model parameters that govern tca cycle enzyme activities by the estimated increase or decrease in reaction fluxes of those enzymes .
we used a freely available software package called bisen ( vanlier et  al . , 2009 ) to implement the tca cycle computer model ( wu et  al . , 2007 ) .
p.  falciparum takes up nutrients to synthesize dna , rna , proteins , and phospholipids .
these processes are not constitutively activated , because the parasite will synthesize them as needed throughout the idc .
we used a validated in - house developed method ( fang et  al . , 2014 ) to simulate time - dependent dna synthesis during the idc . for comprehensiveness
, we introduce the constraint governing the rate of dna synthesis in the metabolic network model : cdatpdatp+cdctpdctp+cdgtpdgtp+cdttpdttpdna+(cdatp+cdctp+cdgtp+cdttp)diphosphate . here ,
datp , dctp , dgtp , and dttp refer to deoxyadenosine , deoxycytidine , deoxyguanosine , and deoxythymidine triphosphates , respectively .
we refer the reader to the original article ( fang et  al . , 2014 ) for further details on this method for simulating time - dependent p.  falciparum macro - molecules during the idc .
the criticality of a gene is determined by optimizing the following objective function : idnacimi , subjected to : smxnvnx1=0 and vlbvvub . here , dna is the gene set containing the metabolites that contribute to dna synthesis ( introduced above ) and ci is the contribution of the ith metabolite ( mi ) .
in addition , s is the stoichiometry matrix , v is a vector containing all metabolic reactions of the p.  falciparum metabolic network model , m is the number of metabolites in the network , n is the number of reactions in the network , and vlb and vub are the upper and lower bounds , respectively , of metabolic reactions . we deemed a reaction as critical if there was no dna synthesis after knocking out the reaction in silico .
chloroquine administration inhibits dna replication in malaria parasites ( cohen and yielding , 1965 , polet and barr , 1968a ) .
2010 ) also show that metabolic genes critical for p.  falciparum dna synthesis are suppressed by chloroquine treatment .
we sought to identify the genes that are substantially altered in response to chloroquine and cause inhibition of p.  falciparum dna synthesis .
obvious candidates are genes that are down - regulated and critical for p.  falciparum dna synthesis .
however , not all such genes will have a pronounced effect on dna synthesis , which depends on the extent to which an individual gene is down - regulated and/or the individual contribution of that gene to the dna synthesis of the parasite .
this suggests that many combinations of such genes may underlie the metabolic response associated with chloroquine treatment .
, 1990 ) , we aimed at identifying the minimal number of genes ( hereinafter referred to as the core genes ) that are sufficient to completely account for the effect of chloroquine on dna synthesis by the parasite .
we identified the core genes based on the bayesian information criteria ( bic ) score ( schwarz , 1978 ) , which is given by , nln(ssen)+kln(n ) , where sse=i=1n(victlvicqn)2max(vctl)min(vctl ) and n denotes the number of time points , k the number of genes being tested , vctl the estimated rate of dna synthesis under control conditions , and vcqn the estimated rate of dna synthesis after chloroquine treatment ( with =1 for the k genes being tested ) .
we implemented an algorithm that starts with a set of size k=1 and increases the size of the test set as long as a given combination of genes for the new set ( size k+1 ) has a bic score lower than the minimal bic score of the previous size .
we exhaustively repeated this process at each chloroquine dose ( 41 , 72 , and 144  nm for pf3d7 , and 43  nm for pfdd2 ) .
we modified our previous approach ( fang et  al . , 2014 ) to simulate the effect of external stress , as depicted in fig .  1 and described in the materials and methods section .
we validated the methodology by demonstrating that reaction fluxes obtained in response to stress were sufficient to predict metabolic phenotypes corresponding to that stress .
2a shows the biochemical pathways incorporated into the computer model used to predict the functional consequences of altered fluxes of tca cycle enzymes in response to previously studied genetic perturbations ( ke et  al . , 2015 ) .
2b shows the predicted alterations in tca cycle metabolites in response to genetic perturbations of tca cycle enzymes ( fig .  2a , red circles ) .
we calculated the altered metabolite concentrations by scaling the control parameters in the tca cycle model ( wu et  al . , 2007 ) by the median of the predicted changes in tca cycle enzyme activities ( fig .  2c ) .
2b , filled bars ) qualitatively captured the trends of relative change in experimentally measured tca cycle metabolites ( fig .  2b , open bars ) .
2b , filled bars ) was not within the margin of error of the experimental measurements ( fig .  2b , open bars ) .
we examined all possible enzymes / transporters ( e.g. , citrate synthase , pyruvate dehydrogenase , aspartate aminotransferase , citrate - malate antiporter ) whose perturbation could yield mitochondrial citrate levels higher than those in the wild - type system . in each scenario , either the predicted citrate levels did not increase or the predicted oxaloacetate levels increased . therefore , based on our model analysis , the source of this increase in the experimentally reported citrate level appears to be non - mitochondrial in nature . in agreement with previous experimental observations ( cohen and yielding , 1965 ,
polet and barr , 1968a ) , our simulations suggested that chloroquine treatment results in dose - dependent inhibition of dna replication ( fig . 
a , blue diamonds ) and 43  nm ( pfdd2 , panel b , blue diamonds ) , the rate of dna synthesis during the ring stage ( the first 18  h after infection ) was less inhibited than it was during the late trophozoite and early schizont phases ( beyond 30  h after infection ) .
these results are in agreement with previous observations ( yayon et  al . , 1983 ) . because many studies suggest that pfdd2 is less sensitive to chloroquine than pf3d7 ( moneriz et  al . , 2011 )
, dna synthesis might have been expected to be more inhibited in pf3d7 than in pfdd2 . in our simulations , however , a similar concentration of chloroquine resulted in greater dna synthesis inhibition in pfdd2 ( fig . 
3b ; blue diamonds versus black circles , paired - sample t - test ; p  < 
0.0001 ) than in pf3d7 ( fig .  3a ; blue diamonds versus black circles , paired - sample t - test ; p  =  0.15 ) .
these differences likely stem from differences in the relative gene expression ratios of pfdd2 and pf3d7 .
notably , the transcriptional profile for pfdd2 was more down - regulated than that for pf3d7 at a similar concentration of chloroquine . in other words , the percentage of down - regulated genes critical for p.  falciparum dna synthesis was greater in pfdd2 than in pf3d7 ( table  s1 ) . together with the highly similar ic50 values reported previously ( hu et  al . , 2010 ) for pf3d7 ( 41  nm ) and
pfdd2 ( 43  nm ) , these results suggest that the pfdd2 strain used here may have diverged from its true ( chloroquine - resistant ) phenotype .
we found that 81 metabolic reactions in our network model are critical for p.  falciparum dna synthesis .
we denoted these reactions as critical because their elimination resulted in abolished dna synthesis ( see section 2.6 ) .
of these 81 reactions , 64 had a known gene association whereas 17 did not .
these genes belonged to pathways involved in glycolysis , pentose phosphate synthesis , pyruvate metabolism , purine metabolism , pyrimidine metabolism , redox metabolism , and mitochondrial metabolism . in the online supplement ( see spreadsheet ) , we provide detailed information , including the reaction name , chemical reaction , gene association , and pathway , of all reactions critical for p.  falciparum dna synthesis .
experimental observations ( cohen and yielding , 1965 , polet and barr , 1968a ) and the present simulations ( fig . 
3 ) indicate that chloroquine treatment inhibits replication of p.  falciparum dna , which in turn suggests that pathways essential for dna synthesis are inhibited .
4 shows the relative change in the amount of dna synthesized by the parasite during the idc period in response to chloroquine treatment .
table  s2 lists the actual amounts of synthesized dna for both strains under each simulated condition .
we observed a dose - dependent decrease in the amount of dna synthesized by the parasite during the idc ( fig .  4 , red markers ) .
we used bic scores ( schwarz , 1978 ) and employed the algorithm described in section 2.7 to identify core genes that caused this net decrease of synthesized dna . preventing only the core genes from being down - regulated by chloroquine was sufficient to restore the amount of dna synthesized to control values ( fig .  4 , green markers ) ;
for these simulations , the reaction expressions of the core genes were set to be equal for control and chloroquine - treated conditions .
a list of the core genes identified for each strain and chloroquine dose ( table  s3 ) suggests that thioredoxin reductase and ribonucleotide reductase substantially affect the amount of synthesized dna .
for pfdd2 , only two genes ( thioredoxin reductase and ribonucleotide reductase ) caused the nearly 50% reduction in the net amount of synthesized dna ( table  s3 and fig .  4 ) .
these genes were not mentioned in the original transcriptome analysis of hu and colleagues , who used a fundamentally different approach to interpret their transcriptome data ( hu et  al . , 2010 ) .
these authors focused on genes that were either co - regulated or showed strongly differential expression ( 8 log2-fold change ) in response to a drug . in contrast , our methodology of coupling transcriptomic changes with a metabolic network model can predict physiological outcomes as they relate to parasite metabolism , regardless of the magnitude of transcriptional change for an individual gene .
transcriptional changes are not proportional to changes in metabolic reaction fluxes ( moxley et  al . , 2009 ) .
our analysis suggests that low - amplitude alterations in these identified genes ( table  s3 ) were sufficient to cause a significant reduction in parasite dna synthesis .
several existing methods integrate whole - genome metabolic networks and transcriptome data ( colijn et  al . , 2009 , song et  al . , 2014 ,
recently , we developed a method that predicts metabolic flux profiles of p.  falciparum during the idc ( fang et  al .
the transcriptome data of p.  falciparum obtained under different stress conditions ( e.g. , drug exposure , genetic mutation , etc . ) contain information about its functional state and can be used to predict qualitative alterations , if any , in metabolic flux profiles during the idc .
however , such data are often collected at a limited number of time points ; for instance , hu and colleagues ( hu et  al . , 2010 ) obtained p.  falciparum transcriptome data at only 8 selected time points during the idc .
currently , there is no approach to integrate limited data , such as those reported previously ( hu et  al . , 2010 ) , with a whole - genome network to predict metabolic flux profiles during the idc .
this is mainly because such an approach requires gene expression data throughout the idc to predict the corresponding metabolic flux profiles .
therefore , we developed an approach that uses gene expression data obtained under a stress condition to predict alterations of metabolic fluxes that correspond to that condition .
we validated the developed method by integrating transcriptome data obtained under a given stress condition and predicting metabolomic data corresponding to that stress ( ke et  al . ,
, we simulated the effect of chloroquine on the metabolic flux profiles of p.  falciparum during the idc .
chloroquine inhibits the formation of hemozoin and leads to heme accumulation , which is toxic to the parasite ( fitch , 1998 ) .
however , the specific mechanism underlying this toxicity is not known ( ginsburg et  al . , 1998 ,
some experiments suggest inhibition of dna replication as a potential mechanism of chloroquine action ( cohen and yielding , 1965 , kwakye - berko and meshnick , 1989 , polet and barr , 1968a ) .
we developed a method that successfully predicted a reduction in the ability of the parasite to synthesize dna in response to chloroquine treatment ( see fig . 
we identified twelve distinct genes that were substantially affected in pf3d7 and pfdd2 after chloroquine treatment ; importantly , preventing their down - regulation in response to chloroquine treatment was sufficient for restoring the p.  falciparum dna synthesis inhibited by chloroquine treatment ( see table  1 ) .
these genes belonged to pathways involved in glycolysis , redox metabolism , and pyrimidine metabolism .
a closer look at these genes reveals a potential mechanism of chloroquine action that unifies previous theories and observations involving heme ( fitch , 1998 ) , dna replication ( polet and barr , 1968a ) , and oxidative stress ( lehane et  al . , 2012 , radfar et  al . ,
5 , suggests that administration of chloroquine has two major metabolic effects : i.it interferes with the function of glyceraldehyde 3-phosphate dehydrogenase ( gapdh ) , which reduces the synthesis of phosphoenolpyruvate , whose carbon fixation activity is essential for anaplerotic tca metabolism and redox balance ( storm et  al .
, 2014).ii.its interference of thioredoxin reductase ( trxr ) reduces conversion of oxidized thioredoxin trx ( s2 ) to reduced thioredoxin trx ( sh2 ) , which is essential for donating electrons to reactions catalyzed by thioredoxin peroxidase ( tpx ) and ribonucleotide reductase ( rnr ) .
tpx is essential for converting h2o2 ( reactive ) to h2o ( harmless ) ( pannala and dash , 2015 ) , and rnr for converting ribonucleoside 5-diphosphates species ( adp , udp , gdp and cdp ) , to deoxy - species ( dadp , dudp , dgdp and dcdp ) which are used for dna synthesis ( munro and silva , 2012 ) .
it interferes with the function of glyceraldehyde 3-phosphate dehydrogenase ( gapdh ) , which reduces the synthesis of phosphoenolpyruvate , whose carbon fixation activity is essential for anaplerotic tca metabolism and redox balance ( storm et  al . , 2014 ) .
its interference of thioredoxin reductase ( trxr ) reduces conversion of oxidized thioredoxin trx ( s2 ) to reduced thioredoxin trx ( sh2 ) , which is essential for donating electrons to reactions catalyzed by thioredoxin peroxidase ( tpx ) and ribonucleotide reductase ( rnr ) .
tpx is essential for converting h2o2 ( reactive ) to h2o ( harmless ) ( pannala and dash , 2015 ) , and rnr for converting ribonucleoside 5-diphosphates species ( adp , udp , gdp and cdp ) , to deoxy - species ( dadp , dudp , dgdp and dcdp ) which are used for dna synthesis ( munro and silva , 2012 ) .
for example , our analysis points to trxr and gapdh as candidate molecules that determine the efficacy of chloroquine . here , we focused on candidates susceptible to a low dose of chloroquine because several pathways will be inhibited at higher concentrations ( see table  s3 ; chloroquine  =  144  nm ) .
in addition , proteins that transport chloroquine out of food vacuoles should decrease the effective concentration of chloroquine . targeting these pathways while administering chloroquine treatment should enhance the efficacy of chloroquine at a given effective concentration .
the hypothesis that chloroquine kills the malarial parasite by dna intercalation is decades old ( meshnick , 1990 ) .
although it was the prime candidate mechanism of chloroquine from the 1950s to the 1970s , it has now largely been abandoned for the past two decades .
our model analyses do not identify a drug target or targets ; however , they suggest that the effect of chloroquine administration leads to inhibition of dna synthesis . according to our model
, the parasite responds to chloroquine by reducing the expression of genes critical for dna synthesis .
although chloroquine could specifically intercalate with the identified core gene set in accord with the dna intercalation theory of chloroquine ( meshnick , 1990 , o'brien et  al . , 1966a
heme is a more plausible molecular inhibitor , given that it accumulates in large amounts within the parasite 's food vacuole and is capable of causing dna damage ( ishikawa et  al .
, 2010 , rahman et  al . , 1997 ) , protein oxidation ( kitatsuji et  al . ,
2016 ) , nadph oxidation ( bodaness , 1983 ) , and apoptosis ( chiabrando et  al . , 2014 ) .
however , there is no direct evidence suggesting that heme causes any of these damages in p.  falciparum .
chloroquine , which was once a highly effective drug , is now abandoned in many malaria endemic regions because of widespread resistance .
knowledge of this mechanism allowed us to identify metabolic pathways that could be exploited to increase chloroquine efficacy .
for example , our model analyses suggest that administration of h2o2 would potentiate parasite killing by chloroquine , consistent with in  vitro experiments ( malhotra et  al . ,
similarly , our analyses suggest the hypothesis that targeting parasite - specific pathways that synthesize phosphoenolpyruvate and reduce thioredoxin would potentiate parasite killing . targeting these pathways ,
while administering amodiaquine ( a chloroquine analogue ) , should result in higher efficacy than administering chloroquine , especially because amodiaquine accumulates more than chloroquine in the parasite s food vacuole and has a mechanism similar to that of chloroquine ( hawley et  al . , 1996 ) .
our aim was to identify and investigate enzymes that contribute to decreased dna replication in response to chloroquine treatment . using a transcriptomics - based approach
, we can identify cellular responses that are specific to chloroquine by expression data taken from chloroquine - treated parasites .
although other antimalarial drugs may activate similar stress pathways , the details of the response will be specific to the particular conditions .
the approach makes no assumptions on the underlying molecular actions that cause or affect the response .
the effect of chloroquine on dna replication is likely to be indirect , being mediated by heme accumulation , which evokes multiple toxic effects in the cell .
we simulated the effect of chloroquine over the entire idc based on the available transcriptome data taken from an eight - hour window during the trophozoite stage .
ideally , for any specific inhibitor , it is desirable to have hourly transcriptional readouts to predict the effect of a systemic perturbation throughout the entire idc . in the case
presented here , the largest effect of chloroquine occurs during the trophozoite stage and the protocol for collecting experimental transcriptomic data was designed to characterize transcriptional changes during this stage .
our use of these data to identify changes in metabolite fluxes throughout the idc relative to untreated controls rests on the assumption that altering the expression of genes not associated with any biological function specific to a particular stage will have minimal physiological effects at other stages .
this is an assumption that is specific to p.  falciparum , because malaria relies on the execution of a highly regimented transcriptional program during the idc , which initiates specific stage - dependent cellular functions .
these assumptions were ultimately corroborated because the observed metabolite phenotype changes are mainly associated with the trophozoite stage and metabolic alterations during previous time points were non - existent or immaterial to dna synthesis .
herein , we investigated the potential mechanism by which chloroquine inhibits dna synthesis in p.  falciparum and kills the parasite .
we first developed an approach ( that requires limited transcriptome data ) to simulate the effect of exogenous stress on p.  falciparum metabolic fluxes during the idc .
we then developed an algorithm employing bayesian information theory to identify potential pathways that inhibit dna synthesis in response to chloroquine treatment .
our analysis suggests that 1 ) chloroquine inhibits dna synthesis via inhibition of redox metabolism , carbon fixation , and anaplerotic tca cycle metabolism ; 2 ) p.  falciparum may maintain dna synthesis after chloroquine treatment by preserving a reduced pool of thioredoxin and by invoking alternate mechanisms for synthesizing phosphoenolpyruvate , and 3 ) a drug or combination of drugs that target redox metabolism , carbon fixation , and heme accumulation could be an effective way to target the parasite .
the method developed here used limited transcriptome data to identify pathways that underlie the inhibition of dna synthesis of p.  falciparum in response to chloroquine treatment .
this suggests that the developed method could also be suitable to simulate and capture the effect of other drugs ( or stress conditions ) on complex metabolic and physiological responses of p.  falciparum , based solely on transcriptome data .
 chloroquine , long the default first - line treatment against malaria , is now abandoned in large parts of the world because of widespread drug - resistance in plasmodium falciparum . in spite of its importance as a cost - effective and efficient drug , a coherent understanding of the cellular mechanisms affected by chloroquine and how they influence 
 the fitness and survival of the parasite remains elusive . 
 here , we used a systems biology approach to integrate genome - scale transcriptomics to map out the effects of chloroquine , identify targeted metabolic pathways , and translate these findings into mechanistic insights . specifically , we first developed a method that integrates transcriptomic and metabolomic data , which we independently validated against a recently published set of such data for krebs - cycle mutants of p.  falciparum . 
 we then used the method to calculate the effect of chloroquine treatment on the metabolic flux profiles of p.  falciparum during the intraerythrocytic developmental cycle . 
 the model predicted dose - dependent inhibition of dna replication , in agreement with earlier experimental results for both drug - sensitive and drug - resistant p.  falciparum strains . 
 our simulations also corroborated experimental findings that suggest differences in chloroquine sensitivity between ring- and schizont - stage p.  falciparum . 
 our analysis also suggests that metabolic fluxes that govern reduced thioredoxin and phosphoenolpyruvate synthesis are significantly decreased and are pivotal to chloroquine - based inhibition of p.  falciparum dna replication . 
 the consequences of impaired phosphoenolpyruvate synthesis and redox metabolism are reduced carbon fixation and increased oxidative stress , respectively , both of which eventually facilitate killing of the parasite . 
 our analysis suggests that a combination of chloroquine ( or an analogue ) and another drug , which inhibits carbon fixation and/or increases oxidative stress , should increase the clearance of p.  falciparum from the host system . 
Introduction
Materials and methods
Results
Discussion
Limitations of the current study
Conclusion