Structural and enzymatic insights into species-specific resistance to schistosome parasite drug therapy

The antischistosomal prodrug oxamniquine is activated by a sulfotransferase (SULT) in the parasitic flatworm Schistosoma mansoni. Of the three main human schistosome species, only S. mansoni is sensitive to oxamniquine therapy despite the presence of SULT orthologs in Schistosoma hematobium and Schistosoma japonicum. The reason for this species-specific drug action has remained a mystery for decades. Here we present the crystal structures of S. hematobium and S. japonicum SULTs, including S. hematobium SULT in complex with oxamniquine. We also examined the activity of the three enzymes in vitro; surprisingly, all three are active toward oxamniquine, yet we observed differences in catalytic efficiency that implicate kinetics as the determinant for species-specific toxicity. These results provide guidance for designing oxamniquine derivatives to treat infection caused by all species of schistosome to combat emerging resistance to current therapy.

Middle East, South America, and the West Indies. Schistosomiasis, the chronic illness resulting from schistosome infection, is classified by the World Health Organization as a neglected tropical disease. Three species of the human blood fluke, Schistosoma mansoni, Schistosoma hematobium, and Schistosoma japonicum, account for almost all cases of human schistosomiasis (2). Symptoms vary according to the infecting species, and chronic infection can damage the liver, spleen, lungs, intestines, and urogenital tract (3). Effective drug treatment is available using the drug praziquantel, currently the recommended drug for schistosome infection. Praziquantel intervention provides enormous health benefits, but the selective pressure the monotherapy applies to the parasite has the potential to drive resistance to the drug. Schistosomes with decreased susceptibility to praziquantel have already been documented in patients in Egypt, Senegal, and Tanzania (4); therefore, combination treatments are needed.
Oxamniquine was the first-line antischistosomal drug in Brazil until the 1980s when it was gradually replaced by praziquantel (5) because of its lower cost and the threat of emerging oxamniquine resistance (6 -9). We recently identified a sulfotransferase (SULT) 4 enzyme in S. mansoni (SmSULT) as the oxamniquine target, leading to the X-ray crystal structure of the enzyme in a complex with oxamniquine and the reaction product 3Ј-phosphoadenosine 5Ј-phosphate (PAP) (10). Oxamniquine is a species-specific treatment effective against S. mansoni but not S. hematobium or S. japonicum despite their SULTs sharing 71 and 58% sequence identity with SmSULT. The reason for the differences in oxamniquine efficacy between species is unknown. Here we describe the enzyme kinetics of SULT from all three sources and present the structures of S. hematobium SULT (ShSULT) in complexes with PAP with and without oxamniquine and of S. japonicum SULT (SjSULT) in complex with PAP. The results provide a framework for redesign of oxamniquine as a pan-species drug treatment to be used in conjunction with praziquantel.

