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.

Awakening the sleeping giant: Epstein-Barr virus reactivation by biological agents

Epstein-Barr virus (EBV) may cause harm in immunocompromised conditions or on stress stimuli. Various chemical agents have been utilized to induce the lytic cycle in EBV-infected cells. However, apart from chemical agents and external stress stimuli, certain infectious agents may reactivate the EBV. In addition, the acute infection of other pathogens may provide suitable conditions for EBV to thrive more and planting the roots for EBV-associated pathologies. Various bacteria such as periodontal pathogens like Aggregatibacter, Helicobacter pylori, etc. have shown to induce EBV reactivation either by triggering host cells directly or indirectly. Viruses such as Human simplex virus-1 (HSV) induce EBV reactivation by HSV US3 kinase while other viruses such as HIV, hepatitis virus, and even novel SARS-CoV-2 have also been reported to cause EBV reactivation. The eukaryotic pathogens such as Plasmodium falciparum and Aspergillus flavus can also reactivate EBV either by surface protein interaction or as an impact of aflatoxin, respectively. To highlight the underexplored niche of EBV reactivation by biological agents, we have comprehensively presented the related information in this review. This may help to shedding the light on the research gaps as well as to unveil yet unexplored mechanisms of EBV reactivation.

Genomic surveillance reveals early detection and transition of delta to omicron lineages of SARS-CoV-2 variants in wastewater treatment plants of Pune, India

The COVID-19 pandemic has emphasized the urgency for rapid public health surveillance methods to detect and monitor the transmission of infectious diseases. The wastewater-based epidemiology (WBE) has emerged as a promising tool for proactive analysis and quantification of infectious pathogens within a population before clinical cases emerge. In the present study, we aimed to assess the trend and dynamics of SARS-CoV-2 variants using a longitudinal approach. Our objective included early detection and monitoring of these variants to enhance our understanding of their prevalence and potential impact. To achieve our goals, we conducted real-time quantitative polymerase chain reaction (RT-qPCR) and Illumina sequencing on 442 wastewater (WW) samples collected from 10 sewage treatment plants (STPs) in Pune city, India, spanning from November 2021 to April 2022. Our comprehensive analysis identified 426 distinct lineages representing 17 highly transmissible variants of SARS-CoV-2. Notably, fragments of Omicron variant were detected in WW samples prior to its first clinical detection in Botswana. Furthermore, we observed highly contagious sub-lineages of the Omicron variant, including BA.1 (~28%), BA.1.X (1.0-72%), BA.2 (1.0-18%), BA.2.X (1.0-97.4%) BA.2.12 (0.8-0.25%), BA.2.38 (0.8-1.0%), BA.2.75 (0.01-0.02%), BA.3 (0.09-6.3%), BA.4 (0.24-0.29%), and XBB (0.01-21.83%), with varying prevalence rates. Overall, the present study demonstrated the practicality of WBE in the early detection of SARS-CoV-2 variants, which could help track future outbreaks of SARS-CoV-2. Such approaches could be implicated in monitoring infectious agents before they appear in clinical cases.

The many paths to artemisinin resistance in Plasmodium falciparum

Emerging resistance against artemisinin (ART) poses a major challenge in controlling malaria. Parasites with mutations in PfKelch13, the major marker for ART resistance, are known to reduce hemoglobin endocytosis, induce unfolded protein response (UPR), elevate phosphatidylinositol-3-phosphate (PI3P) levels, and stimulate autophagy. Nonetheless, PfKelch13-independent resistance is also reported, indicating extensive complementation by reconfiguration in the parasite metabolome and transcriptome. These findings implicate that there may not be a single 'universal identifier' of ART resistance. This review sheds light on the molecular, transcriptional, and metabolic pathways associated with ART resistance, while also highlighting the interplay between cellular heterogeneity, environmental stress, and ART sensitivity.

RNA Polymerase II evolution and adaptations: Insights from Plasmodium and other parasitic protists

The C-terminal domain (CTD) of RNA polymerase II plays a crucial role in regulating transcription dynamics in eukaryotes. The phosphorylation of serine residues within the CTD controls transcription initiation, elongation, and termination. While the CTD is highly conserved across eukaryotes, lower eukaryotes like protists, including Plasmodium, exhibit some differences. In this study, we performed a comparative analysis of CTD in eukaryotic systems to understand why the parasites evolved in this particular manner. The Plasmodium falciparum RPB1 is exceptionally large and feature a gap between the first and second heptad repeats, resulting in fifteen canonical heptad repeats excluding the initial repeat. Analysis of this intervening sequence revealed sub motifs of heptads where two serine residues occupy the first and fourth positions (S1X2X3S4). These motifs lie in the intrinsically disordered region of RPB1, a characteristic feature of the CTD. Interestingly, the S1X2X3S4 sub-motif was also observed in early-divergingeukaryotes like Leishmania major, which lack canonical heptad repeats. Furthermore, eukaryotes across the phylogenetic tree revealed a sigmoid pattern of increasing serine frequency in the CTD, indicating that serine enrichment is a significant step in the evolution of heptad-rich RPB1. Based on these observations and analysis, we proposed an evolutionary model for RNA Polymerase II CTD, encompassing organisms previously deemed exceptions, notably Plasmodium species. Thus, our study provides novel insights into the evolution of the CTD and will prompt further investigations into the differences exhibited by Plasmodium RNA Pol II and determine if they confer a survival advantage to the parasite.

A tale of two waves: Delineating diverse genomic and transmission landscapes driving the COVID-19 pandemic in Pune, India

BACKGROUND

Modern response to pandemics, critical for effective public health measures, is shaped by the availability and integration of diverse epidemiological outbreak data. Tracking variants of concern (VOC) is integral to understanding the evolution of SARS-CoV-2 in space and time, both at the local level and global context. This potentially generates actionable information when integrated with epidemiological outbreak data.

METHODS

A city-wide network of researchers, clinicians, and pathology diagnostic laboratories was formed for genome surveillance of COVID-19 in Pune, India. The genomic landscapes of 10,496 sequenced samples of SARS-CoV-2 driving peaks of infection in Pune between December-2020 to March-2022, were determined. As a modern response to the pandemic, a "band of five" outbreak data analytics approach was used. This integrated the genomic data (Band 1) of the virus through molecular phylogenetics with key outbreak data including sample collection dates and case numbers (Band 2), demographics like age and gender (Band 3-4), and geospatial mapping (Band 5).

RESULTS

The transmission dynamics of VOCs in 10,496 sequenced samples identified B.1.617.2 (Delta) and BA(x) (Omicron formerly known as B.1.1.529) variants as drivers of the second and third peaks of infection in Pune. Spike Protein mutational profiling during pre and post-Omicron VOCs indicated differential rank ordering of high-frequency mutations in specific domains that increased the charge and binding properties of the protein. Time-resolved phylogenetic analysis of Omicron sub-lineages identified a highly divergent BA.1 from Pune in addition to recombinant X lineages, XZ, XQ, and XM.

