NMR analysis demonstrates immunoglobulin G N-glycans are accessible and dynamic.

The N-glycan at Asn297 of the immunoglobulin G Fc fragment modulates cellular responses of the adaptive immune system. However, the underlying mechanism remains undefined, as existing structural data suggest the glycan does not directly engage cell surface receptors. Here we characterize the dynamics of the glycan termini using solution NMR spectroscopy. Contrary to previous conclusions based on X-ray crystallography and limited NMR data, our spin relaxation studies indicate that the termini of both glycan branches are highly dynamic and experience considerable motion in addition to tumbling of the Fc molecule. Relaxation dispersion and temperature-dependent chemical shift perturbations demonstrate exchange of the α1-6Man-linked branch between a protein-bound and a previously unobserved unbound state, suggesting the glycan samples conformational states that can be accessed by glycan-modifying enzymes and possibly glycan recognition domains. These findings suggest a role for Fc-glycan dynamics in Fc-receptor interactions and enzymatic glycan remodeling.

As a primary component in the second phase of the adaptive immune system, Immunoglobulin G (IgG) binds to foreign particles in the blood and targets these invaders for destruction. IgG engaged with an antigen triggers this process by intepan class="Disease">racting with cell surface receptors and the class="Chemical">pan class="Gene">C1q component of complement through its Fc domain1. Specific recognition of IgG then leads to localized inflammation by activating the immune system; this activation threshold must be balanced to avoid hyposensitivity to foreign particles and hypersensitivity to autoantigens. Either extreme results in serious disease, enhanced sensitivity to infection on one hand and debilitating autoimmune disease on the other hand2. There is good reason to believe that the carbohydrates attached to the Fc fragment play a role in maintaining this balance in sensitivity3 and that understanding the conformation and conformational excursions of the glycans could help in the design in appropriate modulators of immune sensitivity. NMR studies presented here take a step toward that understanding.

IgG Fc from healthy pan class="Species">humans has a complex-type biantennary class="Chemical">pan class="Chemical">N-glycan attached to Asn297 of the Cγ2 domains of the dimeric Fc fragment (see Figure 1)3. While the Fc glycans are quite heterogeneous, a high percentage is normally terminated with galactose. Deviations from this pattern correlate well with certain diseases as well as differential effects on inflammatory response4. For example, in those suffering from rheumatoid arthritis (RA), a crippling autoimmune disorder, these glycans are mostly ungalactosylated leaving N-acetylglucosamine residues at the termini of one or both branches of the glycans5. Recent studies also show that the Fc fragment, once enzymatically sialylated in vitro, is potently anti-inflammatory and reduces RA disease severity in a mouse model6,7. This observation has significant implications for possible improvements in intravenous IgG therapy for patients suffering from severe cases of RA by enzymatically adding sialic acids to the normal glycan termini of pooled IgG7. Other correlations between glycan compositions and Fc interactions with its receptors have also been noted6,8–11. In the case of the anti-inflammatory activity of Fc with sialylated glycans, a specific receptor has been postulated12, but, it is not clear how changes in glycan composition modulate interactions with receptors.

Figure 1

Immunoglobulin G contains an N-linked glycan at Asn297

(a) The domain organization of IgG1showing the antigen-binding Fab fragments, the Fc fragment, and the N-glycan. (b) The Fc glycan as remodeled in this work is predominantly a galactose-terminated complex-type, biantennary, fucosylated carbohydrate. Open circles indicate galactose, filled squares N-acetylglucosamine, grey circles mannose, and open triangles fucose residues.

Our current structupan class="Disease">ral knowledge of Fc and its class="Chemical">pan class="Chemical">glycans is actually quite substantial, but this has done surprisingly little to resolve question regarding glycan modulation of receptor interactions. For example, in x-ray diffraction derived models of Fc – receptor complexes13–15 the glycans do not directly interact with the surface of the receptors and glycan modifications seem to induce little change in Fc geometry near the interaction site9. The ends of IgG Fc glycans in galactose terminated cases are well resolved in several x-ray diffraction-derived models and appear to be buried between two Cγ2 domains16,17, leaving their accessibility for any sort of interaction in question. Solution NMR measurements of the same Fc glycoform indicated that the α1–6Mannose (Man) branch was completely immobilized through interactions with the protein surface while the α1–3Man terminus was more dynamic18,19. Immobilization presumably would strongly restrict the α1–6Man branch terminus in interactions with receptors and glycan modifying enzymes such as the sialyltransferases and glycosidases (neuraminidase, β-galactosidase, N-acetylglucosaminidase). It would also indirectly restrict the α1–3Man branch in these same interactions.