Structures of PAP complexes of schistosome sulfotransferases
Initial efforts to express full-length recombinant ShSULT were unsuccessful. Examination of the sequences of the three schistosome enzymes showed that ShSULT contains an additional N-terminal extension of 15 residues compared with SmSULT and SjSULT. Engineering this extension out of the ShSULT expression construct yielded a soluble protein amenable to crystallization. Although this modified ShSULT could be expressed using a construct containing a cleavable N-terminal His 8 tag, little soluble protein was produced from N-terminally tagged constructs for SjSULT, and refolded protein failed to crystallize. The addition of a short Gly-Ser linker plus a C-terminal His 8 tag to full-length SjSULT produced soluble protein that formed crystals. For both enzymes, the crystals were grown in the presence of PAP, yielding the structures of the enzyme-PAP complexes. Soaking the ShSULT crystals with racemic oxamniquine or the individual stereoisomers yielded the structures of the enzyme-PAP-oxamniquine complexes ( Fig. 1 and Table 1).
The overall structures of both ShSULT and SjSULT are similar to that of SmSULT (Protein Data Bank (PDB) code 4MUA). The ShSULT molecule superimposes onto SmSULT with an r.m.s.d. of 0.8 Å over 243 of 253 residues with no significant deviations other than a fully ordered loop (ShSULT residues 76 -80) in the ␣3-loop-␣4 region that is disordered in the SmSULT crystal structure. The N-terminal 17-20 residues are disordered. The SjSULT molecule superimposes onto SmSULT with an r.m.s.d. of 1.0 Å over 232 of 253 residues with one significant deviation in a C-terminal ␣ helix (␣12) caused by its domain swap into a neighboring molecule related by 2-fold crystal symmetry. The swapped ␣12 does not superimpose on the C-terminal ␣ helices of SmSULT or ShSULT associated with ␣11 in a helix-turn-helix motif. Rather, SjSULT ␣12 is displaced and packs against ␣2 (supplemental Fig. 1, A and B) with half of the C-terminal tag (Gly-Ser-His-His-His) visible in the electron density map. The tag region makes van der Waals interactions with the protein surface formed by ␣1-␤1-␣2 secondary structures. We speculate that this arrangement is an artifact caused by the domain swap. Threading the SjSULT C-terminal sequence onto the C-terminal structure of either SmSULT or ShSULT reveals that the SjSULT sequence is compatible with the ShSULT fold (supplemental Fig. 1C) provided small adjustments are made such as selecting common rotamers.
All three schistosome SULTs adopt an ␣/␤ fold characteristic to the PAPS-dependent SULT family of enzymes with a fivestranded parallel ␤ sheet sandwiched between two clusters of five and seven ␣ helices. A search of the Protein Data Bank (11) for structural homologs using the PDBeFold server (12) yields the human integral membrane tyrosylprotein sulfotransferase 2 (TPST2) (13) as the top match to SmSULT (Z score, 7.2; r.m.s.d., 2.8 Å for 157 of 253 residues), ShSULT (Z score, 6.8; r.m.s.d., 2.8 Å for 156 of 246 residues), and SjSULT (Z score, 7.9; r.m.s.d., 2.7 Å for 161 of 253 residues). Fig. 2A presents a secondary structure topology diagram of the schistosome SULT fold indicating elements in common with the mammalian SULTs. Schistosome SULTs have an additional three ␣ helices in common with TPST2 that are absent from the well characterized canonical mammalian SULTs in the SULT1, SULT2, and SULT4 families with known structures. Five ␣ helices in the schistosome SULT do not have a mammalian enzyme counterpart (Fig. 2, A-C).
The PAPS-binding site in the schistosome SULTs is almost identical to that seen in other members of the PAPS-dependent SULT family (14,15). A conserved arginine residue in the ␤1-loop-␣1 bridges the 3Ј-and 5Ј-phosphates of PAP, whereas a conserved serine residue in the same loop forms a hydrogen bond to N7 of the adenine (Fig. 3A). Separate conserved arginine and serine residues from the ␤4-loop-␣6 interact with the negative charge of the 3Ј-phosphate. The lysine residue that forms a hydrogen bond with a sulfate oxygen in PAPS in mammalian SULTs is replaced by a histidine in all three schistosome enzymes (Table 2 and Fig. 3, A and B).

Structures of S. hematobium sulfotransferase with oxamniquine plus PAP
To investigate the interaction of oxamniquine with ShSULT, the ShSULT-PAP crystals were soaked overnight with the individual stereoisomers or racemic oxamniquine. Difference electron density maps (F o Ϫ F c ) clearly indicated the presence of oxamniquine in the predicted substrate-binding region in all cases. In the crystals soaked with (S)-oxamniquine, the drug is in a similar position to that seen in the SmSULT-PAP-(S)oxamniquine complex (PDB code 5BYK) (Fig. 4A). ShSULT Figure 1. A, schematic representation of the ShSULT-PAP-oxamniquine (OXA) complex crystal structure. ␤ strands are colored orange, and ␣ helices are colored slate blue. B, schematic representation of the SjSULT-PAP complex crystal structure. Note that the predicted biological assembly is shown with a chain-swapped C-terminal ␣ helix from the crystal symmetry-related molecule colored magenta. For clarity, the C-terminal Gly-Ser His 8 affinity tag is not shown. C, walleyed stereoview of the superposition of the C␣ traces of ShSULT (dark green) and SjSULT (light orange) onto SmSULT (gray). Table 1 Data collection and refinement statistics UTHSCSA, University of Texas Health Science Center at San Antonio.   Asp-100 accepts a hydrogen bond from the putative reactive hydroxyl group. Ser-166 donates a hydrogen bond to the nitro group, and Asp-153 accepts a hydrogen bond from the amine in the isopropylaminomethyl group of oxamniquine. In the crystals soaked with the (R)-oxamniquine, the drug is in a position that can be described as a possibly unproductive binding mode. It is oriented similarly along the longest axis of the molecule when compared with (S)-oxamniquine, but its hydroxyl group is not in hydrogen-bonding distance to Asp-100. Rather, an intervening solvent molecule is positioned in hydrogen-bonding distance to both Asp-100 and the oxamniquine hydroxyl group (Fig. 4B). In both complexes, the pucker of the oxamniquine piperidine ring differs from that seen in SmSULT.