CONCLUSIONS

The band of five outbreak data analytics approach, which integrates five different types of data, highlights the importance of a strong surveillance system with high-quality meta-data for understanding the spatiotemporal evolution of the SARS-CoV-2 genome in Pune. These findings have important implications for pandemic preparedness and could be critical tools for understanding and responding to future outbreaks.

Chasing SARS-CoV-2 XBB.1.16 Recombinant Lineage in India and the Clinical Profile of XBB.1.16 Cases in Maharashtra, India

Background SARS-CoV-2 has evolved rapidly, resulting in the emergence of lineages with a competitive advantage over one another. Co-infections with different SARS-CoV-2 lineages can give rise to recombinant lineages. To date, the XBB lineage is the most widespread recombinant lineage worldwide, with the recently named XBB.1.16 lineage causing a surge in the number of COVID-19 cases in India. Methodology The present study involved retrieval of SARS-CoV-2 genome sequences from India (between December 1, 2022 and April 8, 2023) through GISAID; sequences were curated, followed by lineage and phylogenetic analysis. Demographic and clinical data from Maharashtra, India were collected telephonically, recorded in Microsoft® Excel, and analyzed using IBM® SPSS statistics, version 29.0.0.0 (241). Results A total of 2,944 sequences were downloaded from the GISAID database, of which 2,856 were included in the study following data curation. The sequences from India were dominated by the XBB.1.16* lineage (36.17%) followed by XBB.2.3* (12.11%) and XBB.1.5* (10.36%). Of the 2,856 cases, 693 were from Maharashtra; 386 of these were included in the clinical study. The clinical features of COVID-19 cases with XBB.1.16* infection (XBB.1.16* cases, 276 in number) showed that 92% of those had a symptomatic disease, with fever (67%), cough (42%), rhinorrhea (33.7%), body ache (14.5%) and fatigue (14.1%) being the most common symptoms. The presence of comorbidity was found in 17.7% of the XBB.1.16* cases. Among the XBB.1.16* cases, 91.7% were vaccinated with at least one dose of vaccine against COVID-19. While 74.3% of XBB.1.16* cases were home-isolated; 25.7% needed hospitalization/institutional quarantine, of these, 33.8% needed oxygen therapy. Out of 276 XBB.1.16* cases, seven (2.5%) cases succumbed to the disease. The majority of XBB.1.16* cases who died belonged to an elderly age group (60 years and above), had underlying comorbid condition/s, and needed supplemental oxygen therapy. The clinical features of COVID-19 cases infected with other co-circulating Omicron variants were similar to XBB.1.16* cases. Conclusion The study reveals that XBB.1.16* lineage has become the most predominant SARS-CoV-2 lineage in India. The study also shows that the clinical features and outcome of XBB.1.16* cases were similar to those of other co-circulating Omicron lineage infected cases in Maharashtra, India.

Genomic analysis of Indian isolates of Plasmodium falciparum: Implications for drug resistance and virulence factors

The emergence of drug resistance to frontline treatments such as Artemisinin-based combination therapy (ACT) is a major obstacle to the control and eradication of malaria. This problem is compounded by the inherent genetic variability of the parasites, as many established markers of resistance do not accurately predict the drug-resistant status. There have been reports of declining effectiveness of ACT in the West Bengal and Northeast regions of India, which have traditionally been areas of drug resistance emergence in the country. Monitoring the genetic makeup of a population can help to identify the potential for drug resistance markers associated with it and evaluate the effectiveness of interventions aimed at reducing the spread of malaria. In this study, we performed whole genome sequencing of 53 isolates of Plasmodium falciparum from West Bengal and compared their genetic makeup to isolates from Southeast Asia (SEA) and Africa. We found that the Indian isolates had a distinct genetic makeup compared to those from SEA and Africa, and were more similar to African isolates, with a high prevalence of mutations associated with antigenic variation genes. The Indian isolates also showed a high prevalence of markers of chloroquine resistance (mutations in Pfcrt) and multidrug resistance (mutations in Pfmdr1), but no known mutations associated with artemisinin resistance in the PfKelch13 gene. Interestingly, we observed a novel L152V mutation in PfKelch13 gene and other novel mutations in genes involved in ubiquitination and vesicular transport that have been reported to support artemisinin resistance in the early stages of ACT resistance in the absence of PfKelch13 polymorphisms. Thus, our study highlights the importance of region-specific genomic surveillance for artemisinin resistance and the need for continued monitoring of resistance to artemisinin and its partner drugs.

Clinical Characteristics of SARS-CoV-2 Omicron Cases in Pune, Maharashtra, India

Background The SARS-CoV-2 Omicron variant, within two months of its detection, replaced the Delta variant to become the dominant circulating variant globally. Therefore, it is essential to understand the characteristics of the disease caused by the variant and its impact on vaccination. Methods A total of 165 confirmed Omicron cases attending a tertiary care hospital in Pune, Maharashtra, between December 2021 to February 2022 were studied. Their demographic, clinical, and immunization history was recorded. Results Among the 165 cases, 7.88% were B.1.1.529 Omicron cases, 25.45% were BA.1 Omicron cases, and 66.67% were BA.2 Omicron cases. Of these 165 patients, 146 (88.48%) were discharged after treatment, 12 (7.27%) died during hospitalization, and seven (4.24%) were brought dead. The presence of one or more comorbid conditions was seen in 15.15%, of which diabetes mellitus and hypertension (28% each) were the most common conditions. Older age (greater than 60 years), an important risk factor for poor outcomes, was present in 9.1% of cases. Among the 165 cases, vaccination with at least one dose of vaccine was found in 80.61% of cases. Out of 165 cases, clinical data was available for 158 cases. Of these 158 cases, 86.71% had symptoms, and 13.29% were asymptomatic. Fever, followed by cough, myalgia, runny nose, and headache, were the most common presenting symptoms. The mean duration of illness was 2.69 days, with 91.14% of cases having the illness for less than five days, and 89.24% of cases had a National Early Warning Score (NEWS) of 1-4, suggesting a good prognosis. In 93.90% of cases, the chest X-ray findings were normal. Of the 158 cases, 92.41% of cases recovered with supportive treatment, and only 7.59% of cases required oxygen therapy. Conclusion The current study shows that the Omicron variant caused mild disease with reduced need for hospital admission and oxygen therapy in India.

Clinical Characteristics and Outcomes of Laboratory-Confirmed SARS-CoV-2 Cases Infected With Omicron Subvariants and the XBB Recombinant Variant

