Short-term metabolic adjustments in Plasmodium falciparum counter hypoxanthine deprivation at the expense of long-term viability

Physiological Indirect Response Model to Omics-Powered Quantitative Systems Pharmacology Model

Over the past several decades, mathematical modeling has been applied to increasingly wider scopes of questions in drug development. Accordingly, the range of modeling tools has also been evolving, as showcased by contributions of Jusko and colleagues: from basic pharmacokinetics/pharmacodynamics (PK/PD) modeling to today's platform-based approach of quantitative systems pharmacology (QSP) modeling. Aimed at understanding the mechanism of action of investigational drugs, QSP models characterize systemic effects by incorporating information about cellular signaling networks, which is often represented by omics data. In this perspective, we share a few examples illustrating approaches for the integration of omics into mechanistic QSP modeling. We briefly overview how the evolution of PK/PD modeling into QSP has been accompanied by an increase in available data and the complexity of mathematical methods that integrate it. We discuss current gaps and challenges of integrating omics data into QSP models and propose several potential areas where integrated QSP and omics modeling may benefit drug development.

Cultivation of Asexual Intraerythrocytic Stages of Plasmodium falciparum

Successfully developed in 1976, the continuous in vitro culture of Plasmodium falciparum has many applications in the field of malaria research. It has become an important experimental model that directly uses a human pathogen responsible for a high prevalence of morbidity and mortality in many parts of the world and is a major source of biological material for immunological, biochemical, molecular, and pharmacological studies. Until present, the basic techniques described by Trager and Jensen and Haynes et al. remain unchanged in many malaria research laboratories. Nonetheless, different factors, including culture media, buffers, serum substitutes and supplements, sources of erythrocytes, and conditions of incubation (especially oxygen concentration), have been modified by different investigators to adapt the original technique in their laboratories or enhance the in vitro growth of the parasites. The possible effects and benefits of these modifications for the continuous cultivation of asexual intraerythrocytic stages of P. falciparum, as well as future challenges in developing a serum-free cultivation system and axenic cultures, are discussed.

Nutrient Limitation Mimics Artemisinin Tolerance in Malaria

Mounting evidence demonstrates that nutritional environment can alter pathogen drug sensitivity. While the rich media used for in vitro culture contains supraphysiological nutrient concentrations, pathogens encounter a relatively restrictive environment in vivo. We assessed the effect of nutrient limitation on the protozoan parasite that causes malaria and demonstrated that short-term growth under physiologically relevant mild nutrient stress (or "metabolic priming") triggers increased tolerance of a potent antimalarial drug. We observed beneficial effects using both short-term survival assays and longer-term proliferation studies, where metabolic priming increases parasite survival to a level previously defined as resistant (>1% survival). We performed these assessments by either decreasing single nutrients that have distinct roles in metabolism or using a media formulation that simulates the human plasma environment. We determined that priming-induced tolerance was restricted to parasites that had newly invaded the host red blood cell, but the effect was not dependent on genetic background. The molecular mechanisms of this intrinsic effect mimic aspects of genetic tolerance, including translational repression and protein export. This finding suggests that regardless of the impact on survival rates, environmental stress could stimulate changes that ultimately directly contribute to drug tolerance. Because metabolic stress is likely to occur more frequently in vivo compared to the stable in vitro environment, priming-induced drug tolerance has ramifications for how in vitro results translate to in vivo studies. Improving our understanding of how pathogens adjust their metabolism to impact survival of current and future drugs is an important avenue of research to slow the evolution of resistance. IMPORTANCE There is a dire need for effective treatments against microbial pathogens. Yet, the continuing emergence of drug resistance necessitates a deeper knowledge of how pathogens respond to treatments. We have long appreciated the contribution of genetic evolution to drug resistance, but transient metabolic changes that arise in response to environmental factors are less recognized. Here, we demonstrate that short-term growth of malaria parasites in a nutrient-limiting environment triggers cellular changes that lead to better survival of drug treatment. We found that these strategies are similar to those employed by drug-tolerant parasites, which suggests that starvation "primes" parasites to survive and potentially evolve resistance. Since the environment of the human host is relatively nutrient restrictive compared to growth conditions in standard laboratory culture, this discovery highlights the important connections among nutrient levels, protective cellular pathways, and resistance evolution.

Metabolic responses in blood-stage malaria parasites associated with increased and decreased sensitivity to PfATP4 inhibitors

BACKGROUND