ShSULT
In the (S)-oxamniquine and (R)-oxamniquine complexes of ShSULT, carbon 3 is on the same side of the piperidine plane as the isopropylaminomethyl group. In both complexes of SmSULT, carbon 3 is on the opposite side of the plane (supplemental Fig. 2).
Oxamniquine is a racemic mixture in the formulation given as treatment to patients, although the S-stereoisomer has higher potency against S. mansoni in vivo (16). When racemic oxamniquine is soaked into crystals of SmSULT, the active site contains only the S-isomer. The ShSULT-PAP-racemic oxamniquine complex contains two molecules per asymmetric unit that superimpose with an r.m.s.d. of 0.5 Å. The active sites in both chains superimpose closely with the exceptions of Tyr-54, Ser-166, and the oxamniquine molecule. In chain B, both Tyr-54 and Ser-166 have similar positions to those of their counterparts in SmSULT, Phe-39 and Thr-157 (Fig. 5A). Oxamniquine also binds in a position similar to that seen in SmSULT, and the electron density clearly shows that only (S)oxamniquine is bound in ShSULT. In chain A, Tyr-54, Ser-166, and oxamniquine are in two conformations (Fig. 5B), and all refine to nearly 0.5 occupancy. For Tyr-54, one position is shifted away from the central cavity, and the other matches the single position adopted in chain B. One position of the side chain of Ser-166 matches its position in chain B, allowing formation of a hydrogen bond with the nitro moiety of oxamniquine. In the other position, the side chain is rotated so that the hydroxyl faces away from the oxamniquine molecule. For oxamniquine, one conformation matches the position of oxamniquine observed in chain B and in the crystal structure of ShSULT complexed with pure (S)-oxamniquine, but the other conformation is slightly rotated around an axis perpendicular to the piperidine plane. To summarize, in chain A, Tyr-54 positioned "in" toward the substrate pocket with the Ser-166 O␥ positioned "out" mimics the S. mansoni-like oxamniquine position, whereas Tyr-54 out and Ser-166 in promote an alternate oxamniquine position. Comparison with the ShSULT-PAP-(R)-oxamniquine and ShSULT-PAP-(S)-oxamniquine structures suggests that both oxamniquine conformations seen in chain A contain only (S)-oxamniquine. The electron density does not support modeling an (R)-oxamniquine as one of the two conformations depicted in Fig. 5B. Figure 2. A, diagram of secondary structure topology shared between the three schistosome SULTs illustrating the variation in fold as compared with mammalian SULTs. The ␤ strands and ␣ helices in yellow are common between the known schistosome and mammalian SULT protein structures, including SULTs 1A1, 1A3, 1B1, 1C1, 1C2, 1C3, 1D1, 1E1, 2A1, 2A3, and 4A1 and TPST2. Additionally, the ␣ helices in cyan are common between schistosome SULT and TPST2, but ␣ helices in red do not share a mammalian enzyme counterpart. N and C termini are labeled. B, diagram of ShSULT, TPST2 (PDB code 3AP1), and SULT1E1 (PDB code 1G3M) for comparison. The color scheme of the secondary structure elements pictured is as in A. Secondary structures of TPST2 and SULT1E1 without a schistosome counterpart are shown in gray. C, superposition of common secondary structure elements (shown in isolation) between schistosome and mammalian SULTs. Mammalian secondary structures are colored as in B, and ShSULT secondary structures are colored dark green.