Background SARS-CoV-2 has evolved to produce new variants causing successive waves of infection. Currently, six variants are being monitored by the World Health Organization that are replacing BA.5. These include BF.7 (BA.5 + R346T in spike), BQ.1 (and BQ.1.1, with BA.5 + R346T, K444T, N460K mutations in spike), BA.2.75 (including BA.2.75.2 and CH.1.1), and XBB (including XBB.1.5). BQ.1 and XBB variants are more immune evasive and have spread quickly throughout the world. Concerning the potential severity of infections caused by these variants, the present study describes the clinical characteristics and outcomes of these major variants in Maharashtra. Methodology A total of 1,141 reverse transcriptase-polymerase chain reaction (RT-PCR)-positive SARS-CoV-2 samples, with a cycle threshold (Ct) value of less than 25, were processed for SARS-CoV-2 whole genome sequencing between July 10, 2022, and January 12, 2023. All corresponding demographic and clinical data were recorded and analyzed using Microsoft® Excel and Epi Info™. Results Out of the 1,141 samples sequenced, BA.2.75* (63.78%) was the predominant Omicron variant, followed by the XBB* (18.88%), BA.2.38* (4.94%), BA.5* (4.06%), BA.2.10* (3.51%), and BQ.1* (1.65%). A total of 540 cases were contacted telephonically, of whom 494 (91.48%) were symptomatic with mild symptoms. Fever (77.73%) was the most common symptom, followed by cold (47.98%), cough (42.31%), and myalgia and fatigue (18.83%). Of the 540 cases, 414 (76.67%) cases recovered at home, and 126 (23.33%) were institutionally quarantined/hospitalized. Among the home-isolated and hospitalized cases, 416 (99.76%) and 108 (87.80%), respectively, recovered with symptomatic treatment, while one (0.24%) and 15 (12.20%), respectively, succumbed to the disease. Out of the 540 cases, 491 (90.93%) were vaccinated with at least one dose of the COVID-19 vaccine, 41 (7.59%) were unvaccinated, and for eight (1.48%) cases, vaccination data was not available. Conclusions The current study indicates that the XBB* variant is causing mild disease in India. However, as XBB* possesses both immune-escape and infectivity-enhancing mutations, it has the potential to spread to other parts of the world rapidly. Further, anti-SARS-CoV-2 vaccination improves survival rates in COVID-19.

Histidinal-Based Potent Antimalarial Agents

Herein we report the synthesis and evaluation of peptide-histidinal conjugated drug scaffolds, which have the potential to target the hemoglobin-degrading proteases falcipain-2/3 from the human malaria parasite. Scaffolds with various substitutions were tested for antimalarial activity, and compounds 8 g, 8 h, and 15 exhibited EC50 values of ∼0.018 μM, ∼0.069 μM, and ∼0.02 μM, respectively. Structure-based docking studies on falcipain-2/3 proteases (PDB:2GHU and PDB:3BWK) revealed that compounds 8 g, 8 h, and 15 interact strongly with binding sites of falcipain-2/3 in a substrate-like manner. In silico ADME studies revealed that the molecules of interest showed no or minimal violations of drug-likeness parameters. Further, phenotypic assays revealed that compound 8 g and its biotinylated version inhibit hemoglobin degradation in the parasite food vacuole. The identification of falcipain-2/3 targeting potent inhibitors of the malaria parasite can serve as a starting point for the development of lead compounds as future antimalarial drug candidates.

Chromodomain Protein Interacts with H3K9me3 and Controls RBC Rosette Formation by Regulating the Expression of a Subset of RIFINs in the Malaria Parasite

Plasmodium falciparum expresses clonally variant proteins on the surface of infected erythrocytes to evade the host immune system. The clonally variant multigene families include var, rifin, and stevor, which express Erythrocyte Membrane Protein 1 (EMP1), Repetitive Interspersed Families of polypeptides (RIFINs), and Sub-telomeric Variable Open Reading frame (STEVOR) proteins, respectively. The rifins are the largest multigene family and are essentially involved in the RBC rosetting, the hallmark of severe malaria. The molecular regulators that control the RIFINs expression in Plasmodium spp. have not been reported so far. This study reports a chromodomain-containing protein (PfCDP) that binds to H3K9me3 modification on P. falciparum chromatin. Conditional deletion of the chromodomain (CD) gene in P. falciparum using an inducible DiCre-LoxP system leads to selective up-regulation of a subset of virulence genes, including rifins, a few var, and stevor genes. Further, we show that PfCDP conditional knockout (PfΔCDP) promotes RBC rosette formation. This study provides the first evidence of an epigenetic regulator mediated control on a subset of RIFINs expression and RBC rosetting by P. falciparum.

Pervasive sequence-level variation in the transcriptome of Plasmodium falciparum

Single-nucleotide variations (SNVs) in RNA, arising from co- and post-transcriptional phenomena including transcription errors and RNA-editing, are well studied in a range of organisms. In the malaria parasite Plasmodium falciparum, stage-specific and non-specific gene-expression variations accompany the parasite's array of developmental and morphological phenotypes over the course of its complex life cycle. However, the extent, rate and effect of sequence-level variation in the parasite's transcriptome are unknown. Here, we report the presence of pervasive, non-specific SNVs in the P. falciparum transcriptome. SNV rates for a gene were correlated to gene length (r[Formula: see text]0.65-0.7) but not to the AT-content of that gene. Global SNV rates for the P. falciparum lines we used, and for publicly available P. vivax and P. falciparum clinical isolate datasets, were of the order of 10-3 per base, ∼10× higher than rates we calculated for bacterial datasets. These variations may reflect an intrinsic transcriptional error rate in the parasite, and RNA editing may be responsible for a subset of them. This seemingly characteristic property of the parasite may have implications for clinical outcomes and the basic biology and evolution of P. falciparum and parasite biology more broadly. We anticipate that our study will prompt further investigations into the exact sources, consequences and possible adaptive roles of these SNVs.

Identification of Co-Existing Mutations and Gene Expression Trends Associated With K13-Mediated Artemisinin Resistance in Plasmodium falciparum

Plasmodium falciparum infects millions and kills thousands of people annually the world over. With the emergence of artemisinin and/or multidrug resistant strains of the pathogen, it has become even more challenging to control and eliminate the disease. Multiomics studies of the parasite have started to provide a glimpse into the confounding genetics and mechanisms of artemisinin resistance and identified mutations in Kelch13 (K13) as a molecular marker of resistance. Over the years, thousands of genomes and transcriptomes of artemisinin-resistant/sensitive isolates have been documented, supplementing the search for new genes/pathways to target artemisinin-resistant isolates. This meta-analysis seeks to recap the genetic landscape and the transcriptional deregulation that demarcate artemisinin resistance in the field. To explore the genetic territory of artemisinin resistance, we use genomic single-nucleotide polymorphism (SNP) datasets from 2,517 isolates from 15 countries from the MalariaGEN Network (The Pf3K project, pilot data release 4, 2015) to dissect the prevalence, geographical distribution, and co-existing patterns of genetic markers associated with/enabling artemisinin resistance. We have identified several mutations which co-exist with the established markers of artemisinin resistance. Interestingly, K13-resistant parasites harbor α-ß hydrolase and putative HECT domain-containing protein genes with the maximum number of SNPs. We have also explored the multiple, publicly available transcriptomic datasets to identify genes from key biological pathways whose consistent deregulation may be contributing to the biology of resistant parasites. Surprisingly, glycolytic and pentose phosphate pathways were consistently downregulated in artemisinin-resistant parasites. Thus, this meta-analysis highlights the genetic and transcriptomic features of resistant parasites to propel further exploratory studies in the community to tackle artemisinin resistance.

