An integrative analysis of small molecule transcriptional responses in the human malaria parasite Plasmodium falciparum

PfHDAC1 is an essential regulator of P. falciparum asexual proliferation and host cell invasion genes with a dynamic genomic occupancy responsive to artemisinin stress

Plasmodium falciparum, the deadly protozoan parasite responsible for malaria, has a tightly regulated gene expression profile closely linked to its intraerythrocytic development cycle. Epigenetic modifiers of the histone acetylation code have been identified as key regulators of the parasite's transcriptome but require further investigation. In this study, we map the genomic distribution of Plasmodium falciparum histone deacetylase 1 (PfHDAC1) across the erythrocytic asexual development cycle and find it has a dynamic occupancy over a wide array of developmentally relevant genes. Overexpression of PfHDAC1 results in a progressive increment in parasite load over consecutive rounds of the asexual infection cycle and is associated with enhanced gene expression of multiple families of host cell invasion factors (merozoite surface proteins, rhoptry proteins, etc.) and with increased merozoite invasion efficiency. With the use of class-specific inhibitors, we demonstrate that PfHDAC1 activity in parasites is crucial for timely intraerythrocytic development. Interestingly, overexpression of PfHDAC1 results in decreased sensitivity to frontline-drug dihydroartemisinin in parasites. Furthermore, we identify that artemisinin exposure can interfere with PfHDAC1 abundance and chromatin occupancy, resulting in enrichment over genes implicated in response/resistance to artemisinin. Finally, we identify that dihydroartemisinin exposure can interrupt the in vitro catalytic deacetylase activity and post-translational phosphorylation of PfHDAC1, aspects that are crucial for its genomic function. Collectively, our results demonstrate PfHDAC1 to be a regulator of critical functions in asexual parasite development and host invasion, which is responsive to artemisinin exposure stress and deterministic of resistance to it.

IMPORTANCE

Malaria is a major public health problem, with the parasite Plasmodium falciparum causing most of the malaria-associated mortality. It is spread by the bite of infected mosquitoes and results in symptoms such as cyclic fever, chills, and headache. However, if left untreated, it can quickly progress to a more severe and life-threatening form. The World Health Organization currently recommends the use of artemisinin combination therapy, and it has worked as a gold standard for many years. Unfortunately, certain countries in southeast Asia and Africa, burdened with a high prevalence of malaria, have reported cases of drug-resistant infections. One of the major problems in controlling malaria is the emergence of artemisinin resistance. Population genomic studies have identified mutations in the Kelch13 gene as a molecular marker for artemisinin resistance. However, several reports thereafter indicated that Kelch13 is not the main mediator but rather hinted at transcriptional deregulation as a major determinant of drug resistance. Earlier, we identified PfGCN5 as a global regulator of stress-responsive genes, which are known to play a central role in artemisinin resistance generation. In this study, we have identified PfHDAC1, a histone deacetylase as a cell cycle regulator, playing an important role in artemisinin resistance generation. Taken together, our study identified key transcriptional regulators that play an important role in artemisinin resistance generation.

Regulation of Axon Guidance by Slit2 and Netrin-1 Signaling in the Lacrimal Gland of Aqp5 Knockout Mice

Purpose

Dry eye disease (DED) is multifactorial and associated with nerve abnormalities. We explored an Aquaporin 5 (AQP5)-deficiency-induced JunB activation mechanism, which causes abnormal lacrimal gland (LG) nerve distribution through Slit2 upregulation and Netrin-1 repression.

Methods

Aqp5 knockout (Aqp5-/-) and wild-type (Aqp5+/+) mice were studied. LGs were permeabilized and stained with neuronal class III β-tubulin, tyrosine hydroxylase (TH), vasoactive intestinal peptide (VIP), and calcitonin gene-related peptide (CGRP). Whole-mount images were acquired through tissue clearing and 3D fluorescence imaging. Mouse primary trigeminal ganglion (TG) neurons were treated with LG extracts and Netrin-1/Slit2 neutralizing antibody. Transcription factor (TF) prediction and chromatin immunoprecipitation-polymerase chain reaction (ChIP-PCR) experiments verified the JunB binding and regulatory effect on Netrin-1 and Slit2.

Results