Insights into schistosome parasite drug resistance Schistosome sulfotransferase sulfation kinetics
The specific activities of all three recombinant enzymes were determined using the common sulfotransferase substrate quercetin (17). With this flavanol, SmSULT and SjSULT have similar specific activities, whereas ShSULT has a higher specific activity by a factor of 0.3 (Table 3). Alanine substitutions of the putative catalytic base in SmSULT (D91A) or ShSULT (D100A) decrease the activity of the enzymes 5-and 11-fold compared with the wild-type enzymes, respectively.
To determine which oxamniquine stereoisomer is preferred by each of the enzymes, k cat /K m values were determined. In preliminary analyses with the racemic mixtures, significant substrate inhibition, frequently observed with sulfotransferases (18 -21), was seen with SmSULT and ShSULT but not SjSULT. This precluded determining reliable k cat and K m values for these two enzymes. For SmSULT, the highest activity with racemic oxamniquine of 9 min Ϫ1 was at 6 -7 M; ShSULT had the highest activity with racemic oxamniquine, 16 min Ϫ1 , at ϳ60 M; the k cat and K m values for SjSULT with the racemic substrate were 1.2 Ϯ 0.2 min Ϫ1 and 50 Ϯ 19 M. It was still possible to determine k cat /K m values in all cases by determining initial rates at low concentrations of each stereoisomer of oxamniquine ( Table 4). The results show significant differences among the abilities of the three enzymes to utilize the two isomers of oxamniquine as substrates. SmSULT does not show a significant preference for either stereoisomer. In contrast, (S)-oxamniquine is a better substrate by an order of magnitude for ShSULT, whereas SjSULT shows about a 3-fold preference for (R)-oxamniquine.

Characterization of the product of oxamniquine sulfation by schistosome SULTs
The sulfated product of the SULT reaction with oxamniquine is proposed to rapidly eliminate sulfate to form a highly    (Fig. 6A). A smaller signal (1) was also detectable at m/z 382.1043. This mass is consistent with the formula C 14 H 21 N 3 O 6 SNa (predicted m/z, 382.1049) and can thus be identified as a sulfated oxamniquine with an associated sodium ion. This ion was not detectable when PAPS or enzyme was omitted from the reaction or when the D91A mutant enzyme was used in a complete reaction (Fig. 6B). The same ion was produced in reactions with recombinant ShSULT (Fig. 6C). Thus, the enzymes from both sources are able to catalyze the transfer of sulfate from PAPS to oxamniquine to produce a sulfated product. No ion corresponding to the proposed ethylene product was detected. The mass spectra from the reactions of both enzymes also show a major species (2) with m/z 322.2240; expansion of the spectra shows that a second smaller ion (3) with m/z 322.1762 is also produced by both enzymes. Neither ion is formed by the inactive mutant enzyme, establishing that they are products of the SULT reaction. The m/z values of the two ions correspond to ions with the formulae C 16 H 28 N 5 O 2 (predicted m/z, 322.2238) and C 16 H 24 N 3 O 4 (predicted m/z, 322.1761), respectively. The formation of both 2 and 3 was time-dependent with the ion intensities being 2-3-fold greater at 6 s than at 1.2 s (supplemental Table 1). The formation of 2 can be attributed to nucleophilic displacement of the sulfate from oxamniquine sulfate by ethanediamine in the buffer (Fig. 6E), whereas 3 can be attributed to similar displacement of the sulfate by acetate. Consistent with this proposal, if the SmSULT reaction was carried out in ammonium acetate buffer without ethanediamine, 2 was not detectable, and significantly more oxamniquine sulfate was formed (Fig. 6D). If the ethanediamine was replaced with Tris as buffer at pH 7.0, 2 was similarly not detectable, but a new peak with m/z 383.2291 was seen, and the amount of oxamniquine sulfate increased, consistent with nucleophilic displacement of sulfate from oxamniquine sulfate by Tris (predicted m/z, 383.2289).

Discussion
The SULTs from the three main species responsible for almost all human blood fluke infections have the hallmark characteristics of PAPS-dependent SULTs, including the ␤ strandloop-␣ helix PAPS-binding motif. Structurally, the schistosome enzymes are closer to the membrane-bound mammalian tyrosine protein sulfotransferases such as TPST2 than to