Nup93 and CTCF modulate spatiotemporal dynamics and function of the HOXA gene locus during differentiation

Nucleoporins regulate nuclear transport and are also involved in DNA damage, repair, cell cycle, chromatin organization, and gene expression. Here, we studied the role of nucleoporin Nup93 and the chromatin organizer CTCF in regulating HOXA expression during differentiation. ChIP sequencing revealed a significant overlap between Nup93 and CTCF peaks. Interestingly, Nup93 and CTCF are associated with the 3' and 5'HOXA genes respectively. Depletions of Nup93 and CTCF antagonistically modulate expression levels of 3'and 5'HOXA genes in undifferentiated NT2/D1 cells. Nup93 also regulates the localization of the HOXA gene locus, which disengages from the nuclear periphery upon Nup93 but not CTCF depletion, consistent with its upregulation. The dynamic association of Nup93 and CTCF with the HOXA locus during differentiation correlates with its spatial positioning and expression. While Nup93 tethers the HOXA locus to the nuclear periphery, CTCF potentially regulates looping of the HOXA gene cluster in a temporal manner. In summary, Nup93 and CTCF complement one another in modulating the spatiotemporal dynamics and function of the HOXA gene locus during differentiation.

Single-Cell RNA Sequencing Reveals Cellular Heterogeneity and Stage Transition under Temperature Stress in Synchronized Plasmodium falciparum Cells

The malaria parasite has a complex life cycle exhibiting phenotypic and morphogenic variations in two different hosts by existing in heterogeneous developmental states. To investigate this cellular heterogeneity of the parasite within the human host, we performed single-cell RNA sequencing of synchronized Plasmodium cells under control and temperature treatment conditions. Using the Malaria Cell Atlas (https://www.sanger.ac.uk/science/tools/mca) as a guide, we identified 9 subtypes of the parasite distributed across known intraerythrocytic stages. Interestingly, temperature treatment results in the upregulation of the AP2-G gene, the master regulator of sexual development in a small subpopulation of the parasites. Moreover, we identified a heterogeneous stress-responsive subpopulation (clusters 5, 6, and 7 [∼10% of the total population]) that exhibits upregulation of stress response pathways under normal growth conditions. We also developed an online exploratory tool that will provide new insights into gene function under normal and temperature stress conditions. Thus, our study reveals important insights into cell-to-cell heterogeneity in the parasite population under temperature treatment that will be instrumental toward a mechanistic understanding of cellular adaptation and population dynamics in Plasmodium falciparum. IMPORTANCE The malaria parasite has a complex life cycle exhibiting phenotypic variations in two different hosts accompanied by cell-to-cell variability that is important for stress tolerance, immune evasion, and drug resistance. To investigate cellular heterogeneity determined by gene expression, we performed single-cell RNA sequencing (scRNA-seq) of about 12,000 synchronized Plasmodium cells under physiologically relevant normal (37°C) and temperature stress (40°C) conditions phenocopying the cyclic bouts of fever experienced during malarial infection. In this study, we found that parasites exhibit transcriptional heterogeneity in an otherwise morphologically synchronized culture. Also, a subset of parasites is continually committed to gametocytogenesis and stress-responsive pathways. These observations have important implications for understanding the mechanisms of drug resistance generation and vaccine development against the malaria parasite.

Dynamic association of the H3K64 trimethylation mark with genes encoding exported proteins in Plasmodium falciparum

Epigenetic modifications have emerged as critical regulators of virulence genes and stage-specific gene expression in Plasmodium falciparum. However, the specific roles of histone core epigenetic modifications in regulating the stage-specific gene expression are not well understood. In this study, we report an unconventional trimethylation at lysine 64 on histone 3 (H3K64me3) and characterize its functional relevance in P. falciparum. We show that PfSET4 and PfSET5 proteins of P. falciparum methylate H3K64 and that they prefer the nucleosome as a substrate over free histone 3 proteins. Structural analysis of PfSET5 revealed that it interacts with the nucleosome as a dimer. The H3K64me3 mark is dynamic, being enriched in the ring and trophozoite stages and drastically reduced in the schizont stages. Stage-specific global chromatin immunoprecipitation –sequencing analysis of the H3K64me3 mark revealed the selective enrichment of this methyl mark on the genes of exported family proteins in the ring and trophozoite stages and a significant reduction of the same in the schizont stages. Collectively, our data identify a novel epigenetic mark that is associated with the subset of genes encoding for exported proteins, which may regulate their expression in different stages of P. falciparum.

Histone acetyltransferase PfGCN5 regulates stress responsive and artemisinin resistance related genes in Plasmodium falciparum

Plasmodium falciparum has evolved resistance to almost all front-line drugs including artemisinin, which threatens malaria control and elimination strategies. Oxidative stress and protein damage responses have emerged as key players in the generation of artemisinin resistance. In this study, we show that PfGCN5, a histone acetyltransferase, binds to the stress-responsive genes in a poised state and regulates their expression under stress conditions. Furthermore, we show that upon artemisinin exposure, genome-wide binding sites for PfGCN5 are increased and it is directly associated with the genes implicated in artemisinin resistance generation like BiP and TRiC chaperone. Interestingly, expression of genes bound by PfGCN5 was found to be upregulated during stress conditions. Moreover, inhibition of PfGCN5 in artemisinin-resistant parasites increases the sensitivity of the parasites to artemisinin treatment indicating its role in drug resistance generation. Together, these findings elucidate the role of PfGCN5 as a global chromatin regulator of stress-responses with a potential role in modulating artemisinin drug resistance and identify PfGCN5 as an important target against artemisinin-resistant parasites.

Analysis of drug resistance marker genes of Plasmodium falciparum after implementation of artemisinin-based combination therapy in Pune district, India

The global emergence and spread of malaria parasites resistant to antimalarial drugs is a major problem in malaria control and elimination. In this study, samples from Pune district were characterized to determine prevalence of molecular markers of resistance to chloroquine (pfcrt codons C72S, M74I, N75E, K76T and pfmdr-1 N86Y, Y184F), pyrimethamine (pfdhfr C50R, N51I, C59R, S108N), sulfadoxine (pfdhps, S436A, A437G, K540E, A581G), and artemisinin (pfkelch13, C580Y, R539T). The pfcrt K76T mutation was found in 78% samples as CVMNT, SVMNT and CVIET haplotype. The pfmdr-1 N86Y and Y184F mutations were found in 54% of samples. The pfdhfr double mutation C59R + S108N was present in 67% of samples, while the pfdhfr triple mutation (N51I + C59R + S108N) was not detected. The pfdhps mutations A437G and K540E were found in 67% of samples. Single mutants of pfdhps were rare, with K540E detected in only 6 patient samples. Similarly, pfdhps A581G was found in 13 of the isolates. The molecular markers associated with artemisinin resistance (mutations in pfkelch13 C580Y, R539T) were not detected in any of the isolates. These results suggest an emerging problem with multidrug-resistant P. falciparum. Though the genotype conventionally associated with artemisinin resistance was not observed, chloroquine-resistant genotype has reached complete fixation in the population. Moreover, the prevalence of mutations in both pfdhfr and pfdhps, with the presence of the quadruple mutant, indicates that continued monitoring is required to assess whether sulfadoxine-pyrimethamine can be used efficiently as a partner drug for artemisinin for the treatment of P. falciparum.

Role of PfGCN5 in nutrient sensing and transcriptional regulation in Plasmodium falciparum

Malaria is a deadly, infectious disease caused by the parasite Plasmodium, leading to millions of deaths worldwide. Plasmodium requires a coordinated pattern of sequential gene expression for surviving in both invertebrate and vertebrate host environments. As parasites largely depend on host resources, they also develop efficient mechanisms to sense and adapt to variable nutrient conditions in the environment and modulate their virulence. Earlier we have shown that PfGCN5, a histone acetyltransferase, binds to the stress-responsive and virulence-related genes in a poised state and regulates their expression under temperature and artemisinin treatment conditions in P. falciparum. In this study, we show upregulation of PfGCN5 upon nutrient stress condition. With the help of chromatin immunoprecipitation coupled high-throughput sequencing (ChIP-seq) and transcriptomic (RNA-sequencing) analyses, we show that PfGCN5 is associated with the genes that are important for the maintenance of parasite cellular homeostasis upon nutrient stress condition. Furthermore, we identified various metabolic enzymes as interacting partners of PfGCN5 by immunoprecipitation coupled with mass spectroscopy, possibly acting as a sensor of nutrient conditions in the environment. We also demonstrated that PfGCN5 interacts and acetylates PfGAPDH in vitro. Collectively, our data provides important insights into transcriptional deregulation upon nutrient stress condition and elucidate the role of PfGCN5 during nutrient stress condition.

Peroxidation of 2-oxindole and barbituric acid derivatives under batch and continuous flow using an eco-friendly ethyl acetate solvent

We have demonstrated the magnetically retrievable Fe(OH)₃Fe₃O₄catalyzed C–H peroxidation of 2-oxindole and barbituric acid derivatives under batch and continuous flow process for the first time.

Ru-Catalyzed dehydrogenative synthesis of antimalarial arylidene oxindoles

Ru(ii)-NHC catalyzes α-olefination of 2-oxindoles using diaryl methanols in the absence of an acceptor. A wide array of symmetrical and unsymmetrical diaryl methanols undergoes dehydrogenative coupling with 2-oxindole selectively to generate various substituted 3-(diphenylmethylene)indolin-2-one derivatives in good yields and produces environmentally benign by-products, H2 and H2O. This methodology was successfully applied for the synthesis of a bioactive drug i.e. TAS-301. The biological activities of the synthesized 3-(diphenylmethylene)indolin-2-one derivatives were screened against the Plasmodium falciparum parasite and found to exhibit a significant activity with IC50 = 2.24 μM.

Isolation and structure elucidation of halymeniaol, a new antimalarial sterol derivative from the red alga Halymenia floresii

A new mono-hydroxy acetylated sterol derivative: 12β-hydroxy-3β, 15α, 16β-triacetoxy-cholest-5-en-7-one (halymeniaol) (1), and cholesterol (2) were isolated from the marine red alga Halymenia floresii. The structure of the compound 1 (halymeniaol) was established from its spectral data, derived from HRMS/MS and 2D NMR. Compound 1 exhibited growth inhibitory activity against chloroquine-resistant Plasmodium falciparum 3D7 strain with an IC₅₀ of 3.0 μM.

Genome-wide identification of novel intergenic enhancer-like elements: implications in the regulation of transcription in Plasmodium falciparum

Background

The molecular mechanisms of transcriptional regulation are poorly understood in Plasmodium falciparum. In addition, most of the genes in Plasmodium falciparum are transcriptionally poised and only a handful of cis-regulatory elements are known to operate in transcriptional regulation. Here, we employed an epigenetic signature based approach to identify significance of previously uncharacterised intergenic regions enriched with histone modification marks leading to discovery of enhancer-like elements.

Results

We found that enhancer-like elements are significantly enriched with H3K4me1, generate unique non-coding bi-directional RNAs and majority of them can function as cis-regulators. Furthermore, functional enhancer reporter assay demonstrates that the enhancer-like elements regulate transcription of target genes in Plasmodium falciparum. Our study also suggests that the Plasmodium genome segregates functionally related genes into discrete housekeeping and pathogenicity/virulence clusters, presumably for robust transcriptional control of virulence/pathogenicity genes.

Conclusions

This report contributes to the understanding of parasite regulatory genomics by identification of enhancer-like elements, defining their epigenetic and transcriptional features and provides a resource of functional cis-regulatory elements that may give insights into the virulence/pathogenicity of Plasmodium falciparum.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4052-4) contains supplementary material, which is available to authorized users.