Three-dimensional tissue and section immunofluorescence showed reduced LG nerves in Aqp5-/- mice, with sympathetic and sensory nerves significantly decreased. Netrin-1 was reduced and Slit2 increased in Aqp5-/- mice LGs. Aqp5+/+ mice LG tissue extracts (TEs) promoted Aqp5-/- TG neurons axon growth, but Netrin-1 neutralizing antibody (NAb) could inhibit that promotion. Aqp5-/- mice LG TEs inhibited Aqp5+/+ TG axon growth, but Slit2 NAb alleviated that inhibition. Furthermore, JunB, a Netrin-1 and Slit2 TF, could bind them and regulate their expression. SR11302, meanwhile, reversed the Netrin-1 and Slit2 shifts caused by AQP5 deficiency.

Conclusions

AQP5 deficiency causes LG nerve abnormalities. Persistent JunB activation, the common denominator for Netrin-1 suppression and Slit2 induction, was found in Aqp5-/- mice LG epithelial cells. This affected sensory and sympathetic nerve fibers' distribution in LGs. Our findings provide insights into preventing, reversing, and treating DED.

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 (Dd2^A211V) 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 Dd2^A211V 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.

In silico Structure Prediction, Molecular Docking, and Dynamic Simulation of Plasmodium falciparum AP2-I Transcription Factor

Plasmodium falciparum Apicomplexan Apetala 2 Invasion (PfAP2-I) transcription factor (TF) is a protein that regulates the expression of a subset of gene families involved in P. falciparum red blood cell (RBC) invasion. Inhibiting PfAP2-I TF with small molecules represents a potential new antimalarial therapeutic target to combat drug resistance, which this study aims to achieve. The 3D model structure of PfAP2-I was predicted ab initio using ROBETTA prediction tool and was validated using Save server 6.0 and MolProbity. Computed Atlas of Surface Topography of proteins (CASTp) 3.0 was used to predict the active sites of the PfAP2-I modeled structure. Pharmacophore modeling of the control ligand and PfAP2-I modeled structure was carried out using the Pharmit server to obtain several compounds used for molecular docking analysis. Molecular docking and postdocking studies were conducted using AutoDock vina and Discovery studio. The designed ligands' toxicity predictions and in silico drug-likeness were performed using the SwissADME predictor and OSIRIS Property Explorer. The modeled protein structure from the ROBETTA showed a validation result of 96.827 for ERRAT, 90.2% of the amino acid residues in the most favored region for the Ramachandran plot, and MolProbity score of 1.30 in the 98th percentile. Five (5) best hit compounds from molecular docking analysis were selected based on their binding affinity (between -8.9 and -11.7 Kcal/mol) to the active site of PfAP2-I and were considered for postdocking studies. For the absorption, distribution, metabolism, elimination, and toxicity (ADMET) properties, compound MCULE-7146940834 had the highest drug score (0.63) and drug-likeness (6.76). MCULE-7146940834 maintained a stable conformation within the flexible protein's active site during simulation. The good, estimated binding energies, drug-likeness, drug score, and molecular dynamics simulation interaction observed for MCULE-7146940834 against PfAP2-I show that MCULE-7146940834 can be considered a lead candidate for PfAP2-I inhibition. Experimental validations should be carried out to ascertain the efficacy of these predicted best hit compounds.

Sources of transcription variation in Plasmodium falciparum

Variation in transcript abundance can contribute to both short-term environmental response and long-term evolutionary adaptation. Most studies are designed to assess differences in mean transcription levels and do not consider other potentially important and confounding sources of transcriptional variation. Detailed quantification of variation sources will improve our ability to detect and identify the mechanisms that contribute to genome-wide transcription changes that underpin adaptive responses. To quantify innate levels of expression variation, we measured mRNA levels for more than 5000 genes in the malaria parasite, Plasmodium falciparum, among clones derived from two parasite strains across biologically and experimentally replicated batches. Using a mixed effects model, we partitioned the total variation among four sources - between strain, within strain, environmental batch effects, and stochastic noise. We found 646 genes with significant variation attributable to at least one of these sources. These genes were categorized by their predominant variation source and further examined using gene ontology enrichment analysis to associate function with each source of variation. Genes with environmental batch effect and within strain transcript variation may contribute to phenotypic plasticity, while genes with between strain variation may contribute to adaptive responses and processes that lead to parasite strain-specific survival under varied conditions.

Chemogenomic Fingerprints Associated with Stage-Specific Gametocytocidal Compound Action against Human Malaria Parasites