Despite these inferences from structupan class="Disease">ral data, class="Chemical">pan class="Chemical">glycan termini have been remodeled in vitro using various enzymes6,18,20,21. Attempts to add sialic acids to branch termini of Fc domains using human sialyltransferase St6GalI21 suggest a strong preference for addition to the more dynamic α1–3Man branch. However, the sialyltransferase prefers the α1–3Man galactose over the α1–6Man galactose even using isolated glycans, with this preference only moderately enhanced using Fc-conjugated glycans. Thus, there is an apparent contradiction between the current structural picture of tight coordination of galactose-terminated glycans and the availability of the glycan termini to various glycan modifying enzymes. This provides motivation for the study presented here using solution nuclear magnetic resonance (NMR) spectroscopy to study dynamic sampling of other glycan conformations.

Solution NMR spectroscopy is an ideal technique for the chapan class="Disease">racterization of the structure and dynamics of conformationally labile systems22. class="Chemical">pan class="Chemical">Polysaccharides themselves are conformationally labile due a lack of short or long range stabilizing interactions. This has often made their study by purely structural methods difficult. Terminal residues of glycans are often unresolved in x-ray crystal structures of glycoproteins either because of conformational heterogeneity or because of enzymatic truncation to improve crystallizability23,24. Even when observed, one must remember that the conformations seen may represent only one low energy conformer of an ensemble of structures accessible in solution. In contrast to x-ray diffraction, all members of a conformationally labile moiety often contribute to observable NMR parameters.

Using NMR, the timescales, and sometimes amplitudes of motion in pan class="Chemical">carbohydrates, can be chaclass="Chemical">pan class="Disease">racterized based on well-defined spin relaxation techniques and theory25,26. Current applications of NMR spin relaxation techniques have been developed primarily for observation of backbone amide or sidechain methyl groups in studies of dynamics in isotopically labeled proteins25–27. Techniques to label, assign, and characterize the dynamics of glycoprotein glycans are less well developed (reviewed in ref24). Here we present methods to incorporate 13C labels into glycoproteins and assign the resonances using a combination of enzymology, mass spectrometry and solution NMR spectroscopy. NMR pulse sequences specific for the resulting labeling schemes are also reported that can be used to address a number of challenges in glycoprotein structural biology. Using these methods, we here report on the dynamics and probable structure of the residues at the termini of Fc glycans.

Results

Isotope labeling of IgG glycans

To selectively enhance the signals from the IgG pan class="Chemical">N-glycans, we added class="Chemical">pan class="Chemical">13C-enriched galactose residues to the glycan termini of intact IgG1 from myeloma cells. A two dimensional 1H-13C correlation experiment (a heteronuclear multiple quantum coherence (HMQC) spectrum) showed two sets of resonances, a stronger set and a weaker set (Figure 2). The origin of the two sets was further investigated by separating the Fc and Fab fragments using proteolysis of the IgG1 molecule followed by preferential retention of the Fc fragment on an affinity column. A similar two-dimensional spectrum of the non-retained fraction (Figure 2) showed the more intense set of resonances to belong to the Fab fragment. Despite the lack of a consensus N-glycosylation site on the constant regions of IgG Fab fragments, ~30% of circulating IgGs contain an N-glycan attached to the variable region of the molecule28, and higher glycan mobility at these glycosylated sites likely allows these to be detected with more intensity. An experiment using purified Fc fragment (Figure 3) also showed the less intense signals observed in the IgG1 spectrum. It is interesting to note that the resonances appeared in the same position when using the intact IgG1 molecule and the isolated fragments; however, lineshapes were sharper and signals more intense in spectra of the smaller fragments. These data suggested the presence of the Fab fragment did not affect the chemical environment of the Fc glycan, allowing studies of Fc fragment resonances in the isolated fragment.

Figure 2

A 2D 13C-1H HMQC correlation spectrum of 13C-Galactose labeled IgG1 and IgG1 Fab fragment

(a) IgG1 and (b) the IgG1 Fab fragment at 50°C. Vertical grey lines indicate resonances at the same 13C and 1H frequencies in both spectra.

Figure 3

2D 13C-HMQC spectra and assignments of 13C-Galactose labeled IgG Fc

(a) uniformly 13C-Galactose labeled Fc fragment and (b) 13C2-Galactose labeled Fc fragment. The extracted lines show 1H linewidth measurements of the H2-C2 crosspeaks extracted from this spectrum. (c) IgG Fc fragment with the α1–6Man-linked Galactose enriched with 13C-Galactose. Filled (hollow) arrows denote a residue on the α1–3Manlinked (α1–6Man-linked) branch of the Fc-conjugated biantennary glycan. Anomeric peaks were overlapped with the residual HOD signal but were observed at lower temperatures.