HOXA repression is mediated by nucleoporin Nup93 assisted by its interactors Nup188 and Nup205

Background

The nuclear pore complex (NPC) mediates nuclear transport of RNA and proteins into and out of the nucleus. Certain nucleoporins have additional functions in chromatin organization and transcription regulation. Nup93 is a scaffold nucleoporin at the nuclear pore complex which is associated with human chromosomes 5, 7 and 16 and with the promoters of the HOXA gene as revealed by ChIP-on-chip studies using tiling microarrays for these chromosomes. However, the functional consequences of the association of Nup93 with HOXA is unknown.

Results

Here, we examined the association of Nup93 with the HOXA gene cluster and its consequences on HOXA gene expression in diploid colorectal cancer cells (DLD1). Nup93 showed a specific enrichment ~1 Kb upstream of the transcription start site of each of the HOXA1, HOXA3 and HOXA5 promoters, respectively. Furthermore, the association of Nup93 with HOXA was assisted by its interacting partners Nup188 and Nup205. The depletion of the Nup93 sub-complex significantly upregulated HOXA gene expression levels. However, expression levels of a control gene locus (GLCCI1) on human chromosome 7 were unaffected. Three-dimensional fluorescence in situ hybridization (3D-FISH) analyses revealed that the depletion of the Nup93 sub-complex (but not Nup98) disengages the HOXA gene locus from the nuclear periphery, suggesting a potential role for Nup93 in tethering and repressing the HOXA gene cluster. Consistently, Nup93 knockdown increased active histone marks (H3K9ac), decreased repressive histone marks (H3K27me3) on the HOXA1 promoter and increased transcription elongation marks (H3K36me3) within the HOXA1 gene. Moreover, the combined depletion of Nup93 and CTCF (a known organizer of HOXA gene cluster) but not Nup93 alone, significantly increased GLCCI1 gene expression levels. Taken together, this suggests a novel role for Nup93 and its interactors in repressing the HOXA gene cluster.