Insights into schistosome parasite drug resistance
cytosolic mammalian sulfotransferases such as SULTE1 (Table 2), although they lack obvious sequences for membrane binding. The schistosome SULTs contain five ␣ helices without a counterpart in mammalian SULT structures; these helices enclose the active site and are likely to be important to the specificity of the schistosome enzymes. The enclosed active sites of the schistosome enzymes would be expected to prevent access of residues on protein surfaces in the absence of significant conformational changes, suggesting that the schistosome enzymes sulfate small molecule substrates.
The oxamniquine-binding sites of all three schistosome enzymes are located within cavities surrounded by ␣ helices 2, 6, 7, 8, 11, and 12. The cavities have the capacity to accommodate substrates of various sizes considering the internal volumes of 984, 985, and 1189 Å 3 for SmSULT, ShSULT, and SjSULT, respectively (note that the cavity volume for SjSULT is based on the arrangement of the domain-swapped C-terminal helix, ␣12, in the crystal structure). For comparison, the calculated molecular volume of oxamniquine is 212 Å 3 . The active sites of the schistosome enzymes are consistent with their utilizing a different set of substrates from the mammalian enzymes (Figs. 2 and 3). The mammalian SULTs have structural elements that would appear to impede productive oxamniquine binding (supplemental Fig. 3), suggesting that oxamniquine is not a substrate. Indeed, a study of metabolism of oxamniquine in several mammals, including humans, showed that the major metabolite excreted in urine has a carboxylate substitution for the hydroxymethyl group, consistent with metabolism by a cytochrome P 450 rather than a SULT (23). Thus, the selective toxicity of oxamniquine toward S. mansoni and not humans is likely attributable to the structural variation between the SULT homologs and the availability of an alternative detoxification pathway in humans.
The mechanism of SULTs is proposed to involve nucleophilic attack of the deprotonated substrate hydroxyl group on the sulfur of PAPS, resulting in sulfuryl transfer from PAPS to the substrate (13, 14, 24 -26). Mutational analyses of conserved active-site residues in TPST2 support catalytic functions for Arg-78, Glu-99, Lys-158, and Ser-285 (13). Superposition of the structures of SmSULT, ShSULT, and SjSULT onto that of TPST2 identifies the counterparts of these residues in the schistosome enzymes (Fig. 2B and Table 2). All three contain an active-site aspartate in place of Glu-99 appropriately positioned to act as the catalytic base to abstract a proton from the hydroxyl group of oxamniquine and other substrates (Fig. 2, A  and B). The finding that substitution of this residue to alanine in SmSULT and ShSULT reduces activity supports this functional assignment. The schistosome enzymes also contain a conserved arginine in the same position as Arg-78 in TPST2. This residue was proposed to donate a proton to the bridging oxygen of the phosphosulfate moiety in PAPS (13), but the computational analysis of Marforio et al. (25) is more consistent with this residue acting as a shuttle for the SO 3 Ϫ as it is transferred from donor to acceptor. The latter role is also more consistent with the high pK a of arginine. The asparagine found in the schistosome enzymes in place of Ser-285 is positioned to play a similar role of stabilizing the transition state by forming a strong hydrogen bond with a non-bridging oxygen on the phosphate (13,25). Lys-158, proposed to stabilize the transition state, is replaced with histidine in the schistosome SULTs. This residue could readily serve as a positive charge to stabilize the transition state if the imidazole ring were protonated.
When racemic oxamniquine is soaked into crystals of SmSULT, the active site contains only the S-isomer, establishing that this stereoisomer binds more tightly to that enzyme. One possible explanation for oxamniquine inefficacy in S. hematobium is that ShSULT preferentially binds (R)-oxamniquine. However, soaking racemic oxamniquine into ShSULT crystals similarly shows (S)-oxamniquine bound. (R)-Oxamniquine does bind but adopts a potentially unproductive binding mode where its hydroxyl group is separated from the catalytic aspartate by an intervening solvent molecule (Fig. 4B). The propensity for this binding mode may correlate to the observation that k cat /K m for (R)-oxamniquine is an order of magnitude lower than for (S)-oxamniquine as an ShSULT substrate.
One of the two molecules in the asymmetric unit in the ShSULT crystals with racemic oxamniquine shows (S)-oxamniquine bound in more than one conformation, whereas in the other subunit the (S)-oxamniquine position is comparable with the oxamniquine position in the (S)-oxamniquine crystal soak of both SmSULT and ShSULT (Figs. 4A and 5A). The combination of side chain rotamers of Tyr-54 and Ser-166 in ShSULT appears to select the position for bound (S)-oxamniquine (Fig.  5B). In SmSULT, the equivalent residue to Ser-166 is Thr-157, and its methyl group may stabilize the side chain where it borders a hydrophobic patch comprising Leu-92, Val-128, Leu-129, and Leu-161. The combination of a stable Thr-157 rotamer and a Phe-39 unable to donate a hydrogen bond to its neighboring asparagine likely imposes a more rigid active site for oxamniquine binding in SmSULT than in ShSULT. The crystal structures of oxamniquine-bound SmSULT and ShSULT also differ in that ShSULT-bound oxamniquine adopts a piperidine ring pucker opposite to SmSULT-bound oxamniquine (supplemental Fig. 2), although it is not clear that this plays a significant role in oxamniquine affinity.
A number of residues in contact with oxamniquine in the ShSULT-PAP-oxamniquine and SmSULT-PAP-oxamniquine crystal structures are different in SjSULT (supplemental Fig. 4). Although SjSULT still has the aspartate that serves as the active-site base (Asp-87), an asparagine (Asn-140) replaces an aspartate seen in the other two enzymes (Asp-153 in ShSULT and Asp-144 in SmSULT). In addition, SjSULT contains Asn-153 in place of Ser-166 (ShSULT) and Thr-157 (SmSULT). It is not clear that the two asparagine substitutions would prevent oxamniquine binding to SjSULT as the asparagine side chain is still capable of forming a hydrogen bond with the oxamniquine nitro and amine groups. The nearby substitution of Val-139 for glycine (SmSULT Gly-143 and ShSULT Gly-152) likely has a more significant impact on binding. The additional bulk of SjSULT Val-139 appears to create steric clash with the isopropylaminomethyl group of oxamniquine according to structural superpositions with SmSULT-PAP-oxamniquine (supplemental Fig. 4) and ShSULT-PAP-oxamniquine complexes.
The continuous-flow mass spectral analyses establish directly that schistosome SULTs catalyze a sulfotransferase reaction with oxamniquine as substrate. The observation of additional products resulting from nucleophilic displacement of sulfate by buffer components is consistent with the oxamniquine sulfate generated by the enzymes being highly susceptible to nucleophilic attack. Oxamniquine has been proposed to exert its schistosomicidal activity upon the decay of the sulfated product of the SULT reaction to a reactive ethylene with toxic alkylating activity within the parasite (17,22). Assays using tritiated oxamniquine show the formation of labeled macromolecules in worm extracts, consistent with activation of oxamniquine by a SULT generating a reactive species (10,17,22). However, the variation in the amount of oxamniquine sulfate when the buffer is changed is more consistent with an S N 2-like mechanism in which a nucleophile displaces the sulfate in a concerted reaction than with an S N 1-like mechanism in which loss of sulfate generates a discrete ethylene intermediate that is subsequently attacked by a nucleophile. Thus, the present data suggest an alternative hypothesis for oxamniquine toxicity in which macromolecules in the parasite participate directly in the nucleophilic displacement of sulfate from activated oxamniquine.
The kinetic analyses confirm the activity of all three schistosome enzymes with oxamniquine as substrate and allow direct comparison of their enzymatic activities with oxamniquine. Based on the k cat /K m values, SmSULT has the highest activity with oxamniquine as substrate and shows no significant preference for either stereoisomer. ShSULT is about one-half as active as SmSULT with (S)-oxamniquine but has much lower activity with the R-isomer. SjSULT is less active with oxamniquine as substrate by about an order of magnitude compared with SmSULT and shows a preference for the R-isomer.
These differences may account for the efficacy of oxamniquine as a drug against S. mansoni infection and its failure against S. hematobium and S. japonicum.
The reduction in ShSULT catalytic efficiency toward oxamniquine may be related to its ability to bind oxamniquine in multiple positions (Figs. 4 and 5) because the active sites of the two enzymes are quite similar. In the case of SjSULT, the substitution of Gly-143 with Val-139 and the substitution of Asp-153 with Asn-140 (supplemental Fig. 4) provide a structural rationale for the significant decrease in catalytic efficiency toward oxamniquine. If this is correct, small modifications in the structure of oxamniquine that avoid the clash with Val-139 may generate a more effective drug.
In summary, the present results provide a rationale for the lack of efficacy of oxamniquine with S. hematobium and S. japonicum and provide structural and kinetic information that will be critical to develop oxamniquine derivatives with efficacy against these two species. We speculate that efficient turnover of oxamniquine in S. mansoni creates a large enough pool of toxic product to kill the blood fluke, whereas oxamniquine is not effective against S. hematobium and S. japonicum due to its diminished turnover. The finding that ShSULT is almost as active as SmSULT with a very similar active site suggests that small structural changes to oxamniquine may be sufficient to develop an improved derivative. In addition, we have identified a possible structural basis for the low enzymatic activity of SjSULT with oxamniquine that will be useful in developing a derivative targeting that enzyme. The structures of all three schistosome SULTs show that their active sites are similar to but not identical to their mammalian counterparts, providing a rationale for the different effects in humans and worms. On the whole, the crystal structures of the schistosome SULTs and the evidence for in vitro oxamniquine binding and turnover in both ShSULT and SjSULT support the initiative for developing an oxamniquine derivative as a pan-species antischistosomal drug.

