Protein Translation Enzyme lysyl-tRNA Synthetase Presents a New Target for Drug Development against Causative Agents of Loiasis and Schistosomiasis

Helminth parasites are an assemblage of two major phyla of nematodes (also known as roundworms) and platyhelminths (also called flatworms). These parasites are a major human health burden, and infections caused by helminths are considered under neglected tropical diseases (NTDs). These infections are typified by limited clinical treatment options and threat of drug resistance. Aminoacyl-tRNA synthetases (aaRSs) are vital enzymes that decode genetic information and enable protein translation. The specific inhibition of pathogen aaRSs bores well for development of next generation anti-parasitics. Here, we have identified and annotated aaRSs and accessory proteins from Loa loa (nematode) and Schistosoma mansoni (flatworm) to provide a glimpse of these protein translation enzymes within these parasites. Using purified parasitic lysyl-tRNA synthetases (KRSs), we developed series of assays that address KRS enzymatic activity, oligomeric states, crystal structure and inhibition profiles. We show that L. loa and S. mansoni KRSs are potently inhibited by the fungal metabolite cladosporin. Our co-crystal structure of Loa loa KRS-cladosporin complex reveals key interacting residues and provides a platform for structure-based drug development. This work hence provides a new direction for both novel target discovery and inhibitor development against eukaryotic pathogens that include L. loa and S. mansoni.

also been reported to synthesize signaling molecules like diadenosine polyphosphates (Ap4A, Ap5A) that modulate variety of cellular functionalities such as DNA replication, gene expression and ion channel regulation to mention a few [18,22,30,31]. Crystal structures and functional analyses of human cytoplasmic KRS have shown that this enzyme can exist in tetrameric and dimeric forms, where the tetrameric form is bound to multi-synthetase complex and is translationally active, and the dimeric form can participate in transcription regulation and may have cytokine-like functions [32,33]. Thus, determining the oligomeric status of KRSs is of key importance in understanding their functionality and mechanism. Previous reports on the malaria parasite Plasmodium falciparum (Pf) KRS showed that this dimeric protein is inhibited by the fungal secondary metabolite cladosporin with high potency [15,22,34]. Cladosporin targets malaria in both blood and liver stages with IC 50 values below 100 nM [34]. This antimalarial effect is highly selective and mammalian cells are protected as assessed by cytotoxicity assays.
In an effort to understand the protein translation components responsible for supplying charged tRNAs for ribosomal protein synthesis within L. loa and S. mansoni, we first cataloged all their aaRSs and associated protein factors. We noted that L. loa and S. mansoni KRSs are present as single copy in both parasites. We discovered that cladosporin is a very potent inhibitor of L. loa and S. mansoni KRSs, and has enzyme inhibitory IC 50 values in low nanomolar ranges (~52 nM and~97 nM respectively). We provide the X-ray co-crystal structure of this drug bound L. loa KRS to demonstrate its binding mechanism and selectivity. Our proof-ofconcept data on pathogen aaRSs predicts that targeting of other schistosome members and trypanosomes should also be feasible using the same chemical scaffold.

Results
The L. loa and S. mansoni aaRSs Usually one aaRS is required per amino acid and thus a complete set of 20 aaRSs is necessary for protein translation in any cellular context when alternate pathways for producing charged tRNAs are not available [25,26]. We thus looked for all sets of aaRSs in L. loa and S. mansoni genomes. L. loa genome encodes 35 putative aaRSs along with 5 accessory proteins-and this set likely fulfills aminoacylation requirements of its two translational chambers of cytoplasm and mitochondria (S1 Table) [35,36]. A careful analysis of the predicted cellular distribution shows that 19 aaRSs are likely to be localized in the cytoplasm with an absence of KRS (S1 Table). Gene structure suggests that the single copy LlKRS absence from the cytoplasm as it contains a predicted mitochondrial N-terminal signal sequence, but its sequence alignment, domain and motif analyses suggest it to be a eukaryotic-type protein, possibly dual localized. Putative aaRSs with specificity for 16 amino acids are present in L. loa mitochondria with 4 (cysteine, glutamine, glycine and threonine specific) missing aaRSs (S1 Table). Amongst these, the glutamine charged tRNA can be provided by indirect pathway involving putative mitochondrial glutamyl-tRNA amidotransferase [35,36]. In comparison to Loa loa, S. mansoni contains a set of 19 cytoplasmic aaRSs with notable absence of glycine-specific aaRS (S2 Table). However, the mitochondrial GlyRS is a eukaryotic-type enzyme and its dual localization is likely, as has been reported in other organisms [13,37]. The S. mansoni mitochondrial aaRSs set is deficient in GlnRS, HisRS and LysRS (S2 Table). Based on data from several laboratories including ours on aaRSs cellular distributions in eukaryotic pathogens, it is likely that the twin mechanisms of dual aaRs localizations and trafficking of charged tRNA across translational compartments are active in L. loa and S. mansoni as well, in order to provide all substrates required for protein translation in both compartments [13,17,37,38].