We, therefore, investigated the Fc pan class="Chemical">glycan more thoroughly using Fc fclass="Chemical">pan class="Disease">ragments prepared from pooled human serum. Figure 3 shows a 13C- HMQC spectrum for a sample of IgG-Fc labeled with uniformly 13C-enriched galactose at the termini of both branches of its glycans. This spectrum clearly showed two types of resonances, one broad set and one sharp (the broad set was less easily observed in the intact IgG molecule of Figure 2 because of the additional line broadening due to its larger size). The difference in the two sets for the Fc sample was more clearly shown in a spectrum of Fc fragments labeled with galactose enriched at just the C2 position (Figure 3). This labeling allowed specific assignments of the C2 carbon resonances and permitted examination of linewidths in cross-sections from rows through the two peaks (see insert in Figure 3).

To obtain bpan class="Disease">ranch specific assignments we preclass="Chemical">pared an Fc fclass="Chemical">pan class="Disease">ragment enriched primarily on the α1–6Man terminus with 13CU-galactose. We first introduced unlabeled galactose residues at the termini of the Fc glycan, which natively are not fully galactosylated on either the α1–6 or α1–3Man branch and are seldom sialylated, then sialylated the glycans using St6Gal1 which placed a sialic acid residue preferentially on the α1–3Man branch21. Exogalactosidase treatment then removed the α1–6Man branch galactose but failed to remove the largely unlabeled α1–3Man galactose as it was protected by the sialic acid. Finally, we galactosylated the glycans with uniformly 13C-labeled galactose on the α1–6Man branch terminus (see Methods). The resulting spectrum is shown in Figure 3. The sharper set of peaks was much reduced in intensity in comparison to Figure 3, allowing assignment of these peaks to the α1–3Man branch and the broad peaks to the α1–6Man branch. This assignment is consistent with previous observations which identified a broad galactose anomeric resonance as the α1–6Man galactose of the Fc-conjugated glycan18.

A comparison of our spectpan class="Disease">ra to the chemical shifts of a free class="Chemical">pan class="Chemical">glycan and free galactose revealed the atom-specific identities of these peaks (ref 29–31, Biological Magnetic Resonance data Bank: www.bmrb.wisc.edu). The chemical shift pattern of the α1–3Man branch resonances was quite similar to that reported in the literature for isolated glycans of similar structure; however, all resonances were shifted slightly downfield in both 1H and 13C dimensions. In contrast, the α1–6Man branch resonances experienced significant deviations when compared to a free glycan. The anomeric resonances of the α1–3 and α1–6Man branches (not shown in Figures 2 and 3) had the least similar chemical shift values, followed by the 2 position; the remaining positions had more similar shifts (Supplementary Results, Supplementary Figure 1).

Dynamics of IgG glycans

NMR resonance linewidths, which are directly related to tpan class="Disease">ransverse relaxation class="Chemical">pan class="Species">rates (Rs), provide a very ready measure of dynamic properties with narrower lines indicating decreased relaxation rates and increased rates of rotational motion. Figure 3 shows that the proton linewidth of the H2 resonance of the α1–6Man branch was nearly twice that of the corresponding resonance from the α1–3Man branch. This immediately suggested the α1–6Man branch to be more immobilized. 13C linewidths can be more quantitatively interpreted because the dominant relaxation mechanisms can often be assumed to originate in the dipole-dipole interaction with directly bonded protons. Previous work on 13C spin relaxation had in fact interpreted linewidth data of the Gal C1s18. This linewidth was very close to that for protein 13C αs and suggested an interaction between the α1–6Man branch and the protein strong enough to completely immobilize the branch for periods longer than the protein tumbling time (~20ns). However, at higher fields the 13C line widths of the α1–6Man 13C2 resonance in our construct were more than 3 times that of the corresponding carbon on the α1–3Man branch (Supplementary Figure 2) and a simple interpretation using an isotropic model gave an effective correlation time of 37ns, a time considerably longer than the hydrodynamic predictions for the Fc fragment.

These pan class="Disease">anomalies were confirmed in a more complete set of spin relaxation measurements (Figure 4 and Table 1). There was, in fact, field dependence to simple R measurements with a smaller relaxation class="Chemical">pan class="Species">rate observed at 14.0T as opposed to 21.1T. However, R measurements, even at this lower field were not consistent with predictions based on an isotropic tumbling model and suggested that if solely dipolar in origin, R would be even less at this field strength. Both the field dependence and these inconsistencies suggested contributions to R from chemical exchange mechanisms that originate from the α1–6Man branch sampling multiple conformational states on the µs to ms timescale and modulating chemical shift in this process. At least some of the states would have to have substantial internal motion to raise the effective R value.

