Abstract
The polypeptide GalNAc-transferases (GalNAc-Ts), that initiate mucin-type O-glycosylation, consist of a catalytic and a lectin domain connected by a flexible linker. In addition to recognizing polypeptide sequence, the GalNAc-Ts exhibit unique long-range N- and/or C-terminal prior glycosylation (GalNAc-O-Ser/Thr) preferences modulated by the lectin domain. Here we report studies on GalNAc-T4 that reveal the origins of its unique N-terminal long-range glycopeptide specificity, which is the opposite of GalNAc-T2. The GalNAc-T4 structure bound to a monoglycopeptide shows that the GalNAc-binding site of its lectin domain is rotated relative to the homologous GalNAc-T2 structure, explaining their different long-range preferences. Kinetics and molecular dynamics simulations on several GalNAc-T2 flexible linker constructs show altered remote prior glycosylation preferences, confirming that the flexible linker dictates the rotation of the lectin domain, thus modulating the GalNAc-Ts' long-range preferences. This work for the first time provides the structural basis for the different remote prior glycosylation preferences of the GalNAc-Ts.
GalNAc transferases’ (GalNAc-Ts) catalytic domains are connected to a lectin domain through a flexible linker. Here the authors present a structural analysis of GalNAc-T4 that implicates the linker region as modulator of the orientations of the lectin domain, which in turn imparts substrate specificity.
Introduction
Polypeptide N-acetylgalactosaminyltransferases (GalNAc-Ts) are higher eukaryotic-retaining glycosyltransferases (GTs) that transfer a GalNAc moiety from uridine diphosphate N-acetylgalactosamine (UDP-GalNAc) onto Ser/Thr residues of proteins1. This family of enzymes, which initiates the most common type of O-glycosylation in metazoans1–5, plays profound roles for the structure, stability, and function of proteins6. GalNAc-O-glycosylation is one of the most abundant types of protein glycosylation in eukaryotes that is chiefly found in the densely clustered heavily glycosylated domains of mucins and mucin-domain-containing glycoproteins, hence the name, mucin-type O-glycosylation. However, it is also clear that many other proteins contain isolated sites of GalNAc-O-glycosylation, where more than 5000 human proteins trafficking the secretory pathway have been identified to date containing one or more mucin-type O-glycans7.
The initiation of mucin-type O-glycosylation by the GalNAc-Ts is one of the most complex regulated types of protein glycosylation. The GalNAc-Ts comprise a family of up to 20 isoenzymes possessing a range of kinetic properties and expression patterns. These enzymes orchestrate with high fidelity, the initial patterns of O-glycosylation on diverse proteins, including the high-density regions in mucins, where 30–50% of the amino acids may be glycosylated1, 8. The GalNAc-Ts are unique among metazoan GT enzymes because in addition to their N-terminal catalytic domain adopting a GT-A fold (family 27 in the CAZy database9), they possess a C-terminal GalNAc-binding lectin domain with a β-trefoil fold4, 10–13 (classified as a carbohydrate-binding module 13 in the CAZy database), which provides additional functions to these enzymes14 (see below for further details). Both domains are linked through a short flexible linker whose motion has been suggested to be responsible for the dynamic conformational landscape of these enzymes4.
GalNAc-Ts can be classified according to their different glycosylation capacity on peptide and glycopeptide substrates. Three distinctive glycosylation modes have been reported: catalytic domain-dependent glycosylation on “naked” peptides (Fig. 1a) and on glycopeptides (Fig. 1b), hereafter termed short-range glycosylation15, 16; and remote lectin domain-dependent glycosylation on glycopeptides (Fig. 1c), termed long-range glycosylation15, 16. Adjacent and proximal sites relative to an existing GalNAc glycosite are glycosylated in a catalytic domain-dependent manner (Fig. 1), with GalNAc-T4 being part of a small number of GalNAc-Ts (including GalNAc-T7 and T10) that will also glycosylate contiguous sites16. Distant sites can be glycosylated in a lectin domain-dependent manner (Fig. 1) with individual GalNAc-Ts having distinct preferences on remote sites located N or C terminus from a prior GalNAc glycosite15, 16. For example, GalNAc-T1/T2/T14 tend to glycosylate preferentially N-terminal sites remote from prior C-terminal GalNAc glycosites in glycopeptides, whereas GalNAc-T3/T4/T6/T12 exhibit the opposite preference (Fig. 1c). Furthermore, some GalNAc-Ts, including GalNAc-T5/T13/T16, appear to exhibit both orientational preferences for long-range glycosylation16. Previous structural work with GalNAc-T2 in complex with “naked” peptides and glycopeptides has provided insight into the catalytic domain-dependent reaction, and also demonstrated how the GalNAc-T2 lectin domain guides catalysis to N-terminal acceptor sites very distant from a prior C-terminal GalNAc glycosite4, 12, 17. However, the molecular basis of how other GalNAc-Ts, such as GalNAc-T3/T4/T6/T12, achieve the opposite long-range glycosylation preferences remains unclear15, 16.
The GalNAc-T4 isoform is of particular interest as it is the only isoform capable of glycosylating two out of the five acceptor sites in the partially glycosylated MUC1 mucin tandem repeat in a lectin domain-dependent manner18. The density of glycosylation together with the structure of O-glycans in the mucin repeat regions are important features for the cancer-glycoforms of MUC1 being exploited for diagnostic as well as therapeutic purposes19. Generally, GalNAc-T4 glycosylates very few isolated glycosylation sites in “naked” peptide acceptors, although it readily glycosylates the important Thr57 in PSGL-1 that is the essential O-glycan required for P-selectin-mediated leukocyte trafficking20, 21. We report herein a multidisciplinary approach on GalNAc-T4, combined with the characterization of different chimeras and mutants of GalNAc-T2, which begins to reveal the molecular basis of the long-range glycosylation preferences of the GalNAc-Ts on glycopeptide acceptors substrates.
Results
Models explaining the long-range glycosylation preferences
Based on previous results, we can infer two possible models (see Fig. 1c) that might provide a plausible explanation for the GalNAc-T3/T4/T6/T12 long-range glycosylation preference in contrast to the one found for enzymes such as GalNAc-T1/T2/T1415, 16. In model 1, the rotation of the lectin domain GalNAc-binding site is required, thereby maintaining the same orientation of the peptide on the catalytic domain. In model 2, the lectin domain GalNAc-binding site would adopt the same orientation as found in GalNAc-T2 while the glycopeptide adopts an inverse orientation on the catalytic domain, an assumption unlikely to happen providing the high similarity of the peptide-binding groove among GalNAc-Ts11 (Fig. 1).
GalNAc-T4-substrate interaction and kinetics studies
To test the two models proposed above, three simplified peptides enriched in glycine and alanine residues were designed, synthetized, and evaluated as substrates of GalNAc-T4 (Table 1). One of these peptides was a “naked” peptide, 1 (denoted –TT– for simplicity), whereas the other two were monoglycopeptides, 2 and 3 (denoted -TT–T*- and -T*–TT-, where T* represents a GalNAc glycosylated Thr), that contain a GalNAc moiety located either at the C or N terminus, respectively. All peptides have two potential Thr acceptor sites and also contain an adjacent PXP motif, which is recognized by most GalNAc-Ts4, 16. Our kinetic studies showed that GalNAc-T4 glycosylated selectively peptide 3 (with a high affinity, Km ~ 40 µM, and catalytic efficiency, ~75 (nmole min)−1) over peptides 1 and 2, both of which were imperceptibly glycosylated over a range of peptide concentrations (Fig. 2a). This is consistent with previous findings16 on GalNAc-T4 and illustrates that the lectin domain guides GalNAc-T4 to preferentially glycosylate C-terminal sites of monoglycopeptides that are remote from an N terminus GalNAc-glycosylated residue. Although, the “naked” peptide 1 was not a substrate, GalNAc-T4 has been reported to readily glycosylate other highly specific peptide sequences20, 21.
Table 1.
Peptide | Sequence |
---|---|
1 (–TT–) | GAGAGAGTTPGPG |
2 (-TT–T*-) | AGAGTTPGPGAGAT*GA |
3 (-T*–TT-) | GAT*GAGAGAGTTPGPG |
T* denotes the Thr-O-GalNAc moiety
To understand why GalNAc-T4-only glycosylated monoglycopeptide 3, we performed binding studies using saturation-transfer difference (STD) NMR experiments on peptide substrates 1–3 and α-methyl-GalNAc in the presence of UDP and MnCl2 (Fig. 2b and Supplementary Fig. 1). The absence of STD-NMR signals for peptide 1 revealed that the “naked” peptide was very poorly recognized by GalNAc-T4 in agreement with the low glycosylation observed during the kinetics assays. For both monoglycopeptides 2 and 3 relatively weak STD-NMR signals were observed for the Ala and Thr methyl protons, suggesting weak recognition of the peptide by GalNAc-T4. However, the O-GalNAc glycosylated Thr of the monoglycopeptides 2 and 3 receive a higher STD enhancement. Remarkably, the sugar protons of the GalNAc moiety of the monoglycopeptides 2 and 3 (region 3.5 a 4.2 ppm of the spectrum together with the singlet at 2 ppm) showed unambiguous and more intense STD-NMR signals (Fig. 2b and Supplementary Fig. 1). In addition, the STD-NMR-derived epitope mapping of the GalNAc moiety is identical for both monoglycopeptide 2 and control α-methyl-GalNAc. In fact, for both compounds, the highest STD response corresponded to the H2 proton of the GalNAc moiety followed by a modest STD of protons H4, H3, and N-acetyl group, suggesting similar lectin-binding modes for these molecules. In contrast, monoglycopeptide 3 presents a somewhat different GalNAc STD-NMR-binding pattern, where H2, H4, and the N-acetyl group of the GalNAc moiety have the highest STD-NMR response closely followed by the H3 proton. This suggests that the monoglycopeptide 3 has a distinct binding mode of interaction with the enzyme, which is different to the one found for monoglycopeptide 2 and α-methyl-GalNAc control. From our kinetic and STD-NMR results, we can infer that the unique binding of the GalNAc residue of monoglycopeptide 3 to the lectin domain, facilitates the correct orientation of the peptide acceptor onto the catalytic domain for optimal C-terminal glycosylation. Thus, even though monoglycopeptide 2 shows GalNAc binding to the lectin domain, its N-terminal acceptor site must be directed away from the catalytic domain, leading to the inability of GalNAc-T4 to glycosylate this substrate. Hence, the correct location of a prior glycosite in a monoglycopeptide, with respect to potential acceptor sites, is essential for an optimal interaction of the peptide within the catalytic domain, which in turn is critical for an efficient catalysis and specificity.
Crystal structure of GalNAc-T4
To further elucidate the molecular basis of GalNAc-T4’s monoglycopeptide glycosylation preference and rule out one of the two proposed models (Fig. 1), we obtained triclinic crystals of GalNAc-T4 that were subsequently soaked with monoglycopeptide 3 and UDP and MnCl2. The resulting crystal allowed us to solve the corresponding structure at a high resolution (1.90 Å) (Supplementary Table 1). Within the asymmetric unit, two independent GalNAc-T4 molecules were present. The crystal structure shows a compact structure with the typical GT-A fold and the lectin domain located at the N- and C-terminal regions, respectively (Fig. 3a). A sequence alignment between GalNAc-T2 and T4 displays a sequence identity of ~40% and a high resemblance between these two enzymes at the secondary structure level (Supplementary Fig. 2). In fact, superposition analysis of both enzymes shows that both catalytic domains are more structurally similar (root-mean-square deviation (RMSD) of 1.23 Å for aligned Cα atoms) than the lectin domains (RMSD of 1.82 Å) (Fig. 3b and Supplementary Fig. 3). A global superposition of both enzymes renders a larger RMSD of 2.12 Å, which suggests that a major shift has occurred between the orientation of the lectin domains (see below for a further discussion; Fig. 3b). As expected and accordingly to model 1, the GalNAc-binding site in GalNAc-T4 is located on the right side of the lectin domain, which requires a rotation of 28° with respect to the homologous site in GalNAc-T2 (Fig. 3b and Supplementary Fig. 4). To calculate the angle of rotation, the Cαs of Asp458, Trp282, and Asp459 were used. Note that the angle of rotation could be also determined using the Cαs of Asp458, Trp286, and Asp459. In this case, the rotation angle was 30°, which is very close to 28°. This is due to Trp282 of GalNAc-T2 and Trp286 of GalNAc-T4 adopt almost identical positions in the catalytic domain (Fig. 3b). This different position of the GalNAc-T4 GalNAc-binding site, consistent with the STD-NMR results, may explain why GalNAc-T4 promotes a distinct type of long-range glycosylation, consisting of the glycosylation of acceptor sites located at the C-terminal from prior GalNAc glycosites in glycopeptides. This also suggests that the location of the GalNAc-binding site of the lectin domain is coupled with the long-range glycosylation preferences of these enzymes.
Lectin domain-binding site of GalNAc-T4
Even though the crystals were soaked with a high concentration of peptide 3, UDP and MnCl2, we could only visualize well-defined density for the GalNAc-O-Thr moiety of monoglycopeptide 3 bound to the lectin α-subdomain GalNAc-binding site (Fig. 3a, c). Note that GalNAc-Ts lectin domains contain three potential carbohydrate-binding subdomains named α, β, and γ; however, in most GalNAc-Ts only one of them is functional15, 16, 18, 22, 23. The almost absence of density for the peptide backbone of monoglycopeptide 3, is consistent with the STD-NMR results (showing no significant STD with the acceptor peptide residues) and with the weak binding affinity of this peptide (Kd estimated in the mM range), determined by surface plasmon resonance (SPR) experiments in the presence of UDP and MnCl2 (Supplementary Fig. 5). The discrepancies between GalNAc-T4’s high-affinity Km obtained from our kinetics studies (Fig. 2a) versus the poor Kd values from the binding studies can be attributed to the fact that the enzyme kinetics was performed with UDP-GalNAc. It is known that this donor substrate stabilizes the flexible loop of the catalytic domain active site in a closed conformation, thus completing the formation of the peptide-binding groove and leading to an active enzyme4, 12. Thus, our observations in the absence of UDP-GalNAc represent the weak binding of peptide substrate to the open flexible loop conformation of the enzyme.
A closer inspection at the lectin domain GalNAc-binding site depicts that the GalNAc moiety was tethered by conserved residues, including Asp459, Asn483, Phe475, and His478, of the lectin α-subdomain GalNAc-binding site (equivalent residues in GalNAc-T2 are Asp458, Asn479, Tyr471, and His474; Fig. 3d). It is important to note that only the lectin α-subdomain GalNAc-binding site has been shown to be functional in both GalNAc-T2 and GalNAc-T4 lectin domains4, 18. Most of the interactions are through hydrogen bonds in both enzymes, although they also share CH–π interactions between the GalNAc moiety and Phe475/Tyr471. The importance of the Asp residue in recognition of the GalNAc moiety and its role in the lectin domain-mediated glycosylation has been previously reported4, 18. Furthermore, the crystal structure explained the GalNAc STD-NMR-derived epitope mapping. In particular, the structure shows that H4 and H3 of the GalNAc moiety are in closer contact with Phe475 and His478, which is further supported by the STD response (92% for H4 and 72% for H3) observed for these protons in monoglycopeptide 3.
Molecular docking and dynamics simulations
To address how the lectin domain guides the delivery of potential peptide acceptor sites to the catalytic domain, we performed molecular docking and dynamics (MD) simulations on GalNAc-T4 in complex with UDP/Mn+2 and monoglycopeptide 3 (see “Methods” section). The calculations suggest that although the peptide bound at the catalytic domain was highly dynamic, the two potential peptide acceptor Thr residues were in close contact to the β-phosphate of UDP during the simulated timeframe (200 ns), thereby accounting for the glycosylation of this peptide (Fig. 3e, f, g and Supplementary Fig. 6 and Supplementary Movie 1). A comparison with the GalNAc-T2 structure in complex with MUC5AC-13 and UDP/Mn+2 clearly provides evidence of the two different orientations of the lectin domains and supports our conclusion on the coupling of the location of the GalNAc-binding site with the observed long-distance glycosylation preference (Fig. 3e, f).
The flexible linker dictates the lectin domain rotation
Previous studies suggested that the flexible linker allowed for the interdomain translational-like motion in these enzymes4. We therefore rationalized that different flexible linkers found in these GTs might also be behind the rotation capacity of the lectin domain (Supplementary Fig. 7). To test this hypothesis, we initially performed 500 ns MD simulations in water solution of native GalNAc-T2 and GalNAc-T4. Simulations clearly showed that the lectin domains of GalNAc-T2 and GalNAc-T4 apo enzymes do not rotate in the considered timescale (Supplementary Movies 2 and 3). This suggests that the orientation of the GalNAc-binding site with respect to the catalytic domain may remain relatively fixed in a fully folded transferase. We then performed MD simulations on chimeras 1 and 2 that corresponded to GalNAc-T2 with two different lengths of the flexible linker from GalNAc-T3 (Table 2). In agreement with our hypothesis, we observed the prompt rotation of the lectin domain of the GalNAc-T2 chimera 2 toward a position similar to that in GalNAc-T4 (Fig. 4a, Supplementary Fig. 8, and Supplementary Movie 4). This stepwise motion was completed in ~30 ns and involved sequential extension of the flexible linker, rotation of the lectin domain around the Z and Y axes, and finally compression of the linker (Fig. 4a). Thus, initially the GalNAc-T3 linker was artificially forced to adopt the compact conformation of the GalNAc-T2 native linker, but early on in the simulation (after ~10 ns), the linker quickly recovered its more native extended conformation and as a consequence, the lectin domain springs and rotates sidewise to finally re-assemble in a structure similar to that of our GalNAc-T4. This motion is likely smoothed by the absence of strong interdomain interactions in GalNAc-T2, which do take in place in GalNAc-T4 (see for instance the persistent salt bridge between Arg397 and Glu487 in Supplementary Movie 5). Thus, the intrinsic conformational preference(s) of a short, unfolded fragment of ~13 residues, was found to determine the relative orientation of two very large protein domains (see also Supplementary Movies 6 and 7; these movies show the intrinsic different conformational preferences of the GalNAc-T2 (folded) and GalNAc-T3 (extended) flexible linkers by performing 500 ns MD simulations of these isolated fragments).
Table 2.
Transferase | Sequence of flexible linkera |
---|---|
GalNAc-T2 | P435ELRVPDHQDIAF447 |
GalNAc-T2-t3Flexible linker (Chimera 1) | P494EVYVPDLNPVIS506 |
GalNAc-T2-t3Flexible linker-AF (Chimera 2) | P494EVYVPDLNPVIS506AF |
GalNAc-T2-t4Flexible linker (Chimera 3) | P433NLHVPEDRPGWH445 |
GalNAc-T2-t3Flexible linker-AF-P503A | P494EVYVPDLNAVIS506AF |
(Chimera 2_P503A) | |
GalNAc-T2 (double mutant, R438A-D444A) | P435ELAVPDHQAIAF447 |
GalNAc-T2 (triple mutant, R438A-D444A-F447A) | P435ELAVPDHQAIAA447 |
GalNAc-T3 | P494EVYVPDLNPVIS506 |
GalNAc-T4 | P433NLHVPEDRPGWH445 |
GalNAc-T4-t2Flexible linker (Chimera 4) | P435ELRVPDHQDI445 |
GalNAc-T4-t2Flexible linker-AF (Chimera 5) | P435ELRVPDHQDIAF447 |
GalNAc-T4-t2Flexible linker-A (Chimera 6) | P435ELRVPDHQDIA446H |
Above are the sequences of the native transferases and linker domain mutants and chimeras utilized in this work. Only the sequences of the linkers that are exchanged are shown. Shortened names in (bold) are referred to in the text and figures
aNote that some of the exchanged linker regions also include part of the succeeding β-strand of GalNAc-T2 and GalNAc-T4. This is due to the boundaries of flexible linkers are not well defined at the C terminus. The additional Pro in the flexible linker of GalNAc-T3/T4, that is not present in GalNAc-T2, is underlined
Kinetics characterization of the chimeras
To validate our MD simulations predictions, we expressed and purified chimeras 1–3 and performed a series of partial and full kinetic analysis against our peptide substrates (Fig. 4b, c) (also note that chimeras4–6 could not be expressed in our hands, see “Methods” section and Table 2). Remarkably, the GalNAc-T2 chimeras rendered glycosylation preferences similar to GalNAc-T3 and GalNAc-T415, 16 (preferring monoglycopeptide 3 over monoglycopeptide 2, by approximately fourfold, which is significantly different from native GalNAc-T2) despite comprising most of the GalNAc-T2 architecture (Fig. 4b, c and Supplementary Fig. 9). Interestingly, GalNAc-T2’s preference for these substrates varied with peptide substrate concentration, ranging from a high preference for monoglycopeptide 2 at low substrate concentrations to nearly equal preferences for both monoglycopeptides 2 and 3 at higher substrate concentrations (Fig. 4c). These differences arise from monoglycopeptide 2’s ~fourfold lower Km compared to monoglycopeptide 3, thus giving an over threefold increase in catalytic efficiency for monoglycopeptide 2 over monoglycopeptide 3 (Fig. 4c and Supplementary Table 2). These findings readily explain the previously reported differences in GalNAc-T2’s long-range glycosylation preferences15, 17. This type of behavior was not observed in the kinetics of GalNAc-T4 or chimera 3 (note that detailed kinetics were not performed on chimera 1 or 2). Further analysis revealed that the changes in specificity for chimeras 1–3 were chiefly due to a decrease in the observed specific activity (chimera 1 and 2), and changes in Vmax and Km (chimera 3) against monoglycopeptide 2 compared to that of monoglycopeptide 3 whose specific activity and kinetic parameters remained relative constant (Fig. 4b, c, Supplementary Table 2, and Supplementary Fig. 9). It is also worth noting that the specific activities (and detailed kinetics when obtained) for chimeras 1–3 toward the “naked” peptide 1 are quite similar to those of GalNAc-T2, indicating that the GalNAc-T3/T4 flexible linkers in these chimeras did not affect the overall architecture of the GalNAc-T2’s active site (Fig. 4b, c, Supplementary Table 2, and Supplementary Fig. 9). Finally, we evaluated which Thr residues of the substrates were preferably glycosylated by Edman amino acid sequencing, finding that the Thr residue adjacent to the PXP motif was glycosylated by chimeras 1–3, as would be expected for the PXP motif16 (Supplementary Fig. 10). Little to no glycosylation was observed at the proximal Thr.
Together, our results suggest that the glycosylation preferences found for the GalNAc-T2 chimeras can be reasonably explained by the rotation of the lectin domain as predicted by the MD simulations of chimera 2. Thus, the lectin domains of these chimeras must begin to adopt an equivalent position to that of GalNAc-T4 to account for their altered long-range glycosylation preferences. We conclude that the flexible linker modulates the rotation of the lectin domain, which subsequently determines the orientation of the functional GalNAc-binding site of the lectin domain with respect to the catalytic domain. This modulation of the lectin domain relative to the catalytic domain likely determines the distinct long-range glycosylation preferences for all of the GalNAc-T isoforms. Our studies below further examine the role of the linker in dictating these long-range preferences.
Site-directed mutagenesis of the GalNAc-T2 and chimera 2
Next, we explored whether specific particular residues within the flexible linkers might be responsible for the lectin domain orientation. A multiple alignment of the flexible linkers (Supplementary Fig. 7) shows that GalNAc-T3/T4/T6/T12, which share the same long-distance glycosylation preferences15, 16, contain an additional Pro that is not present in GalNAc-Ts encompassing the opposed glycosylation preferences such as GalNAc-T1/T2/T14. This Pro in the flexible linkers of GalNAc-T3 (P494EVYVPDLNPVIS506; underlined) and GalNAc-T4 (P433NLHVPEDRPGWH445; underlined) is located at the end of the loop with respect to the flexible linker of GalNAc-T2 (P435ELRVPDHQDIAF447). Considering that Pro residues are abundant in protein turns, we hypothesized that this additional Pro might be, at least in part, responsible of the different orientation found crystallographically for the GalNAc-T4 lectin domain, and computationally for chimera 2. To explore this hypothesis, we mutated this Pro (in the GalNAc-T3 flexible linker sequence of the GalNAc-T2 chimera 2) to Ala-producing chimera 2_P503A with the expectation that it would reverse its preference (Table 2). This chimera, however, showed an identical glycosylation preference as the starting chimera 2 (Fig. 4b and Supplementary Fig 9), rather than reverting to the expected wild-type GalNA-T2 preference. This suggests that more complex events within the flexible linker might take place to explain the interplay between the flexible linker dynamics and its coupling with the lectin domain rotation. Also note that this mutated chimera displays a very low-specific activity compared to GalNAc-T2 (~3–10%) (Supplementary Fig. 9). This may suggest that the architecture of the catalytic site of this chimera may have been altered or perhaps blocked by an altered orientation of its lectin domain.
Next, we evaluated in detail the interactions between residues of the flexible linker in the compact and extended crystal structures of GalNAc-T24. Two major interactions, namely a salt bridge between Arg438 and Asp444, and a CH–π interaction between Pro440 and Phe447, appeared to be important to maintain the compact, folded structure of its flexible linker (Supplementary Fig. 11). We reasoned that these interactions might be important for fixing the GalNAc-T2 lectin domain orientation with respect to its catalytic domain. To test this hypothesis, we generated the double mutant R438A-D444A and the triple mutant R438A-D444A-F447A (Table 2). Our assumption was that the disruption of these interactions might lead to more “flexible” linkers, which in turn could cause some degree of rotation of the lectin domain toward the orientation observed in GalNAc-T4. MD simulations showed that, whereas the lectin domain did not rotate in the appropriate direction in the double mutant R438A-D444A (Supplementary Movie 8), it did undergo complete rotation in the triple mutant R438A-D444A-F447A toward the orientation observed in chimera 2 and GalNAc-T4, although such motion required a much longer simulation time (complete rotation observed after 160–450 ns, Supplementary Movie 9).
These results were validated by their observed substrate preferences (obtained for both mutants; Fig. 4b and Supplementary Fig. 9) and by a kinetic analysis of the triple mutant R438A-D444A-F447A (Fig. 4c and Supplementary Table 2). Whereas the double mutant gave nearly equal monoglycopeptide preferences as observed for GalNAc-T2 (at high substrate concentration), the triple mutant showed a 1.5–2-fold preference for monoglycopeptide 3 over monoglycopeptide 2, resembling the preference of GalNAc-T4 (Fig. 4b,c, Supplementary Table 2, and Supplementary Fig. 9). As with chimeras 1–3, the change in preferences observed with the GalNAc-T2 triple mutant was due to a decrease in Vmax and increase in Km for monoglycopeptide 2 rather than due to a significant change in the kinetic properties of monogycopeptide 3 (Fig. 4b,c, Supplementary Table 2, and Supplementary Fig. 9). The specific activity (and kinetic constants, when obtained) of these linker mutants against peptide 1 are very similar to GalNAc-T2, as well as chimeras 1–3, further supporting the notion that their catalytic domains have not been significantly altered (Supplementary Table 2, Fig. 4b, c, and Supplementary Fig. 9). In summary, the stepwise change in preference (Fig. 4b, c) from monoglycopeptide 2 toward that of monoglycopeptide 3, for native GalNAc-T2, the triple mutant and chimera 3, mostly reflects changes in the kinetic constants against monoglycopeptide 2 rather than those for monoglycopeptide 3.
All together, these results strongly suggest that the flexible linker plays a significant role in directing the long-range glycopeptide specificity of this family of transferases by altering the relative orientations of the catalytic and lectin domains. In addition, our results suggest that specific interactions within the flexible linker are sufficient to modulate this orientation thereby dictating each isoforms long-range glycosylation preferences.
Discussion
Most types of protein O-glycosylation in higher eukaryotes (O-Man, O-Xyl, O-Glc, O-GlcNAc, and O-Fuc) are initiated by only one or two GT isoenzymes1. However, these types of protein O-glycosylation are only found on limited proteins and sequence motifs. On the contrary, Nature has imposed a large family of GalNAc-Ts isoforms to account for mucin-type O-glycosylation. Initially, these GTs were considered redundant isoenzymes but recent findings clearly demonstrate that individual GalNAc-T isoforms serve unique biological functions in health and disease. For example, deficiency in GALNT3 and GALNT2 causes familial tumoral calcinosis24 and dyslipidemia25, respectively. How such transferase specificity may be regulated is best illustrated by GalNAc-T3, which is the only isoform capable of glycosylating Thr178 at a proprotein convertase (PC) processing site (RHTR179↓) of FGF23, an important regulator of phosphate homeostasis22. However, glycosylation of Thr178 by GalNAc-T3 requires prior O-glycosylation of residue Thr171 of FGF23, demonstrating that GalNAc-T3 utilizes its long-range N-terminal preference for prior GalNAc glycosylation to target or amplify the glycosylation of Thr17815. Thus, the combination of local peptide sequence and remote prior glycosylation is used by these transferases to ensure the fidelity of isoform-specific glycosylation. With the recent quantitative differential O-glycoproteomics studies7, it is now clear that many individual GalNAc-T isoforms serve unique non-redundant contributions to glycosylation of the proteome, presumably employing similar strategies as that of GalNAc-T3 and FGF23. The long-range lectin domain-mediated functions of GalNAc-Ts are likely instrumental for the coordinated glycosylation of additional diverse sequence motifs and in particularly high-density acceptor sequence motifs such as in mucins. Therefore, it is likely that the large number of GalNAc-Ts and their lectin domains evolved to recognize diverse peptide sequence motifs as well as to differentially recognize prior sites of GalNAc glycosylation. This would ensure the fidelity of glycosylation in both low- and high-density glycosylated regions by sequential orchestrated glycosylation reactions, coordinated by both the specificity of the catalytic domain as well as the long-range recognition of the lectin domain, that interacts with the prior attached GalNAc residue15, 16. Therefore, to fully understand the biological functions of these enzymes and to elucidate their roles in disease, we must understand the molecular basis, and mechanisms therein, which dictate and modulate their complex long- and short-range substrate specificity.
Here we have addressed the molecular mechanisms behind the unique long-range lectin domain-dependent glycosylation functions of GalNAc-Ts and have provided a model explaining how these transferases differentially recognize such N- or C-terminal long-range prior glycosylation15, 16. We previously reported that the location of the lectin α-subdomain GalNAc-binding site, relative to the catalytic domain, was key to explaining how GalNAc-T2 delivered N-terminal acceptor sites of C-terminal monoglycopeptides to the catalytic domain in a lectin domain-dependent manner4. We also demonstrated that the flexible linker, located ~28–36 Å from the active site of GalNAc-T2, and ~18–25 Å from GalNAc-T2’s lectin α-subdomain was in charge of the conformational heterogeneity of GalNAc-T24. This was clearly due to the translational capacity of the lectin domain with respect to the catalytic domain, which was further modulated in the presence of substrates4. Note also that in all of the GalNAc-T2 structures reported, irrespectively of whether they are compact or extended structures, the lectin domain is positioned such that its α-subdomain GalNAc-binding site is always located on the left side of the lectin domain4, 11, 12. However, it has remained a conundrum how other GalNAc-Ts (i.e., GalNAc-T3, T4, T12, etc.) presented the opposite long-range glycosylation preference of GalNAc-T2 and others. Our multidisciplinary approach presented herein has revealed additional conformational properties of the short flexible linker that provide the molecular basis for the distinct long-range glycosylation preferences found for these enzymes. Our crystal structure of GalNAc-T4 revealed that its lectin domain is positioned such that its lectin α-subdomain GalNAc-binding site is now located on the right side of the lectin domain, in agreement with our proposed Model 1 (Figs. 1c and 3). This is fundamental to explaining why GalNAc-T4 and similar behaving GalNAc-Ts prefer to glycosylate C-terminal acceptor sites of N-terminal glycosylated glycopeptides. Our STD-NMR studies together with the transferase kinetic studies suggested that the orientation of the lectin domain-bound GalNAc moiety was fundamentally important for the correct orientation and presentation of the acceptor residues to the catalytic site, thereby leading to more efficient catalysis. Our MD simulations on multiple chimera transferases together with their kinetic analysis, further provide compelling evidence that the flexible linkers not only allow translational motion to the lectin domain, but also rotational motion. Both types of motions likely differ between the GalNAc-T isoforms and therefore will contribute to the different glycosylation preferences found among these enzymes. All together, our studies point to the dissimilar flexible linkers as major contributors to the distinct long-range glycosylation preferences of these transferases (Supplementary Fig. 7). It is further likely that some flexible linkers will have conformational properties in solution that allow two or more conformational states of the lectin domain relative to the catalytic domain thus accounting for the dual long-range preferences observed for GalNAc-T2 (and its chimeras and mutants) as well as the several GalNAc-Ts that display nearly equal N- and C-terminal long-range glycosylation preferences (i.e., GalNAc-T5, 13,1615). That we did not observe such conformational flexibility in our MD trajectories of GalNAc-T2 (and some chimeras and mutants) may suggest that the activation energy for this interconversion may be very high and/or that our simulations were not sufficiently long enough.
In conclusion, we have provided for the first time the molecular basis for the distinctive long-range glycosylation preferences of the GalNAc-Ts, which is based on a very small flexible linker that provides rotational and translational capacity to the lectin domain. This work further exemplifies how a structural feature, very distant from both the active site and the lectin domain GalNAc-binding site, is capable of tuning the activity and specificity of these biologically important multidomain enzymes.
Methods
Cloning and purification of human GalNAc-T4/chimeras
The DNA sequence encoding amino acid residues 36–578 of the human GalNAc-T4, defined as galnact4, was codon optimized and synthesized by GenSript (USA) for expression in Pichia pastoris. The DNA, containing at the 5′ end a recognition sequence for XhoI and a KEX2 cleavage signal, and at the 3′ end a sequence encoding a histidine tag, a stop codon and a recognition sequence for SacII, was cloned into the pUC57 vector (GenScript). Following digestion with XhoI and SacII, the construct was subcloned into the protein expression vector pPICZαA, resulting in the expression plasmid pPICZαA-galnact4. Subsequently, the plasmid was linearized by SacI-HF and transformed into SMD1168. Transformants were selected and colonies were grown as described before26, 27. The DNA sequences encoding the GalNAc-T4 chimeras were also synthesized by GenScript and cloned in the same above vector. DNA linearization and transformation in SMD1168 were also performed in the same manner described above. Note that none of the GalNAc-T4 chimeras were expressed as a soluble secreted form in P. pastoris supernatants, impeding their biophysical characterization.
Supernatant containing the human GalNAc-T4 was dialyzed against buffer A (20 mM Na2HPO4 pH 7.4, 20 mM imidazol 500 mM NaCl) and loaded into a His-Trap Column (GE Healthcare). The protein was eluted with an imidazol gradient in buffer A from 20 mM to 500 mM. Buffer exchange of GalNAc-T4 into buffer B (25 mM Tris pH 8, 150 mM NaCl) was carried out using a HiPrep 26/10 desalting column (GE Healthcare). GalNAc-T4 was further purified by size exclusion chromatography using a Superdex 75 XK26/60 column (Sigma) previously equilibrated with buffer B. Fractions containing GalNAc-T4 were dialized against buffer C (25 mM Tris pH 8, 1 mM TCEP), concentrated and used for biophysical experiments.
Cloning and purification of GalNAc-T2 chimeras/mutants
The DNA sequences encoding the GalNAc-T2 chimeras were synthesized by GenScript and cloned in the vector pPICZαA vector as previously described12. Site-directed mutagenesis experiments were also performed by GenScript using the previous reported template pPICZαAgalnact2 (K75-Q571)12 and the plasmid encoding the chimera 2. The linearization and transformation of all constructs to SMD1168 were also performed in the same manner described for GalNAc-T4. The purification of all these GalNAc-T2 chimeras and mutants were expressed and purified using the purification protocol of the wild-type enzyme, as described previously12. Purity of enzymes was evaluated by SDS-PAGE coomassie staining and quantification of enzymes was quantified by absorbance at 280 nm using their theoretical extinction coefficients.
Crystallization
Crystals of the GalNAc-T4 were grown by hanging drop experiments at 18 °C by mixing 0.5 μl of protein solution (4 mg/ml GalNAc-T4, 5 mM UDP, 2 mM MnCl2, and 5 mM MUC5AC-3,134 (GTT*PSPVPTTSTT*SAP) in 25 mM Tris pH 7.5, 0.5 mM EDTA, and 1 mM tris(2-carboxyethyl)phosphine (TCEP)) with an equal volume of a reservoir solution (18% PEG3350, 0.1 M ammonium nitrate). Under these conditions, crystals appeared within 2–5 days. Note that MUC5AC-3-13 acted as an additive to improve the diffraction quality and the resolution of the crystals. The crystals were soaked for 30 min with a mix containing 20 mM monoglycopeptide 3 and 20 mM UDP in 25 mM Tris pH 7.5 and 2 mM MnCl2. Then, the crystals were cryo-protected with 25% ethylene glicol, 18% PEG3350, and 0.1 M ammonium nitrate, and frozen in a nitrogen gas stream cooled to 100 K.
Structure determination and refinement
The data were collected in the beamline I03 of Diamond (DLS) at a wavelength and temperature of 0.97 Å and a temperature of 100 K, respectively. The data were processed and scaled using the XDS package28 and CCP429, 30 software. Relevant statistics are given in Supplementary Table 1. The crystal structure was solved by molecular replacement with Phaser29, 30 and using the PDB entry 5AJP as the template that corresponds to the human GalNAc-T2. Initial phases were further improved by cycles of manual model building in Coot31 and refinement with REFMAC532. Once the GalNAc-T4 catalytic domain was unambiguously built and refined, ARP/wARP29, 30 was used to fully build the GalNAc-T4 lectin domain. Again, new rounds of Coot and refinement with REFMAC5 were performed. The crystal structure of GalNAc-T4 from crystals co-crystallized with MUC5AC-3-13 displayed no electron density for the MUC5AC-3-13 diglycopeptide. The final model of crystals soaked with monoglycopeptide 3 and UDP was validated with PROCHECK, model statistics are given in Supplementary Table 1. The asymmetric unit of the triclinic crystal contained two molecules of GalNAc-T4. Despite the long soaking and high concentration of monoglycopeptide 3 and UDP, only clear density was visualized for the presence of GalNAc-O-Thr and GalNAc moieties of monoglycopeptide 3 bound to the lectin α-subdomain GalNAc-binding site (Fig. 3a, c). The Ramachandran plot shows that 96.71, 2.89, and 0.40% of the amino acids are in most favored, allowed, and disallowed regions, respectively.
Surface plasmon resonance experiments
SPR experiments were performed at 25 °C with a Biacore X-100 apparatus (Biacore AB) in 25 mM Tris buffer, 1 mM DTT, 4 mM MnCl2, 100 µM UDP, 0.01% surfactant P20, pH 7.5 (running buffer) at 25 °C. Flow cells (CM5 sensor chip; Biacore) were activated for 7 min by injecting 140 µl of 50 mM N-hydroxysuccinimide (NHS): 200 mM ethyl-3(3dimethylamino) propylcarbodiimide (EDC). GalNAc-T4 was immobilized on a flow cell 2 by injection of a 100 µg/ml protein solution diluted with 10 mM sodium acetate buffer with a flow rate of 10 µl/min for 7 min followed by an injection of 130 µl ethanolamine to block any remaining activated groups on the surface. The level of immobilization reached was about 8000 RU. Flow cell 1, used as reference, was blocked with ethanolamine at the same conditions of flow cell 2 without immobilization of protein. Affinity experiments were made using a series of different concentrations of peptide 3 in the range of 0.01–5 mM with a flow rate of 30 µl/min during 100 s. Each injection was followed by a 100 s injection of running buffer (dissociation phase). No regeneration steps were performed between injections. Response data were collected at real time and analyzed with the Biacore® X-100 Evaluation software and plotted as response shift versus analyte concentration. We could not determine a Kd value due to binding saturation was not achieved.
Synthesis of peptides
Peptides were synthesized via stepwise microwave-assisted solid-phase peptide synthesis on a Liberty Blue synthesizer using the Fmoc strategy on Rink Amide MBHA resin (0.1 mmol). The glycosylated amino acid building blocks (2.0 equiv) were synthesized as described in the literature33 and manually coupled using HBTU, while the other Fmoc amino acids (5.0 equiv) were automatically coupled using oxyma pure/DIC. The O-acetyl groups of (AcO)3GalNAc moiety were removed in a mixture of NH2NH2/MeOH (7:3). The peptides were then released from the resin, and all acid sensitive side-chain protecting groups were simultaneously removed using TFA 95%, TIS 2.5%, H2O 2.5%, followed by precipitation with cold diethyl ether. Finally, they were purified by HPLC using a Phenomenex Luna C18(2) column (10 μ, 250 mm × 21.2 mm) and a dual absorbance detector, with a flow rate of 20 ml/min. Peptides were further subjected to Edman amino acid sequencing on a Shimadzu PPSQ-53A peptide sequencer prior to use.
Peptide 1: HPLC: Rt = 11.53 min (Grad: water 0.1% TFA/acetonitrile (95:5) → (87:13), 13 min, λ = 212 nm). HRMS ESI + (m/z) calcd. for C39H66N14O15 [M + 2 H]2+ 485.2411, found 485.2301.
Monoglycopeptide 2: HPLC: Rt = 11.75 min (Grad: water 0.1% TFA/acetonitrile (95:5) → (85.6:14.4), 15 min, λ = 212 nm). HRMS ESI + (m/z) calcd. for C57H95N18O24 [M + H]+ 1415.6761, found 1415.6739.
Monoglycopeptide 3: HPLC: Rt = 13.72 min (Grad: water 0.1% TFA/acetonitrile (93:7) → (86:14), 15 min, λ = 212 nm). HRMS ESI + (m/z) calcd. for C56H93N18O24 [M + H]+ 1401.6605, found 1401.6546.
NMR experiments
All NMR experiments were recorded on a Bruker Avance 600 MHz spectrometer equipped with a triple channel cryoprobe head. The1H NMR resonances of the peptides 1–3 were completely assigned through standard 2D-TOCSY (30 and 80 ms mixing time) and 2D-NOESY experiments (400 ms mixing time). Solution conditions used for the NMR characterization studies were 1–3 mM (glyco)peptide, 25 mM perdeuterated tris-d11 in 90:10 H2O/D2O, 7.5 mM NaCl, and 1 mM DTT, uncorrected pH 7.4. The assignments were accomplished either at 278 or 298 K. The resonance of 2,2,3,3-tetradeutero-3-trimethylsilylpropionic acid (TSP) was used as a chemical shift reference in the1H NMR experiments (δ TSP = 0 ppm). Peak lists for the 2D-TOCSY and 2D-NOESY spectra were generated by interactive peak picking using the computer aided resonance assignment (CARA) software.
Samples for STD experiments were prepared in perdeuterated 25 mM TRIS-d11 in deuterated water, 7.5 mM NaCl, and 1 mM DTT, uncorrected pH 7.4. STD-NMR experiments were performed at 298 K in the presence of 75 μM UDP, 75 μM MnCl2 with ~880 μM peptide (or GalNAc-O-Me) and 13.5 μM GalNAc-T4 giving a molar ratio of 65:1 peptide:GalNAc-T4.
The STD-NMR spectra were acquired with 1920 transients in a matrix with 64 k data points in t2 in a spectral window of 12335.53 Hz centered at 2819.65 Hz. An excitation sculpting module with gradients was employed to suppress the water proton signals. Selective saturation of the protein resonances (on resonance spectrum) was performed by irradiating at −1 ppm using a series of Eburp2.1000-shaped 90° pulses (50 ms, 1 ms delay between pulses) for a total saturation time of 2.0 s. For the reference spectrum (off resonance), the samples were irradiated at 100 ppm. Proper control experiments were performed with the ligands in the presence and absence of the protein in order to optimize the frequency for protein saturation (−1 ppm) and to ensure that the ligand signals were not affected. However, all the glycopeptides when irradiated at −1 ppm in the absence of protein showed residual saturation on the aliphatic methyl groups in the STD-NMR spectra. This nonspecific saturation was taken into account, by subtraction, when quantifying the STD-NMR data in the presence of the transferase. As well, a blank STD experiment with only the protein was also recorded. The substraction of this protein STD spectrum allowed eliminating the signal background of the protein. In all cases, to accomplish the STD-NMR-derived epitope mapping of each ligand, the STD-NMR total intensities were normalized with respect to the highest STD-NMR response. For Fig. 2b and Supplementary Fig. 1, the STD response of each amino acid corresponds to the average of STD percentages of all amino acid proton resonances that were measured with sufficient accuracy. The signal of the anomeric proton of as well as, the Hα protons of the Ala amino acids of glycopeptides could not be analyzed in the STD-NMR spectra due to their close proximity to the HDO resonance. Proton resonances, from Gly and Pro, appear in the same chemical shift region of the spectrum and were not discriminated.
Molecular docking
Two different docking calculations were conducted to generate the ternary complex.
First, UDP was docked into GalNAc-T4 with the aid of of AutoDock Vina 1.1.234. The predicted binding energies ranged from −6.4 to −5.2 kcal/mol. The Autogrid grid point spacing was set at 0.375 Å, center coordinates of the grid box were 5.4, −2.4, − 7.2 (x, y, z), and number of grid points in xyz was 76, 40, 42, respectively. All allowed torsional bonds were considered rotatable. The 3D structure of the docked ligand to the protein with the lowest-binding energy was used for further calculations (see below).
In a second simulation, glycopeptide T(α-O-GalNAc)GAGAGAGTTPGPG was docked into GalNAc-T4. The Autogrid grid point spacing was set at 0.375 Å, center coordinates of the grid box were 0.0, −9.7, 9.4 (x, y, z), and number of grid points in xyz was 86, 60, 90, respectively. All allowed torsional bonds were considered rotatable, except the bonds involved in the glycosidic linkage, which were fixed to ϕ = 64°, ψ = 110°. The 3D structure of the docked ligand to the protein with lowest-binding energy and with the GalNAc unit located in the lectin domain region, was then used for further calculations.
The combination of both structures was used to build up the ternary complex, which was subjected to MD simulations (see below).
MD simulations
The starting coordinates for the complex between glycopeptide 3, UDP/Mn+2, and GalNAc-T4 were generated combining the X-ray structure of GalNAc-T4 solved in this manuscript and the crystal structure of the activated form of GalNAc-T2 in complex with UDP/Mn+2 and MUC5AC-13 (PDB entry 5AJP). The GalNAc unit of the glycopeptide 3 was located by superimposition to the one solved in the X-ray structure of GalNAc-T4, and UDP was located by superimposition of GalNAc-T2 and T4 proteins. The peptide backbone of compound 3 was situated so that the Thr residues were laying on the proximity of the UDP moiety.
The starting coordinates for chimera 2 were generated by superimposing the X-ray structures of GalNAc-T2 (PDB entry: 5AJP) and GalNAc-T4. Residues 435–445 (PELRVPDHQDI) were manually deleted from GalNAc-T2 and replaced by residues 494–506 (PEVYVPDLNPVIS) of GalNAc-T3 using PyMol (http://www.pymol.org). Residues A507F508 were maintained in GalNAc-T2 to facilitate the linkage between its catalytic and lectin domains. The double (R438A-D444A) and triple mutant (R438A-D444A-F447A) were generated using PyMol. The same protocol was used for other computational chimeras.
Force field parameters for the substrates were generated with the antechamber module of Amber14 using a combination of GLYCAM0635 parameters for the GalNAc unit and the general Amber force field (GAFF) for GDP, with partial charges set to fit the electrostatic potential generated with HF/6-31 G(d) by RESP. The charges are calculated according to the Merz–Singh–Kollman scheme using Gaussian 0936. Each protein was immersed in a truncated octahedral box with a 10 Å buffer of TIP3P water molecules and neutralized by adding explicit counter ions (Na+, Cl−). All subsequent simulations were performed using the ff14SB force field37. A two-stage geometry optimization approach was used. The first stage minimizes only the positions of solvent molecules and ions, and the second stage is an unrestrained minimization of all the atoms in the simulation cell. The systems were then gently heated by incrementing the temperature from 0 to 300 K under a constant pressure of 1 atm and periodic boundary conditions. Harmonic restraints of 30 kcal/mol were applied to the solute, and the Andersen temperature coupling scheme was used to control and equalize the temperature. The time step was kept at 1 fs during the heating stages. Water molecules are treated with the SHAKE algorithm such that the angle between the hydrogen atoms is kept fixed. Long-range electrostatic effects are modeled using the particle-mesh-Ewald method. An 8-Å cutoff was applied to Lennard–Jones and electrostatic interactions. Each system was equilibrated for 2 ns with a 2-fs time step at a constant volume and temperature of 300 K. Production trajectories were then run for additional 100–500 ns under the same simulation conditions.
Transferase assays and kinetics
Specific activity determinations: GalNAc-T glycosylation reactions against (glyco)peptides 1–3 were performed with 75 mM sodium cacodylate, pH 6.5, 1 mM 2-mercaptoethanol, 10 mM MnCl2, 0.25 mM [3H]-radiolabeled UDP-GalNAc (~6 × 108 DPM/µmole, American Radiolabeled Chemicals Inc.), and 1.4 mM (~1.25 mg/ml) of (glyco)peptide (from 4 mM stock in 1 mM TRIS, pH ~7) and varying concentrations (0.01 to 0.5 µM) of transferase (determined by OD280), giving a final reaction volumes of 20 µl in 100 µl capped Eppendorf tubes. Reactions were incubated at 37 °C in a thermostated microplate shaker (Taitec Microincubator M-36) and quenched at the appropriate reaction time by the addition of 20 µl of 250 mM EDTA and placed on ice prior to final workup. Typically, for each transferase, glycosylation reactions were performed at the same time with all three (glyco)peptide substrates (including a no peptide control) using the same batch of transferase and the same UDP-[3H]-GalNAc stock. Before performing time course experiments, the relative activity of each transferase was determined by trial and error against all three peptides to determine the optimal transferase concentration or incubation time that gave no more than 10% peptide glycosylation against the best substrate. Final time course experiments typically consisted of incubation times of 10, 20, and 30, or 30, 60, and 90 or 120 min (depending on activity) and were performed at least twice. Most transferases were also characterized at two different enzyme concentrations (0.01–0.03 and 0.2–0.5 µM) and the results combined. After quenching the reaction and diluting to 4 ml, free UDP and non-hydrolyzed UDP-GalNAc were removed by passage over a Dowex 1 × 8 anion exchange resin (~3 ml column). Total UDP-[3H]-GalNAc utilization (transfer to peptide substrate and transfer to water, i.e, hydrolysis) was determined by difference after scintillation counting (Beckman LS5801 scintillation counter) 1/20 of the sample before and after passage over Dowex 1 × 8. The extent of [3H]-GalNAc transfer to peptide and the extent of hydrolysis was determined by Sephadex G10 gel filtration analysis, typically for the longest reaction time point, as previously described15. Example gel filtration chromatograms are given in Supplementary Fig. 12. Significant hydrolysis was only observed for GalNAc-T4, which was corrected for in determining its specific activity. Specific activity was obtained by calculating at each individual data time point (of the combined time course plots, Fig. 4b), an individual-specific activity, which was then averaged to obtain the overall average specific activity (and standard deviation). Specific activity is reported as nmol of glycopeptide product formed per min, per pmol of enzyme. Substrate activity ratios (i.e., substrate ratios 2/1, 3/1, and 3/2) were also obtained at each individual time point and averaged (with standard deviation). These data are tabulated in Supplementary Fig. 9. Note that the Microsoft Excel least squares fitting function LINEST was used to obtain the linear plots shown in Fig. 4b while fixing the y intercept to zero. Note that the specific activity calculated from the least square slope differed slightly (less than 10%) from that determined from the average of the individual time points.
Detailed enzyme kinetics on GalNAc-T4, GalNAc-T2, the GalNAc-T2 triple mutant and GalNAc-T2 chimera 3: Glycosylation reactions were performed as described above using a fixed enzyme concentration (0.06 µM for GalNAc-T4, and 0.02 µM for GalNAc-T2 and its linker constructs) with varying peptide concentrations (typically 45, 90, 180, 350, 700, 1400, and 2000 μM). Incubation times (10–90 min) were chosen such that no more than 30% of the UDP-GalNAc donor was depleted while typically giving less than 10% peptide glycosylation, after correction for UDP-GalNAc hydrolysis described above. Specific activities were obtained from one to two individual incubation time points performed in two–three separate experiments. These individual-specific activity values were used to calculate the kinetic constants of Km and Vmax using the nonlinear Michaelis–Menten fitting program in GraphPad Prism 7.03.
Determination of site of glycosylation
The determination of substrate glycosylation sites was performed by Edman amino acid sequencing (Applied Biosystems Procise 494 peptide sequencer) of the G10 isolated [3H]-GalNAc glycosylated substrate as previously described15, 16 where each cycle was collected on a fraction collector and scintillation counted (Beckman LS5801 scintillation counter) for [3H]-GalNAc content. Although there is commonly sample-to-sample variability in the [3H]-GalNAc content loaded on the sequencer (due sample losses and different initial UDP-[3H]-GalNAc-specific activities), the observed sites of incorporation were found to be identical between different experiments with the same transferase and substrate. The presence of [3H]-GalNAc lag after a peak of [3H]-GalNAc incorporation is commonly observed in these determinations, which is due to the poor extraction from the sample filter of the glycosylated-PTH residues compared with the standard amino acid PTH derivatives15, 16.
Data availability
The coordinate and structure factor for GalNAc-T4 in complex with monoglycopeptide 3 has been deposited in the Worldwide Protein Data Bank (wwPDB) with the accession code 5NQA. All relevant data are available from the authors on reasonable request.
Electronic supplementary material
Acknowledgements
We thank synchrotron radiation sources DLS (Oxford) and in particular beamline I03 (experiment number MX10121-7). We thank ARAID, MEC (CTQ2013-44367-C2-2-P, BFU2016-75633-P, CTQ2015-67727-R, CTQ2015-70524-R, and RYC-2013-14706), the National Institutes of Health (GM113534, and instrument grant GM113534-01S), the Danish National Research Foundation (DNRF107), the FCT-Portugal (UID/Multi/04378/2013 and PTNMR Project No 022161), and the DGA (B89) for the financial support. I.C. thanks Universidad de La Rioja for the FPI grant. F.M. thanks FCT-Portugal for IF Investigator. E.L.-N. acknowledges her postdoctoral EMBO fellowship ALTF 1553-2015 co-funded by the European Commission (LTFCOFUND2013, GA-2013-609409) and Marie Curie Actions. H.C. and J.J.-B. thank EU for the TOLLerant project. The research leading to these results has also received funding from the FP7 (2007–2013) under BioStruct-X (grant agreement No. 283570 and BIOSTRUCTX_5186). We also thank BIFI (Memento cluster) and CESGA for computer support.
Author contributions
R.H.-G. designed the crystallization construct and solved the crystal structure. E.L.-N., M.R., and R.H.-G. purified the enzymes, crystallized the complex, and refined the crystal structure. I.C., F.C., and J.M.P. synthetized the glycopeptides. F.C. and G.J.-O. performed the MD simulations. H.Co., A.D., J.J.-B., and F.M. performed and analyzed the NMR experiments. T.A.G. and E.J.P.D. performed the kinetic studies together with the Edman amino acid sequencing. R.H.-G. wrote the article with the main contribution of T.A.G., F.C., G.J.-O., H.Cl., and F.M. All authors read and approved the final manuscript.
Competing interests
The authors declare no competing financial interests.
Footnotes
Matilde de las Rivas, Erandi Lira-Navarrete and Earnest James Paul Daniel contributed equally to this work.
Electronic supplementary material
Supplementary Information accompanies this paper at 10.1038/s41467-017-02006-0.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
- 1.Bennett EP, et al. Control of mucin-type O-glycosylation: a classification of the polypeptide GalNAc-transferase gene family. Glycobiology. 2012;22:736–756. doi: 10.1093/glycob/cwr182. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 2.Gill DJ, Clausen H, Bard F. Location, location, location: new insights into O-GalNAc protein glycosylation. Trends Cell Biol. 2011;21:149–158. doi: 10.1016/j.tcb.2010.11.004. [DOI] [PubMed] [Google Scholar]
- 3.Hurtado-Guerrero R. Recent structural and mechanistic insights into protein O-GalNAc glycosylation. Biochem. Soc. Trans. 2016;44:61–67. doi: 10.1042/BST20150178. [DOI] [PubMed] [Google Scholar]
- 4.Lira-Navarrete E, et al. Dynamic interplay between catalytic and lectin domains of GalNAc-transferases modulates protein O-glycosylation. Nat. Commun. 2015;6:6937. doi: 10.1038/ncomms7937. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 5.Tabak LA. The role of mucin-type O-glycans in eukaryotic development. Semin. Cell Dev. Biol. 2010;21:616–621. doi: 10.1016/j.semcdb.2010.02.001. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 6.Schjoldager KT, Clausen H. Site-specific protein O-glycosylation modulates proprotein processing - deciphering specific functions of the large polypeptide GalNAc-transferase gene family. Biochim. Biophys. Acta. 2012;1820:2079–2094. doi: 10.1016/j.bbagen.2012.09.014. [DOI] [PubMed] [Google Scholar]
- 7.Schjoldager KT, et al. Deconstruction of O-glycosylation–GalNAc-T isoforms direct distinct subsets of the O-glycoproteome. EMBO Rep. 2015;16:1713–1722. doi: 10.15252/embr.201540796. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 8.Hanson, R. L. & Hollingsworth, M. A. Functional consequences of differential O-glycosylation of MUC1, MUC4, and MUC16 (downstream effects on signaling). Biomolecules6, E34 (2016). [DOI] [PMC free article] [PubMed]
- 9.Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014;42:D490–D495. doi: 10.1093/nar/gkt1178. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 10.Fritz TA, Hurley JH, Trinh LB, Shiloach J, Tabak LA. The beginnings of mucin biosynthesis: the crystal structure of UDP-GalNAc:polypeptide alpha-N-acetylgalactosaminyltransferase-T1. Proc. Natl Acad. Sci. USA. 2004;101:15307–15312. doi: 10.1073/pnas.0405657101. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 11.Fritz TA, Raman J, Tabak LA. Dynamic association between the catalytic and lectin domains of human UDP-GalNAc:polypeptide alpha-N-acetylgalactosaminyltransferase-2. J. Biol. Chem. 2006;281:8613–8619. doi: 10.1074/jbc.M513590200. [DOI] [PubMed] [Google Scholar]
- 12.Lira-Navarrete E, et al. Substrate-guided front-face reaction revealed by combined structural snapshots and metadynamics for the polypeptide N-acetylgalactosaminyltransferase 2. Angew. Chem. Int. Ed. Engl. 2014;53:8206–8210. doi: 10.1002/anie.201402781. [DOI] [PubMed] [Google Scholar]
- 13.Kubota T, et al. Structural basis of carbohydrate transfer activity by human UDP-GalNAc: polypeptide alpha-N-acetylgalactosaminyltransferase (pp-GalNAc-T10) J. Mol. Biol. 2006;359:708–727. doi: 10.1016/j.jmb.2006.03.061. [DOI] [PubMed] [Google Scholar]
- 14.Wandall HH, et al. The lectin domains of polypeptide GalNAc-transferases exhibit carbohydrate-binding specificity for GalNAc: lectin binding to GalNAc-glycopeptide substrates is required for high density GalNAc-O-glycosylation. Glycobiology. 2007;17:374–387. doi: 10.1093/glycob/cwl082. [DOI] [PubMed] [Google Scholar]
- 15.Gerken TA, et al. The lectin domain of the polypeptide GalNAc transferase family of glycosyltransferases (ppGalNAc Ts) acts as a switch directing glycopeptide substrate glycosylation in an N- or C-terminal direction, further controlling mucin type O-glycosylation. J. Biol. Chem. 2013;288:19900–19914. doi: 10.1074/jbc.M113.477877. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 16.Revoredo L, et al. Mucin-type O-glycosylation is controlled by short- and long-range glycopeptide substrate recognition that varies among members of the polypeptide GalNAc transferase family. Glycobiology. 2016;26:360–376. doi: 10.1093/glycob/cwv108. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 17.Raman J, et al. The catalytic and lectin domains of UDP-GalNAc:polypeptide alpha-N-Acetylgalactosaminyltransferase function in concert to direct glycosylation site selection. J. Biol. Chem. 2008;283:22942–22951. doi: 10.1074/jbc.M803387200. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 18.Hassan H, et al. The lectin domain of UDP-N-acetyl-D-galactosamine: polypeptide N-acetylgalactosaminyltransferase-T4 directs its glycopeptide specificities. J. Biol. Chem. 2000;275:38197–38205. doi: 10.1074/jbc.M005783200. [DOI] [PubMed] [Google Scholar]
- 19.Posey AD, Jr, et al. Engineered CAR T cells targeting the cancer-associated Tn-Glycoform of the membrane mucin MUC1 control adenocarcinoma. Immunity. 2016;44:1444–1454. doi: 10.1016/j.immuni.2016.05.014. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 20.McEver RP, Cummings RD. Role of PSGL-1 binding to selectins in leukocyte recruitment. J. Clin. Invest. 1997;100:S97–S103. doi: 10.1172/JCI119556. [DOI] [PubMed] [Google Scholar]
- 21.Bennett EP, et al. Cloning of a human UDP-N-acetyl-alpha-D-Galactosamine:polypeptide N-acetylgalactosaminyltransferase that complements other GalNAc-transferases in complete O-glycosylation of the MUC1 tandem repeat. J. Biol. Chem. 1998;273:30472–30481. doi: 10.1074/jbc.273.46.30472. [DOI] [PubMed] [Google Scholar]
- 22.Kato K, et al. Polypeptide GalNAc-transferase T3 and familial tumoral calcinosis. Secretion of fibroblast growth factor 23 requires O-glycosylation. J. Biol. Chem. 2006;281:18370–18377. doi: 10.1074/jbc.M602469200. [DOI] [PubMed] [Google Scholar]
- 23.Pedersen JW, et al. Lectin domains of polypeptide GalNAc transferases exhibit glycopeptide binding specificity. J. Biol. Chem. 2011;286:32684–32696. doi: 10.1074/jbc.M111.273722. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 24.Topaz O, et al. Mutations in GALNT3, encoding a protein involved in O-linked glycosylation, cause familial tumoral calcinosis. Nat. Genet. 2004;36:579–581. doi: 10.1038/ng1358. [DOI] [PubMed] [Google Scholar]
- 25.Khetarpal SA, et al. Loss of function of GALNT2 lowers high-density lipoproteins in humans, nonhuman primates, and rodents. Cell Metab. 2016;24:234–245. doi: 10.1016/j.cmet.2016.07.012. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 26.Lira-Navarrete E, et al. Structural insights into the mechanism of protein O-fucosylation. PLoS ONE. 2011;6:e25365. doi: 10.1371/journal.pone.0025365. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 27.Valero-Gonzalez J, et al. A proactive role of water molecules in acceptor recognition by protein O-fucosyltransferase 2. Nat. Chem. Biol. 2016;12:240–246. doi: 10.1038/nchembio.2019. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 28.Kabsch W. Xds. Acta Crystallogr. D Biol. Crystallogr. 2010;66:125–132. doi: 10.1107/S0907444909047337. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 29.Winn MD, et al. Overview of the CCP4 suite and current developments. Acta Crystallogr. D Biol. Crystallogr. 2011;67:235–242. doi: 10.1107/S0907444910045749. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 30.Collaborative Computational Project, Number 4. The CCP4 suite programs for protein crystallography. Acta Crystallogr. D Biol. Crystallogr.50, 760–763 (1994). [DOI] [PubMed]
- 31.Emsley P, Cowtan K. Coot: model-building tools for molecular graphics. Acta Crystallogr. D Biol. Crystallogr. 2004;60:2126–2132. doi: 10.1107/S0907444904019158. [DOI] [PubMed] [Google Scholar]
- 32.Murshudov GN, et al. REFMAC5 for the refinement of macromolecular crystal structures. Acta Crystallogr. D Biol. Crystallogr. 2011;67:355–367. doi: 10.1107/S0907444911001314. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 33.Plattner C, Hofener M, Sewald N. One-pot azidochlorination of glycals. Org. Lett. 2011;13:545–547. doi: 10.1021/ol102750h. [DOI] [PubMed] [Google Scholar]
- 34.Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2010;31:455–461. doi: 10.1002/jcc.21334. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 35.Kirschner KN, et al. GLYCAM06: a generalizable biomolecular force field. Carbohydrates. J. Comput. Chem. 2008;29:622–655. doi: 10.1002/jcc.20820. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 36.Frisch, M. J. et al. Gaussian 09, Revision D.01 (Gaussian, Inc., Wallingford, 2009).
- 37.Hornak V, et al. Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins. 2006;65:712–725. doi: 10.1002/prot.21123. [DOI] [PMC free article] [PubMed] [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.
Supplementary Materials
Data Availability Statement
The coordinate and structure factor for GalNAc-T4 in complex with monoglycopeptide 3 has been deposited in the Worldwide Protein Data Bank (wwPDB) with the accession code 5NQA. All relevant data are available from the authors on reasonable request.