L. loa and S. mansoni KRSs are dimeric and enzymatically active
LlKRS and SmKRS enzymes that contained the aminoacylation and anticodon binding domains were produced recombinantly in E. coli (Fig 1A). Our localization predictions and comparative sequence analysis hinted that the LlKRS contained a mitochondrial signal peptide , while the SmKRSs was predicted to be a cytoplasmic enzyme (Fig 1A). To assess the oligomer status of purified proteins, we performed gel permeation chromatography (GPC) experiments using a calibrated column with known standards [22]. Both wild type worm KRSs eluted at sizes corresponding to expected dimers in our GPC experiment, unlike the human counterpart that purportedly shows a tetrameric form (Fig 1B) [32]. Purified worm parasite proteins were used for enzyme assays using SmtRNA Lys and were found to be enzymatically active (Fig 1C and 1D).

Cladosporin binding and inhibition of L. loa and S. mansoni KRSs
Cladosporin is a 3,4-dihydro-6,8-dihydroxy-3-(6-methyl-tetrahydro-2H-pyran-2-yl) compound that mimics the adenosine part of ATP (Fig 2A). To test the activity inhibition and IC 50 values of wild type SmKRS and LlKRS, enzyme assays were performed in presence of cladosporin. Human-like LlKRS V329Q/S346T mutant protein was also produced by taking cues from previous reports and structural data analysis in this work (see next sections). The drug showed concentration-dependent enzymatic inhibition and revealed IC 50 values of 52 nM and 97 nM for wild type LlKRS and SmKRS respectively, while a significantly higher IC 50 value of 1370 nM was observed for the human-like LlKRS mutant protein ( Fig 2B). We also performed protein thermal shift assays to determine the binding of cladosporin to human-like LlKRS and wild type worm KRSs in presence of L-lysine and (a) either no ligand, or (b) with non-hydrolysable ATP analogue (adenosine 5'-(β, γ-imido) triphosphate (AMPPNP)), or (c) with cladosporin in equal micromolar amounts. Data indicated that AMPPNP in a 10:1 molar ratio to KRSs was able to induce a small shift of~0.3°C,~0.4°C and~0.4°C for LlKRS, SmKRS and for human-like LlKRS respectively indicating very weak binding ( Fig 2C). As expected, cladosporin when used in ten-fold higher molar concentration (20uM cladosporin) relative to KRSs (2 μM) was able to induce substantial shifts of 15.5°C and 10.8°C in both LlKRS and SmKRS respectively, indicating high affinity interactions with parasitic KRSs (Fig 2C). On the other hand, significantly smaller shift of 5.3 and 7.2°C were observed when human-like LlKRS (2 μM) was incubated with ten fold (20 μM) or even hundred fold higher (200 μM) concentrations of cladosporin, indicating much poorer binding of cladosporin to the mutant LlKRS ( Fig  2C). To determine the binding affinity of cladosporin, we performed ITC experiments and discovered K d values of 45.2 ± 8.4 nM and 62.8 ± 7.8 nM for LlKRS and SmKRS respectively ( Fig  2D). Similar changes in binding enthalpy (ΔH) and entropy (TΔS) factors indicated conserved mechanism of worm KRS-cladosporin complexation (Table 1). Together, our enzyme inhibition, TSA and ITC data demonstrate strong affinity of cladosporin for these worm KRSs, and validate the potential of cladosporin to selectively bind parasitic KRSs over human counterpart where its affinity is relatively poor (3.3 μM) [15,39].

Structure of LlKRS-CLD-K complex
To understand the atomic basis of KRS-cladosporin binding, we solved the crystal structure L. loa KRS in complex with cladosporin (CLD) and L-lysine (K). Crystal packing analysis showed two dimers of LlKRS in the crystallographic asymmetric unit, validating our GPC results on recombinant LlKRS (Figs 3A and 1B, Table 2). LlKRS folds into a canonical eukaryotic KRS and contains N-terminal OB fold anticodon binding domain and a C-terminal catalytic domain ( Fig 3B). The signature motifs 1, 2 and 3 present in the catalytic domain are also  conserved ( Fig 3B). Cladosporin fits into the ATP binding site in LlKRS and interacts with most of the residues that accommodate adenosine moiety of ATP (Fig 3C and 3D). The isocoumarin ring of cladosporin is stabilized mainly by π-π staking with Phe344, T-stacking with His-340 and hydrogen bondings with Asn341 backbone and Glu334 ( Fig 3D). In addition, guanidine group from Arg563 and Arg332 also stabilize the isocoumarin moiety ( Fig 3D). Gly560 provides hydrophobic support to tetrahydropyran ring (THP) whereas the Ser346 lends suitable space for its methyl moeity. The L-lysine binds in the inner region of active site pocket and is accommodated by series of hydrogen bondings with protein atoms (Fig 3E).

Comparisons with human KRS
A comprehensive sequence alignment of cladosporin-sensitive pathogen KRSs like PfKRS, SmKRS, LlKRS was used to map drug-binding residues in the active sites of these KRSs (Fig 4A  and 4B). HsKRS (PDB: 4YCU) and LlKRS share~66% overall sequence identity and show r.m. s.d. of 1.88 Å in their cladosporin-bound forms for 463 C α atoms. The earlier reported human tetramer enzyme (PDB: 3BJU) and our observed LlKRS dimer show differences in the amino acid sequences and topology of tetramer interface regions 1 and 2 ( Fig 4A and 4B) [32]. Despite the comparatively conserved eukaryotic insertion 1 in LlKRS, the sequence and structural differences in tetramer interface region within it appear to have endowed only dimeric conformation to the worm KRSs (Fig 4A and 4B). Interestingly, LlKRS is sequence-wise and in architectural terms (r.m.s.d.) closer to HsKRS than the PfKRS (54% sequence identity and 2.22 Å r. m.s.d. with LlKRS for 455 C α atoms), and yet possesses cladosporin sensitivity like PfKRS. Availability of recent cladosporin-bound HsKRS structure provided an opportunity for us to compare the human and worm enzymes. The active site region and binding mechanism of cladosporin for both LlKRS and HsKRS is remarkably similar, with only difference of Ser346/ Thr337 (Ll/Hs) and distant Val329/Gln321 (Ll/Hs) substitutions ( Fig 5A). It is clear that the residues Ser346 and Val329 provide extra space for accommodating the methyl moiety of THP ring, and they thus contribute to species selectivity [15,34,39]. PfKRS is currently the beststudied model for understanding cladosporin-binding mechanism, and in addition to two selective residues, many other structural aspects of malaria KRS that contribute to selectivity have become apparent during recent structural investigations by our group and from others [15,22,39]. To understand the LlKRS cladosporin binding and specificity, we analyzed it in backdrop of known Plasmodium and human KRS structures [15,22,39]. In PfKRS, binding of cladosporin induces a loop movement (near motif II) of approximately~2.4 Å towards the active site, and rearrangements of His338 (LlHis331), Phe342 (LlPhe335) and Arg559 (LlArg553) occur to accommodate isocoumarin moiety of cladosporin (Fig 5B and S1 Movie).
(C) Protein thermal shift profile of LlKRS, SmKRS and human-like LlKRS (all three at 2 μM) in presence of cladosporin, or AMPPNP or without these two (but with L-lysine). The plot shows measured derivative Tm and data as plotted against fluorescence (arbitrary units) in y-axis and temperature in x-axis. Human-like LlKRS thermal shift at two concentrations of cladosporin 20 μM (in red) and 200 μM (in black) is shown. Mean data from triplicates are presented. (D) Binding data using ITC. Cladosporin was titrated into the protein samples and K d was determined using Microcal origin software. These events coincide with formation of disulfide bond in a disordered loop region of PfKRS ( Fig 5B). The L-lysine induces an inward mobility in the active site roof region and also stabilizes a disordered loop ( Fig 5B). All four major transitions are present in the PfKRS-CLD-K complex, and the recent crystal structure of HsKRS-CLD-K also shows a structural state similar to PfKRS (r.m.s.d. 1.43 Å for 490 Cα atoms) [15,39] (Fig 5C). The L-lysine induced changes have recently been proposed to be the main factor driving cladosporin species selectivity ( Fig  5B) [39]. To address this, we compared cladosporin and L-lysine bound LlKRS structure to the already available HsKRS and PfKRS cladosporin-bound structures. We found that a helix in Llysine-induced mobile body is disordered in LlKRS (Fig 5D). Additionally, the stable helix (in case of HsKRS) or the disulfide stabilized loop (in case of PfKRS) are also absent in LlKRS and that this region is disordered (Fig 5D). These structural observations hence support the observation that most likely the conserved pair of Val329 and Ser346 found in pathogen KRSs drive cladosporin selectivity.

Discussion
Helminths represent some of the most prevalent neglected tropical disease parasites, and schistosomiasis likely ranks just below malaria as a cause of misery in context of public health [2,7,40]. The currently available monodrug treatment of schistosomiasis using praziquantel poses threat of possible drug resistance [2,7,8,40]. On the other hand, loiasis is prevalent in rainforest and low socioeconomic regions, and has gained prominence in recent years due to adverse effects of drug treatments during co-endemicity with other filarial pathogens [3,6]. Research efforts directed at understanding vital cellular processes such as protein translation machinery can hugely benefit drug discovery initiatives, especially given the promise of utility in context of other infectious diseases like malaria. This is especially of benefit to neglected tropical disease research, where efforts to develop drugs needs to be cost effective. Prompted by these concerns, we sought to dissect worm aaRSs that are responsible for protein translation in these organisms. In this report, we have provided a comprehensive overview of the aaRS distributions in genomes of L. loa and S. mansoni. In both these organisms, it is likely that aaRSs fulfill translational requirements in two cellular compartments of mitochondria and cytoplasm by evolutionarily successful mechanisms of aaRS dual localization, indirect aminoacylation pathways and trafficking of charged tRNAs [13,17,37,38]. Further, the presence of single copy KRSs in both pathogens presents a lucrative opportunity to target this enzyme so as to dismantle protein synthesis process in two translational compartments simultaneously. We additionally found that other members of Schistosoma genus like S. japonicum (Gene-Bank: CAX83109.1) and S. haematobium (GeneBank: KGB31491.1) also possess the conserved cladosporin-sensitive motif and hence could be targeted by cladosporin via their KRSs ( Fig  6A). Interestingly though, sequence analyses show that L. loa and O. volvulus could be specific targets amongst pathogenic filarial nematodes (OVOC_0000240101-mRNA-1) (Fig 6A). Sequence differences in L. loa, W. bancrofti (GeneBank: EJW79634.1) and Brugia malayi (Gen-eBank: CDP92701.1) KRSs suggest that the latter two might be less sensitive to cladosporin, hence providing an opportunity to selectively target L. loa and O. volvulus during co-infections ( Fig 6A). Poor bioavailability is the main limitation in development of cladosporin as a lead  Cladosporin binding induces closing-in of the loop that contains motif 2, with rotameric adjustments in motif 2 residues Phe342, His338 and Arg559 (S1 Movie). This is accompanied by disulfide bond formation in the disordered loop (blue). L-molecule [34]. Chemical synthesis protocols for cladosporin are now available which can aid in structure-guided rational synthesis of more drug-like cladosporin derivatives [41]. Apart from loiasis and schistosomiasis, suitable derivatization of cladosporin for better ADME (absorption, distribution, metabolism, and excretion) properties will be highly valuable for drug development against host of parasitic infections including malaria and feasibly trypanosomosis (given the conservation in active site residues that recognize cladosporin, Fig 6B and 6C) [15]. Hence, cladosporin-based small molecular libraries could be very good starting points for cell-based and phenotypic screening against a number of eukaryotic pathogens. The presented data here therefore provide new avenues for novel drug development against parasitic worm diseases and highlights numerous potential aaRS targets that can now be exploited.

Identification and annotation of aaRSs from L. loa and S. mansoni
The aaRSs and accessory proteins were identified using HMM-search tool in the HMMER web server (http://www.ebi.ac.uk/Tools/hmmer/) by restricting the taxonomy against L. loa and S. mansoni and with an E-value cut-off of 0.01. Pfam IDs of aaRSs and accessory proteins were used in the HMMER based search. Each hit was verified further by sequence, domain and motif analyses using SMART [42], CD-search [43] and superfamily servers [44]. Localizations were predicted using online servers MitoProt (mitochondrial localization-http://ihg.gsf.de/ ihg/mitoprot.html) and NucPred (nuclear localization) [45]. Available mitochondrial localization prediction softwares are trained on non-helminthes organisms and thus sequence alignments with respective validated mitochondrial aaRSs (from NCBI) were also used to identify prokaryotic/mitochondrial type aaRSs. Mitochondrial/prokaryotic type aaRSs found in our analysis, but without the predicted mitochondrial signal sequence were also assigned as putatively mitochondrial. Splice variants from single gene and any atypical aaRSs found in Uniprot [46] were verified via EnsemblMetazoa transcript database (http://metazoa.ensembl.org/index. html). Protein sequences for related worm parasites were obtained from LlKRS or SmKRS protein blast while O. volvulus KRS (OVOC_0000240101-mRNA-1) sequence was obtained from LlKRS protein blast in http://www.sanger.ac.uk/.

Molecular cloning, expression and purification
LlKRS and SmKRS protein sequences were aligned against PfKRS PDB sequence (4PG3) to identify N-terminal sequences of unknown functions or a possible signal sequence, i.e. 1 to 27 for SmKRS and 1 to 77 for the LlKRS. Gene sequences encoding protein residues 28-634 of SmKRS and 78-599 of LlKRS were designed for expression in Escherichia coli and subcloned into pETM-41 vector using NcoI and KpnI restriction sites. Human-like LlKRS V329Q / S346T mutant was created by site directed mutagenesis in two positions V329 to Q and S346 to T. Cloning, expression and purification of the human-like LlKRS protein was performed as for the wild type LlKRS. Protein expression for wild type KRSs and human like LlKRS was induced by adding 0.5 mM isopropyl -d-1-thiogalactopyranoside (IPTG) to cells grown till OD 600 of 0.6-0.8 at 37°C. These cells were grown for 20 h post-induction at 18°C. Bacterial cells were lysine binding further induces a closing-in of mobile element present at roof of active site pocket and stabilization of the loop residues 580-590. The final PfKRS-CLD-K complex with all four major transitions is shown in orange. (C) HsKRS-CLD-K complex (cyan) overall conformation is similar to PfKRS-CLD-K complex (orange). The disulfide-stabilized loop is in an ordered helix in HsKRS. (D) LlKRS differs from previous KRS-cladosporin structures in that the incoming mobile roof and the disulfide regions are both disordered. Hence, cladosporin selectivity for P. falciparum, L. loa and perhaps S. mansoni KRSs is likely driven by the conserved residues Ser and Val that distally line active site pockets in these pathogen KRSs.
doi:10.1371/journal.pntd.0005084.g005 harvested by centrifugation at 5000 g for 15 min and the bacterial pellets were suspended in a buffer consisting of 50 mM Tris-HCl (pH 8.0), 200 mM NaCl, 10 mM beta-mercaptoethanol (βMe), 15%(v/v) glycerol, 0.1 mg ml -1 lysozyme and 1 mM phenylmethylsulfonyl fluoride (PMSF). Cells were lysed by sonication and cleared by centrifugation at 20,000 g for 45 min. The cleared supernatants were applied onto amylose beads (GE Healthcare). All three MBP-KRS fusion proteins were eluted with 10 mM maltose in 50 mM Tris-HCl (pH 8.0), 200 mM NaCl, 10 mM βMe, 1 mM DTT, 0.5 mM EDTA. The MBP tag was removed by incubation of eluted pure fractions with TEV protease at 293 K for 24 h. Wild type and mutant cleaved worm KRSs were concentrated using a 10 kDa cutoff Centricon centrifugal device (Millipore) and were purified by gel-filtration chromatography on a GE HiLoad 60/600 Superdex column. Pure fractions were checked by SDS-PAGE and pooled for crystallization. Before crystallization, the wild type LlKRS was concentrated to 10 mg/ml (A280, extinction coefficient-46760 M -1 cm -1 ) and stored at -80°C. To determine the oligomer status, high molecular weight calibration standards (GE Healthcare) and purified proteins were run on GPC column using protein buffer mobile phase at flow rate of 0.5 ml/ min. Molecular weights of LlKRS and SmKRS were estimated from their elution profiles by plotting log molecular weight (X-axis) against partition coefficient (K av , Y-axis) for known standards.

Aminoacylation assays
SmtRNA Lys was synthesized using an in-vitro transcription method described earlier with minor modifications [47]. A double stranded DNA sequence (TCAGTAGCTG AGTGGATAAT GCGA TGGCGT TTTAAGCGAA CGGTACTGGG TTCGAGTCCCAGAGTGAACCA) encoding cytoplasmic SmtRNA Lys (GeneDB: SmtRNA_01463_Lys_TTT.1.1) containing 5' T7 RNA polymerase promoter, CCA sequence at 3' (in italics) and 2'-O methyl substitution in last two nucleotides of antisense strand was purchased (Sigma). This sequence was transcribed using T7 quick high yield RNA synthesis kit (NEB) at 37°C for 16 h according to manufacturer's instructions. DNA template was removed by DNase (10U/ml) treatment for 3 h in ice followed by addition of EDTA (50 mM). Transcripts were extracted using phenol/CH 3 Cl/isoamyl alcohol (25:24:1) and ethanol precipitation method and reconstituted in 20 mM HEPES, 5 mM EDTA. Samples were further purified using DEAE column (binding buffer 100 mM HEPES-KOH, 12 mM MgCl 2 , 200 mM NaCl (pH 7.5) and elution buffer 100 mM HEPES-KOH, 12 mM MgCl 2 , 800 mM NaCl (pH 7.5)). Fractions containing tRNA were ethanol precipitated, their quality was checked on SDS-urea PAGE and samples were resuspended in 5 mM HEPES-KOH, 1 mM EDTA for storage at -80°C at a concentration of 50 μM. Aminoacylation and enzyme inhibition assays for both KRSs were performed as described elsewhere [17,48]. SmtRNA Lys was refolded prior to enzyme assays by heating at 70°C for 10 minutes followed by addition of 10 mM MgCl 2 and slow cooling to room temperature. Aminoacylation buffer for both worm KRSs contained 30 mM HEPES (pH 7.5), 150 mM NaCl, 30 mM KCl, 50 mM MgCl 2 , 1mM DTT, 100 μM ATP, 500 μM L-lysine, 18 μM SmtRNA Lys , 2 U/ml E. coli inorganic pyrophosphatase (NEB) and 500 nM recombinant SmKRS or LlKRS protein at 37°C. Reaction at different time points was stopped by addition of 40 mM EDTA followed by transfer to ice. Recombinant maltose binding protein (MBP) was used as a control. Cladosporin inhibition assays were performed using inhibitor concentrations in log values ranging from 0.01 nM to 10 μM in the assay buffer.

Thermal shift assays
Protein melt curve assays for both worm KRSs were performed as reported earlier [49]. Both KRSs were diluted in buffer containing 20 mM Tris (pH 8.0), 200 mM NaCl, 5 mM MgCl 2 , 1 mM L-lysine and 2 X SYPRO orange dye (Life Technologies). Total of 20 μM of each ligand AMPPNP (Sigma) and cladosporin (gifted by Bart Staker, SSGCID) was added separately to 2 μM KRS proteins and incubated in ice for 10 min. Ligand bound and unbound samples of both KRSs were heated from 20 to 96°C at a rate of 1°C min -1 and fluorescence signals were monitored by StepOnePlus quantitative real-time PCR system (Life Technologies). Humanlike LlKRS (2 μM) was tested in presence of 20 μM and 200 μM cladosporin concentrations to demonstrate weak binding even with higher inhibitor concentrations. Each curve was an average of three measurements and data were analyzed using thermal shift software (Life technologies). Samples without the addition of ligands were used to determine thermal melting profile of apo proteins. Cladosporin alone and AMPPNP alone in assay buffers, along with no protein controls were used and flat lines were observed for these fluorescence readings at all temperatures. Derivative Tm (melting temperature) was used for analysis.

Isothermal Titration Calorimetry experiments (ITC)
ITC experiments were conducted at 30°C in the MicroCal ITC-200 apparatus (GE Healthcare) and results were fitted into graph using Microcal origin software. Both parasite KRSs were prepared in PBS (phosphate-buffered saline) pH 7.4 with 1 mM L-Lys and 2 mM MgCl 2 , and cladosporin was solubilized in the same buffer. Cladosporin at concentrations of 240 μM and 360 μM was titrated into 21 μM and 22.5 μM protein concentrations of SmKRS and LlKRS respectively. For LlKRS, titrations consisted of 0.4 μl of first injection followed by 39 injections of 1 μl with 150 s intervals between injections. For SmKRS, same titrations were performed with 120 s intervals. Titration of cladosporin in buffer alone was performed to determine the change in enthalpy caused by ligand dilution and then subtracted as background from the actual ligand-binding experiments. Limited protein precipitation with LlKRS during our multiple ITC experiments (in various conditions for n-value optimization trials) was observed and possibly contributed to the lower n-value, which nonetheless is verified as value of 1 based on the crystal structure of the enzyme-drug complex we present here.

Crystallization, data collection and structure determination
Crystallization trials for both SmKRS and LlKRS were performed and crystals of LlKRS were obtained at 20°C by the hanging-drop vapor-diffusion method in 1:1 ratio of LlKRS (10 mg ml -1 ) and mother liquor 20%(w/v) PEG 3350, 200 mM magnesium acetate and 10mM spermidine. Thin, plate-shaped crystal clusters soaked in cryoprotectant 20% glycerol were directly mounted in cooled nitrogen gas at 100 K. X-ray diffraction data were collected using a MAR CCD detector on beamline BM14 of the European Synchrotron Radiation Facility, Grenoble, France. A total of 150 images were collected with 1 min exposure and 1 oscillation per frame. The diffraction images were processed and scaled with HKL-2000 program suite [50]. The structure was solved using phaser-MR with HsKRS as template (66% sequence identity, PDB: 3BJU) [51]. Initial model was built by AutoBuild in PHENIX [51] followed by multiple rounds of manual building using Coot [52] in combination with phenix.refine refinement cycles in PHENIX [51]. All structural superimpositions and preparation of figures was conducted using Chimera [53] and PyMol (http://www.pymol.org). Efforts to crystallize SmKRS, though challenging, are ongoing.  Table. Putative S. mansoni aaRSs and their predicted localizations. Uniprot IDs of the predicted aaRSs in S. mansoni are shown. Single gene variants are shown in italics. (N) denotes the predicted nuclear localization. Putative cytoplasmic phenylalanyl-tRNA synthetase is a heterodimer and the subunits are denoted as (α) and (β) alongside their gene IDs. Two subunits of glutamyl-tRNA amidotransferase are denoted as (A) and (B). (DOCX) S1 Movie. Morph movie for cladosporin (shown in yellow) binding to PfKRS (shown in blue) is shown. In unbound state, Phe342 conformation disallows cladosporin stacking. Additionally, other active site residues are in an unfavorable non-binding conformation. Cladosporin selects for a conformation suitable for stacking and possibly induces rotameric adjustments in active site residues that together stabilize binding. (AVI)