Structural and functional analysis of the anti-malarial drug target prolyl-tRNA synthetase

Aminoacyl tRNA Synthetases: Implications of Structural Biology in Drug Development against Trypanosomatid Parasites

The ensemble of aminoacyl tRNA synthetases is regarded as a key component of the protein translation machinery. With the progressive increase in structure-based studies on tRNA synthetase-ligand complexes, the detailed picture of these enzymes is becoming clear. Having known their critical role in deciphering the genetic code in a living system, they have always been chosen as one of the important targets for development of antimicrobial drugs. Later on, the role of aminoacyl tRNA synthetases (aaRSs) on the survivability of trypanosomatids has also been validated. It became evident through several gene knockout studies that targeting even one of these enzymes affected parasitic growth drastically. Such successful studies have inspired researchers to search for inhibitors that could specifically target trypanosomal aaRSs, and their never-ending efforts have provided fruitful results. Taking all such studies into consideration, these macromolecules of prime importance deserve further investigation for the development of drugs that cure spectrum of infections caused by trypanosomatids. In this review, we have compiled advancements of over a decade that have taken place in the pursuit of devising drugs by using trypanosomatid aaRSs as a major target of interest. Several of these inhibitors work on an exemplary low concentration range without posing any threat to the mammalian cells which is a very critical aspect of the drug discovery process. Advancements have been made in terms of using structural biology as an important tool to analyze the architecture of the trypanosomatids aaRSs and concoction of inhibitors with augmented specificities toward their targets. Some of the inhibitors that have been tested on other parasites successfully but their efficacy has so far not been validated against these trypanosomatids have also been appended.

Computational screening of FDA-approved drugs to identify potential TgDHFR, TgPRS, and TgCDPK1 proteins inhibitors against Toxoplasma gondii

Toxoplasma gondii (T. gondii) is one of the most successful parasites in the world, because about a third of the world's population is seropositive for toxoplasmosis. Treatment regimens for toxoplasmosis have remained unchanged for the past 20 years, and no new drugs have been introduced to the market recently. This study, performed molecular docking to identify interactions of FDA-approved drugs with essential residues in the active site of proteins of T. gondii Dihydrofolate Reductase (TgDHFR), Prolyl-tRNA Synthetase (TgPRS), and Calcium-Dependent Protein Kinase 1 (TgCDPK1). Each protein was docked with 2100 FDA-approved drugs using AutoDock Vina. Also, the Pharmit software was used to generate pharmacophore models based on the TgDHFR complexed with TRC-2533, TgPRS in complex with halofuginone, and TgCDPK1 in complex with a bumped kinase inhibitor, RM-1-132. Molecular dynamics (MD) simulation was also performed for 100 ns to verify the stability of interaction in drug-protein complexes. Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) analysis evaluated the binding energy of selected complexes. Ezetimibe, Raloxifene, Sulfasalazine, Triamterene, and Zafirlukast drugs against the TgDHFR protein, Cromolyn, Cefexim, and Lactulose drugs against the TgPRS protein, and Pentaprazole, Betamethasone, and Bromocriptine drugs against TgCDPK1 protein showed the best results. These drugs had the lowest energy-based docking scores and also stable interactions based on MD analyses with TgDHFR, TgPRS, and TgCDPK1 drug targets that can be introduced as possible drugs for laboratory investigations to treat T. gondii parasite infection.

Structural characterization of glutamyl-tRNA synthetase (GluRS) from Plasmodium falciparum

Aminoacyl-tRNA synthetases (aaRSs) are essential enzymes in protein translation machinery that provide the charged tRNAs needed for protein synthesis. Over the past decades, aaRSs have been studied as anti-parasitic, anti-bacterial, and anti-fungal drug targets. This study focused on the cytoplasmic glutamyl-tRNA synthetase (GluRS) from Plasmodium falciparum, which belongs to class Ib in aaRSs. GluRS unlike most other aaRSs requires tRNA to activate its cognate amino acid substrate L-Glutamate (L-Glu), and fails to form an intermediate adenylate complex in the absence of tRNA. The crystal structures of the Apo, ATP, and ADP-bound forms of Plasmodium falciparum glutamyl-tRNA synthetase (PfGluRS) were solved at 2.1 Å, 2.2 Å, and 2.8 Å respectively. The structural comparison of the Apo- and ATP-bound holo-forms of PfGluRS showed considerable conformational changes in the loop regions around the ATP-binding pocket of the enzyme. Biophysical characterization of the PfGluRS showed binding of the enzyme substrates L-Gluand ATP.. The sequence and structural conservation were evident across GluRS compared to other species. The structural dissection of the PfGluRS gives insight into the critical residues involved in the binding of ATP substrate, which can be harvested to develop new antimalarial drugs.

Targeting prolyl-tRNA synthetase via a series of ATP-mimetics to accelerate drug discovery against toxoplasmosis

The prolyl-tRNA synthetase (PRS) is a validated drug target for febrifugine and its synthetic analog halofuginone (HFG) against multiple apicomplexan parasites including Plasmodium falciparum and Toxoplasma gondii. Here, a novel ATP-mimetic centered on 1-(pyridin-4-yl) pyrrolidin-2-one (PPL) scaffold has been validated to bind to Toxoplasma gondii PRS and kill toxoplasma parasites. PPL series exhibited potent inhibition at the cellular (T. gondii parasites) and enzymatic (TgPRS) levels compared to the human counterparts. Cell-based chemical mutagenesis was employed to determine the mechanism of action via a forward genetic screen. Tg-resistant parasites were analyzed with wild-type strain by RNA-seq to identify mutations in the coding sequence conferring drug resistance by computational analysis of variants. DNA sequencing established two mutations, T477A and T592S, proximal to terminals of the PPL scaffold and not directly in the ATP, tRNA, or L-pro sites, as supported by the structural data from high-resolution crystal structures of drug-bound enzyme complexes. These data provide an avenue for structure-based activity enhancement of this chemical series as anti-infectives.