Conclusions

This study reveals that the nucleoporin Nup93 assisted by its interactors Nup188 and Nup205 mediates the repression of HOXA gene expression.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0106-0) contains supplementary material, which is available to authorized users.

Plasmodium falciparum epigenome: A distinct dynamic epigenetic regulation of gene expression

Histone modification profiles are predictive of gene expression and most of the knowledge gained is acquired through studies done in higher eukaryotes. However, genome-wide studies involving Plasmodium falciparum, the causative agent of malaria, have been rather few, at lower resolution (mostly using ChIP-on-chip), and covering limited number of histone modifications. In our recent study [1], we have performed extensive genome-wide analyses of multiple histone modifications including the active (H3K4me2, H3K4me3, H3K9ac, H3K14ac, H3K27ac and H4ac), inactive (H3K9me3 and H3K27me3), elongation (H3K79me3) and regulatory element (H3K4me1) in a stage-specific manner. Furthermore, we used a ligation-based method suitable for sequencing homopolymeric stretches as seen in P. falciparum for next-generation sequencing library amplification [2], enabling highly quantitative analysis of the extremely AT-rich P. falciparum genome. Our recently published study suggests that transcription regulation by virtue of poised chromatin and differential histone modifications is unique to P. falciparum [1]. Here we describe the experiments, quality controls and chromatin immunoprecipitation-sequencing data analysis of our associated study published in Epigenetics and Chromatin [1]. Stage-specific ChIP-sequencing data for histone modifications is submitted to Gene Expression Omnibus (GEO) database under the accession number GSE63369.

A comprehensive epigenome map of Plasmodium falciparum reveals unique mechanisms of transcriptional regulation and identifies H3K36me2 as a global mark of gene suppression

Background

Role of epigenetic mechanisms towards regulation of the complex life cycle/pathogenesis of Plasmodium falciparum, the causative agent of malaria, has been poorly understood. To elucidate stage-specific epigenetic regulation, we performed genome-wide mapping of multiple histone modifications of P. falciparum. Further to understand the differences in transcription regulation in P. falciparum and its host, human, we compared their histone modification profiles.

Results

Our comprehensive comparative analysis suggests distinct mode of transcriptional regulation in malaria parasite by virtue of poised genes and differential histone modifications. Furthermore, analysis of histone modification profiles predicted 562 genes producing anti-sense RNAs and 335 genes having bidirectional promoter activity, which raises the intriguing possibility of RNA-mediated regulation of transcription in P. falciparum. Interestingly, we found that H3K36me2 acts as a global repressive mark and gene regulation is fine tuned by the ratio of activation marks to H3K36me2 in P. falciparum. This novel mechanism of gene regulation is supported by the fact that knockout of SET genes (responsible for H3K36 methylation) leads to up-regulation of genes with highest occupancy of H3K36me2 in wild-type P. falciparum. Moreover, virulence (var) genes are mostly poised and marked by a unique set of activation (H4ac) and repression (H3K9me3) marks, which are mutually exclusive to other Plasmodium housekeeping genes.

Conclusions

Our study reveals unique plasticity in the epigenetic regulation in P. falciparum which can influence parasite virulence and pathogenicity. The observed differences in the histone code and transcriptional regulation in P. falciparum and its host will open new avenues for epigenetic drug development against malaria parasite.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-015-0029-1) contains supplementary material, which is available to authorized users.

Camello, a novel family of Histone Acetyltransferases that acetylate histone H4 and is essential for zebrafish development

In this study, we have investigated genome-wide occurrence of Histone Acetyltransferases (HATs) in genomes of Mus musculus and Danio rerio on the basis of presence of HAT domain. Our study identified a group of proteins that lacks characteristic features of known HAT families, relatively smaller in size and has no other associated domains. Most of the proteins in this unclassified group are Camello proteins, which are not yet known and classified as functional HATs. Our in vitro and in vivo analysis revealed that Camello family proteins are active HATs and exhibit specificity towards histone H4. Interestingly, Camello proteins are among the first identified HATs showing perinuclear localization. Moreover, Camello proteins are evolutionarily conserved in all chordates and are observed for the first time in cnidarians in phylogeny. Furthermore, knockdown of Camello protein (CMLO3) in zebrafish embryos exhibited defects in axis elongation and head formation. Thus, our study identified a novel family of active HATs that is specific for histone H4 acetylation, exhibits perinuclear localization and is essential for zebrafish development.

H3K9 and H3K14 acetylation co-occur at many gene regulatory elements, while H3K14ac marks a subset of inactive inducible promoters in mouse embryonic stem cells

Background

Transcription regulation in pluripotent embryonic stem (ES) cells is a complex process that involves multitude of regulatory layers, one of which is post-translational modification of histones. Acetylation of specific lysine residues of histones plays a key role in regulating gene expression.

Results

Here we have investigated the genome-wide occurrence of two histone marks, acetylation of histone H3K9 and K14 (H3K9ac and H3K14ac), in mouse embryonic stem (mES) cells. Genome-wide H3K9ac and H3K14ac show very high correlation between each other as well as with other histone marks (such as H3K4me3) suggesting a coordinated regulation of active histone marks. Moreover, the levels of H3K9ac and H3K14ac directly correlate with the CpG content of the promoters attesting the importance of sequences underlying the specifically modified nucleosomes. Our data provide evidence that H3K9ac and H3K14ac are also present over the previously described bivalent promoters, along with H3K4me3 and H3K27me3. Furthermore, like H3K27ac, H3K9ac and H3K14ac can also differentiate active enhancers from inactive ones. Although, H3K9ac and H3K14ac, a hallmark of gene activation exhibit remarkable correlation over active and bivalent promoters as well as distal regulatory elements, a subset of inactive promoters is selectively enriched for H3K14ac.

Conclusions

Our study suggests that chromatin modifications, such as H3K9ac and H3K14ac, are part of the active promoter state, are present over bivalent promoters and active enhancers and that the extent of H3K9 and H3K14 acetylation could be driven by cis regulatory elements such as CpG content at promoters. Our study also suggests that a subset of inactive promoters is selectively and specifically enriched for H3K14ac. This observation suggests that histone acetyl transferases (HATs) prime inactive genes by H3K14ac for stimuli dependent activation. In conclusion our study demonstrates a wider role for H3K9ac and H3K14ac in gene regulation than originally thought.

SAGA and ATAC histone acetyl transferase complexes regulate distinct sets of genes and ATAC defines a class of p300-independent enhancers

Histone acetyltransferase (HAT) complexes are coactivators that are important for transcriptional activation by modifying chromatin. Metazoan SAGA and ATAC are distinct multisubunits complexes that share the same catalytic HAT subunit (GCN5 or PCAF). Here, we show that these human HAT complexes are targeted to different genomic loci representing functionally distinct regulatory elements both at broadly expressed and tissue-specific genes. While SAGA can principally be found at promoters, ATAC is recruited to promoters and enhancers, yet only its enhancer binding is cell-type specific. Furthermore, we show that ATAC functions at a set of enhancers that are not bound by p300, revealing a class of enhancers not yet identified. These findings demonstrate important functional differences between SAGA and ATAC coactivator complexes at the level of the genome and define a role for the ATAC complex in the regulation of a set of enhancers.