Kinase-focused inhibitors previously revealed compounds with differential activity against different stages of Plasmodium falciparum gametocytes. MMV666810, a 2-aminopyrazine, is more active on late-stage gametocytes, while a pyrazolopyridine, MMV674850, preferentially targets early-stage gametocytes. Here, we probe the biological mechanisms underpinning this differential stage-specific killing using in-depth transcriptome fingerprinting. Compound-specific chemogenomic profiles were observed with MMV674850 treatment associated with biological processes shared between asexual blood stage parasites and early-stage gametocytes but not late-stage gametocytes. MMV666810 has a distinct profile with clustered gene sets associated primarily with late-stage gametocyte development, including Ca²⁺-dependent protein kinases (CDPK1 and 5) and serine/threonine protein kinases (FIKK). Chemogenomic profiling therefore highlights essential processes in late-stage gametocytes, on a transcriptional level. This information is important to prioritize compounds that preferentially compromise late-stage gametocytes for further development as transmission-blocking antimalarials.

Machine Learning Uses Chemo-Transcriptomic Profiles to Stratify Antimalarial Compounds With Similar Mode of Action

The rapid development of antimalarial resistance motivates the continued search for novel compounds with a mode of action (MoA) different to current antimalarials. Phenotypic screening has delivered thousands of promising hit compounds without prior knowledge of the compounds’ exact target or MoA. Whilst the latter is not initially required to progress a compound in a medicinal chemistry program, identifying the MoA early can accelerate hit prioritization, hit-to-lead optimization and preclinical combination studies in malaria research. The effects of drug treatment on a cell can be observed on systems level in changes in the transcriptome, proteome and metabolome. Machine learning (ML) algorithms are powerful tools able to deconvolute such complex chemically-induced transcriptional signatures to identify pathways on which a compound act and in this manner provide an indication of the MoA of a compound. In this study, we assessed different ML approaches for their ability to stratify antimalarial compounds based on varied chemically-induced transcriptional responses. We developed a rational gene selection approach that could identify predictive features for MoA to train and generate ML models. The best performing model could stratify compounds with similar MoA with a classification accuracy of 76.6 ± 6.4%. Moreover, only a limited set of 50 biomarkers was required to stratify compounds with similar MoA and define chemo-transcriptomic fingerprints for each compound. These fingerprints were unique for each compound and compounds with similar targets/MoA clustered together. The ML model was specific and sensitive enough to group new compounds into MoAs associated with their predicted target and was robust enough to be extended to also generate chemo-transcriptomic fingerprints for additional life cycle stages like immature gametocytes. This work therefore contributes a new strategy to rapidly, specifically and sensitively indicate the MoA of compounds based on chemo-transcriptomic fingerprints and holds promise to accelerate antimalarial drug discovery programs.

Cellular thermal shift assay for the identification of drug-target interactions in the Plasmodium falciparum proteome

Despite decades of research, little is known about the cellular targets and the mode of action of the vast majority of antimalarial drugs. We recently demonstrated that the cellular thermal shift assay (CETSA) protocol in its two variants: the melt curve and the isothermal dose-response, represents a comprehensive strategy for the identification of antimalarial drug targets. CETSA enables proteome-wide target screening for unmodified antimalarial compounds with undetermined mechanisms of action, providing quantitative evidence about direct drug-protein interactions. The experimental workflow involves treatment of P. falciparum-infected erythrocytes with a compound of interest, heat exposure to denature proteins, soluble protein isolation, enzymatic digestion, peptide labeling with tandem mass tags, offline fractionation, and liquid chromatography-tandem mass spectrometry analysis. Methodological optimizations necessary for the analysis of this intracellular parasite are discussed, including enrichment of parasitized cells and hemoglobin depletion strategies to overcome high hemoglobin abundance in the host red blood cells. We outline an effective data processing workflow using the mineCETSA R package, which enables prioritization of drug-target candidates for follow-up studies. The entire protocol can be completed within 2 weeks.

Simultaneous genome-wide gene expression and transcript isoform profiling in the human malaria parasite

Gene expression DNA microarrays have been vital for characterizing whole-genome transcriptional profiles. Nevertheless, their effectiveness relies heavily on the accuracy of genome sequences, the annotation of gene structures, and the sequence-dependent performance of individual probes. Currently available gene expression arrays for the malaria parasite Plasmodium falciparum rely on an average of 2 probes per gene, usually positioned near the 3′ end of genes; consequently, existing designs are prone to measurement bias and cannot capture complexities such as the occurrence of transcript isoforms arising from alternative splicing or alternative start/ stop sites. Here, we describe two novel gene expression arrays with exon-focused probes designed with an average of 12 and 20 probes spanning each gene. This high probe density minimizes signal noise inherent in probe-to-probe sequence-dependent hybridization intensity. We demonstrate that these exon arrays accurately profile genome-wide expression, comparing favorably to currently available arrays and RNA-seq profiling, and can detect alternatively spliced transcript isoforms as well as non-coding RNAs (ncRNAs). Of the 964 candidate alternate splicing events from published RNA-seq studies, 162 are confirmed using the exon array. Furthermore, the exon array predicted 330 previously unidentified alternate splicing events. Gene expression microarrays continue to offer a cost-effective alternative to RNA-seq for the simultaneous monitoring of gene expression and alternative splicing events. Microarrays may even be preferred in some cases due to their affordability and the rapid turn-around of results when hundreds of samples are required for fine-scale systems biology investigations, including the monitoring of the networks of gene co-expression in the emergence of drug resistance.