Exploration of aminoacyl-tRNA synthetases from eukaryotic parasites for drug development

Parasitic diseases result in considerable human morbidity and mortality. The continuous emergence and spread of new drug-resistant parasite strains is an obstacle to controlling and eliminating many parasitic diseases. Aminoacyl-tRNA synthetases (aaRSs) are ubiquitous enzymes essential for protein synthesis. The design and development of diverse small molecule, drug-like inhibitors against parasite-encoded and expressed aaRSs have validated this enzyme family as druggable. In this work, we have compiled the progress to date towards establishing the druggability of aaRSs in terms of their biochemical characterization, validation as targets, inhibitor development, and structural interpretation from parasites responsible for malaria (Plasmodium), lymphatic filariasis (Brugia,Wuchereria bancrofti), giardiasis (Giardia), toxoplasmosis (Toxoplasma gondii), leishmaniasis (Leishmania), cryptosporidiosis (Cryptosporidium), and trypanosomiasis (Trypanosoma). This work thus provides a robust framework for the systematic dissection of aaRSs from these pathogens and will facilitate the cross-usage of potential inhibitors to jump-start anti-parasite drug development.

Anti-Phytophthora Activity of Halofuginone and the Corresponding Mode of Action

Febrifugine, a natural alkaloid, exhibits specific anti-phytophthora activity; however, its mode of action is unclear. In this study, halofuginone, a synthetic derivative of febrifugine, showed significantly higher anti-phytophthora activities than those of febrifugine and the commercial drug metalaxyl against Phytophthora sojae, Phytophthora capsici, and Phytophthora infestans with effective concentration for 50% inhibition (EC₅₀) values of 0.665, 0.673, and 0.178 μg/mL, respectively. Proline could alleviate the growth inhibition of halofuginone on P. capsici, implying that halofuginone might target prolyl-tRNA synthetase (PcPRS). The anti-phytophthora mechanism of halofuginone was then investigated by molecular docking, fluorescence titration, and enzymatic inhibition assays. The results revealed that halofuginone could bind to PcPRS and shared a similar binding site with the substrate proline. Point mutations at Glu316 and Arg345 led to 24.5 and 16.1% decreases in the enzymatic activity of PcPRS but 816.742- and 459.557-fold increases in the resistance to halofuginone, respectively. The results further confirmed that halofuginone was a competitive inhibitor of proline against PcPRS, and Glu316 and Arg345 played important roles in the binding of halofuginone and proline. Taken together, the results indicated that halofuginone is an alternative anti-phytophthora drug candidate and that PcPRS represents a potential target for the development of new pesticides.

Genomic analysis of single nucleotide polymorphisms in malaria parasite drug targets

Malaria is a life-threatening parasitic disease caused by members of the genus Plasmodium. The development and spread of drug-resistant strains of Plasmodium parasites represent a major challenge to malaria control and elimination programmes. Evaluating genetic polymorphism in a drug target improves our understanding of drug resistance and facilitates drug design. Approximately 450 and 19 whole-genome assemblies of Plasmodium falciparum and Plasmodium vivax, respectively, are currently available, and numerous sequence variations have been found due to the presence of single nucleotide polymorphism (SNP). In the study reported here, we analysed global SNPs in the malaria parasite aminoacyl-tRNA synthetases (aaRSs). Our analysis revealed 3182 unique SNPs in the 20 cytoplasmic P. falciparum aaRSs. Structural mapping of SNPs onto the three-dimensional inhibitor-bound complexes of the three advanced drug targets within aaRSs revealed a remarkably low mutation frequency in the crucial aminoacylation domains, low overall occurrence of mutations across samples and high conservation in drug/substrate binding regions. In contrast to aaRSs, dihydropteroate synthase (DHPS), also a malaria drug target, showed high occurrences of drug resistance-causing mutations. Our results show that it is pivotal to screen potent malaria drug targets against global SNP profiles to assess genetic variances to ensure success in designing drugs against validated targets and tackle drug resistance early on.

Inhibition of Plasmodium falciparum phenylalanine tRNA synthetase provides opportunity for antimalarial drug development

Bicyclic azetidine compounds possess antimalarial activity via targeting of the cytoplasmic Plasmodium falciparum (Pf) protein translation enzyme phenylalanine-tRNA synthetase (cFRS). These drugs kill parasites both in vitro and in vivo, including the blood, liver, and transmission developmental stages. Here we present the co-crystal structure of PfcFRS with a potent inhibitor, the bicyclic azetidine BRD7929. Our studies reveal high-affinity binding of BRD7929 with PfcFRS along with exquisite specificity compared with the human enzyme, leading in turn to potent and selective inhibition of the parasite enzyme. Our co-crystal structure shows that BRD7929 binds in the active site in the α subunit of PfcFRS, where it occupies the amino acid site, an auxiliary site, and partially the ATP site. This structural snapshot of inhibitor-bound PfcFRS thus provides a platform for the structure-guided optimization of novel antimalarial compounds.

Double drugging of prolyl-tRNA synthetase provides a new paradigm for anti-infective drug development

Toxoplasmosis is caused by Toxoplasma gondii and in immunocompromised patients it may lead to seizures, encephalitis or death. The conserved enzyme prolyl-tRNA synthetase (PRS) is a validated druggable target in Toxoplasma gondii but the traditional ‘single target–single drug’ approach has its caveats. Here, we describe two potent inhibitors namely halofuginone (HFG) and a novel ATP mimetic (L95) that bind to Toxoplasma gondii PRS simultaneously at different neighbouring sites to cover all three of the enzyme substrate subsites. HFG and L95 act as one triple-site inhibitor in tandem and form an unusual ternary complex wherein HFG occupies the 3’-end of tRNA and the L-proline (L-pro) binding sites while L95 occupies the ATP pocket. These inhibitors exhibit nanomolar IC₅₀ and EC₅₀ values independently, and when given together reveal an additive mode of action in parasite inhibition assays. This work validates a novel approach and lays a structural framework for further drug development based on simultaneous targeting of multiple pockets to inhibit druggable proteins.