ATAC and Mediator coactivators form a stable complex and regulate a set of non-coding RNA genes

The Ada-Two-A-containing (ATAC) histone acetyltransferase and Mediator coactivator complexes regulate independent and distinct steps during transcription initiation and elongation. Here, we report the identification of a new stable molecular assembly formed between the ATAC and Mediator complexes in mouse embryonic stem cells. Moreover, we identify leucine zipper motif-containing protein 1 as a subunit of this meta-coactivator complex (MECO). Finally, we demonstrate that the MECO regulates a subset of RNA polymerase II-transcribed non-coding RNA genes. Our findings establish that transcription coactivator complexes can form stable subcomplexes to facilitate their combined actions on specific target genes.

Deciphering the key residues in Plasmodium falciparum beta-ketoacyl acyl carrier protein reductase responsible for interactions with Plasmodium falciparum acyl carrier protein

The type II fatty acid synthase (FAS) pathway of Plasmodium falciparum is a validated unique target for developing novel antimalarials, due to its intrinsic differences from the typeI pathway operating in humans. beta-Ketoacyl acyl carrier protein (ACP) reductase (FabG) performs the NADPH-dependent reduction of beta-ketoacyl-ACP to beta-hydroxyacyl-ACP, the first reductive step in the elongation cycle of fatty acid biosynthesis. In this article, we report intensive studies on the direct interactions of Plasmodium FabG and Plasmodium ACP in solution, in the presence and absence of its cofactor, NADPH, by monitoring the change in intrinsic fluorescence of P.falciparum FabG (PfFabG) and by surface plasmon resonance. To address the issue of the importance of the residues involved in strong, specific and stoichiometric binding of PfFabG to P.falciparum ACP (PfACP), we mutated Arg187, Arg190 and Arg230 of PfFabG. The activities of the mutants were assessed using both an ACP-dependent and an ACP-independent assay. The affinities of all the PfFabG mutants for acetoacetyl-ACP (the physiological substrate) were reduced to different extents as compared to wild-type PfFabG, but were equally active in biochemical assays with the substrate analog acetoacetyl-CoA. Kinetic analysis and studies of direct binding between PfFabG and PfACP confirmed the identification of Arg187 and Arg230 as critical residues for the PfFabG-PfACP interactions. Our studies thus reveal the significance of the positively charged/hydrophobic patch located adjacent to the active site cavities of PfFabG for interactions with PfACP.

A unique and differential effect of denaturants on cofactor mediated activation of Plasmodium falciparum β-ketoacyl-ACP reductase

The urea and guanidinium chloride (GdmCl) induced unfolding of FabG, a beta-ketoacyl-ACP reductase of Plasmodium falciparum, was examined in detail using intrinsic fluorescence of FabG, UV-circular dichroism (CD), spectrophotometric enzyme activity measurements, glutaraldehyde cross-linking, and size exclusion chromatography. The equilibrium unfolding of FabG by urea is a multistep process as compared with a two-state process by GdmCl. FabG is fully unfolded at 6.0M urea and 4.0M GdmCl. Approximately 90% of the enzyme activity could be recovered on dialyzing the denaturants, showing that denaturation by both urea and GdmCl is reversible. We found two states in the reversible unfolding process of FabG in presence of NADPH; one is an activity-enhanced state and the other, an inactive state in case of equilibrium unfolding with urea. On the contrary, in presence of NADPH, there is no stabilization of FabG in case of equilibrium unfolding with GdmCl. We hypothesize that the hydrogen-bonding network may be reorganized by the denaturant in the activity-enhanced state formed in presence of 1.0M urea, by interrupting the association between dimer-dimer interface and help in accommodating the larger substrate in the substrate binding tunnel thus, increasing the activity. Furthermore, binding of the active site organizer, NADPH leads to compaction of the FabG in presence of urea, as evident by acrylamide quenching. We have shown here for the first time, the detailed inactivation kinetics of FabG, which have not been evaluated in the past from any of the FabG family of enzymes from any of the other sources. These findings provide impetus for exploring the influences of ligands on the structure-activity relationship of Plasmodium beta-ketoacyl-ACP reductase.

Mass spectrometry-based systems approach for identification of inhibitors of Plasmodium falciparum fatty acid synthase

The emergence of strains of Plasmodium falciparum resistant to the commonly used antimalarials warrants the development of new antimalarial agents. The discovery of type II fatty acid synthase (FAS) in Plasmodium distinct from the FAS in its human host (type I FAS) opened up new avenues for the development of novel antimalarials. The process of fatty acid synthesis takes place by iterative elongation of butyryl-acyl carrier protein (butyryl-ACP) by two carbon units, with the successive action of four enzymes constituting the elongation module of FAS until the desired acyl length is obtained. The study of the fatty acid synthesis machinery of the parasite inside the red blood cell culture has always been a challenging task. Here, we report the in vitro reconstitution of the elongation module of the FAS of malaria parasite involving all four enzymes, FabB/F (beta-ketoacyl-ACP synthase), FabG (beta-ketoacyl-ACP reductase), FabZ (beta-ketoacyl-ACP dehydratase), and FabI (enoyl-ACP reductase), and its analysis by matrix-assisted laser desorption-time of flight mass spectrometry (MALDI-TOF MS). That this in vitro systems approach completely mimics the in vivo machinery is confirmed by the distribution of acyl products. Using known inhibitors of the enzymes of the elongation module, cerulenin, triclosan, NAS-21/91, and (-)-catechin gallate, we demonstrate that accumulation of intermediates resulting from the inhibition of any of the enzymes can be unambiguously followed by MALDI-TOF MS. Thus, this work not only offers a powerful tool for easier and faster throughput screening of inhibitors but also allows for the study of the biochemical properties of the FAS pathway of the malaria parasite.

Discovery of a rhodanine class of compounds as inhibitors of Plasmodium falciparum enoyl-acyl carrier protein reductase

Enoyl acyl carrier protein (ACP) reductase, one of the enzymes of the type II fatty acid biosynthesis pathway, has been established as a promising target for the development of new drugs for malaria. Here we present the discovery of a rhodanine (2-thioxothiazolidin-4-one) class of compounds as inhibitors of this enzyme using a combined approach of rational selection of compounds for screening, analogue search, docking studies, and lead optimization. The most potent inhibitor exhibits an IC(50) of 35.6 nM against Plasmodium falciparum enoyl ACP reductase (PfENR) and inhibits growth of the parasite in red blood cell cultures at an IC(50) value of 750 nM. Many more compounds of this class were found to inhibit PfENR at low nanomolar to low micromolar concentrations, expanding the scope for developing new antimalarial drugs. The structure-activity relationship of these rhodanine compounds is discussed.