Cloning, expression, and purification
DNA encoding the ShSULT gene was subcloned into the pAG8H vector derived from pKM260 (27), which contains an inducible lacZ promoter and an N-terminal His 8 tag that can be released from the target protein via a tobacco etch virus (TEV) protease cleavage site. The plasmid was transformed into Escherichia coli strain BL21 pLysS (Promega, Madison, WI), and the cells were grown at 37°C. When the A 600 reached 0.7, the temperature was decreased to 18°C, and expression was induced by the addition of isopropyl ␤-D-thiogalactoside to a final concentration of 1 mM. Cells from 6 liters of cell culture were washed and resuspended in 50 ml of 50 mM Tris, pH 8.0, 500 mM NaCl (column buffer) containing 250 l of Sigma protease inhibitor mixture (Sigma-Aldrich)/liter of cultured cells and lysed by sonication on ice. The full-length ShSULT localized in inclusion bodies, and little soluble protein was expressed. Sequence analysis indicated that ShSULT contained an N-terminal extension compared with SmSULT; thus, to attempt to alleviate the solubility problem, a more compact ShSULT construct was designed with 16 residues truncated at the N terminus. The truncated ShSULT was successfully overexpressed and localized in the soluble fraction. The clarified supernatant was loaded onto a Ni 2ϩ -nitrilotriacetic acid affinity chromatography column (GE Healthcare), washed with 5 volumes of column buffer, and eluted using a 10 -500 mM imidazole gradient. Fractions identified as ShSULT by SDS-PAGE were pooled, and His-tagged TEV protease was added in a protein to protease ratio of 15:1. The resulting solution was dialyzed overnight at 4°C against 50 mM Tris, pH 8.0, 100 mM NaCl, 2 mM dithiothreitol. The dialysate was passed over the Ni 2ϩ -nitrilotriacetic acid column again to remove the His tag and the His-tagged TEV protease, and the cleaved target protein flowed through. The sample was loaded onto a GE Healthcare prepacked Q anionexchange column and eluted with a 0.1-1.0 M NaCl gradient. Fractions containing ShSULT as identified by SDS-PAGE were pooled and dialyzed overnight at 4°C in 25 mM Tris, pH 8.5, 50 mM NaCl, 2 mM tris(carboxyethyl)phosphine. The product was ϳ98% pure as estimated by SDS-PAGE. When attempts to produce sufficient soluble N-terminal His 8 -tagged SjSULT failed, DNA encoding the SjSULT gene was subcloned into the pNC8H vector derived from pKM260, which contains an inducible lacZ promoter and a non-cleavable C-terminal His 8 tag with a two-residue Gly-Ser linker. SjSULT was prepared as described above for ShSULT omitting the affinity tag cleavage step. Purified ShSULT and SjSULT proteins were concentrated to 10 mg/ml for crystallization trials. PAP was added in a 4:1 stoichiometric ratio over protein, and the sample was incubated for 1 h on ice immediately prior to crystallization screening. The D100A mutation was introduced into the pAG8H-ShSULT plasmid using the Phusion (Thermo Fisher Scientific, Waltham, MA) site-directed mutagenesis protocol, and the protein was prepared as described above. The D91A mutation was introduced into the pAG8H-SmSULT vector in the same way, and wild-type SmSULT and SmSULT D91A proteins were prepared as described previously (10).