Among infectious diseases, parasitic diseases are a major cause of death and morbidity. Toxoplasmosis is caused by an infection of the apicomplexan parasite Toxoplasma gondii. In immunocompromised patients Toxoplasmosis may lead to seizures, encephalitis or death. Novel therapeutics for human parasites are constantly needed. In recent years, the aminoacyl-tRNA synthetase (aaRS) enzyme family has been validated as a drug target for several parasitic infections. The Toxoplasma gondii prolyl-tRNA synthetase inhibitor halofuginone (HFG) has been validated earlier but here we show that an ATP-mimic called L95 is a potent inhibitor and can bind to the target enzyme in the presence of HFG. Thus, the two inhibitors described in this study simultaneously occupy all three natural substrate (ATP, L-amino acid and 3’-end of tRNA) binding pockets and thereby inhibit the enzyme leading to parasite death. This unprecedented double drugging of a pathogen enzyme may delay resistance mutation generation and this approach opens the path to multi-drugging of validated parasite proteins.

Drug targeting of aminoacyl-tRNA synthetases in Anopheles species and Aedes aegypti that cause malaria and dengue

BACKGROUND

Mosquito-borne diseases have a devastating impact on human civilization. A few species of Anopheles mosquitoes are responsible for malaria transmission, and while there has been a reduction in malaria-related deaths worldwide, growing insecticide resistance is a cause for concern. Aedes mosquitoes are known vectors of viral infections, including dengue, yellow fever, chikungunya, and Zika. Aminoacyl-tRNA synthetases (aaRSs) are key players in protein synthesis and are potent anti-infective drug targets. The structure-function activity relationship of aaRSs in mosquitoes (in particular, Anopheles and Aedes spp.) remains unexplored.

METHODS

