Biochemical and functional characterization of Plasmodium falciparum GTP cyclohydrolase I

Functional profiling of the Toxoplasma genome during acute mouse infection

Within a host, pathogens encounter a diverse and changing landscape of cell types, nutrients, and immune responses. Examining host-pathogen interactions in animal models can therefore reveal aspects of infection absent from cell culture. We use CRISPR-based screens to functionally profile the entire genome of the model apicomplexan parasite Toxoplasma gondii during mouse infection. Barcoded gRNAs were used to track mutant parasite lineages, enabling detection of bottlenecks and mapping of population structures. We uncovered over 300 genes that modulate parasite fitness in mice with previously unknown roles in infection. These candidates span multiple axes of host-parasite interaction, including determinants of tropism, host organelle remodeling, and metabolic rewiring. We mechanistically characterized three novel candidates, including GTP cyclohydrolase I, against which a small-molecule inhibitor could be repurposed as an antiparasitic compound. This compound exhibited antiparasitic activity against T. gondii and Plasmodium falciparum, the most lethal agent of malaria. Taken together, we present the first complete survey of an apicomplexan genome during infection of an animal host, and point to novel interfaces of host-parasite interaction that may offer new avenues for treatment.

GTP Cyclohydrolase I as a Potential Drug Target: New Insights into Its Allosteric Modulation via Normal Mode Analysis

Guanosine triphosphate (GTP) cyclohydrolase I (GCH1) catalyzes the conversion of GTP into dihydroneopterin triphosphate (DHNP). DHNP is the first intermediate of the folate de novo biosynthesis pathway in prokaryotic and lower eukaryotic microorganisms and the tetrahydrobiopterin (BH4) biosynthesis pathway in higher eukaryotes. The de novo folate biosynthesis provides essential cofactors for DNA replication, cell division, and synthesis of key amino acids in rapidly replicating pathogen cells, such as Plasmodium falciparum (P. falciparum), a causative agent of malaria. In eukaryotes, the product of the BH4 biosynthesis pathway is essential for the production of nitric oxide and several neurotransmitter precursors. An increased copy number of the malaria parasite P. falciparum GCH1 gene has been reported to influence antimalarial antifolate drug resistance evolution, whereas mutations in the human GCH1 are associated with neuropathic and inflammatory pain disorders. Thus, GCH1 stands as an important and attractive drug target for developing therapeutics. The GCH1 intrinsic dynamics that modulate its activity remains unclear, and key sites that exert allosteric effects across the structure are yet to be elucidated. This study employed the anisotropic network model to analyze the intrinsic motions of the GCH1 structure alone and in complex with its regulatory partner protein. We showed that the GCH1 tunnel-gating mechanism is regulated by a global shear motion and an outward expansion of the central five-helix bundle. We further identified hotspot residues within sites of structural significance for the GCH1 intrinsic allosteric modulation. The obtained results can provide a solid starting point to design novel antineuropathic treatments for humans and novel antimalarial drugs against the malaria parasite P. falciparum GCH1 enzyme.

Polymorphism of Antifolate Drug Resistance in Plasmodium vivax From Local Residents and Migrant Workers Returned From the China-Myanmar Border

Drug-resistant Plasmodium vivax malaria impedes efforts to control, eliminate, and ultimately eradicate malaria in Southeast Asia. P. vivax resistance to antifolate drugs derives from point mutations in specific parasite genes, including the dihydropteroate synthase (pvdhps), dihydrofolate reductase (pvdhfr), and GTP cyclohydrolase I (pvgch1) genes. This study aims to investigate the prevalence and spread of drug resistance markers in P. vivax populating the China-Myanmar border. Blood samples were collected from symptomatic patients with acute P. vivax infection. Samples with single-clone P. vivax infections were sequenced for pvdhps and pvdhfr genes and genotyped for 6 flanking microsatellite markers. Copy number variation in the pvgch1 gene was also examined. Polymorphisms were observed in six different codons of the pvdhps gene (382, 383, 512, 549, 553, and 571) and six different codons of the pvdhfr gene (13, 57, 58, 61, 99, 117) in two study sites. The quadruple mutant haplotypes 57I/L/58R/61M/117T of pvdhfr gene were the most common (comprising 76% of cases in Myitsone and 43.7% of case in Laiza). The double mutant haplotype 383G/553G of pvdhps gene was also prevalent at each site (40.8% and 31%). Microsatellites flanking the pvdhfr gene differentiated clinical samples from wild type and quadruple mutant genotypes (F_ST= 0.259-0.3036), as would be expected for a locus undergoing positive selection. The lack of copy number variation of pvgch1 suggests that SP-resistant P. vivax may harbor alternative mechanisms to secure sufficient folate.