Knockout of the peroxiredoxin 5 homologue PFAOP does not affect the artemisinin susceptibility of Plasmodium falciparum

Artemisinins are the current mainstay of malaria chemotherapy. Their exact mode of action is an ongoing matter of debate, and several factors have recently been reported to affect an early stage of artemisinin resistance of the most important human malaria parasite Plasmodium falciparum. Here, we identified a locus on chromosome 7 that affects the artemisinin susceptibility of P. falciparum in a quantitative trait locus analysis of a genetic cross between strains 7G8 and GB4. This locus includes the peroxiredoxin gene PFAOP. However, steady-state kinetic data with recombinant PfAOP do not support a direct interaction between this peroxidase and the endoperoxide artemisinin. Furthermore, neither the overexpression nor the deletion of the encoding gene affected the IC₅₀ values for artemisinin or the oxidants diamide and tert-butyl hydroperoxide. Thus, PfAOP is dispensable for blood stage parasite survival, and the correlation between the artemisinin susceptibility and chromosome 7 is probably based on another gene within the identified locus.

Muddled mechanisms: recent progress towards antimalarial target identification

In the past decade, malaria rates have plummeted as a result of aggressive infection control measures and the adoption of artemisinin-based combination therapies (ACTs). However, a potential crisis looms ahead. Treatment failures to standard antimalarial regimens have been reported in Southeast Asia, and devastating consequences are expected if resistance spreads to the African continent. To prevent a potential public health emergency, the antimalarial arsenal must contain therapeutics with novel mechanisms of action (MOA). An impressive number of high-throughput screening (HTS) campaigns have since been launched, identifying thousands of compounds with activity against one of the causative agents of malaria, Plasmodium falciparum. Now begins the difficult task of target identification, for which studies are often tedious, labor intensive, and difficult to interpret. In this review, we highlight approaches that have been instrumental in tackling the challenges of target assignment and elucidation of the MOA for hit compounds. Studies that apply these innovative techniques to antimalarial target identification are described, as well as the impact of the data in the field.

Discovering New Transmission-Blocking Antimalarial Compounds: Challenges and Opportunities

The ability to target human-mosquito parasite transmission challenges global malaria elimination. However, it is not obvious what a transmission-blocking drug will look like; should it target only parasite transmission stages; be combined with a partner drug killing the pathogenic asexual stages; or kill both the sexual and asexual blood stages, preferably displaying polypharmacology? The development of transmission-blocking antimalarials requires objective analyses of the current strategies. Here, pertinent issues and questions regarding the target candidate profile of a transmission-blocking compound, and its role in malaria elimination strategies, are highlighted and novel perspectives proposed. The essential role of a test cascade that integrates screening and validation strategies to identify next-generation transmission-blocking antimalarials is emphasised.

PfCRT and PfMDR1 modulate interactions of artemisinin derivatives and ion channel blockers

Treatment of the symptomatic asexual stage of Plasmodium falciparum relies almost exclusively on artemisinin (ART) combination therapies (ACTs) in endemic regions. ACTs combine ART or its derivative with a long-acting partner drug to maximize efficacy during the typical three-day regimen. Both laboratory and clinical studies have previously demonstrated that the common drug resistance determinants P. falciparum chloroquine resistance transporter (PfCRT) and multidrug resistance transporter (PfMDR1) can modulate the susceptibility to many current antimalarial drugs and chemical compounds. Here we investigated the parasite responses to dihydroartemisinin (DHA) and various Ca²⁺ and Na⁺ channel blockers and showed positively correlated responses between DHA and several channel blockers, suggesting potential shared transport pathways or mode of action. Additionally, we demonstrated that PfCRT and PfMDR1 could also significantly modulate the pharmacodynamic interactions of the compounds and that the interactions were influenced by the parasite genetic backgrounds. These results provide important information for better understanding of drug resistance and for assessing the overall impact of drug resistance markers on parasite response to ACTs.