Figure 4

13C spin relaxation measurements of Galactose resonances

(a) R1 and (b) R2 measurements of 13C6. (c) R1 and (d) R1ρ measurements of 13C2. The data are fit with an equation describing a single exponential decay; using the rates shown in Table 1. Where the rates are significantly different, a solid (dashed) line denotes the α1–3Man-linked (α1–6Man-linked) fits. R1ρ measurements were obtained at an applied field strength of 5000Hz.

Table 1

13C Gal-Fc Relaxation Measurements

	α1–3Man-linked	α1–6Man-linked
	Relaxation Rates	Relaxation Rates
21.1T	Galactose C6	Galactose C6
R1	2.8 ± 0.1 s−1	3.1 ± 0.4 s−1
R2	48.4 ± 5.5 s−1	59.7 ± 9.2 s−1
21.1T	Galactose C2	Galactose C2
R1	1.2 ± 0.1 s−1	1.2 ± 0.2 s−1
R1ρ	20.1 ± 1.0 s−1	56.1 ± 10 s−1
R2 linewidth	40 s−1	150 s−1
14.0T
R1ρ	14.7 ± 1.1 s−1	48.0 ± 5.4 s−1
R2 linewidth	28 s−1	107 s−1
R2	29 s−1	n.d.

n.d.- not determined due to low signal / noise

Fortunately, there are good experiments to explicitly probe the existence of chemical exchange contributions. One involves the use of relaxation dispersion experiments in which tpan class="Disease">ransverse relaxation class="Chemical">pan class="Species">rates are measured with Carr-Purcell-Meiboom-Gill pulse sequences having variable delays between 180° pulses (equivalently having variable pulsing rates (νCPMG))26,27. In these experiments, resonances that experience multiple chemical environments on a timescale near νCPMG become less intense at slow pulsing rates and more intense at higher rates. Results from these experiments at 50°C with a field strength of 21.1T or 18.8T are presented in Figure 5. At low pulsing rates the measurements approach predictions from line widths. At higher pulsing rates chemical shifts are rapidly refocused and, in the limit of very rapid pulsing the effects of chemical exchange, are removed. It is significant that the R value for the α1–6Man branch fell off with pulsing rate, while the R value for the α1–3Man branch did not. The change in R in the presence of chemical exchange between two states (state A and state B) was well described by Eq (1) in which the contribution of chemical exchange to the observed R is dependent upon the rate of exchange, the time between pulses, and the scaling factor φ. Fitting this equation to data for the α1–6Man branch yielded an effective exchange constant of 5300 ± 1700 s−1. The scaling factor φ in principle can yield information about the nature of the exchange process by describing populations of states and chemical shift changes on moving between these states. But these factors are not easily separated without additional information.

Figure 5

Relaxation dispersion and temperature-dependent chemical shift measurements show evidence of two states