Spiroindolone and pyrazoleamide antimalarial compounds target Plasmodium falciparum P-type ATPase (PfATP4) and induce disruption of intracellular Na+ homeostasis. Recently, a PfATP4 mutation was discovered that confers resistance to a pyrazoleamide while increasing sensitivity to a spiroindolone. Transcriptomic and metabolic adaptations that underlie this seemingly contradictory response of P. falciparum to sublethal concentrations of each compound were examined to understand the different cellular accommodation to PfATP4 disruptions.

METHODS

A genetically engineered P. falciparum Dd2 strain (Dd2A211V) carrying an Ala211Val (A211V) mutation in PfATP4 was used to identify metabolic adaptations associated with the mutation that results in decreased sensitivity to PA21A092 (a pyrazoleamide) and increased sensitivity to KAE609 (a spiroindolone). First, sublethal doses of PA21A092 and KAE609 causing substantial reduction (30-70%) in Dd2A211V parasite replication were identified. Then, at this sublethal dose of PA21A092 (or KAE609), metabolomic and transcriptomic data were collected during the first intraerythrocytic developmental cycle. Finally, the time-resolved data were integrated with a whole-genome metabolic network model of P. falciparum to characterize antimalarial-induced physiological adaptations.

RESULTS

Sublethal treatment with PA21A092 caused significant (p < 0.001) alterations in the abundances of 91 Plasmodium gene transcripts, whereas only 21 transcripts were significantly altered due to sublethal treatment with KAE609. In the metabolomic data, a substantial alteration (≥ fourfold) in the abundances of carbohydrate metabolites in the presence of either compound was found. The estimated rates of macromolecule syntheses between the two antimalarial-treated conditions were also comparable, except for the rate of lipid synthesis. A closer examination of parasite metabolism in the presence of either compound indicated statistically significant differences in enzymatic activities associated with synthesis of phosphatidylcholine, phosphatidylserine, and phosphatidylinositol.

CONCLUSION

The results of this study suggest that malaria parasites activate protein kinases via phospholipid-dependent signalling in response to the ionic perturbation induced by the Na+ homeostasis disruptor PA21A092. Therefore, targeted disruption of phospholipid signalling in PA21A092-resistant parasites could be a means to block the emergence of resistance to PA21A092.

Metabolic changes accompanying the loss of fumarate hydratase and malate-quinone oxidoreductase in the asexual blood stage of Plasmodium falciparum

In the glucose-rich milieu of red blood cells, asexually replicating malarial parasites mainly rely on glycolysis for ATP production, with limited carbon flux through the mitochondrial tricarboxylic acid (TCA) cycle. By contrast, gametocytes and mosquito-stage parasites exhibit an increased dependence on the TCA cycle and oxidative phosphorylation for more economical energy generation. Prior genetic studies supported these stage-specific metabolic preferences by revealing that six of eight TCA cycle enzymes are completely dispensable during the asexual blood stages of Plasmodium falciparum, with only fumarate hydratase (FH) and malate-quinone oxidoreductase (MQO) being refractory to deletion. Several hypotheses have been put forth to explain the possible essentiality of FH and MQO, including their participation in a malate shuttle between the mitochondrial matrix and the cytosol. However, using newer genetic techniques like CRISPR and dimerizable Cre, we were able to generate deletion strains of FH and MQO in P. falciparum. We employed metabolomic analyses to characterize a double knockout mutant of FH and MQO (ΔFM) and identified changes in purine salvage and urea cycle metabolism that may help to limit fumarate accumulation. Correspondingly, we found that the ΔFM mutant was more sensitive to exogenous fumarate, which is known to cause toxicity by modifying and inactivating proteins and metabolites. Overall, our data indicate that P. falciparum is able to adequately compensate for the loss of FH and MQO, rendering them unsuitable targets for drug development.

Metabolic adjustments of blood-stage Plasmodium falciparum in response to sublethal pyrazoleamide exposure

Due to the recurring loss of antimalarial drugs to resistance, there is a need for novel targets, drugs, and combination therapies to ensure the availability of current and future countermeasures. Pyrazoleamides belong to a novel class of antimalarial drugs that disrupt sodium ion homeostasis, although the exact consequences of this disruption in Plasmodium falciparum remain under investigation. In vitro experiments demonstrated that parasites carrying mutations in the metabolic enzyme PfATP4 develop resistance to pyrazoleamide compounds. However, the underlying mechanisms that allow mutant parasites to evade pyrazoleamide treatment are unclear. Here, we first performed experiments to identify the sublethal dose of a pyrazoleamide compound (PA21A092) that caused a significant reduction in growth over one intraerythrocytic developmental cycle (IDC). At this drug concentration, we collected transcriptomic and metabolomic data at multiple time points during the IDC to quantify gene- and metabolite-level alterations in the treated parasites. To probe the effects of pyrazoleamide treatment on parasite metabolism, we coupled the time-resolved omics data with a metabolic network model of P. falciparum. We found that the drug-treated parasites adjusted carbohydrate metabolism to enhance synthesis of myoinositol-a precursor for phosphatidylinositol biosynthesis. This metabolic adaptation caused a decrease in metabolite flux through the pentose phosphate pathway, causing a decreased rate of RNA synthesis and an increase in oxidative stress. Our model analyses suggest that downstream consequences of enhanced myoinositol synthesis may underlie adjustments that could lead to resistance emergence in P. falciparum exposed to a sublethal dose of a pyrazoleamide drug.