Probing the Structural Dynamics of the Plasmodium falciparum Tunneling-Fold Enzyme 6-Pyruvoyl Tetrahydropterin Synthase to Reveal Allosteric Drug Targeting Sites

The de novo folate synthesis pathway is a well-established drug target in the treatment of many infectious diseases. Antimalarial antifolate drugs have proven to be effective against malaria, however, rapid drug resistance has emerged on the two primary targeted enzymes: dihydrofolate reductase and dihydroptoreate synthase. The need to identify alternative antifolate drugs and novel metabolic targets is of imminent importance. The 6-pyruvol tetrahydropterin synthase (PTPS) enzyme belongs to the tunneling fold protein superfamily which is characterized by a distinct central tunnel/cavity. The enzyme catalyzes the second reaction step of the parasite’s de novo folate synthesis pathway and is responsible for the conversion of 7,8-dihydroneopterin to 6-pyruvoyl-tetrahydropterin. In this study, we examine the structural dynamics of Plasmodium falciparum PTPS using the anisotropic network model, to elucidate the collective motions that drive the function of the enzyme and identify potential sites for allosteric modulation of its binding properties. Based on our modal analysis, we identified key sites in the N-terminal domains and central helices which control the accessibility to the active site. Notably, the N-terminal domains were shown to regulate the open-to-closed transition of the tunnel, via a distinctive wringing motion that deformed the core of the protein. We, further, combined the dynamic analysis with motif discovery which revealed highly conserved motifs that are unique to the Plasmodium species and are located in the N-terminal domains and central helices. This provides essential structural information for the efficient design of drugs such as allosteric modulators that would have high specificity and low toxicity as they do not target the PTPS active site that is highly conserved in humans.

AMBER force field parameters for the Zn (II) ions of the tunneling-fold enzymes GTP cyclohydrolase I and 6-pyruvoyl tetrahydropterin synthase

The folate biosynthesis pathway is an essential pathway for cell growth and survival. Folate derivatives serve as a source of the one-carbon units in several intracellular metabolic reactions. Rapidly dividing cells rely heavily on the availability of folate derivatives for their proliferation. As a result, drugs targeting this pathway have shown to be effective against tumor cells and pathogens, but drug resistance against the available antifolate drugs emerged quickly. Therefore, there is a need to develop new treatment strategies and identify alternative metabolic targets. The two de novo folate biosynthesis pathway enzymes, GTP cyclohydrolase I (GCH1) and 6-pyruvoyl tetrahydropterin synthase (PTPS), can provide an alternative strategy to overcome the drug resistance that emerged in the two primary targeted enzymes dihydrofolate reductase and dihydropteroate synthase. Both GCH1 and PTPS enzymes contain Zn²⁺ ions in their active sites, and to accurately study their dynamic behaviors using all-atom molecular dynamics (MD) simulations, appropriate parameters that can describe their metal sites should be developed and validated. In this study, force field parameters of the GCH1 and PTPS metal centers were generated using quantum mechanics (QM) calculations and then validated through MD simulations to ensure their accuracy in describing and maintaining the Zn²⁺ ion coordination environment. The derived force field parameters will provide accurate and reliable MD simulations involving these two enzymes for future in-silico identification of drug candidates against the GCH1 and PTPS enzymes. Communicated by Ramaswamy H. Sarma.

GTP cyclohydrolase I activity from Rickettsia monacensis strain Humboldt, a rickettsial endosymbiont of Ixodes pacificus

The complete folate biosynthesis pathway exists in the genome of a rickettsial endosymbiont of Ixodes pacificus, Rickettsia monacensis strain Humboldt (formerly known as Rickettsia species phylotype G021). Recently, our lab demonstrated that the folA gene of strain Humboldt, the final gene in the folate biosynthesis pathway, encodes a functional dihydrofolate reductase enzyme. In this study, we report R. monacensis strain Humboldt has a functional GTP cyclohydrolase I (GCH1), an enzyme required for the hydrolysis of GTP to form 7,8-dihydroneopterin triphosphate in the folate biosynthesis pathway. The GCH1 gene of R. monacensis, folE, share homology with the folE gene of R. monacensis strain IrR/Munich, with a nucleotide sequence identity of 99%. Amino acid alignment and comparative protein structure modeling have shown that the FolE protein of R. monacensis has a conserved core subunit of GCH1 from the T-fold structural superfamily. All amino acid residues, including conserved GTP binding sites and zinc binding sites, are preserved in the FolE protein of R. monacensis. A recombinant GST-FolE protein from R. monacensis was overexpressed in Escherichia coli, purified by affinity chromatography, and assayed for enzyme activity in vitro. The in vitro enzymatic assay described in this study accorded the recombinant GCH1 enzyme of R. monacensis with a specific activity of 0.81 U/mg. Our data suggest folate genes of R. monacensis strain Humboldt have the potential to produce biochemically active enzymes for de novo folate synthesis, addressing the physioecological underpinnings behind tick-Rickettsia symbioses.