We employed computational techniques to identify aaRSs from five different mosquito species (Anopheles culicifacies, Anopheles stephensi, Anopheles gambiae, Anopheles minimus, and Aedes aegypti). The VectorBase database ( https://vectorbase.org/vectorbase/app ) and web-based tools were utilized to predict the subcellular localizations (TargetP-2.0, UniProt, DeepLoc-1.0), physicochemical characteristics (ProtParam), and domain arrangements (PfAM, InterPro) of the aaRSs. Structural models for prolyl (PRS)-, and phenylalanyl (FRS)-tRNA synthetases-were generated using the I-TASSER and Phyre protein modeling servers.

RESULTS

Among the vector species, a total of 37 (An. gambiae), 37 (An. culicifacies), 37 (An. stephensi), 37 (An. minimus), and 35 (Ae. aegypti) different aaRSs were characterized within their respective mosquito genomes. Sequence identity amongst the aaRSs from the four Anopheles spp. was > 80% and in Ae. aegypti was > 50%.

CONCLUSIONS

Structural analysis of two important aminoacyl-tRNA synthetases [prolyl (PRS) and phenylanalyl (FRS)] of Anopheles spp. suggests structural and sequence similarity with potential antimalarial inhibitor [halofuginone (HF) and bicyclic azetidine (BRD1369)] binding sites. This suggests the potential for repurposing of these inhibitors against the studied Anopheles spp. and Ae. aegypti.

A sequence-based method for predicting extant fold switchers that undergo α-helix ↔ β-strand transitions

Extant fold‐switching proteins remodel their secondary structures and change their functions in response to cellular stimuli, regulating biological processes and affecting human health. Despite their biological importance, these proteins remain understudied. Predictive methods are needed to expedite the process of discovering and characterizing more of these shapeshifting proteins. Most previous approaches require a solved structure or all‐atom simulations, greatly constraining their use. Here, we propose a high‐throughput sequence‐based method for predicting extant fold switchers that transition from α‐helix in one conformation to β‐strand in the other. This method leverages two previous observations: (a) α‐helix ↔ β‐strand prediction discrepancies from JPred4 are a robust predictor of fold switching, and (b) the fold‐switching regions (FSRs) of some extant fold switchers have different secondary structure propensities when expressed by themselves (isolated FSRs) than when expressed within the context of their parent protein (contextualized FSRs). Combining these two observations, we ran JPred4 on 99‐fold‐switching proteins and found strong correspondence between predicted and experimentally observed α‐helix ↔ β‐strand discrepancies. To test the overall robustness of this finding, we randomly selected regions of proteins not expected to switch folds (single‐fold proteins) and found significantly fewer predicted α‐helix ↔ β‐strand discrepancies. Combining these discrepancies with the overall percentage of predicted secondary structure, we developed a classifier to identify extant fold switchers (Matthews correlation coefficient of .71). Although this classifier had a high false‐negative rate (7/17), its false‐positive rate was very low (2/136), suggesting that it can be used to predict a subset of extant fold switchers from a multitude of available genomic sequences.

Structural analyses of the malaria parasite aminoacyl-tRNA synthetases provide new avenues for antimalarial drug discovery

Malaria is a parasitic illness caused by the genus Plasmodium from the apicomplexan phylum. Five plasmodial species of P. falciparum (Pf), P. knowlesi, P. malariae, P. ovale, and P. vivax (Pv) are responsible for causing malaria in humans. According to the World Malaria Report 2020, there were 229 million cases and ~ 0.04 million deaths of which 67% were in children below 5 years of age. While more than 3 billion people are at risk of malaria infection globally, antimalarial drugs are their only option for treatment. Antimalarial drug resistance keeps arising periodically and thus threatens the main line of malaria treatment, emphasizing the need to find new alternatives. The availability of whole genomes of P. falciparum and P. vivax has allowed targeting their unexplored plasmodial enzymes for inhibitor development with a focus on multistage targets that are crucial for parasite viability in both the blood and liver stages. Over the past decades, aminoacyl-tRNA synthetases (aaRSs) have been explored as anti-bacterial and anti-fungal drug targets, and more recently (since 2009) aaRSs are also the focus of antimalarial drug targeting. Here, we dissect the structure-based knowledge of the most advanced three aaRSs-lysyl- (KRS), prolyl- (PRS), and phenylalanyl- (FRS) synthetases in terms of development of antimalarial drugs. These examples showcase the promising potential of this family of enzymes to provide druggable targets that stall protein synthesis upon inhibition and thereby kill malaria parasites selectively.

In vitro effects of febrifugine on Schistosoma mansoni adult worms

Background

Reports on the antischistosomal effect of several antimalarial drugs such as artesunate, mefloquine, and amodiaquine suggest that febrifugine, which exerts an antimalarial effect, can also be expected to possess antischistosomal potential. The present study investigates the antischistosomal effects of febrifugine.

Methods

In experiment 1, Schistosoma mansoni adult worm pairs were incubated in a medium alone as a control or supplemented with febrifugine at 0.05, 0.1, 0.2, and 0.5 μg/ml for 14 days. The morphology of the worms and the egg production of the female worms were observed simultaneously. In experiment 2, the incubation was conducted as in experiment 1, except that the febrifugine concentrations were reduced to 0.005, 0.01, and 0.02 μg/ml. In addition, S. mansoni adult worms were incubated with either 0.5 μg/ml febrifugine or none as a control for 5 days and stained with neutral red dye.

Results

Febrifugine significantly reduced the survival of S. mansoni male and female worms at concentrations of 0.02-0.5 μg/ml following incubation for 14 days and remarkably inhibited the daily egg output of the female worms. The non-treated male and female worms remained morphologically normal within the period of 14 days, whereas male and female worms treated with febrifugine at different concentrations gradually twisted and subsequently died. In addition, S. mansoni adult worms were incubated with either 0.5 μg/ml febrifugine or none as a control for 5 days and stained with neutral red dye. Non-treated male worms were morphologically normal and stained dark red with neutral red, while febrifugine-treated male worms appeared similar to those in the control group and were stained at a slightly lower level of dark red than the non-treated male worms. Non-treated female worms were morphologically normal, and their intestinal tract and vitellaria were stained deep red and dark red, respectively. In contrast, febrifugine-treated female worms were morphologically damaged, and their intestinal tract and vitellaria remained mostly unstained and stained dark red, respectively.

Conclusion

Febrifugine exerts potent antischistosomal effects and can be expected to contribute to the development of a novel antischistosomal drug.

Identification of Selective Novel Hits against Plasmodium falciparum Prolyl tRNA Synthetase Active Site and a Predicted Allosteric Site Using in silico Approaches

Recently, there has been increased interest in aminoacyl tRNA synthetases (aaRSs) as potential malarial drug targets. These enzymes play a key role in protein translation by the addition of amino acids to their cognate tRNA. The aaRSs are present in all Plasmodium life cycle stages, and thus present an attractive malarial drug target. Prolyl tRNA synthetase is a class II aaRS that functions in charging tRNA with proline. Various inhibitors against Plasmodium falciparum ProRS (PfProRS) active site have been designed. However, none have gone through clinical trials as they have been found to be highly toxic to human cells. Recently, a possible allosteric site was reported in PfProRS with two possible allosteric modulators: glyburide and TCMDC-124506. In this study, we sought to identify novel selective inhibitors targeting PfProRS active site and possible novel allosteric modulators of this enzyme. To achieve this, virtual screening of South African natural compounds against PfProRS and the human homologue was carried out using AutoDock Vina. The modulation of protein motions by ligand binding was studied by molecular dynamics (MD) using the GROningen MAchine for Chemical Simulations (GROMACS) tool. To further analyse the protein global motions and energetic changes upon ligand binding, principal component analysis (PCA), and free energy landscape (FEL) calculations were performed. Further, to understand the effect of ligand binding on the protein communication, dynamic residue network (DRN) analysis of the MD trajectories was carried out using the MD-TASK tool. A total of ten potential natural hit compounds were identified with strong binding energy scores. Binding of ligands to the protein caused observable global and residue level changes. Dynamic residue network calculations showed increase in betweenness centrality (BC) metric of residues at the allosteric site implying these residues are important in protein communication. A loop region at the catalytic domain between residues 300 and 350 and the anticodon binding domain showed significant contributions to both PC1 and PC2. Large motions were observed at a loop in the Z-domain between residues 697 and 710 which was also in agreement with RMSF calculations that showed increase in flexibility of residues in this region. Residues in this loop region are implicated in ATP binding and thus a change in dynamics may affect ATP binding affinity. Free energy landscape (FEL) calculations showed that the holo protein (protein-ADN complex) and PfProRS-SANC184 complexes were stable, as shown by the low energy with very few intermediates and hardly distinguishable low energy barriers. In addition, FEL results agreed with backbone RMSD distribution plots where stable complexes showed a normal RMSD distribution while unstable complexes had multimodal RMSD distribution. The betweenness centrality metric showed a loss of functional importance of key ATP binding site residues upon allosteric ligand binding. The deep basins in average L observed at the allosteric region imply that there is high accessibility of residues at this region. To further analyse BC and average L metrics data, we calculated the ΔBC and ΔL values by taking each value in the holo protein BC or L matrix less the corresponding value in the ligand-bound complex BC or L matrix. Interestingly, in allosteric complexes, residues located in a loop region implicated in ATP binding had negative ΔL values while in orthosteric complexes these residues had positive ΔL values. An increase in contact frequency between residues Ser263, Thr267, Tyr285, and Leu707 at the allosteric site and residues Thr397, Pro398, Thr402, and Gln395 at the ATP binding TXE loop was observed. In summary, this study identified five potential orthosteric inhibitors and five allosteric modulators against PfProRS. Allosteric modulators changed ATP binding site dynamics, as shown by RMSF, PCA, and DRN calculations. Changes in dynamics of the ATP binding site and increased contact frequency between residues at the proposed allosteric site and the ATP binding site may explain how allosteric modulators distort the ATP binding site and thus might inhibit PfProRS. The scaffolds of the identified hits in the study can be used as a starting point for antimalarial inhibitor development with low human cytotoxicity.

Conformational heterogeneity in apo and drug-bound structures of Toxoplasma gondii prolyl-tRNA synthetase

Prolyl-tRNA synthetase (PRS) is a member of the aminoacyl-tRNA synthetase family that drives protein translation in cells. The apicomplexan PRSs are validated targets of febrifugine (FF) and its halogenated derivative halofuginone (HF). PRSs are of great interest for drug development against Plasmodium falciparum and Toxoplasma gondii. In this study, structures of apo and FF-bound T. gondii (TgPRS) are revealed and the dynamic nature of the conformational changes that occur upon FF binding is unraveled. In addition, this study highlights significant conformational plasticity within two different crystal structures of apo PRSs but not within drug-bound PRSs. The apo PRSs exist in multi-conformational states and manifest pseudo-dimeric structures. In contrast, when FF is bound the PRS dimer adopts a highly symmetrical architecture. It is shown that TgPRS does not display extant fold switching, in contrast to P. falciparum PRS, despite having over 65% sequence identity. Finally, structure-comparison analyses suggest the utility of r.m.s.d. per residue (r.m.s.d.^/res) as a robust tool to detect structural alterations even when the r.m.s.d. is low. Apo TgPRS reveals FF/HF-induced rigidity and this work has implications for drug-design studies that rely on the apo structures of target proteins.

Aminoacyl tRNA synthetases as potential drug targets of the Panthera pathogen Babesia

Background

A century ago, pantheras were abundant across Asia. Illegal hunting and trading along with loss of habitat have resulted in the designation of Panthera as a genus of endangered species. In addition to the onslaught from humans, pantheras are also susceptible to outbreaks of several infectious diseases, including babesiosis. The latter is a hemoprotozoan disease whose causative agents are the eukaryotic parasites of the apicomplexan genus Babesia. Babesiosis affects a varied range of animals including humans (Homo sapiens), bovines (e.g. Bos taurus), pantheras (e.g. Panthera tigris, P. leo, P. pardus) and equines. Babesia spp. are transmitted by the tick vector Ixodes scapularis or ticks of domestic animals, namely Rhipicephalus (Boophilus) microplus and R. (B.) decoloratus. At the level of protein translation within these organisms, the conserved aminoacyl tRNA synthetase (aaRS) family offers an opportunity to identify the sequence and structural differences in the host (Panthera) and parasites (Babesia spp.) in order to exploit these for drug targeting Babesia spp.

Methods

Using computational tools we investigated the genomes of Babesia spp. and Panthera tigris so as to annotate their aaRSs. The sequences were analysed and their subcellular localizations were predicted using Target P1.1, SignalP 3.0, TMHMM v.2.0 and Deeploc 1.0 web servers. Structure-based analysis of the aaRSs from P. tigris and its protozoan pathogens Babesia spp. was performed using Phyre2 and chimera.

Results

We identified 33 (B. bovis), 34 (B. microti), 33 (B. bigemina) and 33 (P. tigris) aaRSs in these respective organisms. Poor sequence identity (~ 20–50%) between aaRSs from Babesia spp. and P. tigris was observed and this merits future experiments to validate new drug targets against Babesia spp.

Conclusions

Overall this work provides a foundation for experimental investigation of druggable aaRSs from Babesia sp. in an effort to control Babesiosis in Panthera.

Inaccurate secondary structure predictions often indicate protein fold switching

Although most proteins conform to the classical one‐structure/one‐function paradigm, an increasing number of proteins with dual structures and functions have been discovered. In response to cellular stimuli, such proteins undergo structural changes sufficiently dramatic to remodel even their secondary structures and domain organization. This “fold‐switching” capability fosters protein multi‐functionality, enabling cells to establish tight control over various biochemical processes. Accurate predictions of fold‐switching proteins could both suggest underlying mechanisms for uncharacterized biological processes and reveal potential drug targets. Recently, we developed a prediction method for fold‐switching proteins using structure‐based thermodynamic calculations and discrepancies between predicted and experimentally determined protein secondary structure (Porter and Looger, Proc Natl Acad Sci U S A 2018; 115:5968–5973). Here we seek to leverage the negative information found in these secondary structure prediction discrepancies. To do this, we quantified secondary structure prediction accuracies of 192 known fold‐switching regions (FSRs) within solved protein structures found in the Protein Data Bank (PDB). We find that the secondary structure prediction accuracies for these FSRs vary widely. Inaccurate secondary structure predictions are strongly associated with fold‐switching proteins compared to equally long segments of non‐fold‐switching proteins selected at random. These inaccurate predictions are enriched in helix‐to‐strand and strand‐to‐coil discrepancies. Finally, we find that most proteins with inaccurate secondary structure predictions are underrepresented in the PDB compared with their alternatively folded cognates, suggesting that unequal representation of fold‐switching conformers within the PDB could be an important cause of inaccurate secondary structure predictions. These results demonstrate that inconsistent secondary structure predictions can serve as a useful preliminary marker of fold switching.

Genomic analyses of aminoacyl tRNA synthetases from human-infecting helminths

BACKGROUND

Helminth infections affect ~ 60% of the human population that lives in tropical and subtropical regions worldwide. These infections result in diseases like schistosomiasis, lymphatic filariasis, river blindness and echinococcosis. Here we provide a comprehensive computational analysis of the aminoacyl tRNA synthetase (aaRS) enzyme family from 27 human-infecting helminths. Our analyses support the idea that several helminth aaRSs can be targeted for drug repurposing or for development of new drugs. For experimental validation, we focused on Onchocerciasis (also known as "river blindness"), a filarial vector-borne disease that is prevalent in Africa and Latin America. We show that halofuginone (HF) can act as a potent inhibitor of Onchocerca volvulus prolyl tRNA synthetase (OvPRS).

RESULTS

The conserved enzyme family of aaRSs has been validated as druggable targets in numerous eukaryotic parasites. We thus embarked on assessing aaRSs from the genomes of 27 helminths that cause infections in humans. In order to delineate the distribution of aaRSs per genome we utilized Hidden Markov Models of aaRS catalytic domains to identify all orthologues. We note that Fasciola hepatica genome encodes the highest number of aaRS-like proteins (69) whereas Taenia asiatica has the lowest count (32). The number of genes for any particular aaRS-like protein varies from 1 to 8 in these 27 studied helminths. Sequence alignments of helminth-encoded lysyl, prolyl, leucyl and threonyl tRNA synthetases suggest that various known aaRS inhibitors like Cladosporin, Halofuginone, Benzoborale and Borrelidin may be of utility against helminths. The recombinantly expressed Onchocerca volvulus PRS was used as proof of concept for targeting aaRS with drug-like molecules like HF.

CONCLUSIONS

Systematic analysis of unique subdomains within helminth aaRSs reveals the presence of a number of non-canonical domains like PAC3, Utp-14, Pex2_Pex12 fused to catalytic domains in the predicted helminth aaRSs. We have established a platform for biochemical validation of a large number of helminth aaRSs that can be targeted using available inhibitors to jump-start drug repurposing against human helminths.

A Weak Spot in Multiple Protozoan Parasites

Eukaryotic protozoan parasites, including the etiological agents of malaria, toxoplasmosis, and leishmaniasis, collectively cause significant mortality in humans. In a recent issue of Structure,Jain et al. (2017) identify a set of quinazolinone-based derivatives targeting the parasitic prolyl-tRNA synthetase enzyme as promising drugs for the clearance of diverse parasites.

Aminoacyl tRNA synthetases as malarial drug targets: a comparative bioinformatics study

BACKGROUND

Treatment of parasitic diseases has been challenging due to evolution of drug resistant parasites, and thus there is need to identify new class of drugs and drug targets. Protein translation is important for survival of malarial parasite, Plasmodium, and the pathway is present in all of its life cycle stages. Aminoacyl tRNA synthetases are primary enzymes in protein translation as they catalyse amino acid addition to the cognate tRNA. This study sought to understand differences between Plasmodium and human aminoacyl tRNA synthetases through bioinformatics analysis.

METHODS

Plasmodium berghei, Plasmodium falciparum, Plasmodium fragile, Plasmodium knowlesi, Plasmodium malariae, Plasmodium ovale, Plasmodium vivax, Plasmodium yoelii and human aminoacyl tRNA synthetase sequences were retrieved from UniProt database and grouped into 20 families based on amino acid specificity. These families were further divided into two classes. Both families and classes were analysed. Motif discovery was carried out using the MEME software, sequence identity calculation was done using an in-house Python script, multiple sequence alignments were performed using PROMALS3D and TCOFFEE tools, and phylogenetic tree calculations were performed using MEGA vs 7.0 tool. Possible alternative binding sites were predicted using FTMap webserver and SiteMap tool.

RESULTS

Motif discovery revealed Plasmodium-specific motifs while phylogenetic tree calculations showed that Plasmodium proteins have different evolutionary history to the human homologues. Human aaRSs sequences showed low sequence identity (below 40%) compared to Plasmodium sequences. Prediction of alternative binding sites revealed potential druggable sites in PfArgRS, PfMetRS and PfProRS at regions that are weakly conserved when compared to the human homologues. Multiple sequence analysis, motif discovery, pairwise sequence identity calculations and phylogenetic tree analysis showed significant differences between parasite and human aaRSs proteins despite functional and structural conservation. These differences may provide a basis for further exploration of Plasmodium aminoacyl tRNA synthetases as potential drug targets.

CONCLUSION

This study showed that, despite, functional and structural conservation, Plasmodium aaRSs have key differences from the human homologues. These differences in Plasmodium aaRSs can be targeted to develop anti-malarial drugs with less toxicity to the host.

Extant fold-switching proteins are widespread

It is commonly thought that each globular protein has a single 3D structure, or fold, that fosters its function. In contrast, recent studies have identified several fold-switching proteins whose secondary structures can be remodeled in response to cellular stimuli. Although thought to be rare, we found 96 literature-validated fold-switching proteins by exhaustively searching the database of protein structures [Protein Data Bank (PDB)]. Characterizing these proteins led us to hypothesize that their abundance may be underrepresented in the PDB. Thus, we developed a computational method that identifies fold-switching proteins and used it to estimate that 0.5–4% of PDB proteins switch folds. These results suggest that proteins switch folds with significant frequency, which has implications for cell biology, genomics, and human health.

A central tenet of biology is that globular proteins have a unique 3D structure under physiological conditions. Recent work has challenged this notion by demonstrating that some proteins switch folds, a process that involves remodeling of secondary structure in response to a few mutations (evolved fold switchers) or cellular stimuli (extant fold switchers). To date, extant fold switchers have been viewed as rare byproducts of evolution, but their frequency has been neither quantified nor estimated. By systematically and exhaustively searching the Protein Data Bank (PDB), we found ∼100 extant fold-switching proteins. Furthermore, we gathered multiple lines of evidence suggesting that these proteins are widespread in nature. Based on these lines of evidence, we hypothesized that the frequency of extant fold-switching proteins may be underrepresented by the structures in the PDB. Thus, we sought to identify other putative extant fold switchers with only one solved conformation. To do this, we identified two characteristic features of our ∼100 extant fold-switching proteins, incorrect secondary structure predictions and likely independent folding cooperativity, and searched the PDB for other proteins with similar features. Reassuringly, this method identified dozens of other proteins in the literature with indication of a structural change but only one solved conformation in the PDB. Thus, we used it to estimate that 0.5–4% of PDB proteins switch folds. These results demonstrate that extant fold-switching proteins are likely more common than the PDB reflects, which has implications for cell biology, genomics, and human health.

Chemical synthesis of febrifugine and analogues

The quinazolinone-containing 2,3-disubstituted piperidines febrifugine and isofebrifugine have been the subject of significant research efforts since their occurrence in Dichroa febrifuga and their anti-malarial actions were first described in the late 1940s. Subsequently they have also been shown to be present in other plants belonging to the hydrangea family and various analogues of febrifugine have been prepared in attempts to tune biological properties. The most notable analogue is termed halofuginone and a substantial body of work now demonstrates that this compound possesses potent human disease relevant activities. This review focuses on the literature associated with efforts dedicated towards uncovering the structures of febrifugine and isofebrifugine, the development of practical methods for their synthesis and the syntheses of structural analogues.

Drug targeting of one or more aminoacyl-tRNA synthetase in the malaria parasite Plasmodium falciparum

Malaria remains a major infectious disease and, despite incidence reduction, it threatens resurgence in drug-resistant forms. Antimalarial drugs remain the mainstay of therapeutic options and hence there is a constant need to identify and validate new druggable targets. Plasmodium falciparum aminoacyl-tRNA synthetases (Pf-aaRSs) drive protein translation and are potent targets for development of next-generation antimalarials. Here, we detail advances made in structural-biology-based investigations in Pf-aaRSs and discuss their distribution of druggable pockets. This review establishes a platform for systematic experimental dissection of malarial parasite aaRSs as a new focus for sustained drug development efforts against malaria.

Targeting Prolyl-tRNA Synthetase to Accelerate Drug Discovery against Malaria, Leishmaniasis, Toxoplasmosis, Cryptosporidiosis, and Coccidiosis

Developing anti-parasitic lead compounds that act on key vulnerabilities are necessary for new anti-infectives. Malaria, leishmaniasis, toxoplasmosis, cryptosporidiosis and coccidiosis together kill >500,000 humans annually. Their causative parasites Plasmodium, Leishmania, Toxoplasma, Cryptosporidium and Eimeria display high conservation in many housekeeping genes, suggesting that these parasites can be attacked by targeting invariant essential proteins. Here, we describe selective and potent inhibition of prolyl-tRNA synthetases (PRSs) from the above parasites using a series of quinazolinone-scaffold compounds. Our PRS-drug co-crystal structures reveal remarkable active site plasticity that accommodates diversely substituted compounds, an enzymatic feature that can be leveraged for refining drug-like properties of quinazolinones on a per parasite basis. A compound we termed In-5 exhibited a unique double conformation, enhanced drug-like properties, and cleared malaria in mice. It thus represents a new lead for optimization. Collectively, our data offer insights into the structure-guided optimization of quinazolinone-based compounds for drug development against multiple human eukaryotic pathogens.

Structure-Based Targeting of Orthologous Pathogen Proteins Accelerates Antiparasitic Drug Discovery

Parasitic diseases caused by eukaryotic pathogens impose significant health and economic burden worldwide. The level of research funding available for many parasitic diseases is insufficient in relation to their adverse social and economic impact. In this article, we discuss that extant 3D structural data on protein-inhibitor complexes can be harnessed to accelerate drug discovery against many related pathogens. Assessment of sequence conservation within drug/inhibitor-binding residues in enzyme-inhibitor complexes can be leveraged to predict and validate both new lead compounds and their molecular targets in multiple parasitic diseases. Hence, structure-based targeting of orthologous pathogen proteins accelerates the discovery of new antiparasitic drugs. This approach offers significant benefits for jumpstarting the discovery of new lead compounds and their molecular targets in diverse human, livestock, and plant pathogens.

Genome-Wide Association Studies of Drug-Resistance Determinants

Population genetic strategies that leverage association, selection, and linkage have identified drug-resistant loci. However, challenges and limitations persist in identifying drug-resistance loci in malaria. In this review we discuss the genetic basis of drug resistance and the use of genome-wide association studies, complemented by selection and linkage studies, to identify and understand mechanisms of drug resistance and response. We also discuss the implications of nongenetic mechanisms of drug resistance recently reported in the literature, and present models of the interplay between nongenetic and genetic processes that contribute to the emergence of drug resistance. Throughout, we examine artemisinin resistance as an example to emphasize challenges in identifying phenotypes suitable for population genetic studies as well as complications due to multiple-factor drug resistance.

Dimerization of Arginyl-tRNA Synthetase by Free Heme Drives Its Inactivation in Plasmodium falciparum

Excess cellular heme is toxic, and malaria parasites regulate its levels during hemoglobin digestion. Aminoacyl-tRNA synthetases are ubiquitous enzymes, and of these, arginyl-tRNA synthetase (RRS) is unique as its enzymatic product of charged tRNA is required for protein synthesis and degradation. We show that Plasmodium falciparum arginyl-tRNA synthetase (PfRRS) is an active, cytosolic, and monomeric enzyme. Its high-resolution crystal structure highlights critical structural differences with the human enzyme. We further show that hemin binds to and inhibits the aminoacylation activity of PfRRS. Hemin induces a dimeric form of PfRRS that is thus rendered enzymatically dead as it is unable to recognize its cognate tRNA(arg). Excessive hemin in chloroquine-treated malaria parasites results in significantly reduced charged tRNA(arg) levels, thus suggesting deceleration of protein synthesis. These data together suggest that the inhibition of Plasmodium falciparum arginyl-tRNA synthetase can now be synergized with existing antimalarials for more potent drug cocktails against malaria parasites.

Recent advances in the biology and drug targeting of malaria parasite aminoacyl-tRNA synthetases

Escalating drug resistance in malaria parasites and lack of vaccine entails the discovery of novel drug targets and inhibitor molecules. The multi-component protein translation machinery is a rich source of such drug targets. Malaria parasites contain three translational compartments: the cytoplasm, apicoplast and mitochondrion, of which the latter two are of the prokaryotic type. Recent explorations by many groups into the malaria parasite protein translation enzymes, aminoacyl-tRNA synthetases (aaRSs), have yielded many promising inhibitors. The understanding of the biology of this unique set of 36 enzymes has become much clearer in recent times. Current review discusses the advances made in understanding of crucial aaRSs from Plasmodium and also the specific inhibitors found against malaria aaRSs.

The cytoplasmic prolyl-tRNA synthetase of the malaria parasite is a dual-stage target of febrifugine and its analogs

The emergence of drug resistance is a major limitation of current antimalarials. The discovery of new druggable targets and pathways including those that are critical for multiple life cycle stages of the malaria parasite is a major goal for developing next-generation antimalarial drugs. Using an integrated chemogenomics approach that combined drug resistance selection, whole-genome sequencing, and an orthogonal yeast model, we demonstrate that the cytoplasmic prolyl-tRNA (transfer RNA) synthetase (PfcPRS) of the malaria parasite Plasmodium falciparum is a biochemical and functional target of febrifugine and its synthetic derivative halofuginone. Febrifugine is the active principle of a traditional Chinese herbal remedy for malaria. We show that treatment with febrifugine derivatives activated the amino acid starvation response in both P. falciparum and a transgenic yeast strain expressing PfcPRS. We further demonstrate in the Plasmodium berghei mouse model of malaria that halofuginol, a new halofuginone analog that we developed, is active against both liver and asexual blood stages of the malaria parasite. Halofuginol, unlike halofuginone and febrifugine, is well tolerated at efficacious doses and represents a promising lead for the development of dual-stage next-generation antimalarials.

Three-dimensional structures in the design of therapeutics targeting parasitic protozoa: reflections on the past, present and future

Parasitic protozoa cause a range of diseases which threaten billions of human beings. They are responsible for tremendous mortality and morbidity in the least-developed areas of the world. Presented here is an overview of the evolution over the last three to four decades of structure-guided design of inhibitors, leads and drug candidates aiming at targets from parasitic protozoa. Target selection is a crucial and multi-faceted aspect of structure-guided drug design. The major impact of advances in molecular biology, genome sequencing and high-throughput screening is touched upon. The most advanced crystallographic techniques, including XFEL, have already been applied to structure determinations of drug targets from parasitic protozoa. Even cryo-electron microscopy is contributing to our understanding of the mode of binding of inhibitors to parasite ribosomes. A number of projects have been selected to illustrate how structural information has assisted in arriving at promising compounds that are currently being evaluated by pharmacological, pharmacodynamic and safety tests to assess their suitability as pharmaceutical agents. Structure-guided approaches are also applied to incorporate properties into compounds such that they are less likely to become the victim of resistance mechanisms. A great increase in the number of novel antiparasitic compounds will be needed in the future. These should then be combined into various multi-compound therapeutics to circumvent the diverse resistance mechanisms that render single-compound, or even multi-compound, drugs ineffective. The future should also see (i) an increase in the number of projects with a tight integration of structural biology, medicinal chemistry, parasitology and pharmaceutical sciences; (ii) the education of more `medicinal structural biologists' who are familiar with the properties that compounds need to have for a high probability of success in the later steps of the drug-development process; and (iii) the expansion of drug-development capabilities in middle- and low-income countries.

Structure of Prolyl-tRNA Synthetase-Halofuginone Complex Provides Basis for Development of Drugs against Malaria and Toxoplasmosis

The Chinese herb Dichroa febrifuga has traditionally treated malaria-associated fever. Its active component febrifugine (FF) and derivatives such as halofuginone (HF) are potent anti-malarials. Here, we show that FF-based derivatives arrest parasite growth by direct interaction with and inhibition of the protein translation enzyme prolyl-tRNA synthetase (PRS). Dual administration of inhibitors that target different tRNA synthetases suggests high utility of these drug targets. We reveal the ternary complex structure of PRS-HF and adenosine 5'-(β,γ-imido)triphosphate where the latter facilitates HF integration into the PRS active site. Structural analyses also highlight spaces within the PRS architecture for HF derivatization of its quinazolinone, but not piperidine, moiety. We also show a remarkable ability of HF to kill the related human parasite Toxoplasma gondii, suggesting wider HF efficacy against parasitic PRSs. Hence, our cell-, enzyme-, and structure-based data on FF-based inhibitors strengthen the case for their inclusion in anti-malarial and anti-toxoplasmosis drug development efforts.

Halofuginone - the multifaceted molecule

Halofuginone is an analog of febrifugine-an alkaloid originally isolated from the plant Dichroa febrifuga. During recent years, halofuginone has attracted much attention because of its wide range of beneficial biological activities, which encompass malaria, cancer, and fibrosis-related and autoimmune diseases. At present two modes of halofuginone actions have been described: (1) Inhibition of Smad3 phosphorylation downstream of the TGFβ signaling pathway results in inhibition of fibroblasts-to-myofibroblasts transition and fibrosis. (2) Inhibition of prolyl-tRNA synthetase (ProRS) activity in the blood stage of malaria and inhibition of Th17 cell differentiation thereby inhibiting inflammation and the autoimmune reaction by activation of the amino acid starvation and integrated stress responses. This review deals with the history and origin of this natural product, its synthesis, its known modes of action, and it's various biological activities in pre-clinical animal models and in humans.