Inter-study and time-dependent variability of metabolite abundance in cultured red blood cells

BACKGROUND

Cultured human red blood cells (RBCs) provide a powerful ex vivo assay platform to study blood-stage malaria infection and propagation. In recent years, high-resolution metabolomic methods have quantified hundreds of metabolites from parasite-infected RBC cultures under a variety of perturbations. In this context, the corresponding control samples of the uninfected culture systems can also be used to examine the effects of these perturbations on RBC metabolism itself and their dependence on blood donors (inter-study variations).

METHODS

Time-course datasets from five independent studies were generated and analysed, maintaining uninfected RBCs (uRBC) at 2% haematocrit for 48 h under conditions originally designed for parasite cultures. Using identical experimental protocols, quadruplicate samples were collected at six time points, and global metabolomics were employed on the pellet fraction of the uRBC cultures. In total, ~ 500 metabolites were examined across each dataset to quantify inter-study variability in RBC metabolism, and metabolic network modelling augmented the analyses to characterize the metabolic state and fluxes of the RBCs.

RESULTS

To minimize inter-study variations unrelated to RBC metabolism, an internal standard metabolite (phosphatidylethanolamine C18:0/20:4) was identified with minimal variation in abundance over time and across all the samples of each dataset to normalize the data. Although the bulk of the normalized data showed a high degree of inter-study consistency, changes and variations in metabolite levels from individual donors were noted. Thus, a total of 24 metabolites were associated with significant variation in the 48-h culture time window, with the largest variations involving metabolites in glycolysis and synthesis of glutathione. Metabolic network analysis was used to identify the production of superoxide radicals in cultured RBCs as countered by the activity of glutathione oxidoreductase and synthesis of reducing equivalents via the pentose phosphate pathway. Peptide degradation occurred at a rate that is comparable with central carbon fluxes, consistent with active degradation of methaemoglobin, processes also commonly associated with storage lesions in RBCs.

CONCLUSIONS

The bulk of the data showed high inter-study consistency. The collected data, quantification of an expected abundance variation of RBC metabolites, and characterization of a subset of highly variable metabolites in the RBCs will help in identifying non-specific changes in metabolic abundances that may obscure accurate metabolomic profiling of Plasmodium falciparum and other blood-borne pathogens.

Metabolic Survival Adaptations of Plasmodium falciparum Exposed to Sublethal Doses of Fosmidomycin

The malaria parasite Plasmodium falciparum contains the apicoplast organelle that synthesizes isoprenoids, which are metabolites necessary for posttranslational modification of Plasmodium proteins. We used fosmidomycin, an antibiotic that inhibits isoprenoid biosynthesis, to identify mechanisms that underlie the development of the parasite's adaptation to the drug at sublethal concentrations. We first determined a concentration of fosmidomycin that reduced parasite growth by ∼50% over one intraerythrocytic developmental cycle (IDC). At this dose, we maintained synchronous parasite cultures for one full IDC and collected metabolomic and transcriptomic data at multiple time points to capture global and stage-specific alterations. We integrated the data with a genome-scale metabolic model of P. falciparum to characterize the metabolic adaptations of the parasite in response to fosmidomycin treatment. Our simulations showed that, in treated parasites, the synthesis of purine-based nucleotides increased, whereas the synthesis of phosphatidylcholine during the trophozoite and schizont stages decreased. Specifically, the increased polyamine synthesis led to increased nucleotide synthesis, while the reduced methyl-group cycling led to reduced phospholipid synthesis and methyltransferase activities. These results indicate that fosmidomycin-treated parasites compensate for the loss of prenylation modifications by directly altering processes that affect nucleotide synthesis and ribosomal biogenesis to control the rate of RNA translation during the IDC. This also suggests that combination therapies with antibiotics that target the compensatory response of the parasite, such as nucleotide synthesis or ribosomal biogenesis, may be more effective than treating the parasite with fosmidomycin alone.