Crystallization, structure determination, and refinement
Automated screening for crystallization was carried out using the sitting-drop vapor-diffusion method with an Art Robbins Instruments (Sunnyvale, CA) Phoenix system in the X-ray Crystallography Core Laboratory at the University of Texas Health Science Center at San Antonio. ShSULT-PAP crystals were obtained from the Molecular Dimensions (Suffolk, UK) Morpheus screen condition number 1-12. ShSULT-PAPoxamniquine complex crystals were prepared by soaking PAPcontaining crystals overnight in the mother liquor saturated with racemic oxamniquine or the individual (R)-and (S)-oxamniquine enantiomers prepared as described previously (16). SjSULT-PAP crystals were obtained from the Rigaku (Tokyo, Japan) Precipitant Synergy 2 screen condition number 151. The SjSULT-PAP crystals grew as polycrystalline clusters with an epitaxially twinned core branching into single crystal rods, and only an isolated rod produced a useful single crystal diffraction pattern. Attempts to soak oxamniquine into isolated SjSULT-PAP crystals returned damaged crystals and unproductive X-ray diffraction. All crystals were mounted in undersized nylon loops with excess mother liquor wicked off and flash-cooled in liquid nitrogen prior to data collection. The structures of ShSULT and SjSULT were determined by the molecular replacement method implemented in Phaser (28) using the SmSULT structure (PDB code 4MUA) as the search model. Coordinates for all models were refined using PHENIX (29), including simulated annealing with torsion angle dynamics, and alternated with manual rebuilding using Coot (30). The occupancies for oxamniquine coordinates were held at 1.0 during refinement with the exception of the racemic oxamniquine soaked into ShSULT crystals where occupancies for alternate conformations were refined independently. Data were collected using our home source Rigaku MicroMax 007-HF equipped with R-AXIS-HTC detectors and at the Advanced Photon Source Northeastern Collaborative Access Team (NE-CAT) beamlines 24-ID-C and 24-ID-E and integrated and scaled using XDS (31). Data collection and refinement statistics are shown in Table 1. Structural pocket volume calculations were performed using the Computed Atlas of Surface Topography of Proteins (CASTp) server (32). Figures were generated using PyMOL (Schrödinger, LLC), TopDraw (33), and Chem-Draw (PerkinElmer Life Sciences).

