If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Department of Biochemistry and Molecular Biology, University of Georgia, Athens, Georgia, USACenter for Tropical and Emerging Global Diseases, University of Georgia, Athens, Georgia, USA
Department of Biochemistry and Molecular Biology, University of Georgia, Athens, Georgia, USACenter for Tropical and Emerging Global Diseases, University of Georgia, Athens, Georgia, USA
Department of Biochemistry and Molecular Biology, University of Georgia, Athens, Georgia, USAComplex Carbohydrate Research Center, University of Georgia, Athens, Georgia, USA
Department of Biochemistry and Molecular Biology, University of Georgia, Athens, Georgia, USAComplex Carbohydrate Research Center, University of Georgia, Athens, Georgia, USA
Department of Biochemistry and Molecular Biology, University of Georgia, Athens, Georgia, USACenter for Tropical and Emerging Global Diseases, University of Georgia, Athens, Georgia, USAComplex Carbohydrate Research Center, University of Georgia, Athens, Georgia, USA
Skp1, a subunit of E3 Skp1/Cullin-1/F-box protein ubiquitin ligases, is modified by a prolyl hydroxylase that mediates O2 regulation of the social amoeba Dictyostelium and the parasite Toxoplasma gondii. The full effect of hydroxylation requires modification of the hydroxyproline by a pentasaccharide that, in Dictyostelium, influences Skp1 structure to favor assembly of Skp1/F-box protein subcomplexes. In Toxoplasma, the presence of a contrasting penultimate sugar assembled by a different glycosyltransferase enables testing of the conformational control model. To define the final sugar and its linkage, here we identified the glycosyltransferase that completes the glycan and found that it is closely related to glycogenin, an enzyme that may prime glycogen synthesis in yeast and animals. However, the Toxoplasma enzyme catalyzes formation of a Galα1,3Glcα linkage rather than the Glcα1,4Glcα linkage formed by glycogenin. Kinetic and crystallographic experiments showed that the glycosyltransferase Gat1 is specific for Skp1 in Toxoplasma and also in another protist, the crop pathogen Pythium ultimum. The fifth sugar is important for glycan function as indicated by the slow-growth phenotype of gat1Δ parasites. Computational analyses indicated that, despite the sequence difference, the Toxoplasma glycan still assumes an ordered conformation that controls Skp1 structure and revealed the importance of nonpolar packing interactions of the fifth sugar. The substitution of glycosyltransferases in Toxoplasma and Pythium by an unrelated bifunctional enzyme that assembles a distinct but structurally compatible glycan in Dictyostelium is a remarkable case of convergent evolution, which emphasizes the importance of the terminal α-galactose and establishes the phylogenetic breadth of Skp1 glycoregulation.
A prominent mechanism of O2 sensing in metazoa involves an O2-dependent prolyl 4-hydroxylase (PHD2) that generates degrons on hypoxia-inducible factor-α (HIFα) that are recognized by the von Hippel Lindau subunit of the E3(VBC) ubiquitin ligase, leading to its polyubiquitination and degradation in the 26S proteasome (
). However, protists lack HIFα, and protist PhyA enzymes instead hydroxylate Skp1, a subunit of the Skp1/Cullin-1/F-box protein (SCF) family of E3 ubiquitin ligases (
) that are related to the E3(VBC) ubiquitin ligases. In the social amoeba Dictyostelium discoideum, Skp1 prolyl hydroxylation is involved in mediating the O2 checkpoint for fruiting body formation (
). Remarkably, the hydroxyproline (Hyp) residue is sequentially modified by a series of five glycosyltransferase (GT) reactions leading to the assembly of a linear pentasaccharide of recently defined structure (
), and radiotracer and biochemical complementation studies indicate that Skp1 is the sole target of these GTs. Skp1 serves as an adaptor that links the F-box protein (FBP) and Cullin-1 subunits of the SCF complex through mostly independent binding events. Interactome studies indicate that glycosylation promotes accumulation of Skp1/FBP subcomplexes in vitro and in vivo (
), thereby potentially activating the respective SCF complexes and contributing to the degradation of all of the substrates that are recognized by the cellular repertoire of FBPs. Current evidence indicates that the glycan assumes a relatively rigid structure that organizes the intrinsically disordered subsite-2 region involved in FBP recognition. By increasing the fraction of time that the C-terminal region is folded into α-helix-8 and extended to form subsite-2, the glycan is thought to promote the interaction of Skp1 with FBPs in the cell (
). However, the importance of the exact glycan sequence and the significance of the terminal (fifth) sugar deserve further attention because there is limited precedent for glycan-imposed order on local protein structure (
). Although chronic infections are usually clinically benign, reactivation, as can occur in immunosuppressed individuals, can lead to serious diseases of the central nervous system and other organs. Moreover, human fetuses are subject to serious birth defects in the case of acute maternal infection. Previous studies documented that the T. gondii prolyl hydroxylase (phyA) can complement a disruption of phyA in D. discoideum and contribute to O2-dependent growth of T. gondii in a fibroblast monolayer growth assay (
). Furthermore, coding sequences for the first three GT activities are present, although, surprisingly, T. gondii still assembles a pentasaccharide glycan on its Skp1 (
) and, as reported here, another enzyme that mediates addition of the final sugar. Analytical studies revealed that the fourth sugar is different from that of D. discoideum, but the fifth has remained unknown. These findings have raised important questions, including why are the GTs and glycans different between these two protists, what constrains their evolution, and are they compatible with the conformational regulation model posited for D. discoideum?
Pythium ultimum is an oomycete plant pathogen that resides in the stramenopile branch of the larger TSAR group of protozoans to which the apicomplexans also belong. P. ultimum is the agent for root rot disease in agriculturally important crops (
). We recently showed P. ultimum PhyA and the first GT Gnt1 constitute a bifunctional protein that is active toward Skp1A but not the second Skp1B encoded by its genome (
). Thus O2-dependent post-translational regulation of Skp1 appears to be widespread, and P. ultimum and protists from other lineages have evidently evolved a second Skp1 that avoids O2-dependent glycoregulation. We show here that the terminal glycosylation of P. ultimum Skp1 is like that of T. gondii rather than that of D. discoideum.
By identifying the Gat1 GT that mediates the addition of the final sugar to the T. gondii, we have been able to determine the complete glycan sequence and to show that Gat1 is specific for Skp1, regulates Skp1, and contributes to optimal growth in cell culture. Despite a distinct structure, the glycan mediates a conserved conformational effect on Skp1. We propose that the selection pressure for a 5-sugar chain was so strong that, in lineages where the Gat1 progenitor was evolving to acquire a new function, the later-evolving amoebozoa acquired a novel GT mechanism to ensure assembly of a glycan that was structurally similar enough to still conformationally control Skp1.
Results
The fifth and final Skp1 sugar is an α-linked Gal and depends on TgGat1
Previous studies described the mechanism of assembly of the first four sugars on TgSkp1-Pro-154 (Fig. 1A), but the left the identity of the final sugar (other than its being a hexose) unresolved (
). The corresponding sugar in Dictyostelium (Fig. 1B) is a 3-linked αGal that is susceptible to removal with green coffee bean α-galactosidase. Similar treatment of tryptic peptides from TgSkp1, isolated from the standard type 1 strain RH by immunoprecipitation, resulted in complete conversion of the pentasaccharide form of the glycopeptide to the tetrasaccharide form, indicating that the terminal Hex is an αGal residue (Fig. 1C).
Figure 1Gat1 is required for terminal α-galactosylation of Skp1 in parasites.A and B, schematic of Skp1 glycosylation pathway in Toxoplasma (
). The identity of Gat1 as the final enzyme in the Toxoplasma pathway and the nature of the final sugar are reported here. C, dependence of the terminal sugar on Gat1 and its characterization. Skp1 was immunoprecipitated from type 1 RH and type 2 ME49 strains and their genetic derivatives, and peptides generated by trypsinization were analyzed by nano-LC-MS. In addition, Skp1 from RH was treated with green coffee bean α-gal after trypsinization. The values represent the levels of pentasaccharide-peptide and tetrasaccharide-peptide levels detected, after normalization to all detected modification states of the peptide; note that the values have only partial relative meaning because of the low and varied detection efficiency of glycopeptides. As indicated, the detection threshold was <0.0005. The open circle for the terminal sugar of the pentasaccharide indicates that it was only known as a hexose at the time of the experiment. See Table S2 and Fig. S5 for primary data. Similar results were obtained in independent samples from RH and Δgat1/RH (not shown).
), and one of these, referred to as Gat1, was predicted to be the missing Skp1 GT on account of the phylogenetic co-distribution of its gene in protists that possess Glt1-like genes (
). To test the dependence of Skp1 glycosylation on Tggat1, the gene was disrupted using CRISPR/Cas9 in the RH strain, yielding gat1Δ as described in Fig. S4A. PCR studies confirmed replacement with a dhfr cassette, and enzyme assays (see below) showed a loss of enzyme activity. The resulting clone produced a version of Skp1 in which the 5-sugar glycopeptide was no longer detectable, but ions corresponding to the 4-sugar glycopeptide were detected at similar abundance (Fig. 1C). Similar results were obtained in the type 2 ME49 strain (Fig. S4A and Fig. 1C) and in strain RHΔΔ in which gat1 was disrupted by homologous recombination using a different selection marker, referred to as Δgat1-1 (Table 1 and Fig. S3B; data not shown). The similar results obtained by different genetic methods in distinct genetic backgrounds indicate that Tggat1 is required for the addition of the terminal sugar, but whether this was a direct or indirect effect was unclear.
). In two-dimensional fibroblast monolayers, parasites lyse out and invade neighboring cells to repeat the cycle, resulting in plaques whose area is a measure of efficiency of these cellular processes. Past studies showed that plaque size growth is compromised by mutational blockade of Skp1 hydroxylation and earlier steps of the glycosylation pathway (
). Similarly, disruption of gat1 also resulted in modestly smaller plaques in both RH or RHΔΔ backgrounds (Fig. 2, A and B) and in RHΔΔ in which gat1 was disrupted using homologous recombination without CRISPR-Cas9 (Fig. 2C). To determine whether the effects were specific to the genetic lesion at the gat1 locus, the RH and RHΔΔ KO strains generated using CRISPR/Cas9 were modified again by CRISPR/Cas9 to introduce single copies of epitope-tagged versions of the gat1 coding locus, downstream of an endogenous gat1 promoter cassette or a tubulin promoter cassette, respectively, into the uprt locus. The expected insertions were confirmed using PCR (Figs. S3C and S4B). As a result, TgGat1-3×HA could be detected by Western blotting of tachyzoites in the RHΔΔ background (Fig. S4D), and, as discussed below, the complemented RHΔΔ strain restored Skp1 glycosylation according to a biochemical complementation test. Although TgGat1-Ty expressed under its own promoter cassette in the RH background was not detected (not shown), enzyme activity was partially restored in the TgGat1-Ty strain (Fig. S4C). Both strains exhibited larger plaque sizes than their respective KO parents (Fig. 2, A and B), confirming that the effect on growth was due to the original loss of Gat1.
Figure 2Parasite growth depends on Gat1. Parasites were plated at clonal density on two-dimensional monolayers of HFFs and allowed to invade, proliferate, lyse out, and reinfect neighboring fibroblasts. After 5.5 days, cultures were fixed, stained, and analyzed for the areas occupied by lysed fibroblasts (plaques). Data from each of three independent trials, which were each normalized to the parental strain, were merged for presentation. A, comparison of the type 1 RH strain before and after gat1 replacement using CRISPR/Cas9 and complementation with gat1 under control of its own promoter cassette in the uprt locus. B, comparison of RHΔΔ, gat1-2Δ/RHΔΔ, and the latter complemented with gat1 under control of a tubulin promoter in the uprt locus. C, comparison of gat1-1Δ/RHΔΔ, prepared by homologous recombination. Significance of differences in plaque areas between parasite strains was assessed by Student's t test. CI, confidence interval.
An evolutionary analysis was conducted to gain further insight into the function of TgGat1. Based on searches of genomic databases using BLASTP, TgGat1 belongs to the CAZy GT8 family. The top-scoring hits, with Expect values of <10−32, were found only in protists that contain Toxoplasma PgtA-like and Glt1-like sequences (Fig. S7) and lack D. discoideum AgtA-like sequences, suggesting a common function. The most similar sequences, in searches seeded with the putative catalytic domain, belong to glycogenin, with Expect values of ≥E−27. All other homologous sequences had Expect values of ≥10−22. Glycogenin is a dimeric α4-glucosyltransferase that can prime the synthesis of glycogen in the cytoplasm of yeast and animals by a mechanism that involves autoglucosylation. Glycogenin appears not to be involved in starch formation and is an evolutionarily recent addition to the glycogen biosynthesis pathways that occur in the absence of glycogenin activity in bacteria and many unicellular eukaryotes (
Developmentally regulated biosynthesis of carbohydrate and storage polysaccharide during differentiation and tissue cyst formation in Toxoplasma gondii.
), an α1,4Glc polymer with α1,6-linked branches that resembles glycogen, in its cytoplasmic compartment. T. gondii amylopectin is assembled by a UDP-Glc–based metabolism that is related to the floridean starch of the red alga Cyanidioschyzon merolae and, to a lesser extent, to that of glycogen storing animals and fungi. Homologs of glycogen synthase (
Developmentally regulated biosynthesis of carbohydrate and storage polysaccharide during differentiation and tissue cyst formation in Toxoplasma gondii.
). Thus, Gat1 might have a function in regulating T. gondii amylopectin formation.
Glycogenin consists of a CAZy GT8 family catalytic domain plus a C-terminal glycogen synthase–binding domain separated by a linker, whereas Gat1 consists only of a single catalytic domain (Fig. 3A). TgGat1 is predicted to be a 345-amino acid protein encoded by a single exon gene in the Type I GT1 strain (TGGT1_310400). It is 34% identical to rabbit (Oryctolagus cuniculus) glycogenin over the catalytic domain but includes a poorly conserved 90-amino acid sequence that interrupts the catalytic domain (Fig. 3A). This region is likely unstructured based on secondary structure prediction by the XtalPred server (
). The Gat1-like sequence from P. ultimum (Uniprot K3WC47) lacks this sequence, so it was analyzed for comparison. PuGat1 is predicted to be a 266-amino acid protein encoded by a two-exon gene, annotated as PYU1_G002535-201 (transcript ID PYU1_T002538).
Figure 3Gat1 is closely related to glycogenin sequences.A, domain analysis of Gat1 from T. gondii and P. ultimum compared with human glycogenin-1. B, the evolutionary history of the sequence of the Gat1 catalytic domain was inferred by using a maximum likelihood method. The tree with the highest log likelihood (−13,279.56) is shown. Gat1 and Gat1-like sequences are colored green, glycogenin and glycogenin-like sequences are in red, and characterized and other selected other CAZy GT8 sequences are in black, or in purple if predicted to reside in the secretory pathway rather than the cytoplasm. The percentage of trees in which the associated taxa clustered together is shown at each branch. Branch lengths are measured by the number of substitutions per site. See Figs. S6–S8 for alignments.
To further evaluate the evolutionary relationship of these putative GTs, their catalytic domains and those of the most closely related or known sequences from the CAZy GT8 family (Fig. S7) were aligned (Fig. S8) and analyzed by a maximum likelihood method (Fig. 3B). The results suggest that Gat1 and glycogenin evolved separately from a common ancestor. Although the last common ancestor was not resolved, Gat1 was potentially the predecessor to glycogenin owing to its presence in other primitive unicellular eukaryotes that bear no evidence of glycogenin-like sequences and because there is currently no evidence for the existence of close homologs of Gat1 and glycogenin within the same clade that would suggest that they evolved as paralogs of a gene duplication. However, the possibility that either one or the other product of an ancestral gene duplication was always lost in every extant derivative cannot be excluded. Gat1 and glycogenin each possess unique conserved sequence motifs (Fig. S6) that potentially support functional differences.
Gat1 is a terminal Skp1 α-galactosyltransferase
To determine whether TgGat1 can directly modify Skp1, the predicted full-length protein was expressed as a His6-tagged conjugate in Escherichia coli, purified on a TALON resin, and treated with TEV protease, leaving an N-terminal GlyHis-dipeptide stub before the start Met (Fig. 4A). The presumptive ortholog from P. ultimum was prepared similarly. A screen for UDP-sugar hydrolysis activity of TgGat1 yielded, after extended reaction times, only UDP-Gal and UDP-Glc as candidate substrates from a panel of six common UDP-sugars (Fig. S9A). A quantitative comparison showed ∼7-fold greater activity toward UDP-Gal than UDP-Glc (Fig. 4B).
Figure 4Gat1 preferentially galactosylates the Skp1 glycan in vitro.A, recombinantly expressed and purified preparations of TgGat1 and PuGat1 were analyzed by SDS-PAGE and staining with Coomassie blue. B, temporal dependence of UDP-Gal and UDP-Glc hydrolysis. The averages and S.D. of three technical replicates are shown. A similar profile was observed with a different enzyme concentration. See Fig. S9E for a trial with higher enzyme concentrations. C, transferase activity utilizing 8 μm UDP-Gal or UDP-Glc toward 20 mm Glcα1,4Glcα-pNP (maltose-pNP) for TgGat1 and PuGat1. The averages and S.D. of two technical replicates are shown; similar profiles were in two independent assays with a different TgGat1 preparation. D, UDP-Gal and UDP-Glc concentration dependence of TgGat1 transferase activity toward 20 mm maltose-pNP. The averages and S.D. of two technical replicates are shown, and an independent trial with TgGat1 and PuGat1 against UDP-Gal is shown in Fig. S9F. E, maltose-pNP concentration dependence of TgGat1 and PuGat1 transferase activity from 20 μm UDP-Gal. The averages and S.D. of two technical replicates are shown. F, relative Gal-transferase activity of TgGat1 and PuGat1 toward different acceptors. The averages and S.D. of three technical replicates are shown. Similar results were obtained in three independent trials. G, effect of UDP-Glc concentration on the Gal-transferase activity of TgGat1. Reactions were incubated for 1 h. The averages and S.D. of two technical replicates are shown. H, Gal-transferase activity of TgGaT1 toward varied concentrations of GlFGaGn-Skp1, in the presence of 40 μm UDP-Gal (1 µCi) after a 1-h incubation. Data from independent preparations of TgSkp1 are colored in different shades. FGaGn-Skp1 is included for comparison. Error bars, S.D. of duplicate measurements. Inset, Western blots of the Skp1 preparations used, where FGaGn-Skp1, which is recognized specifically by pAb UOK104, is largely converted in a 3.5-h reaction using Glt1 and UDP-Glc to GlFGaGn-Skp1, which is recognized only by the pan-specific pAb UOK75. I, reactions with synthetic oligosaccharides conjugated to pNP were conducted in parallel using the same conditions. J, biochemical complementation to detect Gat1 substrates. Desalted S100 extracts of RH and gat1Δ/RH were reacted with recombinant Gat1 in the presence of UDP-[3H]Gal, and the product of the reaction was separated on an SDS-polyacrylamide gel, which was sliced into 40 bands for liquid scintillation counting. The migration position of Skp1 is marked with an arrow. See Fig. S9 (H and I) for trials using different strains.
The ability of Gat1 to transfer Gal or Glc to another sugar, rather than water, was tested using a substrate analog for glycogenin, Glcα1,4Glcα1-pNP (maltose-pNP), which mimics the terminal disaccharide of glycogen and starch and has a terminal αGlc, as found on the Skp1 tetrasaccharide. Although Gat1 from either T. gondii or P. ultimum could modify maltose-pNP using either sugar nucleotide (Fig. 4C), the enzymes strongly preferred UDP-Gal (Fig. 4D). Furthermore, TgGat1 activity was not saturated by UDP-Glc at 0.5 mm (Fig. 4D), whereas TgGat1 and PuGat1 exhibited apparent Km values for UDP-Gal in the range of 30-15 μm (Fig. S9F) These values were greater than the 2-4 μm values reported for rabbit and yeast glycogenins for UDP-Glc (
). TgGat1 and PuGat1 exhibited higher apparent Km values for maltose-pNP in the range of 16-43 mm (Fig. 4D), which were greater than the 4 mm value reported for rabbit glycogenin (
Extending the acceptor to three sugars or decreasing it to one resulted in less activity, but either anomer of Glc-pNP was acceptable (Fig. 4F). A similar pattern was observed for the P. ultimum and T. gondii enzymes. The enzymes were specific for terminal Glc acceptors as activity toward Gal-pNP was not detected. In comparison, pNP-Skp1 tetrasaccharide (GlFGaGn-pNP), which mimics the natural acceptor on Skp1, was a superior acceptor substrate (Fig. 4F) with a Km of 1.5 mm (Fig. S9G), and the truncated trisaccharide FGaGn-pNP was inactive, indicating that Glc was the position of attachment. Thus, the Gat1 enzymes preferred their native tetrasaccharide acceptor substrate and UDP-Gal as a donor but tolerated, with low efficiency, the preferred substrates of glycogenin, UDP-Glc, and α4-linked oligomers Glc.
The importance of Skp1 as context for the acceptor glycan was examined using GlFGaGn-Skp1, which was prepared by reaction of FGaGn-Skp1 with UDP-Glc and Glt1, resulting in loss of the trisaccharide epitope (Fig. 4H, inset). In a comparison of acceptor concentration dependence, TgGat1 was about 33 times more active toward GlFGaGn-Skp1 than free GlFGaGn-pNP (Fig. 4, H and I), based on ∼33 times less activity (dpm incorporated) at 0.001 times the substrate concentration (substrate concentrations were both in the linear response range). The reaction with GlFGaGn-Skp1 did not approach saturation at the highest concentration tested, 6 μm. Thus, the TgGat1 reaction was much more efficient when the tetrasaccharide was associated with its native substrate Skp1 and thus consistent with Gat1 being directly responsible for modifying Skp1 in the cell.
A characteristic of glycogenin is its ability to modify the HO group of a Tyr side chain near its active site with αGlc and then to repeatedly modify the 4-position of the Glc with another αGlc and repeat the process up to 8-12 sugars. Thus, when isolated as a recombinant protein expressed in UDP-Glc–positive E. coli, glycogenin is partly glucosylated (
). However, neither TgGat1 nor PuGat1 prepared in a similar manner were found to be glycosylated, based on an exact mass measurement using nLC/MS (Fig. S10, A–E). Furthermore, following incubation with either UDP-Gal or UDP-Glc, no change in SDS-PAGE mobility (Fig. S10D) or exact mass were observed (Fig. S10E). Thus, no evidence for autoglycosylation activity of Gat1 from either species could be detected.
Skp1 is the only detectable substrate of Gat1 in parasite extracts
GlFGaGn-Skp1 is a substrate for Gat1, but are there others? This was addressed by complementing extracts of gat1Δ parasites with recombinant TgGat1 in the presence of UDP-[3H]Gal and measuring incorporation of 3H after display of the proteome on a one-dimensional SDS-polyacrylamide gel. A high level of incorporation of 3H that depended on the addition of enzyme was observed at the position of Skp1 (Fig. 4J), as expected, but negligible dpm were detected elsewhere in the gel. Furthermore, negligible dpm were incorporated into Skp1 in RH parental cells, indicating that little GlFGaGn-Skp1 accumulates in WT cells. Similar results were observed in studies of gat1Δ clones in the RHΔΔ and Me49 backgrounds (Fig. S9, H and I). Finally, complementation of the gat1Δ clone in RHΔΔ with Gat1 expressed under the tubulin promoter resulted in the absence of measurable incorporation into Skp1, confirming specificity for Gat1 expression per se.
Although Gat1 was unable to serve as its own GT acceptor in the manner of glycogenin, its ability to modify α4Glc oligomers in vitro, albeit with low efficiency, raised the possibility that it may affect amylopectin (starch) biosynthesis in cells by applying αGal residues to its nonreducing termini. Starch normally accumulates to substantial levels in bradyzoites, a slow-growing form of the parasite that is induced by stress. Because induction of bradyzoite differentiation in cell culture is more efficient in type 2 strains, Me49 and gat1Δ/Me49 cells were induced by pH upshift and examined for differentiation by labeling of bradyzoite cyst walls with FITC-DBA-lectin and for starch with the periodic acid/Schiff's base reagent. As shown in Fig. 5, ME49 bradyzoites accumulated substantial levels of starch relative to tachyzoites, and no difference in the pattern or level was ascertained in gat1Δ cells using this qualitative assessment. Thus, Gat1 does not appear to affect starch synthesis, which is consistent with the absence of terminal Gal residues in a sugar composition analysis of T. gondii starch (
Amylopectin biogenesis and characterization in the protozoan parasite Toxoplasma gondii, the intracellular development of which is restricted in the HepG2 cell line.
Figure 5Starch appears unaffected in gat1Δ parasites. To promote normal starch accumulation, rapidly proliferating tachyzoites (A) of the type II strain Me49 (RFP-expressing) and its gat1Δ derivative were induced to differentiate as slow-proliferating bradyzoite cysts (B) in human foreskin fibroblasts. Cultures were fixed and stained with periodic acid/Schiff's base to reveal starch as a purple adduct. Arrow in (A), parasitophorous vacuole containing dozens to hundreds of tachyzoites within a fibroblast. Arrow in (B), cyst containing dozens of slow-growing bradyzoites, as confirmed by labeling of the cyst wall with FITC-DBA lectin (not shown). Scale bars, 50 μm. Two independent trials yielded similar results.
) was modified with TgGat1 using UDP-[U-13C]Gal as the donor substrate. The pentasaccharide reaction product (∼30% conversion) was analyzed together with the tetrasaccharide starting material by NMR. The previous assignment of the chemical shifts of the tetrasaccharide (
Sequence determination of oligosaccharides and regular polysaccharides using NMR spectroscopy and a novel Web-based version of the computer program CASPER.
) and confirmed by analysis of the two-dimensional COSY, TOCSY, and HMBC spectra (Fig. 6). One-dimensional 1H NMR spectra reveal the presence of the mixture of tetra- and pentasaccharides, followed by a downfield shift in the 13C-Glc-H1 peaks upon linkage to the terminal 13C-Gal (Fig. 6A). The HMBC 1H–13C correlation spectrum shows the (Fig. 6B, center panel) connection from Gal-H1 to Glc-H3, consistent with the downfield peak in the 1H–13C HSQC spectrum (Fig. 6B, top panel), and the proton resonances in the HSQC-TOCSY spectrum (Fig. 6B, bottom panel), establishing the glycosidic linkage between the terminal αGal and underlying αGlc as 1→3. Finally, a 1H–1H COSY spectrum confirms the assignments of the underlying Glc-H1, H2, and H3 (Fig. 6C). Consistent results were obtained when the tetrasaccharide was modified in the presence of PuGat1 (data not shown). Taken together, our NMR analyses are most consistent with the glycan structure, Galα1,3Glcα1,3Fucα1,2Galβ1,3GlcNAcα-, indicating that TgGat1 is a retaining UDP-Gal:glucoside α1,3-galactosyltransferase. Not only does Gat1 transfer the α-anomer of a different sugar compared with glycogenin, it attaches it to a different position (4- not 3-) of the acceptor αGlc residue.
Fig. 6Gat1 assembles a Galα1,3Glc linkage on the Skp1 tetrasaccharide. Shown is NMR analysis of the TgSkp1 pentasaccharide. [1-13C]Glcα1,3Fucα1,2Gal-β1,3GlcNAcα1-pNP was partially (30%) modified by TgGat1 or PuGat1 in the presence of UDP-[U-13C]Gal. A, one-dimensional 600-MHz 1H NMR spectrum. Magnification shows the region of 13C-anomeric carbons in the mixture of modified and unmodified tetrasaccharide. A cartoon diagram of the TgSkp1 pentasaccharide attached to pNP (Galα1,3Glcα1,3Fucα1,2Gal-β1,3GlcNAcα1-pNP) is shown at the top using CFG glycan symbols (
). B, 1H–13C HSQC, HMBC, and HSQC-TOCSY spectra demonstrating anomeric carbon to ring proton correlations. The Gal-H1/C1 doublet peaks were too weak to observe and are indicated by boxes. C, 1H–1H COSY and 1H–13C HMBC spectra. Identical results were obtained using [1-13C]GlFGaGn-pNP modified by PuGat1 (not shown). See Fig. S11 for a listing of chemical shift values.
Comparison of crystal structures explains the catalytic differences between Gat1 and glycogenin
To further probe the relationship between Gat1 and glycogenin, we compared their structures by X-ray crystallography. Attempts to crystallize TgGat1 were unsuccessful, even after deletion of its unconserved insert (Fig. 3A). However, PuGat1, which lacks this insert, was co-crystallized in the presence of Mn2+ and UDP.
The crystal structure of PuGat1 in complex with Mn2+ and UDP was solved using single-wavelength anomalous dispersion phasing of a Pt2+ derivative, and the resolution was extended to 1.76 Å using a native data set (Table S3). The asymmetric unit contains a single chain of PuGat1 with unambiguous electron density for the nucleotide and Mn2+ ion (Fig. 7A). The overall structure of PuGat1 reveals a canonical GT-A fold (
) consisting of eight α-helices and eight β-strands. The N terminus (residues 1-8) and two loops (residues 80-96 and 242-244) are disordered and were not modeled. The structure is similar to glycogenin-1 from the rabbit O. cuniculus (Oc-glycogenin-1), which superimposes 213 corresponding Cα atoms with a root mean square deviation of 3.3 Å despite a sequence identity of only 34% (Fig. 7B). The application of crystallographic symmetry shows that PuGat1 forms the same dimer described (
) analysis, the PuGat1 dimer interface buries 1,090 Å2 with a favorable p value of 0.107, which suggests that the dimer contact is stable. Sedimentation velocity analysis of 3.5 μm PuGat1 reveals a c(s) distribution consisting of single species at 4.0 S, which corresponds to the predicted value of 4.2 S for a dimer (Fig. 7C). The slightly slower sedimentation indicates that the enzyme in solution is less compact than that observed in the crystal structure. PuGat1 was dimeric even at 0.3 μm (Fig. S12), suggesting that it forms a dimer with an affinity >2-fold higher than that of Oc-glycogenin-1 (which has a reported Kd of 0.85 μm) (
). Based on gel filtration and preliminary sedimentation velocity experiments (not shown), TgGat1 is also a dimer.
Figure 7Gat1 is structurally related to glycogenin.A, asymmetric unit of crystallized Gat1 from P. ultimum in complex with UDP and Mn2+. α-Helices are shown in red, β-strands in yellow, and loops in green, and the active site with bound ligand is boxed. B, PuGat1 forms a homodimer that closely superimposes on Oc-glycogenin-1 (PDB entry 1LL2). Cylinders, α-helices; arrows, β-sheets; red ellipse, 2-fold symmetry axis perpendicular to the page. C, sedimentation velocity data modeled as a continuous c(s) distribution (normalized to 1.0) yields an S value for 3.5 μm PuGat1 that is close to the predicted value for a stable dimer in solution. Fig. S12 shows that the dimer is stable down to at least 0.3 μm. D, UDP is coordinated in nearly identical fashion to that of glycogenin-1, based on the difference density map (Fo − Fc) that was contoured at 5σ, calculated after omitting UDP and Mn2+ and subjecting the structure to simulated annealing. Octahedral coordination of Mn2+ is satisfied by the DXD motif, His-231, and UDP. The comparison with glycogenin-1 is illustrated in Fig. S13.
) and a conserved His residue coordinate the Mn2+ ion using the Oδ2 atom of Asp-117, both the Oδ1 and Oδ2 atoms of Asp-119, and the Nε2 atom of His-231 (Fig. 7D). The Mn2+ ion is also coordinated by the oxygen atoms from the α- and β-phosphates of UDP. Comparing the PuGat1 and Oc-glycogenin-1 active sites shows that all of the interactions with the nucleotide are conserved, with the exception of the interactions with N3 and O4 of the uracil ring (Fig. S13). Other changes in Gat1 active site include a Leu to Ser substitution at residue 233, which can potentially remove a packing interaction with the donor sugar in PuGat1 (Fig. 8). The hydroxyl of the substituted Ser-233 forms a hydrogen bond with the adjacent side chain of Gln-206. A water molecule (W509) replaces the Leu side chain and forms a hydrogen bond with Asn-149.
Figure 8Active-site geometry explains Gat1's preference for UDP-Gal rather than UDP-Glc. Comparison of the sugar-binding pockets of PuGat1 and Oc-glycogenin-1 displayed as wall-eyed stereo view. A, the Glc moiety is modeled based on the Oc-glycogenin-1 crystal structure with its intact sugar nucleotide. B, the Gal moiety is modeled by flipping the stereochemistry of Glc at the C4′ position. PuGat1 and glycogenin-1 side chains are represented by green and gray sticks, respectively, the yellow dashes represent hydrogen bonds, and water molecules are represented by blue spheres.
Comparing the crystal structure of UDP-Glc bound Oc-glycogenin-1 with PuGat1 immediately suggests a reason for why these enzymes have different donor specificities (Fig. 8). Superimposing the Oc-glycogenin-1:UDP-Glc structure onto PuGat1 shows that UDP-Glc would displace water 499 coordinated by Thr-180 (Fig. 8A). This would leave the Thr-180 hydroxyl group unsatisfied, and the unfavorable burial of a polar group likely explains why UDP-Glc is a poor donor (Fig. 8A). In contrast, we modeled in UDP-Gal by flipping the stereochemistry at the C4 position. The O4 atom of Gal would be ideally positioned to satisfy the Thr-180 hydroxyl group. Concurrently, Asp-176, whose Cα atom underwent a 2.3-Å shift relative to its location in glycogenin, would be in position to receive a hydrogen bond from the O4 atom of the Gal (Fig. 8B). Altogether, Gat1's sugar donor preference for UDP-Gal is likely due to formation of favorable hydrogen bonds with the O4 atom of Gal in contrast to the burial of the Thr-180 hydroxyl group when binding UDP-Glc.
Computational modeling predicts the specificity of Gat1 toward the Skp1 tetrasaccharide
To address the basis of Gat1's preference for the GlFGaGn- glycan, the lowest-energy conformation of the reducing form of the glycan generated by GLYCAM-web (RRID:SCR_018260) was docked using AutoDock Vina. A plausible docking mode was selected based on the requirement that the C′3-hydroxyl group must be oriented toward the anomeric carbon of the donor sugar to serve as the nucleophile for the addition of the Gal and that the glycan does not clash with the other subunit of the dimer. Of 100 docking simulations, only the top-scoring pose with a binding energy score of −5.7 kcal/mol satisfied the selection requirement. In this pose, the glycan adopted an alignment in a groove formed by Gat1 dimerization (Fig. 9, A and B). The glycan is stabilized by hydrogen bond contributions from the side chains or peptide backbones of residues Asp-141, Phe-143, and Ser-234 from subunit A and residues Leu-212, Lys-216, Asn-217, and Tyr-220 from subunit B (Fig. 9C). In addition, nonpolar interactions against the faces of the sugar moieties are provided by residues Thr-208, Leu-212, Phe-143, and Trp-221 from subunit A and residue Tyr-220 from subunit B (Fig. 9D). Phe-143, Tyr-220, Trp-221, and Ser-234 are uniquely conserved in Gat1 proteins relative to glycogenins (Fig. S6). The packing interactions from conserved hydrophobic residues are likely the major contributors in terms of binding energy. The extensive electrostatic and packing complementarity, which could not be achieved using the same approach using the α4Glc-tetramer recognized by glycogenin, can explain the distinct preference of Gat1 for the Skp1 tetrasaccharide acceptor substrate.
Figure 9Computational docking explains specificity of Gat1 for the Skp1-tetrasaccharide.A, top docking pose of the Skp1-tetrasaccharide in the PuGat1 active site and a groove formed by the dimer. B, 90° turn of the image shown in A. C, hydrogen-bonding interactions of the glycan. D, hydrophobic packing of the sugar faces and fucose-methyl (in sticks) with nonpolar surfaces (sticks/dots) of Gat1 subunit A (green) and B (cyan).
The Toxoplasma glycan influences Skp1 helix-8 extension via sugar-protein contacts
The glycan-protein contacts described previously for D. discoideum Skp1 correlated with the extension of helix-8, which was interpreted to provide better access by FBPs. To address whether this mechanism is conserved for T. gondii, despite the difference of the fourth sugar, the computational studies were repeated on glycosylated TgSkp1. Energy-minimized structures of the glycans from DdSkp1 and TgSkp1 revealed only a difference in the position of O4 of the fourth sugar (due to the Glc/Gal configurational inversion) (Fig. S14A). Six all-atom MD simulations were performed without coordinate constraints for 250 ns each to allow greater sampling of conformational space. Three of the simulations began with a 50-ns pre-equilibration of the glycan with Cα atom constraints on the polypeptide, whereas the others proceeded directly. As before, the simulations did not converge on a common structure, so a combined time-resolved linear regression analysis of the six simulations was conducted to identify correlations between helix extension and calculated polar and nonpolar interaction energies with each moiety of the glycan, and the results are shown in Table 2.
Table 2MMGBSA-derived per-residue energies between the protein and glycan that exhibit a strong correlation with helix extension (distance in Fig. 10A) according to a linear regression analysis of 48 bins from the six MD simulations (Fig. S14E)
A representative frame from a simulation with a strong correlation between the extension of helix-8 (dotted green line) and glycan contacts with the polypeptide chain is shown in Fig. 10 (A–C). Glycan-protein contacts involve both polar hydrogen bonds and nonpolar van der Waals interactions between sugar and amino acid residues. The strongest correlated polar interactions (Table 2) correspond to three hydrogen bonds between sugars and amino acids within the unstructured peptide chain located between helix-7 and -8 (residues 147-152), as depicted in Figs. 9C and 10B. Notably, the 4-OH of αGlc (the fourth sugar), which is epimeric to the D. discoideum fourth sugar, is oriented toward the solvent, thus not affecting protein contacts. Nonpolar interactions also occur within this region, with the terminal αGal moiety showing the highest correlation (Table 2). The hydrophobic face of the terminal αGal can pack against a nonpolar pocket consisting of planar faces of peptide backbone (Figs. 9E and 10D). Specifically, its C3–C5 hydrogens can pack with the backbone of residues 147-149, and a hydrogen on C6 can be buried within residues 143 and 144 (Fig. 10D). Both the polar and nonpolar interactions identified here involve amino acids that are conserved across both organisms, with the exception of Val-149, which is Lys in DdSkp1 (Fig. S14D). This substitution would not be expected to affect Gal burial because the majority of interactions are with the protein backbone at this position (not shown). The previously noted hydrogen bonds 1 and 3 to helix-7, involving Glu-140 (Glu-129 in DdSkp1), were also observed in this study (Fig. S14C) but were poorly correlated (r2 = 0.0 versus 0.63 average for the other three) with helix-8 extension (Fig. S14E). The energetic analyses applied to the T. gondii variant confirm its ability to form an organized structure that can still similarly influence local Skp1 polypeptide organization. The new studies emphasize the interaction of the glycan with the interhelix region rather than helix-7 and contributions of packing of the outer sugars, including the terminal αGal. As previously hypothesized, these effects on the conformational ensemble of TgSkp1 have the potential to make it more receptive to FBP binding.
Figure 10Packing of the glycan with Skp1 can explain F-box–binding site conformation.T. gondii GaGlFGaGn-Skp1 was subjected to six 250-ns all-atom molecular dynamics simulations. A, frame representative of the glycan-protein interaction and associated helix-8 extension, from a simulation (Equil-1; see Fig. S14E) in which the glycan was pre-equilibrated for 50 ns prior to the start of the simulation. The dotted green line refers to the distance from the C terminus to the center of mass of residues 1-136 and ranged from 18 to 61 Å. B, zoom-in of A depicting the glycan (carbon atoms in green) and amino acids (carbon atoms in gray) described in Table 2. Dotted black lines depict hydrogen bonds contributing to the polar energies described in Table 2. C, the back side of B. D and E, packing of terminal sugars against the polypeptide. Carbon atoms of the peptide are in orange. D, glycan and peptide represented as sticks. E, as in D, with glycan and peptide represented by spheres.
Glycosylation mediates much of the effect of O2-dependent prolyl hydroxylation of Skp1 in D. discoideum and T. gondii, and the present work emphasizes the importance of the fifth and final sugar of the glycan. The T. gondii glycan adopts a similar constrained conformation as described for D. discoideum despite the sequence difference (Fig. S14, A and B), and new all-atom molecular dynamics simulations showed a conserved cis-interaction with the back side of the intrinsically disordered region of Skp1 that comprises subsite-2 of the FBP-binding pocket (Fig. S14C). A deeper time-resolved energetics analysis of polar and nonpolar contacts showed a strong dynamic correlation between burial of the peripheral (especially the fifth) sugars in the turn from helix-7 to the loop connecting helix-8, and the extension of helix-8 which exposes subsite-2 for FBP docking (Fig. 10 and Table 2). This analysis de-emphasized the previous correlation of hydrogen bonds with helix-7 in favor of contacts with loop residues between helix-7 and helix-8 (Fig. S14E). This quantitative analysis of the trajectories extends critical features of the prior study in D. discoideum (
) by emphasizing the role of nonpolar packing of the glycan terminus in helix-8 extension.
The fifth sugar is attached by Gat1, an unusual GT that resides in the cytoplasmic compartment rather than the secretory pathway of the cell and appears to be dedicated to only Skp1. Although confined to the protist kingdom, it is nevertheless descended from a widely distributed lineage of sugar nucleotide–dependent CAZy GT8 GTs and may be the evolutionary predecessor of glycogenin, a GT that modulates glycogen formation in the cytoplasm of yeast, fungi, and animals. Although Gat1s from T. gondii and P. ultimum catalyze the transfer of αGal from UDP-Gal to nonreducing terminal Glc acceptors (Fig. 4 and Fig. S9), their ability to utilize UDP-Glc at low efficiency could be rationalized by the organization of the active site in the crystal structure (Fig. 8) and anticipates this evolutionary transition. The ability of UDP-Glc to inhibit Gat1's galactosyltransferase activity (Fig. 4G) might be regulatory in cells. The recently reported capability of glycogenin to transfer Gal at certain steps (
) lends further support to the evolutionary relationship. Although Gat1 is able to modify Glc in multiple contexts, it prefers Glc at the terminus of the native Skp1 tetrasaccharide over Glc in a native starch- or glycogen-like glycan (Fig. 4 and Fig. S9). Furthermore, Gat1 was substantially more reactive when the glycan was attached to Skp1 (Fig. 4, H and I), indicating that the apoprotein contributes to increased activity. The increase in activity may be due to the recognition of Skp1 protein by Gat1, but this hypothesis remains to be tested.
Computational docking showed how the Skp1 glycan alone in its calculated lowest-energy state can naturally fit within a groove along the dimer interface of PuGat1 (Fig. 9A), in a manner that cannot be achieved with an α4-glucan as in starch or glycogen in its lowest-energy state (not shown). The orientation of the acceptor Glc of the Skp1 glycan in this docking mode also supports the formation of the α1,3 linkage that was determined by NMR analysis of the product of the reaction of both PuGat1 and TgGat1 with a synthetic version of the Skp1 tetrasaccharide as an acceptor (Fig. 6). Furthermore, a high degree of conservation between the two enzymes is supported by their very similar activity and specificity characteristics (Fig. 4 and Fig. S9). The described binding mode is a possible explanation for the greater activity of the free tetrasaccharide as an acceptor substrate compared with the α4-linked glucans. Overall, the data suggest that recognition of Skp1 as a substrate includes both active-site preference for the specific glycan and separate determinants on the polypeptide. Owing to the limited computational scope of this study, however, further analysis is needed to fully define the basis for Skp1 recognition.
The significance of the terminal sugar is echoed by the slow growth phenotype of gat1Δ parasites in the monolayer plaque assay. This was attributed to the absence of Gat1 because similar results were obtained in independent knockouts in three different genetic backgrounds, and the defect was corrected by genetic complementation under its own promoter or a strong tubulin promoter in the uprt locus (Fig. 2). Because no other Gat1 targets were detected by biochemical complementation of gat1Δ extracts (Fig. 4J and Fig. S9 (H and I)), and no effect on starch formation was observed (Fig. 5), it is likely that failure to fully glycosylate Skp1 is responsible for the growth defect. This interpretation is consistent with similar effects of knocking out earlier GTs in the pathway and their selectivity for Skp1 in vitro (
). Although we cannot completely exclude some other function of Gat1, we note that at 345 amino acids, nearly all of its sequence appears devoted to the enzymatic domain (Fig. 3A). Studies are under way to identify SCF substrates and how they are affected by Skp1 modification and influence parasite growth and fitness.
The importance of the fifth sugar is reinforced by the observation that the amoebozoa, which lack both Glt1 and Gat1 (Fig. 1, A and B), evolved a new unrelated enzyme, AgtA, to fulfill the role. AgtA is a dual-function GT that applies both the fourth and fifth sugars, each an α3-linked Gal. This raises the interesting possibility that AgtA evolved to compensate for the unexplained loss of Glt1 and Gat1, with recursive addition of the same sugar being an accessible evolutionary pathway to recover the pentasaccharide. Gat1 sequences are found in alveolates (including T. gondii), stramenopiles (including P. ultimum), and rhizaria (which together with telonemids are known as the TSAR group) and archaeplastids (which together with the TSAR group are known as bikonts), but not in the unikonts that include the amoebozoa, fungi, and animals (
). Whereas the fate of Glt1 in the unicots is unknown, its disappearance potentially freed Gat1 to evolve into glycogenin, an α4-glucosyltransferase that, in vitro, can autoglucosylate itself to prime glycogen synthesis in fungi and metazoa. The enzymes share catalytic properties (see above) and perfectly conserve structural similarities, including the same dimer interface (Figure 7, Figure 8, Figure 9), although further studies are needed to support the model of lineal descent. Thus, it is interesting to speculate whether the proposed transition conserved a cellular function. Although glycogenin acquired a C-terminal glycogen synthase–binding domain (Fig. 3A) (
) that seems consistent with its in vitro capacity to initiate glycogen assembly, data now exist indicating that glycogenin is not required for glycogen synthesis in yeast and mice (
), although glycogen levels are affected by an unknown mechanism. The trans-glycosyltransferase activity for this enzyme lineage in bikonts protists raises the possibility of an unanticipated activity in unicots that could help explain the glycogenin-KO findings in yeast and mice.
Experimental procedures
Maintenance of host cells and parasite manipulations
Cultures of human foreskin fibroblast (HFF; ATCC SCRC-1041) or hTERT HFF (BJ-5ta, ATCC CRL-4001) were maintained Dulbecco's modified Eagle's medium supplemented with 10% (v/v) fetal bovine serum, 2 mm l-glutamine, and 100 units/ml penicillin/streptomycin (Corning) at 37 °C in a humidified CO2 (5%) incubator. Type 1 RH (
) strains of T. gondii were cultured on HFF or hTERT HFF monolayers in the same medium as described for host cells except that 1% (v/v) fetal bovine serum was used where stated and cloned by limiting dilution in 96-well plates. Parasites were kept in media without drug for plaque assays, which were performed as before (
). Transfections were conducted by electroporation using a Bio-Rad Gene Pulser Xcell at 1.5 kV and 25 microfarads with 2-mm electroporation cuvettes (VWR International) in Cytomix buffer (10 mm KH2PO4/K2HPO4 (pH 7.6), 120 mm KCl, 0.15 mm CaCl2, 5 mm MgCl2, 25 mm HEPES, 2 mm EDTA).
Disruption and complementation of gat1
TGGT1_310400 (ToxoDB, RRID:SCR_013453), referred to as gat1 (Fig. S1), was disrupted by two independent approaches. In the first method, a disruption DNA was prepared from the vector pmini.GFP.ht (a gift from Dr. Gustavo Arrizabalaga) in which the hxgprt gene is flanked by multiple cloning sites as described (
). 5′- and 3′-flanking targeting sequences were PCR-amplified from strain RHΔΔ with primer pairs Fa and Ra and pairs Fb and Rb, respectively (Table S1). The 5′-fragment was digested with ApaI and XhoI and cloned into similarly digested pmini.GFP.ht. The resulting plasmid was digested with XbaI and NotI and ligated to the similarly digested 3′-flanking DNA. The resulting vector was linearized with PacI and electroporated into strain RHΔΔ, selected under 25 μg/ml mycophenolic acid and 25 μg/ml xanthine, and cloned by limiting dilution. Genomic DNA was screened by PCR to identify Tggat1 disruption clones (Fig. S3A), using primers listed in Table S1.
The second approach was based on a double-CRISPR/Cas9 method as detailed previously (
), with minor modifications. To generate the dual guide (DG) plasmid, a fragment of p2 containing the guide RNA gat1-63 expression cassette was PCR-amplified using primers plasmid 3 FOR and plasmid 3 REV (Table S1), digested with NsiI, and ligated into the NsiI site of a dephosphorylated p3 containing guide RNA gat1-968. The type 1 RH and type 2 ME49 strains were co-transfected with pDG-Gat1 (10 μg) and a dihydroxyfolate reductase (dhfr) amplicon (1 μg) by electroporation (Fig. S4A). CRISPR/Cas9-mediated disruption in RHΔΔ parasites was done similarly except that pDG-Gat1 was co-transfected with a dhfr amplicon containing 45-bp homology arms targeting gat1 (Fig. S3B). Gat1Δ parasites were subsequently selected in 1 μm pyrimethamine (Sigma). The expected replacement of nucleotides 63-968 (relative to A of the ATG start codon) with the dhfr cassette was confirmed by PCR using primers listed in Table S1.
To complement gat1Δ RH parasites, a Tggat1 DNA fragment consisting of the coding sequence of gat1 plus ∼1.2 kb each of 5′- and 3′-flanking DNA sequences was generated by PCR from RH genomic DNA using primers Fa and Rb, which contained ApaI and NotI restriction sites, respectively. After treatment with ApaI and NotI, the PCR product was cloned into similarly digested pmini.GFP.ht plasmid in place of its hxgprt cassette, to generate pmini.Tggat1. The plasmid was transformed into E. coli Top10 cells and purified by using a Monarch miniprep kit (New England Biolabs). A Ty tag DNA sequence was inserted at the 3′-end of the gat1 coding sequence using a Q5 site-directed mutagenesis kit (New England Biolabs) and Fn and Rn primers, yielding pmini.Tggat1.Ty. The sequence was confirmed using primers Fl and Rl (Table S1). RH gat1Δ clone 8 was complemented by co-electroporation with a PCR amplicon from pmini.Tggat1.Ty (1 μg) and a sgUPRT CRISPR/Cas9 plasmid (10 μg) targeting the uprt locus using the guide DNA sequence 5′-ggcgtctcgattgtgagagc (
) (Fig. S4B). Transformants were selected with 10 μm fluorodeoxyuridine (Sigma), and drug-resistant clones were screened by PCR with primers listed in Table S1.
To complement gat1 in RHΔΔ parasites, the UPRT Vha1 cDNA shuttle vector containing TgVhaI cDNA (
) was modified to generate a Gat1-HA complementation plasmid using the New England Biolabs HiFi Builder method. The vector backbone, containing 5′- and 3′-flanking uprt targeting sequences and a Tg-tubulin promoter and 3×HA sequence, was PCR-amplified from the shuttle vector using primers Ft and Rt. The coding sequence of TgGat1 was PCR-amplified from pmini.Tggat1 plasmid using primers Fu and Ru, which had 18-21 nucleotides complementary to the terminal ends of the vector (Fig. S3C). The gel-purified PCR fragments were incubated with HiFi DNA assembly enzyme mix (New England Biolabs) and transformed into E. coli Top10 cells, yielding pUPRTgat13×HA. The Tggat1 sequence was confirmed using primers Fl and Rl. Complementation in a gat1Δ clone derived from RHΔΔ was done similarly except that a gat1-3×HA PCR amplicon with 5′ and 3′ uprt homology arms was used in the transfection.
Bradyzoite induction
ME49-RFP tachyzoites were differentiated to bradyzoites using alkaline pH (
). HFF monolayers were preincubated with sodium bicarbonate–free RPMI (Corning) containing 50 mm HEPES-NaOH (pH 8.1) for 24 h and then infected with tachyzoites and maintained at ambient atmosphere at 37 °C with medium replacement every 24 h. Differentiation was monitored by labeling with Dolichos biflorus agglutinin (DBA) (
). Infected HFF monolayers formed on 25-mm coverslips were washed with PBS (Corning), fixed with 4% paraformaldehyde in PBS for 10 min, washed with PBS, and permeabilized with 1% Triton X-100 (Bio-Rad) in PBS for 10 min, all at room temperature. Samples were blocked with 3% (w/v) BSA in PBS for 1 h, incubated for 2 h with 5 μg/ml FITC-DBA lectin (Vector Laboratories Inc.) in 1% BSA in PBS, and washed with PBS. The coverslips were mounted with ProLong Gold antifade reagent (Invitrogen) on glass slides. The slides were imaged by phase-contrast and fluorescence microscopy on a Zeiss Axioskop 2 Mot Plus.
Periodic acid staining
To assess amylopectin levels, parasite-infected HFF monolayers on 25-mm glass coverslips were washed with plain PBS, fixed with ice-cold MeOH for 5 min, washed with PBS, and incubated with 1% periodic acid in deionized H2O in the dark for 10 min (
). The coverslips were washed with deionized H2O, incubated with Schiff reagent for 15 min, washed once with deionized H2O, and finally rinsed with running tap water for 10 min. The stained coverslips were dehydrated by sequential immersion in 70% (v/v), 80, 90, and 100% EtOH, mounted on glass slides using Permount mounting medium (Fisher Scientific), and imaged on an EVOS XL Core microscope (Invitrogen).
Expression and purification of recombinant TgGat1 and PuGat1
The single exon-coding sequence of Tggat1 cDNA was amplified by PCR from RH genomic DNA using primers, Gat1 Fw and Gat1 Rv (Table S1), cloned into PCR4-TOPO TA (Invitrogen), and transformed into E. coli Top 10 cells. The plasmid was double-digested with BamHI and NheI to yield the gat1 fragment that was cloned into similarly digested pET15-TEVi plasmid (Invitrogen), resulting in the original 346-amino acid coding sequence of Gat1 extended at its N terminus with a His6 tag and TEV protease cleavage site (MGSSHHHHHHSSGRENLYFQGH-). A similar N-terminal modification of rabbit glycogenin did not significantly alter its enzymatic activity (
). The predicted coding sequence for Pugat1 was inferred from PYU1_G002535-201 (EnsemblProtists, RRID:SCR_013154, P. ultimum genome). The coding sequence was codon-optimized for E. coli expression (Fig. S2), chemically synthesized by Norclone Biotech (London, Ontario, Canada), and inserted in the NdeI and XhoI sites of pET15b-TEV. The expressed protein was extended at its N terminus with MGSSHHHHHHSSGENLYFQGH-.
TgGat1 or PuGat1 were expressed in and purified from E. coli BL21-Gold cells as described previously for TgGlt1 (
), through the purification on a 5-ml Co2+ TALON resin column. The eluted protein was dialyzed in 50 mm Tris-HCl (pH 8.0), 300 mm NaCl, 1 mm EDTA, and 2 mm β-mercaptoethanol, followed by 25 mm Tris-HCl (pH 8.0), 200 mm NaCl, 1 mm β-mercaptoethanol, and 5 mm MnCl2. The sample was treated with 2 μm His6-TEV protease, 5 μm TCEP in the same buffer overnight at 20 °C and reapplied to another Co2+ TALON column. The flow-through fraction was dialyzed in 25 mm Tris-HCl (pH 8.0), 50 mm NaCl, 1 mm β-mercaptoethanol, 2 mm MnCl2. The sample was concentrated by centrifugal ultrafiltration, and aliquots were stored at −80 °C. As indicated, the preparations were further purified by gel filtration on a Superdex200 column in the same buffer.
SDS-PAGE and Western blotting
Samples were suspended in diluted with Laemmli sample buffer and typically electrophoresed on a 4–12% gradient SDS-polyacrylamide gel (NuPAGE Novex, Invitrogen). Gels were either stained with Coomassie Blue or transferred to a nitrocellulose membrane using an iBlot system (Invitrogen). Blots were typically blocked in 5% nonfat dry milk in TBS and probed with a 1:1,000-fold dilution of the antibody of interest in the milk solution, followed by secondary probing with a 1:10,000-fold dilution of Alexa 680–labeled goat anti-rabbit IgG secondary antibody (Invitrogen). Blots were imaged on a LI-COR Odyssey IR scanner and analyzed in Adobe Photoshop with no contrast enhancement. For measuring incorporation of radioactivity, 1-mm thick 7–20% acrylamide gels were prepared manually, as detailed (
). Briefly, frozen pelleted tachyzoites (1 × 108) were resuspended in 8 m urea in 50 mm HEPES-NaOH (pH 7.4), incubated on ice for 1 h and at 50 °C for 5 min, and diluted 8-fold in IP buffer (0.2% Nonidet P-40 in 50 mm HEPES-NaOH, pH 7.4). The lysates were centrifuged at 21,000 × g for 20 min at 4 °C, and 100 μl of the supernatant (5 × 107 cells) were incubated with 5 μl of UOK75 rabbit antibody (affinity-purified anti-TgSkp1) coupled to protein A/G magnetic agarose beads (Pierce, 78609) for 1 h at 4 °C. Beads were captured in a DynaMag-2 magnet (Life Technologies) according to the manufacturer's directions and washed three times with 50 mm HEPES-NaOH (pH 7.4); three times with 10 mm Tris-HCl (pH 7.4), 50 mm NaCl; and once with 50 mm NaCl. Bound Skp1 was eluted twice with 60 μl of 133 mm triethylamine (sequencing grade; Pierce, 25108) by incubating for 15 min at RT and immediately neutralized with 40 μl of 0.2 m acetic acid. The eluted material was pooled and dried under vacuum, reconstituted in 8 m urea, 10 mm Tris-HCl (pH 7.4), reduced in 10 mm DTT for 40 min at RT, and alkylated in 50 mm chloroacetamide for 30 min at RT. Samples were then diluted to 2 m urea with 10 mm Tris-HCl (pH 7.4) and digested in 10 μg/ml Trypsin Gold (mass spectrometry grade; Promega, V5280) overnight at RT. Excess trypsin was quenched by the addition of 1% (v/v) TFA (Pierce, 28904) on ice for 15 min and centrifuged at 1,800 × g for 15 min at 4 °C. The supernatants were adsorbed to C18 pipette tips (Bond Elut OMIX C18, Agilent, A7003100) and eluted in 50% acetonitrile (ACN, OptimaTM LC–MS grade; Fisher, A955-4), 0.1% formic acid (FA, LC–MS grade; Pierce, 28905). Eluted peptides were vacuum-dried and reconstituted in 40 μl of 5% ACN, 0.05% TFA, and 4 μl were analyzed by nano-LC-MS/MS.
Treatment of TgSkp1 peptides with α-gal
Trypsinates from above were centrifuged at 1,800 × g for 15 min at 4 °C. The supernatants were dried under vacuum, resuspended in 100 mm sodium citrate phosphate buffer (pH 6.0), and treated with 3.6 milliunits of green coffee bean α-gal (Calbiochem) for 18 h at 37 °C. An additional 3.6 milliunits of α-gal was added for 8 h. After treatment, peptides were processed as above.
MS of TgSkp1 peptides
Reconstituted peptides were loaded onto an Acclaim PepMap C18 trap column (300 μm, 100 Å) in 2% (v/v) ACN, 0.05% (v/v) TFA at 5 μl/min, eluted onto and from an Acclaim PepMap RSLC C18 column (75 μm × 150 mm, 2 μm, 100 Å) with a linear gradient consisting of 4–90% solvent B (solvent A: 0.1% FA; solvent B: 90% ACN, 0.08% (v/v) FA) over 180 min, at a flow rate of 300 nl/min with an Ultimate 3000 RSLCnano UHPLC system, into the ion source of an Orbitrap QE+ mass spectrometer (Thermo Fisher Scientific). The spray voltage was set to 1.9 kV, and the heated capillary temperature was set to 280 °C. Full MS scans were acquired from m/z 350 to 2,000 at 70,000 resolution, and MS2 scans following higher energy collision-induced dissociation (
) were collected for the top 10 most intense ions, with a 30-s dynamic exclusion. The acquired raw spectra were analyzed using Sequest HT (Proteome Discoverer 2.2, Thermo Fisher Scientific) with a full MS peptide tolerance of 10 ppm and MS2 peptide fragment tolerance of 0.02 Da and filtered to generate a 1% target decoy peptide-spectrum match false discovery rate for protein assignments. All known glycoforms for TgSkp1-specific glycopeptides were manually searched for and verified. Searches were performed against the T. gondii (strain ATCC 50853/GT1) proteome (Uniprot proteome ID UP000005641, downloaded May 18, 2018; 8,450 entries) plus a list of common keratin, immunoglobulin-derived, and trypsin contaminants. Carbamidomethylation of Cys was set as a fixed/static modification, and oxidation of Met, deamidation of Asn and Gln residues, and acetylation of protein N termini were set as variable/dynamic modifications. Searches were performed with trypsin cleavage specificity, allowing two missed cleavage events and a minimum peptide length of 6 residues. All known glycoforms for TgSkp1-specific glycopeptides were manually searched for and verified. The raw files were uploaded to the Figshare server with ID 10.6084/m9.figshare.12272882 (Skp1 glycopeptides raw data, Fig. S5).
Enzyme assays
Sugar nucleotide hydrolysis
Recombinant TgGat1 (0.625-2.5 μm) was incubated with a given UDP-sugar (50 μm) (Promega) in 50 mm HEPES-NaOH (pH 7.4), 2 mm MnCl2, 5 mm DTT, in a final volume of 20 μl, for 1-16 h at 37 °C. The UDP generated was quantitated using the UDP-Glo assay (Promega) as described (
Glycosyltransferase activity toward small glycosides
In the standard reaction, TgGat1 or PuGat1 was incubated with 2 mm synthetic glycosides (pNP-α-galactoside, pNP-β-galactoside, pNP-α-glucoside, pNP-β-glucoside, pNP-α-maltoside, chloro-4-nitrophenyl-α-maltotrioside, FGaGn-pNP, GlFGaGn-pNP), 8 μm UDP-Gal (unlabeled), 0.17 μm UDP-[3H]Gal (15.6 μCi/nmol; American Radiolabeled Chemicals), 50 mm HEPES-NaOH (pH 7.0), 2 mm MnCl2, 5 mm DTT, in a final volume of 30 μl, for 1 h at 37 °C. Pilot studies indicated a pH optimum of 7.0, with 50% activity at pH 8.0 and 75% activity at pH 6.0 for TgGat1. Salt dependence studies showed maximal activity with no added NaCl or KCl and 35% activity at a 800 mm concentration of either salt. Activity showed a ∼6-fold preference for MnCl2 over MgCl2, with activity maximal at 2 mm MnCl2. The enzyme was essentially inactive in NiCl2, CoCl2, and CaCl2 (Fig. S9, A–D). For kinetic studies, concentrations and times were varied as indicated, and kinetic parameters were analyzed according to the Michaelis–Menten model based on the least squares fitting method in GraphPad Prism software. Reactions were stopped by the addition of 1 ml of 1 mm ice-cold Na-EDTA (pH 8.0), and incorporation of radioactivity into pNP-glycosides was analyzed by capture and release from a Sep-Pak C18 cartridge and scintillation counting (
) was incubated with 1.3 nmol (88 μg) of Tg-His6-Glt1, 4 nmol of UDP-Glc, 1.2 units of alkaline phosphatase (Promega), 50 mm HEPES-NaOH (pH 8.0), 5 mm DTT, 2 mm MnCl2, 2 mm MgCl2, 2 mg/ml BSA in a final volume of 121 μl for 3.5 h at 37 °C. The reaction was initiated by the addition of UDP-Glc and terminated by freezing at −80°C. Reaction progress was monitored by Western blotting with pAb UOK104, which is specific for FGGn-Skp1 from either D. discoideum or T. gondii (
), followed by probing with pAb UOK75, which is pan-specific for all TgSkp1 isoforms, for normalization. Approximately 85% of total Skp1 was modified. GlFGaGn-Skp1 was purified from Glt1 on a mini-QAE column using a Pharmacia Biotech SMART System. ∼19 μg of GlFGaGn-Skp1 was applied to a mini-QAE column pre-equilibrated with 50 mm Tris-HCl (pH 7.8), 5 mm MgCl2, 0.1 mm EDTA (buffer A) and eluted with a gradient from 0% A to 100% buffer B (50 mm Tris-HCl (pH 7.8), 5 mm MgCl2, 0.1 mm EDTA, 300 mm NaCl) in 40 min at a flow rate of 240 μl/min. GlFGaGn-Skp1 fractions were identified based on A280 and Western blotting probing with pAb UOK75.
TgGat1 (0.127 μm) was incubated with GlFGaGn-Skp1 or FGaGn-Skp1 (0.9-3.6 μm), 3.2 μm UDP-[3H]Gal (15.6 μCi/nmol), 40 μm UDP-Gal (unlabeled), 0.2% (v/v) Tween 20, 50 mm HEPES-NaOH (pH 7.0), 2 mm MnCl2, 5 mm DTT, in a final volume of 20 μl, for 1 h at 37 °C. The reaction was stopped by the addition of 4× Laemmli sample buffer, 1 m DTT (50 mm final concentration), and 2 μg of soybean trypsin inhibitor and boiled for 3 min SDS-PAGE, and incorporation of radioactivity into Skp1 was performed as described above.
To detect Gat1 activity in cells, a cytosolic extract was prepared by hypotonic lysis, ultracentrifugation at 100,000 × g for 1 h, and desalting as described previously (
). 25 μg of desalted S100 protein was incubated with 10-50 nmol of Tg-GlFGaGn-Skp1 and 1.0 μCi of UDP-[3H]Gal (15.6 μCi/nmol) for 5 h, and incorporation of radioactivity into protein was assayed by SDS-PAGE and scintillation counting as described above.
To search for Gat1 substrates, cytosolic S100 fractions (180 μg of protein) were incubated with 0.13 μm TgGat1, 2.0 μCi of UDP-[3H]Gal (15.6 μCi/nmol) in a final volume of 60 μl containing 50 mm HEPES-NaOH (pH 7.0), 2 mm MnCl2, 5 mm DTT, 30 mm NaF, 0.2% Tween 20, at 37 °C for 1 h, supplemented with 1.7 μg of Tg-GlFGaGN-Skp1 as indicated. Incorporation of radioactivity was monitored by the SDS-PAGE assay as described above.
MS of Skp1
To evaluate its glycosylation status, recombinant PuGat1 or TgGat1 (purified by gel filtration) was incubated with or without UDP-Gal or UDP-Glc in the absence of added acceptor substrate and diluted to 50 ng/μl Skp1 with 2% acetonitrile, 0.05% (v/v) TFA. 250-500 ng of protein (5-10 μl) was injected into an Acclaim PepMap C4 trap cartridge (300 μm × 5 mm) equilibrated with 0.05% TFA, 2% acetonitrile, ramped up with an increasing gradient to 0.1% formic acid, 25% acetonitrile, and introduced into an Acclaim PepMap analytical C4 column (75 μm × 15 cm, 5-μm pore size) maintained at 35 °C in an Ultimate 3000 RSLC system coupled to a QE+ Orbitrap mass spectrometer (Thermo Scientific). After equilibrating the analytical column in 98% LC–MS Buffer A (water, 0.1% formic acid) for 10 min and a 6-min ramp up to 27% LC–MS Buffer B (90% (v/v) acetonitrile, 0.1% formic acid), separation was achieved using a linear gradient from 27 to 98% Buffer B over 20 min at a flow rate of 300 nl/min. The column was regenerated after each run by maintaining it at 98% Buffer B for 5 min. The effluent was introduced into the mass spectrometer by nanospray ionization in positive ion mode via a stainless-steel emitter with spray voltage set to 1.9 k, capillary temperature set at 250 °C, and probe heater temperature set at 350 °C. The MS method consisted of collecting full ITMS (MS1) scans (400-2,000 m/z) at 140,000 resolution in intact protein mode (default gas P set to 0.2). PuGat1 species eluting between 17.5 and 21.5 min and TgGat1 species eluting between 18.5 and 22.5 min (∼60–80% acetonitrile) were processed with Xcalibur Xtract deconvolution software to generate monoisotopic masses from the multicharged, protonated ion series. Because TgGat1 MS spectra were not isotopically resolved, masses were extracted after MS spectra deconvolution using the ReSpect algorithm in the BioPharma Finder suite (Thermo Scientific), with a 20 ppm deconvolution mass tolerance and 25 ppm protein sequence matching mass tolerance. For consistency, PuGat1 was also deconvoluted and re-extracted using the ReSpect algorithm and the same conditions. The data raw files were uploaded to the Figshare server with ID 10.6084/m9.figshare.12272909 (recombinant Gat1 intact protein raw files, Fig. S10).
Structure determination of PuGat1
A PuGat1:UDP:Mn2+ complex in 50 mm HEPES-NaOH, pH 7.4, 75 mm NaCl, 2 mm DTT, 5 mm UDP, and 5 mm MnCl2 was crystallized at 20 °C using a hanging-drop vapor diffusion method over a reservoir containing 8–12% (w/v) PEG4000, 0.4 m ammonium sulfate, and 0.1 m sodium acetate at pH 4.0. Crystals were obtained overnight and were transferred to a reservoir solution containing 15% (v/v) of a cryoprotectant mixture (1:1:1, ethylene glycol/DMSO/glycerol) and flash-cooled with liquid N2. The complex crystallized in space group P42212 and diffracted to 1.76 Å (Table S3). X-ray data were collected remotely at the SER-CAT 22-BM beamline at the Argonne National Laboratory using a Fast Rayonix 300HS detector and processed using XDS (
), with 5% of the data omitted for cross-validation.
PuGat1:UDP:Mn2+ crystals were soaked with platinum cyanide for heavy-atom phasing, cryoprotected, and frozen as above. PuGat1:UDP:Pt2+ crystals were isomorphous to PuGat1:UDP:Mn2+ crystals and diffracted to 2.1 Å (Table S3). PuGat1:UDP:Mn2+ crystals were alternatively soaked in UO2 (not shown). The crystals, which diffracted to 2.4 Å in the same space group, lacked UDP density, suggesting displacement by UO2.
The crystal structure of PuGat1:Pt2+ was solved using single-wavelength anomalous dispersion. The data were obtained at a wavelength of 1.85 Å for maximum anomalous signal. A single Pt2+ site was located using PHENIX (
), and the resulting phases had an acceptable figure of merit of 0.31. The model was subjected to iterative cycles of refinement and yielded a final model with Rwork/Rfree of 0.21/0.24 (Table S3). The structure of PuGat1:UDP:Mn2+ was solved using rigid-body refinement of PuGat1:Pt2+. The resulting model was subjected to iterative cycles of refinement and yielded a final model with Rwork/Rfree of 0.18/0.21 (Table S3). Images were rendered in PyMOL (
). A grid box with dimensions 26 × 26 × 34 Å was placed over the ligand-binding site based on where the acceptor is bound in the glycogenin/glucan complex. The ligand was kept rigid, because AutoDock Vina is not parameterized specific to glycan torsion angles. 100 binding modes were calculated with the lowest binding energy scored at −4.4 kcal/mol and the highest binding energy scored at −5.7 kcal/mol.
Sedimentation velocity studies
PuGat1 was further purified on a Superdex S200 gel filtration column (GE Healthcare) equilibrated with 20 mm potassium phosphate (pH 7.4), 50 mm KCl, 0.5 mm tris(2-carboxyethyl)phosphine. Protein concentration was calculated from A280 measured in an Agilent 8453 UV-visible spectrophotometer, based on a molar absorptivity (ε280) of 60,390 m−1 cm−1, which was calculated from the PuGat1 sequence using ProtParam (
). Samples were diluted to 0.3-11 μm, loaded into 12-mm double-sector Epon centerpieces equipped with quartz windows, and equilibrated for 2 h at 20 °C in an An60 Ti rotor. Sedimentation velocity data were collected using an Optima XLA analytical ultracentrifuge (Beckman Coulter) at a rotor speed of 50,000 rpm at 20 °C. Data were recorded at 280 nm for protein samples at 3.5–11 μm and at 230/220 nm for samples at 0.3–1.5 μm, in radial step sizes of 0.003 cm. SEDNTERP (
Computer-aided interpretation of analytical sedimentation data for proteins.
in: Harding S.E. Rowe A.J. Horton J. Analytical Ultracentrifugation in Biochemistry and Polymer Science. Royal Society of Chemistry,
Cambridge, UK1992: 90-125
) was used to model the partial specific volume of PuGat1 (0.73818 ml/g) and the density (1.0034 g/ml) and viscosity (0.0100757 poise) of the buffer. Using SEDFIT (
), data were modeled as continuous c(s) distributions and were fit using baseline, meniscus, frictional coefficient, and systematic time-invariant and radial-invariant noise. Predicted sedimentation coefficient (s) values for the PuGat1 monomer (2.8 S) and dimer (4.2 S) were calculated using HYDROPRO (
). Molecular dynamics simulations were performed as described previously. Briefly, MD simulations were performed with the pmemd.cuda version of AMBER14 (
) in a truncated octahedral box with 15-Å distance from the solute to the end of the unit cell. Electrostatic interactions were treated with the particle mesh–Ewald algorithm, and a cut-off for nonbonded interactions was set to 8 Å (
). SHAKE was employed to constrain hydrogen-containing bonds, enabling an integration time step of 2 fs. Restraints were imposed in specific situations and were enforced with a 10-kcal/mol Å2 energy barrier in each case. Each minimization step consisted of 1,000 cycles of the steepest descent method (1,000 cycles), followed by 24,000 cycles using the conjugate gradient approach. The systems were heated to 300 K under NVT conditions over 60 ps, employing the Berendsen thermostat with a coupling time constant of 1 ps. The subsequent simulations were performed under NPT conditions. A torsion term that corrects 4(trans)-hydroxyproline residue (Hyp) ring puckering was included in simulations of the O-linked residue type based on previous studies that indicate that the ring is primarily exo when glycosylated (
). This torsion term has been adopted in GLYCAM06 (version K).
A 50-ns simulation of the protein was performed with Cα cartesian constraints on all amino acids except those generated by Chimera. The fully glycosylated isoform was created by adding the TgSkp1 pentasaccharide to the exo-pucker conformation of hydroxyproline (residue 154). Six independent simulations were performed. Three ran for 250 ns directly, whereas the other three began with an additional 50 ns in which the protein was restrained to allow the glycan time to adapt to the protein conformation.
Computational analysis
Structural images were created with Visual Molecular Dynamics (
). The structure depicted in Fig. 7 was created by identifying the frame from equil-1 that consisted of the lowest root mean square deviation to the average structure as calculated by cpptraj (
). The cpptraj program was also used to distribute the latter 200 ns of the six simulations into 48 bins containing 250 frames each for analysis. Per-residue MMGBSA energies were calculated with MMPBSA.py.MPI with igb = 2 and idecomp = 3. A Bash script was used to calculate the correlation coefficients (
Proteins related to TgGat1 were searched for using a BLASTP (version 2.4.0) search seeded with the full-length TgGat1 protein sequence against the NCBI nonredundant database (December 2016). The evolutionary relationship of Gat1-like sequences was investigated by using a maximum likelihood method (
). Catalytic domains from 43 CAZy GT8 sequences were selected based on their relatedness to Gat1, glycogenin, or known function and consisted of 196 positions. Sequence alignments were manually curated in BioEdit (version 7.2.5). An initial tree(s) for the heuristic search was obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using a JTT model and then selecting the topology with superior log likelihood value. A discrete γ distribution was used to model evolutionary rate differences among sites (5 categories (+G, parameter = 1.1608)). The rate variation model allowed for some sites to be evolutionarily invariable ([+I], 1.02% of sites).
Data availability
All atomic coordinates and structure factors were deposited in the RCSB PDB under accession codes 6MW5 (PuGat1 Pt2+ derivative) and 6MW8 (PuGat1 native). The raw data files for MS were uploaded to the Figshare server with ID 10.6084/m9.figshare.12272882 (Skp1 glycopeptides raw data, Fig. S5) and ID 10.6084/m9.figshare.12272909 (recombinant Gat1 intact protein raw files, Fig. S10).
Acknowledgments
We thank Kentuan Hicks, Nathan Beattie, and Safal Shrestha for assistance.
Developmentally regulated biosynthesis of carbohydrate and storage polysaccharide during differentiation and tissue cyst formation in Toxoplasma gondii.
Amylopectin biogenesis and characterization in the protozoan parasite Toxoplasma gondii, the intracellular development of which is restricted in the HepG2 cell line.
Sequence determination of oligosaccharides and regular polysaccharides using NMR spectroscopy and a novel Web-based version of the computer program CASPER.
Computer-aided interpretation of analytical sedimentation data for proteins.
in: Harding S.E. Rowe A.J. Horton J. Analytical Ultracentrifugation in Biochemistry and Polymer Science. Royal Society of Chemistry,
Cambridge, UK1992: 90-125
Author contributions—M. M., H. W. K., D. T., K. R., H. v. d. W., R. J. W., Z. A. W., and C. M. W. conceptualization; M. M., H. W. K., D. T., M. O. S., E. G.-P., K. R., and C. M. W. formal analysis; M. M. validation; M. M., H. W. K., D. T., M. O. S., E. G.-P., K. R., P. Z., N. G. D., H. v. d. W., H. T. I., and J. N. G. investigation; M. M., H. W. K., D. T., M. O. S., E. G.-P., K. R., P. Z., N. G. D., H. v. d. W., H. T. I., J. N. G., L. W., R. J. W., Z. A. W., and C. M. W. methodology; M. M., H. W. K., D. T., M. O. S., E. G.-P., K. R., and C. M. W. writing-original draft; M. M., H. W. K., D. T., R. J. W., Z. A. W., and C. M. W. writing-review and editing; H. W. K., M. O. S., E. G.-P., P. Z., and C. M. W. data curation; D. T., E. G.-P., and C. M. W. visualization; L. W., R. J. W., Z. A. W., and C. M. W. supervision; C. M. W. resources; C. M. W. funding acquisition; C. M. W. project administration.
Funding and additional information—M. M. was supported in part by National Institutes of Health (NIH) Grant T32-AI060546; H. W. K. was supported in part by NIH Grant T32-GM107004; and K. H. was supported by National Science Foundation REU Grant DBI-1426834. This project was supported in part by NIH Grant RO1-GM084383 (to C. M. W. and Ira Blader), Grant #14-140 from the Mizutani Foundation for Glycoscience (to C. M. W. and Ira Blader), NIH Grant P41-GM103490 (to L. W., senior investigator), NIH Grant 8P41-GM103390 (Resource for Integrated Glycotechnology to J. H. Prestegard), NIH Grant R01-GM114298 (to Z. A. W.), and NIH Grant S10-OD021762 (to John Rose). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Conflict of interest—The authors declare that they have no conflicts of interest with the contents of this article.
Present addresses: David Thieker, Dept. of Biochemistry and Biophysics, University of North Carolina School of Medicine, Chapel Hill, North Carolina, USA.
Present addresses: M. Osman Sheikh, Amicus Therapeutics, Philadelphia, Pennsylvania, USA.
Present addresses: Kazi Rahman, HIV Dynamics and Replication Program, Center for Cancer Research, National Cancer Institute, Frederick, Maryland, USA.
Present addresses: H. Travis Ichikawa, New Materials Institute, University of Georgia, Athens, Georgia, USA.
Abbreviations—The abbreviations used are:
SCF
Skp1/Cullin-1/F-box protein subcomplex of E3 Cullin-RING-1 ubiquitin ligases
50341-5&id=gr1.jpg){kind=link}
50341-5&id=gr2.jpg){kind=link}
50341-5&id=gr3.jpg){kind=link}
50341-5&id=gr4.jpg){kind=link}
50341-5&id=gr5.jpg){kind=link}
50341-5&id=gr6.jpg){kind=link}
50341-5&id=gr7.jpg){kind=link}
50341-5&id=gr8.jpg){kind=link}
50341-5&id=gr9.jpg){kind=link}
50341-5&id=gr10.jpg){kind=link}