Current progress in antimalarial pharmacotherapy and multi-target drug discovery

Discovery and development of antimalarial drugs have long been dominated by single-target therapy. Continuous effort has been made to explore and identify different targets in malaria parasite crucial for the malaria treatment. The single-target drug therapy was initially successful, but it was later supplanted by combination therapy with multiple drugs to overcome drug resistance. Emergence of resistant strains even against the combination therapy has warranted a review of current antimalarial pharmacotherapy. This has led to the development of the new concept of covalent biotherapy, in which two or more pharmacophores are chemically bound to produce hybrid antimalarial drugs with multi-target functionalities. Herein, the review initially details the current pharmacotherapy for malaria as well as the conventional and novel targets of importance identified in the malaria parasite. Then, the rationale of multi-targeted therapy for malaria, approaches taken to develop the multi-target antimalarial hybrids, and the examples of hybrid molecules are comprehensively enumerated and discussed.

PM2.5-induced extensive DNA methylation changes in the heart of zebrafish embryos and the protective effect of folic acid

We previously found that folic acid (FA) attenuated cardiac defects in zebrafish embryos exposed to extractable organic matter (EOM) from PM2.5, but the underlining mechanisms remain to be elucidated. Since DNA methylation is crucial to cardiac development, we hypothesized that EOM-induced aberrant DNA methylation changes could be diminished by FA supplementation. In this study, zebrafish embryos were exposed to EOM in the absence or presence of FA. Genomic-wide DNA methylation analysis identified both DNA hypo- and hyper-methylation changes in CCGG sites in zebrafish embryos exposed to EOM, which were attenuated by FA supplementation. We identified a total of 316 genes with extensive DNA methylation changes in EOM samples but little or no DNA methylation changes in EOM plus FA samples. The genes were involved in critical cellular processes and signaling pathways important for embryo development. In addition, the EOM-decreased SAM/SAH ratio was counteracted by FA supplementation. Furthermore, FA attenuated the EOM-induced changes in the expression of genes involved in the regulation of DNA methylation and in folate biosynthesis. In conclusion, our data suggest that FA supplementation protected zebrafish embryos from the cardiac developmental toxicity of PM2.5 by alleviating EOM-induced DNA methylation changes.

Overexpression of plasmepsin II and plasmepsin III does not directly cause reduction in Plasmodium falciparum sensitivity to artesunate, chloroquine and piperaquine

Artemisinin derivatives and their partner drugs in artemisinin combination therapies (ACTs) have played a pivotal role in global malaria mortality reduction during the last two decades. The loss of artemisinin efficacy due to evolving drug-resistant parasites could become a serious global health threat. Dihydroartemisinin-piperaquine is a well tolerated and generally highly effective ACT. The implementation of a partner drug in ACTs is critical in the control of emerging artemisinin resistance. Even though artemisinin is highly effective in parasite clearance, it is labile in the human body. A partner drug is necessary for killing the remaining parasites when the pulses of artemisinin have ceased. A population of Plasmodium falciparum parasites in Cambodia and adjacent countries has become resistant to piperaquine. Increased copy number of the genes encoding the haemoglobinases Plasmepsin II and Plasmepsin III has been linked with piperaquine resistance by genome-wide association studies and in clinical trials, leading to the use of increased plasmepsin II/plasmepsin III copy number as a molecular marker for piperaquine resistance. Here we demonstrate that overexpression of plasmepsin II and plasmepsin III in the 3D7 genetic background failed to change the susceptibility of P. falciparum to artemisinin, chloroquine and piperaquine by both a standard dose-response analysis and a piperaquine survival assay. Whilst plasmepsin copy number polymorphism is currently implemented as a molecular surveillance resistance marker, further studies to discover the molecular basis of piperaquine resistance and potential epistatic interactions are needed.