Sulfotransferase assay
Standard assay conditions were 10 M quercetin, 0.05 M inositol monophosphatase 3, 50 mM MOPS, pH 7.0, 15 mM MgCl 2 at 37°C. Protein concentrations in the assays were 0.2 M for SmSULT, ShSULT, and SjSULT and 1.2 M for SmSULT D91A and ShSULT D100A. Assays were started by adding 25 l of 100 M PAPS in MOPS/MgCl 2 buffer to 25 l of the reaction mixture. For standard assays, 10-l aliquots were taken at 5-min intervals and quenched with 6 l of malachite green reagent A from the Universal Sulfotransferase Activity kit (R&D Systems Inc., Minneapolis, MN). Water (20 l) and 6 l of malachite green reagent B were then added; after 20 min, the absorbance at 620 nm was determined. To obtain the k cat /K m values for oxamniquine stereoisomers, the concentration of each was varied over a concentration range where the rate was linearly dependent on the substrate concentration (1-5 M for SmSULT and 10 -50 M for ShSULT and SjSULT). Aliquots were taken from the reactions each minute for 4 min to determine the initial rate of product formation. Protein concentrations in the assays were 0.1-0.2 M for SmSULT and ShSULT and 0.8 -1.2 M for SjSULT. All assays were corrected for background hydrolysis of PAPS based on controls from which enzyme was omitted. The respective k cat /K m values were determined from weighted linear fits of the initial rates as a function of substrate concentration. Data from multiple experiments were combined to obtain the values in Table 4.

Continuous-flow mass spectrometry
The products of the reaction of SmSULT and ShSULT with oxamniquine were detected by continuous-flow mass spectrometry using a Thermo Scientific LTQ OrbiTrap Discovery interfaced with a New Objectives (Woburn, MA) PicoView nanospray source and a Peltier-controlled nanovolume continuous-flow mixer (Eksigent, Dublin, CA) (34). A solution of 20 M SULT in buffer in one line was continuously mixed with an equal volume of 40 M oxamniquine, 100 M PAPS, 10% meth-anol, 5 mM ammonium acetate, pH 7.0, from a second line at 25°C and injected into the mass spectrometer after 1.2 or 6 s. Reactions were monitored in positive-ion spectral mode from m/z 110 to 600.