Relaxation dispersion experiments (a reveal motion on a µs timescale for the α1–6Man-linked branch but not the α1–3Man-linked branch. When fitted to a two state model and data collected at three different magnetic fields, an exchange rate (kEX) of 5300 ± 1700 s−1 was estimated by fitting equation (1). (b) The chemical shift of the α1–6Man-linked 13C2 resonance approaches a saturation point at high temperature and permits the estimation of chemical shift values for each of two states using Eq (3). Chemical shift asymptotes for the high temperature state A and low temperature state B are shown with dashed lines.

The effect of chemical exchange pan class="Species">rates on R can also be investigated with R measurements. In low and modeclass="Chemical">pan class="Species">rate locking fields, R should approach the non-exchange contribution to R (R). R for both the α1–6Man branch and the α1–3Man branch did in fact approach the rapid pulsing limit of R from the CPMG experiments (Figure 5 and Table 1). The large deviations from R based on the linewidth for the α1–6Man branch was clearly consistent with large chemical exchange contributions. There was a smaller but significant deviation from R based on the linewidth for the α1–3Man branch. This was not seen in the CPMG experiments suggesting that there may be some small exchange contributions at a faster rate not accessible with CPMG experiments for this branch.

Bound-state properties of IgG glycans

Ideally one would like to chapan class="Disease">racterize the conformations involved in exchange processes more thoroughly. Tclass="Chemical">pan class="Disease">raditional NMR techniques involving the measurement of distance-dependent NOEs or scalar couplings were unsuccessful, likely due to the dynamic nature of these glycans (data not shown) and the tendency of glycan interactions to be mediated by hydrogen bonds from exchangeable hydroxyl protons rather than the protons attached to carbons. Dynamic sampling of multiple states reduces NOEs involving carbon-bound protons from any one state, and under solvent conditions used here, un-observable deuterons replaced hydroxyl protons. While in principle, the relaxation dispersion experiments provided some structural information in terms of a product of fractional populations and chemical shift deviations squared, additional information was needed to separate these factors.

One can in fact use the tempepan class="Species">rature dependence of the α1–6class="Chemical">pan class="Species">Man galactose chemical shifts, along with some assumptions about the nature of exchange to provide some of this information. The Fc fragment is thermostable and permits measurements across a wide temperature range with the first of three unfolding events occurring with a Tm >70°C32. Our observations indicated the chemical shifts of the Fc α1–6Man galactose residues were quite sensitive to temperature. Among these, the 13C2 resonance showed the greatest displacement in the 5–50°C range with a chemical shift value that approached a plateau of 75.3 ± 0.2 ppm at the high temperature limit (Figure 5b). The plateau value, representing state A, was reasonably consistent with chemical shift values for the same site in the free glycan. While this analysis focused on the C2 carbon, it is significant that several other resonances showed similar trends and moved from a perturbed resonance position at lower temperature toward a free glycan position at higher temperatures.

Using the more precisely determined chemical shift of state A (Figure 5 and Table 2) in conjunction with the observed chemical shift at 50°C (Figure 3) and the CPMG relaxation dispersion data (Figure 5) we extpan class="Disease">racted a value for the population of state A (class="Chemical">pan class="Chemical">pA) of 0.8 using Eq (4). The value of pA may then be used to deconvolute φ and derive a value for Δω using the relationship in Eq 2. The state B chemical shift of 76.4 ±0.5 ppm was determined using the resultant value for Δω of 1500 ± 300 rad s−1 along with the state A chemical shift.

Table 2

Kinetic and Thermodynamic Parameters

kEX	5300 ± 1700 s−1
φEX 21.1T	380,000 ± 60,000 rad2 s−2
R2o 21.1T	26 ± 13 s−1
φEX 18.8T	300,000 ± 100,000 rad2 s−2
R2o 18.8T	39 ± 7 s−1
ΔH	9.2 ± 0.6 kcal mol−1
ΔS	0.032 ± 0.002 kcal mol−1 K−1
δA	75.3 ± 0.2 ppm
δB1	76.4 ± 0.5 ppm
δB2	76.4 ± 0.5 ppm
Δω1	1600 rad s−1
Δω2	1500 ± 300 rad s−1

values derived from chemical shift vs. T fit

values derived from CPMG measurements

The tempepan class="Species">rature-dependent chemical shift data alone may be analyzed to provide an independent measure of population changes as well as the state B chemical shift, as long as one can access a reasonable class="Chemical">pan class="Disease">range of fractional population (0.2> pApB > 0.1) and the chemical shift difference (Δω) is relatively large. The data in Figure 5,fit with Eq (3), show the chemical shift at the low (state B) temperature limit and gave ΔH and ΔS values (Table 2). Despite the limited fractional population range, these values were in good agreement with those determined using the CPMG data. The data indicated that populations were approximately equal at 15°C and moved to higher populations of state A (the state with shifts similar to a free glycan) at temperatures above 15°C. At 37° state A was approximately 70% populated. The trend was consistent with a large positive ΔS value and a significant distribution of sampled conformations in state A.

Assuming that we have properly removed effects of chemical exchange, it is of some interest to assess the internal dynamics of the pan class="Chemical">galactose residues at the ends of the two bclass="Chemical">pan class="Disease">ranches. For the α1–6Man branch this can be done by determining the order parameter (S) and the internal motion correlation time (τe) that characterize the motion of the galactose 13C2-H bond vectors relative to the Fc polypeptide in the un-bound state. This can be done using the Lipari-Szabo formalism33,34 and R and R values corrected for contributions of the bound state assuming isotropic tumbling of Fc at its hydrodynamic limit. Order parameters range from 0 to 1, with 0 corresponding to unrestrained excursions and 1 corresponding to a completely immobilized state. The value for S of the 13C2-H bond vector was 0.3, suggesting this vector is allowed considerable motional excursions relative to the Fc polypeptide. The value of τe was 2 ns. This was considerably shorter than the correlation time characterizing Fc tumbling indicating a high degree of flexibility. Motions of the α1–3Man branch can be analyzed in a similar fashion, but without having to correct R and R values for a bound state. The order parameters and internal correlation time were 0.4 and 2 ns suggesting similar degrees of conformational sampling and levels of flexibility for both branches when in an unbound state.

Discussion

Prior X-pan class="Disease">ray structuclass="Chemical">pan class="Disease">ral data and NMR dynamics analysis had clearly led to the impression that glycans attached to the Fc fragment of IgG are both occluded within a cavity between the halves of the dimeric structure and that the α1–6Man branch of the glycan was firmly anchored to the protein surface. IgG structures are in fact one of a small number of glycoprotein structures showing discrete electron density from extended parts of attached glycans. The existence of a protein-bound state of this glycan branch has been accepted beginning with the first high-resolution structure of the Fc fragment17. Prior studies of the Fc glycan using solution NMR concluded that the α1–6Man branch terminus was firmly bound to the protein surface based on one-dimensional 1H spectroscopy measurements of relaxation rates, or the use of linewidths from one-dimensional 13C spectroscopy to estimate these rates18,19. The examination here, based on multiple field studies and an extended set of spin relaxation experiments, leads to a somewhat different picture that may be more consistent with other biochemical data. Our prior work on sialylation of human IgG Fc, as well as other previous work on IgG Fc, shows that at least the galactose terminated α1–3Man branch is readily accessible for enzymatic addition of sialic acid11,21,35. Our work suggests rates of addition are attenuated by approximately a factor of five by attachment of the glycans to the protein21. It was also shown that the preferential addition to the α1–3Man branch is not solely a result of attachment of the glycan to Fc, but is inherent in enzyme specificity observed even when using free glycan as the substrate. There have been some mutational studies supporting specific interactions between residues on the proteins and the glycan. Interestingly, the proposed polypeptide interactions, notably through Phe241 and Phe243, inhibit but do not prevent glycan processing32,36.

Based on data presented here the α1–3pan class="Species">Man bclass="Chemical">pan class="Disease">ranch showed substantial levels of internal motion and chemical shift patterns that well approximate those of the free glycan (Supplementary Figure 1), though all resonances were shifted slightly downfield. Both suggested nothing more than minimal interactions of at least the terminal galactose with the protein surface. This could be consistent with multiple conformers that populate the cavity in the Fc crystal structure. Terminal residues on this branch do not usually give sufficient electron density to allow definitive modeling. However, the data presented here, as well as the previous enzymatic data are more consistent with conformational excursions that fully expose the terminus of the α1–3Man branch.

Data for the α1–6pan class="Species">Man terminus of the Fc class="Chemical">pan class="Chemical">glycan also showed substantial mobility with excursions between at least two distinct states. The time scales of these excursions (~5000 s−1) and the populations of the states would be adequate to allow access to enzymatic modification with reductions in efficiency entirely consistent with observation. The chemical shifts estimated for the two states of the α1–6Man branch of the glycan support a transition between something similar to that observed in crystal structures (pdb - 1l6x16, 1dn237) and one that is more fully exposed. The dominant, high temperature, state again has chemical shifts approaching those of a free glycan. This would be consistent with a more exposed structure.

We can pan class="Species">rationalize chemical shift changes based on existing crystal structures. The high resolution Fc structures with low tempeclass="Chemical">pan class="Species">rature (B) factors for the α1–6Man linked terminus (pdb -1l6x, 1dn2) have glycosidic torsion angles ϕ and ψ for the α1–6Man linked galactose-N-acetylglucosamine unit of ~65° and ~11°, respectively, which are close to the values of 55° and 6°, respectively, for a calculated energy-minimized structure and 45° and 15°, respectively for the crystal structure of methyl-β-d-lactoside38. This observation is consistent with the general conservation of the α1–6Man linked terminus structure between a protein coordinated and an uncoordinated state, suggesting the large chemical shift deviations characterizing the low-temperature state of the glycan are due to interactions with the protein surface rather than changes in internal conformation of the glycan. Resonances at the 1, 2 and 3 positions of the galactose ring experienced the greatest deviation in chemical shifts (Supplementary Figure 1). This is consistent with the proximity of this portion of the carbohydrate to a surface cavity observed in high resolution Fc models (pdb 1l6x, 1dn2). The 4,5 and 6 positions, however, saw less chemical shift deviations and are located away from the cavity. Therefore, we interpret the minor, low temperature state (B) of the α1–6Man galactose residue as one bound to the surface of the Cγ2 domain in a manner quite consistent with crystal structures.

The precise origin of the chemical shift perturbations in state B are difficult to deduce. There are both residues carrying aromatic sidechains and charged side chains in the region of binding. The aromatic residues are known to restrict pan class="Chemical">glycan processing from mutational studies32,36, however, the shifts of the C2 and class="Chemical">pan class="Chemical">H2 resonances, which are in opposite directions are not consistent with remote group effects from the phenyl rings; these would be expected to produce shifts in the same direction. Electrostatic effects are more likely; an electric field from a properly positioned positively-charged residue, for example, would be expected to pull electron density onto the H2 proton and away from the C2 carbon producing shifts in the opposite direction. Lys246 is in a position where it could produce such effects and is part of the galactose binding pocket.

In accordance with the spin relaxation data and chemical shift estimates we propose a much more dynamic model of the pan class="Chemical">glycans on the Fc fclass="Chemical">pan class="Disease">ragment of IgG as depicted in Figure 6. State B is represented by the X-ray structure of the Fc fragment (pdb-1l6x) with the α1–3Man branch built out in an arbitrary way beyond the branching Man residue. The α1–6Man branch remains fixed in the bound state as seen in the crystal structure. State A (the free state) is represented by structure A. Here only the linkage of the glycan to the N of Asn297 is as in the X-ray structure. The conformation of both branches are extended but in different conformations in the two structures to indicate additional motion of both branches. The location of the dynamic α1–3Man terminus in states B and A, as well as the α1–6Man terminus in A, is purely speculative; however, a location proximal to the face of the Fc fragment as shown in Figure 6 would expose a larger portion of the glycan if the relative location of the two Cγ2 domains as observed from x-ray diffraction data is retained between all states.

Figure 6

Models for Fc glycan dynamics and accessibility showing exposed glycan conformations are possible

The terminal galactose residues of the glycan are labeled as “α1–3”Man-linked or “α1–6”Man-linked. These hypothetical models show accessible locations of the glycans in states A and B based on X-ray structures of the Fc fragment (pdb-1l6x). State B is characterized by the α1–6Man-linked branch coordinated by the surface of the polypeptide, and both branches of state A are highly dynamic.

These models are consistent with observations in pan class="Species">many respects. The class="Chemical">pan class="Chemical">glycan is shielded by the Fc polypeptide accounting for some reduction in rates of modification by glycan modifying enzymes, but still accessible, especially in state A. Glycan – polypeptide interactions, notably through Phe241 and Phe243, inhibit but do not prevent glycan processing32,36. In the models the glycan in the B state is shielded by these residues as highlighted in Figure 6 potentially accounting for the observations.

The labeling, assignment and pulse techniques reported here surmount the pan class="Species">many challenges facing protein-conjugated class="Chemical">pan class="Chemical">glycan studies and are applicable to a wide range of glycoprotein systems, though further methodological improvements will be required for more complex systems. These methods will be central to studying the association of the Fc fragment with cell-surface receptors and characterizing the role the glycans play in these interactions.

Methods

Materials

All enzymes and chemicals were purchased from Sigma-Aldrich (St. Louis, MO), unless otherwise noted. pan class="Chemical">13CU-galactose, class="Chemical">pan class="Chemical">13C2-galactose and deuterium oxide (D2O) were from Cambridge Isotopes (Cambridge, MA), the Fc fragment of IgG and IgG1 was from Athens Research and Technologies (Athens, GA).

Remodeling the IgG1 and IgG Fc glycans

The Fc fpan class="Disease">ragment and IgG1 class="Chemical">pan class="Chemical">glycans were remodeled with 13C-galactose and Nacetylneuraminic acid as previously described21. IgG1 proteolysis was performed using 2 mg / mL papain activated for 10 min at 37°C in a buffer containing 1 mM EDTA, 10 mM cysteine and 50 mM sodium phosphate, pH 7.0. This mixture was washed in a 10kDa-cutoff centrifugal filter (Millipore) to exchange the buffer with 25mM sodium phosphate, 100 mM potassium chloride, pH 7.0. Activated papain (0.1mg) was incubated with 10 mg IgG1 for 4 h at 37°C. Fc and Fab fragments were separated using a Protein A-sepharose column (GE Healthcare) to retain the Fc fragment. Following a complete digestion as judged by sodium dodecyl sulfate poly acrylamide gel electrophoresis, the Fab fragments were pooled from the column flow through plus wash fractions and the Fc fragment was recovered in fractions eluted with a buffer containing 0.1M glycine, pH 3.0.

Fc fpan class="Disease">ragment (10 mg) with class="Chemical">pan class="Chemical">13CU-galactose enrichment on the α1–6Man branch of the glycan was prepared by first removing existing galactose and N-acetylneuraminic acid residues with 80 units β-galactosidase and 250 units sialidase (New England Biolabs) in a 3 mL reaction containing 50 mM Tris, 200 mM sodium chloride and 0.05% sodium azide, pH 8.0. The Fc fragment was then purified over a Protein A-sepharose column, dialyzed against a buffer containing 10 mM MOPS, pH 7.2, then lyophilized. Galactosylation and sialylation were completed as described21 using 12CU sugar nucleotides. Following the four day sialylation reaction, the protein was exchanged into a buffer containing 50 mM sodium citrate, 100 mM sodium chloride and 80 units β-galactosidase, pH 6.0. After an overnight incubation at 37°C, Fc was again purified over a Protein A –sepharose column to remove the exoglycosidase. Following dialysis and lyophilization, Fc was incubated with β1–4galactosyl transferase and UDP-13CU-galactose for two days at 37°C. The completion of the reactions was confirmed by MALDI-MS analysis of the released glycans21.

NMR spectroscopy

NMR spectroscopy was performed on spectrometers opepan class="Species">rating at 21.1T, 18.8T or 14.0T. The 21.1T instrument was equipped with a Varian VNMRS console and 5 mm cryogenically-cooled probe with a cold class="Chemical">pan class="Chemical">carbon preamp suited for direct 13C observation. One 14.0T machine had a similar arrangement with a 3 mm cryogenically-cooled probe. The other 14.0T and the 18.8T instruments were equipped with a Varian Inova console and a 5 mm cryogenically-cooled probe. For each experiment, the temperature was calibrated with an ethylene glycol standard under conditions approximating data acquisition. Once the instrument temperature had equilibrated and the shimming had stabilized acquisition was initiated. DSS was included as an internal reference in each sample. Additionally, the DSS methyl linewidth was measured before and at the end of each experiment (< 1 Hz) to establish field heterogeneity contributions to glycan resonance linewidths.

NMR samples consisted of 2–10 mg of protein in 100% pan class="Chemical">D2O buffer containing 25 mM class="Chemical">pan class="Chemical">sodium phosphate, 100 mM potassium chloride and 1 mM DSS, pH 7.0. Samples were degassed using a vacuum pump until no bubbles were formed. CD spectroscopy and 1H NMR confirmed the samples were folded (data not shown).

A standard pan class="Chemical">13C heteronuclear single-quantum coherence (HSQC) experiment (Varian Bioclass="Chemical">pack) was modified to measure the spin relaxation class="Chemical">pan class="Species">rates of galactose 13C6 by specifically transferring magnetization to C6 and decoupling C5 from C6 during R2 measurement (see Supplementary Methods, Supplementary Figure 3). This is possible at high magnetic fields because the C5 and C6 chemical shifts are sufficiently disparate (~14 ppm) to be independently refocused using shaped pulses. C6 is a desirable observation site in systems with large R2s because the presence of two 1Hs cuts the delay required for transferring H-C antiphase magnetization into pure in-phase carbon magnetization (and vice-versa) in half. It does, however, carry complications in interpretation due to the possible presence of internal motions about the C5-C6 bond.

Ring pan class="Chemical">carbons do not carry the complication of additional bond rotations. Hence, C2 was also a target for spin relaxation measurements. However, the Fc α1–6class="Chemical">pan class="Species">Man galactose ring resonances are especially sensitive to pulse sequence length; in fact, they are not observed in a sensitivity-enhanced carbon HSQC experiment, owing to additional refocusing delays before acquisition when compared to a standard phase-cycled HSQC experiment (data not shown). For this reason, combined with the larger R2s, 13C2 measurements were made, primarily, with direct carbon excitation / observation experiments using the same relaxation elements as those shown in Supplementary Figure 3. CPMG-type relaxation dispersion experiments were 1D 13C2 excite/observe or 13C2 excite/1H observe using methods described27.

Relaxation dispersion curve fitting

Relaxation dispersion data were fit with: where R0 is the relaxation pan class="Species">rate in the absence of chemical exchange, k is the class="Chemical">pan class="Species">rate of exchange, τ is the delay between refocusing pulses and Δωis the frequency difference of states A and B in angular units25.

Chemical shift extraction from temperature dependencies

An equation to describe the chemical shift of a resonance in fast exchange at a given tempepan class="Species">rature in terms of enthalpy and entropy was derived by assuming two-state exchange between states A and B and tempeclass="Chemical">pan class="Species">rature independent entropy and enthalpy values as shown in the Supplementary Methods section: Where δX is the observed chemical shift, δA and δB are the chemical shifts of states A and B, respectively, R is the gas constant, H is enthalpy, S is entropy and T is temperature.

Calculation of pA using CPMG and chemical shift data

The population of the state A may be estimated using φ and δ: where ω is the Larmor frequency of the nucleus in MHz. A derivation of this equation is presented in Supplementary Methods.