Amplification of GTP-cyclohydrolase 1 gene in Plasmodium falciparum isolates with the quadruple mutant of dihydrofolate reductase and dihydropteroate synthase genes in Ghana

Sulfadoxine-pyrimethamine (SP) is used as malaria chemoprophylaxis for pregnant women and children in Ghana. Plasmodium falciparum resistance to SP is linked to mutations in the dihydropteroate synthase gene (pfdhps), dihydrofolate reductase gene (pfdhfr) and amplification of GTP cyclohydrolase 1 (pfgch1) gene. The pfgch1 duplication is associated with pfdhfr L164, a crucial mutant for high level pyrimethamine resistance which is rare in Ghana. The presence of amplified pfgch1 in Ghanaian isolates could be an indicator of the evolution of the L164 mutant. This study therefore determined the pfgch1 copy number variations and SP resistance mutations in clinical isolates from Ghana. One hundred and ninety-two (192) blood samples collected from children aged ≤14 years with uncomplicated malaria in 2013-14 and 2015-16 were used. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to detect the pfgch1 copy number and nested PCR-Sanger sequencing used to detect mutations in pfdhps and pfdhfr genes. Twelve parasites (6.3%) harbored double copies of the pfgch1 gene out of the 192 samples. Of the 12, 75% had the pfdhfr I51-R59-N108, 92% had the pfdhps G437 mutant, 8% had the pfdhps E540 and 67% had the SP resistance haplotype IRNG. No L164 was detected in samples with amplified pfgch1. The rare T108 mutant associated with cycloguanil resistance showed predominance (60%) over N108 in the 2015-16 isolates. The observation of parasites with increased copy number of pfgch1 gene is indicative of the future evolution of the rare quadruple pfdhfr mutant, I51-R59-N108-L164, in Ghanaian parasites. Mutant pfdhps isolates also had increased gch1 copy number suggestive that it may also facilitate sulphadoxine resistance. The selection of parasites with pfgch1 gene amplification will enhance the sustenance and persistence of parasites with SP resistance in the country. Policy makers need to begin the search for a replacement chemoprophylaxis drug for malaria vulnerable groups in Ghana.

How to combat emerging artemisinin resistance: Lessons from "The Three Little Pigs"

It is rare to come across an Aesop’s fable in respectable journals. It might catch scientists outside the malaria field by surprise to learn that the famous story of “The Boy Who Cried Wolf” has been repeatedly compared to the threat from artemisinin-resistant malaria parasites, including the two latest reports on the rise of a specific haplotype in Cambodia and Thailand, sensationally dubbed “Super Malaria” by the media [1, 2]. The comparison to a children’s tale should not negate the fact that malaria drug resistance is one of the most pressing threats to the global public health community. Here, the findings leading to this contentious discourse will be delineated in order to provide a perspective. Possible solutions will be presented to stimulate further research and discussion to solve one of the greatest public health challenges of our lifetime.

Folate biosynthesis pathway: mechanisms and insights into drug design for infectious diseases

Folate pathway is a key target for the development of new drugs against infectious diseases since the discovery of sulfa drugs and trimethoprim. The knowledge about this pathway has increased in the last years and the catalytic mechanism and structures of all enzymes of the pathway are fairly understood. In addition, differences among enzymes from prokaryotes and eukaryotes could be used for the design of specific inhibitors. In this review, we show a panorama of progress that has been achieved within the folate pathway obtained in the last years. We explored the structure and mechanism of enzymes, several genetic features, strategies, and approaches used in the design of new inhibitors that have been used as targets in pathogen chemotherapy.

A Genetic Variant in FIGN Gene Reduces the Risk of Congenital Heart Disease in Han Chinese Populations