15-Deoxyspergualin Primarily Targets the Trafficking of Apicoplast Proteins in Plasmodium falciparum

15-Deoxyspergualin, an immunosuppressant with tumoricidal and antimalarial properties, has been implicated in the inhibition of a diverse array of cellular processes including polyamine synthesis and protein synthesis. Endeavoring to identify the mechanism of antimalarial action of this molecule, we examined its effect on Plasmodium falciparum protein synthesis, polyamine biosynthesis, and transport. 15-Deoxyspergualin stalled protein synthesis in P. falciparum through Hsp70 sequestration and subsequent phosphorylation of the eukaryotic initiation factor eIF2alpha. However, protein synthesis inhibition as well as polyamine depletion were invoked only by high micromolar concentrations of 15-deoxyspergualin, in contrast to the submicromolar concentrations sufficient to inhibit parasite growth. Further investigations demonstrated that 15-deoxyspergualin in the malaria parasite primarily targets the hitherto underexplored process of trafficking of nucleus-encoded proteins to the apicoplast. Our finding that 15-deoxyspergualin kills the malaria parasite by interfering with targeting of nucleus-encoded proteins to the apicoplast not only exposes a chink in the armor of the malaria parasite, but also reveals new realms in our endeavors to study this intriguing biological process.

Conformational stability and thermodynamic characterization of homotetrameric Plasmodium falciparum β-ketoacyl-ACP reductase

The conformational stability of the homotetrameric Plasmodium falciparum beta-ketoacyl-ACP reductase (FabG) was determined by guanidinium chloride-induced isothermal and thermal denaturation. The reversible unfolding transitions were monitored by intrinsic fluorescence, circular dichroism (CD) spectroscopy and by measuring the enzyme activity of FabG. The denaturation profiles were analyzed to obtain the thermodynamic parameters associated with unfolding of the protein. The data confirm the simple A(4) <--> 4A model of unfolding, based on the corroboration of CD data by fluorescence transition and similar Delta G estimation for denaturation curves obtained at four different concentration of the FabG. Denaturation is well described by the linear extrapolation model for denaturant-protein interactions. In addition, the conformational stability (Delta G(s)) as well as the Delta C(p) for the protein unfolding is quite high, 22.68 kcal/mole and 5.83 kcal/(mole K), respectively, which may be a reflection of the relatively large size of the tetrameric molecule (Mr 120, 000) and a large buried hydrophobic core in the folded protein. This study provides a prototype for determining conformational stability of other members of the short-chain alcohol dehydrogenase/reductase superfamily of proteins to which PfFabG belongs.

Inhibitors of nonhousekeeping functions of the apicoplast defy delayed death in Plasmodium falciparum

Targeting of apicoplast replication and protein synthesis in the apicomplexan Toxoplasma gondii has conventionally been associated with the typical "delayed death" phenotype, characterized by the death of parasites only in the generation following drug intervention. We demonstrate that antibiotics like clindamycin, chloramphenicol, and tetracycline, inhibitors of prokaryotic protein synthesis, invoke the delayed death phenotype in Plasmodium falciparum, too, as evident from a specific reduction of apicoplast genome copy number. Interestingly, however, molecules like triclosan, cerulenin, fops, and NAS-91, inhibitors of the recently discovered fatty acid synthesis pathway, and succinyl acetone, an inhibitor of heme biosynthesis that operates in the apicoplast of the parasite, display rapid and striking parasiticidal effects. Our results draw a clear distinction between apicoplast functions per se and the apicoplast as the site of metabolic pathways, which are required for parasite survival, and thus subserve the development of novel antimalarial therapy.

Analyses of co-operative transitions in Plasmodium falciparum beta-ketoacyl acyl carrier protein reductase upon co-factor and acyl carrier protein binding

The type II fatty acid synthase pathway of Plasmodium falciparum is a validated unique target for developing novel antimalarials because of its intrinsic differences from the type I pathway operating in humans. beta-Ketoacyl-acyl carrier protein reductase is the only enzyme of this pathway that has no isoforms and thus selective inhibitors can be developed for this player of the pathway. We report here intensive studies on the direct interactions of Plasmodiumbeta-ketoacyl-acyl carrier protein reductase with its cofactor, NADPH, acyl carrier protein, acetoacetyl-coenzyme A and other ligands in solution, by monitoring the intrinsic fluorescence (lambdamax 334 nM) of the protein as a result of its lone tryptophan, as well as the fluorescence of NADPH (lambdamax 450 nM) upon binding to the enzyme. Binding of the reduced cofactor makes the enzyme catalytically efficient, as it increases the binding affinity of the substrate, acetoacetyl-coenzyme A, by 16-fold. The binding affinity of acyl carrier protein to the enzyme also increases by approximately threefold upon NADPH binding. Plasmodiumbeta-ketoacyl-acyl carrier protein reductase exhibits negative, homotropic co-operative binding for NADPH, which is enhanced in the presence of acyl carrier protein. Acyl carrier protein increases the accessibility of NADPH to beta-ketoacyl-acyl carrier protein reductase, as evident from the increase in the accessibility of the tryptophan of beta-ketoacyl-acyl carrier protein reductase to acrylamide, from 81 to 98%. In the presence of NADP+, the reaction proceeds in the reverse direction (Ka=23.17 microM-1). These findings provide impetus for exploring the influence of ligands on the structure-activity relationship of Plasmodiumbeta-ketoacyl-acyl carrier protein reductase.

Production and purification of refolded recombinant Plasmodium falciparum β-ketoacyl-ACP reductase from inclusion bodies

A recombinant form of Plasmodium falciparum beta-ketoacyl-ACP reductase (PfFabG) was overexpressed in Escherichia coli BL-21 codon plus (DE3). The resulting insoluble inclusion bodies were separated from cellular debris by extensive washing with buffer containing 0.05% Tween 20 and solubilized by homogenization with 8 M urea. Attempts to refold PfFabG from solubilized inclusion bodies employing Rotofor (separation based on different pIs of proteins in a mixture) followed by Ni(2+) or cation exchange chromatography were not successful either by bringing down the urea concentration instantaneously, stepwise, or by dialysis. Denatured PfFabG was therefore initially purified by cation exchange chromatography and was then correctly refolded at a final concentration of 100-200 microg/ml in a 20 mM Na-acetate buffer, pH 5.3, with 300 mM NaCl, 10% glycerol, and 0.05% Tween 20. The protein was found to be properly folded only in the presence of the cofactor NADPH and salt at a concentration 300 mM by drop dilution method at 2-8 degrees C for 12 h. The purified final product was >98% pure by denaturing gel electrophoresis. The purified protein was biologically active in a standard enzymatic assay using acetoacetyl-CoA as a substrate. The enzyme was found to be stable up to fourth day of purification and glycerol was found to stabilize enzyme activity for several weeks, during storage. This effort paves the way for elucidation of the structure-function correlations for PfFabG as well as exploration of the enzyme for developing inhibitors against it for combating malaria.