Metabolic alterations in the erythrocyte during blood-stage development of the malaria parasite

Background

Human blood cells (erythrocytes) serve as hosts for the malaria parasite Plasmodium falciparum during its 48-h intraerythrocytic developmental cycle (IDC). Established in vitro protocols allow for the study of host–parasite interactions during this phase and, in particular, high-resolution metabolomics can provide a window into host–parasite interactions that support parasite development.

Methods

Uninfected and parasite-infected erythrocyte cultures were maintained at 2% haematocrit for the duration of the IDC, while parasitaemia was maintained at 7% in the infected cultures. The parasite-infected cultures were synchronized to obtain stage-dependent information of parasite development during the IDC. Samples were collected in quadruplicate at six time points from the uninfected and parasite-infected cultures and global metabolomics was used to analyse cell fractions of these cultures.

Results

In uninfected and parasite-infected cultures during the IDC, 501 intracellular metabolites, including 223 lipid metabolites, were successfully quantified. Of these, 19 distinct metabolites were present only in the parasite-infected culture, 10 of which increased to twofold in abundance during the IDC. This work quantified approximately five times the metabolites measured in previous studies of similar research scope, which allowed for more detailed analyses. Enrichment in lipid metabolism pathways exhibited a time-dependent association with different classes of lipids during the IDC. Specifically, enrichment occurred in sphingolipids at the earlier stages, and subsequently in lysophospholipid and phospholipid metabolites at the intermediate and end stages of the IDC, respectively. In addition, there was an accumulation of 18-, 20-, and 22-carbon polyunsaturated fatty acids, which produce eicosanoids and promote gametocytogenesis in infected erythrocyte cultures.

Conclusions

The current study revealed a number of heretofore unidentified metabolic components of the host–parasite system, which the parasite may exploit in a time-dependent manner to grow over the course of its development in the blood stage. Notably, the analyses identified components, such as precursors of immunomodulatory molecules, stage-dependent lipid dynamics, and metabolites, unique to parasite-infected cultures. These conclusions are reinforced by the metabolic alterations that were characterized during the IDC, which were in close agreement with those known from previous studies of blood-stage infection.

A mevalonate bypass system facilitates elucidation of plastid biology in malaria parasites

Malaria parasites rely on a plastid organelle for survival during the blood stages of infection. However, the entire organelle is dispensable as long as the isoprenoid precursor, isopentenyl pyrophosphate (IPP), is supplemented in the culture medium. We engineered parasites to produce isoprenoid precursors from a mevalonate-dependent pathway, creating a parasite line that replicates normally after the loss of the apicoplast organelle. We show that carbon-labeled mevalonate is specifically incorporated into isoprenoid products, opening new avenues for researching this essential class of metabolites in malaria parasites. We also show that essential apicoplast proteins, such as the enzyme target of the drug fosmidomycin, can be deleted in this mevalonate bypass parasite line, providing a new method to determine the roles of other important apicoplast-resident proteins. Several antibacterial drugs kill malaria parasites by targeting basic processes, such as transcription, in the organelle. We used metabolomic and transcriptomic methods to characterize parasite metabolism after azithromycin treatment triggered loss of the apicoplast and found that parasite metabolism and the production of apicoplast proteins is largely unaltered. These results provide insight into the effects of apicoplast-disrupting drugs, several of which have been used to treat malaria infections in humans. Overall, the mevalonate bypass system provides a way to probe essential aspects of apicoplast biology and study the effects of drugs that target apicoplast processes.

Malaria parasites rely on an organelle called the apicoplast for growth and survival. Antimalarial drugs such as azithromycin inhibit basic processes in the apicoplast and result in the disruption of the organelle. Surprisingly, addition of a single metabolite, isopentenyl pyrophosphate (IPP), allows the parasites to survive in culture after disruption of the apicoplast. Unfortunately, using IPP to study this phenomenon has several limitations: IPP is prohibitively expensive, has to be used at high concentrations, and has a half-life less than 5 hours. To address these problems, we engineered parasites to express four enzymes from an alternative pathway capable of producing IPP in the parasites. We validated this new system and used it to metabolically label essential metabolites, to delete an essential apicoplast protein, and to characterize the state of apicoplast-disrupted parasites. A key finding from these studies comes from transcriptomic and metabolomic analysis of parasites treated with the drug azithromycin. We found that apicoplast disruption results in few changes in parasite metabolism. In particular, the expression of hundreds of nuclear-encoded apicoplast proteins are not affected by disruption of the apicoplast organelle, making it likely that apicoplast metabolic pathways and processes are still functional in apicoplast-disrupted parasites.