Congenital heart disease (CHD) is one of the most common birth anomalies worldwide. Folate deficiency is an independent risk factor for CHD. Genome-wide association studies (GWAS) revealed that human folate level could be significantly influenced by fidgetin (FIGN), methylenetetrahydrofolate reductase (MTHFR), prickle homolog 2 (PRICKLE2), synaptotagmin 9 (SYT9), gamma-aminobutyric acid B receptor 2 (GABBR2), and alkaline phosphatase (ALPL) genes. The association between the above-mentioned six variants and CHD was examined in the two independent case-control studies in a total of 868 CHD patients and 931 healthy controls. Our results showed that the G > C (rs2119289) variant in intron 4 of FIGN led to a significant reduction of CHD susceptibility in both the separate and combined case-control studies (allele distribution P < 0.001, genotype distribution P < 0.001). Specifically, by analyzing the combined samples, we observed that the risks of CHD in individuals carrying the heterozygous G/C and homozygous C/C genotypes were reduced by 45% (adjusted OR 0.55, 95% CI 0.47-0.67) and 66% (adjusted OR 0.34, 95% CI 0.23-0.50), respectively, in comparison with individuals carrying the wild-type G/G genotype. Our findings have demonstrated that the C allele of variant rs2119289 of FIGN gene is an important genetic marker for decreased CHD risk. Considering that the rs2119289 of FIGN gene is related to the appropriate folate level, FIGN might play an important role in CHD by upregulating plasma folate concentration during embryo heart development. This work provides a new insight for risk assessment of CHD.

A pivot mutation impedes reverse evolution across an adaptive landscape for drug resistance in Plasmodium vivax

BACKGROUND

The study of reverse evolution from resistant to susceptible phenotypes can reveal constraints on biological evolution, a topic for which evolutionary theory has relatively few general principles. The public health catastrophe of antimicrobial resistance in malaria has brought these constraints on evolution into a practical realm, with one proposed solution: withdrawing anti-malarial medication use in high resistance settings, built on the assumption that reverse evolution occurs readily enough that populations of pathogens may revert to their susceptible states. While past studies have suggested limits to reverse evolution, there have been few attempts to properly dissect its mechanistic constraints.

METHODS

Growth rates were determined from empirical data on the growth and resistance from a set of combinatorially complete set of mutants of a resistance protein (dihydrofolate reductase) in Plasmodium vivax, to construct reverse evolution trajectories. The fitness effects of individual mutations were calculated as a function of drug environment, revealing the magnitude of epistatic interactions between mutations and genetic backgrounds. Evolution across the landscape was simulated in two settings: starting from the population fixed for the quadruple mutant, and from a polymorphic population evenly distributed between double mutants.

RESULTS

A single mutation of large effect (S117N) serves as a pivot point for evolution to high resistance regions of the landscape. Through epistatic interactions with other mutations, this pivot creates an epistatic ratchet against reverse evolution towards the wild type ancestor, even in environments where the wild type is the most fit of all genotypes. This pivot mutation underlies the directional bias in evolution across the landscape, where evolution towards the ancestor is precluded across all examined drug concentrations from various starting points in the landscape.

CONCLUSIONS

The presence of pivot mutations can dictate dynamics of evolution across adaptive landscape through epistatic interactions within a protein, leaving a population trapped on local fitness peaks in an adaptive landscape, unable to locate ancestral genotypes. This irreversibility suggests that the structure of an adaptive landscape for a resistance protein should be understood before considering resistance management strategies. This proposed mechanism for constraints on reverse evolution corroborates evidence from the field indicating that phenotypic reversal often occurs via compensatory mutation at sites independent of those associated with the forward evolution of resistance. Because of this, molecular methods that identify resistance patterns via single SNPs in resistance-associated markers might be missing signals for resistance and compensatory mutation throughout the genome. In these settings, whole genome sequencing efforts should be used to identify resistance patterns, and will likely reveal a more complicated genomic signature for resistance and susceptibility, especially in settings where anti-malarial medications have been used intermittently. Lastly, the findings suggest that, given their role in dictating the dynamics of evolution across the landscape, pivot mutations might serve as future targets for therapy.

Computer-Aided Drug Discovery Approaches against the Tropical Infectious Diseases Malaria, Tuberculosis, Trypanosomiasis, and Leishmaniasis

Despite the tremendous improvement in overall global health heralded by the adoption of the Millennium Declaration in the year 2000, tropical infections remain a major health problem in the developing world. Recent estimates indicate that the major tropical infectious diseases, namely, malaria, tuberculosis, trypanosomiasis, and leishmaniasis, account for more than 2.2 million deaths and a loss of approximately 85 million disability-adjusted life years annually. The crucial role of chemotherapy in curtailing the deleterious health and economic impacts of these infections has invigorated the search for new drugs against tropical infectious diseases. The research efforts have involved increased application of computational technologies in mainstream drug discovery programs at the hit identification, hit-to-lead, and lead optimization stages. This review highlights various computer-aided drug discovery approaches that have been utilized in efforts to identify novel antimalarial, antitubercular, antitrypanosomal, and antileishmanial agents. The focus is largely on developments over the past 5 years (2010-2014).