Consequences of glycan truncation on Fc structural integrity

Effective characterization of protein-based therapeutic candidates such as monoclonal antibodies (mAbs) is important to facilitate their successful progression from early discovery and development stages to marketing approval. One challenge relevant to biopharmaceutical development is, understanding how the stability of a protein is affected by the presence of an attached oligosaccharide, termed a glycan. To explore the utility of molecular dynamics simulations as a complementary technique to currently available experimental methods, the Fc fragment was employed as a model system to improve our understanding of protein stabilization by glycan attachment. Long molecular dynamics simulations were performed on three Fc glycoform variants modeled using the crystal structure of a human IgG1 mAb. Two of these three glycoform variants have their glycan carbohydrates partially or completely removed. Structural differences among the glycoform variants during simulations suggest that glycan truncation and/or removal can cause quaternary structural deformation of the Fc as a result of the loss or disruption of a significant number of inter-glycan contacts that are not formed in the human IgG1 crystal structure, but do form during simulations described here. Glycan truncation/removal can also increase the tertiary structural deformation of CH2 domains, demonstrating the importance of specific carbohydrates toward stabilizing individual CH2 domains. At elevated temperatures, glycan truncation can also differentially affect structural deformation in locations (Helix-1 and Helix-2) that are far from the oligosaccharide attachment point. Deformation of these helices, which form part of the FcRn, could affect binding if these regions are unable to refold after temperature normalization. During elevated temperature simulations of the deglycosylated variant, CH2 domains collapsed onto CH3 domains. Observations from these glycan truncation/removal simulations have improved our understanding on how glycan composition can affect mAb stability.


Introduction
Proteins-based therapeutics are an important class of biopharmaceuticals that have substantially increased the number of treatable diseases. 1,2 Among the biopharmaceuticals currently in development, monoclonal antibodies (mAbs) make up the largest portion. Antibodies perform multiple functions-each antibody is composed of two Fab regions capable of binding antigen and one Fc region, which modulates antibody clearance and activation of antigen destruction through effector function pathways such as antibody-dependent cell-mediated cytotoxicity (ADCC) and complement-dependent cytotoxicity (CDC). 3 The Fc region consists of four structural domains, two C H 2 and two C H 3 domains (Fig. 1). In antibodies of the IgG class, each C H 2 domain in the Fc region has a complex biantennary oligosaccharide, termed a glycan (Fig. 2), attached at Asn297 (in Kabat numbering scheme 4 ). Due to the presence of different terminal sugars, the composition of these glycans can be Effective characterization of protein-based therapeutic candidates such as monoclonal antibodies (mAbs) is important to facilitate their successful progression from early discovery and development stages to marketing approval. One challenge relevant to biopharmaceutical development is, understanding how the stability of a protein is affected by the presence of an attached oligosaccharide, termed a glycan. To explore the utility of molecular dynamics simulations as a complementary technique to currently available experimental methods, the Fc fragment was employed as a model system to improve our understanding of protein stabilization by glycan attachment. Long molecular dynamics simulations were performed on three Fc glycoform variants modeled using the crystal structure of a human IgG1 mAb. Two of these three glycoform variants have their glycan carbohydrates partially or completely removed. Structural differences among the glycoform variants during simulations suggest that glycan truncation and/or removal can cause quaternary structural deformation of the Fc as a result of the loss or disruption of a significant number of inter-glycan contacts that are not formed in the human IgG1 crystal structure, but do form during simulations described here. Glycan truncation/removal can also increase the tertiary structural deformation of C H 2 domains, demonstrating the importance of specific carbohydrates toward stabilizing individual C H 2 domains. At elevated temperatures, glycan truncation can also differentially affect structural deformation in locations (Helix-1 and Helix-2) that are far from the oligosaccharide attachment point. Deformation of these helices, which form part of the FcRn, could affect binding if these regions are unable to refold after temperature normalization. During elevated temperature simulations of the deglycosylated variant, C H 2 domains collapsed onto C H 3 domains. Observations from these glycan truncation/removal simulations have improved our understanding on how glycan composition can affect mAb stability. Abbreviations: mAb, monoclonal antibody; RMSD, root mean square deviation; RMSF, root mean square fluctuation; Rg, radius of gyration; HpASA, hydrophobic accessible surface area; GAL, galactose; FUC, fucose, MAN, mannose; NAG, N-acetyl glucosamine one of several forms G2F, G0F, and G2FS2 (Fig. 2). Microheterogeneity in Fc glycan composition varies with species and expression systems used to produce candidate mAbs during research and commercial manufacturing. 5,6 Glycan composition is also known to influence the activation of specific effector function pathways differentially. 7,8 For example, truncating the G2F glycoform by degalactosylation to G0F can decrease CDC activation without affecting ADCC activation. 9,10 Effector function activity can also be modulated by defucosylation, which has been shown to increase ADCC activity. 11 These discoveries have inspired the development of different engineering approaches to enhance/optimize antibody effector function activity by altering the Fc-bound glycan profile or by exchanging amino acids in the protein backbone. 12,13 Glycan modification techniques such as these have led to the first market-approved glycoengineered monoclonal antibody, 14 and may be useful in optimizing effector function activities of therapeutic Fc-fusion proteins as well. 15 www.landesbioscience.com mAbs 905 REPORT

REPORT
In a separate study by Zheng et al., full deglycosylation was shown to expand the Fc hydrodynamic radius, increasing its susceptibility to proteolytic cleavage by papain. 20 Other studies have shown that partial and full deglycosylation can also increase the propensity of IgG1 Fc fragments to aggregate. 21 Therefore, understanding the mechanism that underlies the stability differences of fully and partially deglycosylated IgG Fc fragments is important to the development of therapeutic mAbs. 22 Effective characterization of protein-based therapeutics such as mAbs is important to facilitate their translation from discovery to commercially-available drugs. Achieving this goal partly depends on the sensitivity of analytical tools to provide chemical/physical information. Due to the size and complexity of fulllength mAbs, however, characterization of their higher order structure with currently available experimental methodologies is challenging. The tools available for characterization, including differential scanning calorimetry, analytical ultracentrifugation, circular dichroism, isothermal calorimetry, can provide a global view of mAb structure/stability, but cannot determine specific protein sequence or structural sites that can affect experimental observations. On the other hand, X-ray crystallography can provide detailed structural information, but structures are relatively static and do not provide dynamic information beyond crystallographic B-factors. 23 In principle, protein solution dynamics can be determined using nuclear magnetic resonance (NMR), but intact and IgG fragments are too large for routine NMR 3-D structural analyses. Moreover, NMR experiments require samples with highly concentrated protein, which can distort antibody structure and lead to aggregation. 24 Therefore, it is important for biopharmaceutical development to seek new analytical tools to detect protein structural changes with a high degree of sensitivity.
Atomistic and coarse-grained molecular dynamics (MD) simulations can play an important role in fulfilling this need. [25][26][27][28][29] MD can assess both the global and local conformational behavior of mAbs in solution under a variety of different conditions. Our recent analyses of a murine IgG2a using molecular dynamics provided new insights into mAb dynamics by comparing structural change differences upon thermal stress and deglycosylation. 30 Studies such as these and others illustrate how MD can be an important addition to the set of tools available for characterization of protein structural changes that may occur during the shipping and storage of biopharmaceuticals. 31,32 In this work, explicit water simulations were performed on the following Fc glycoform variants DEGLYCO, M3N2, and G2F, which were modeled using the crystal structure of a human mAb, PDB entry 1HZH (Table 1; Fig. 2). Each system was simulated at 300 K and 375 K for 200 ns in two independent production runs of 100 ns each. Total simulation time was 400 ns per Fc glycoform variant, and the overall total simulated time for this study was 1.2 μs. Simulations were evaluated for differences in domain ensemble averages, FcR (I-III) binding site closure, local structure perturbations, and sugar dynamics/ conformational preferences. Clustering glycan-glycan conformations in simulation trajectories revealed a preferred glycan-glycan conformation that is significantly different than the one in the crystal structure (1HZH). Moreover, the preferred glycan-glycan From a biopharmaceutical development perspective, truncating or removing glycans attached to therapeutic mAb candidates can also have undesirable consequences such as reduced Fc thermostability and increased aggregation. Consequently, it is important to ensure a consistent glycan profile across manufacturing batches or to match the glycan profile of an innovator product when developing a biosimilar version. 16 Therefore, much effort has been devoted to the development of rapid high through-put methods for glycan profile determination. 17,18 In a set of experiments demonstrating the effect glycan composition can have on Fc effector function activity, the thermostability of C H 2 domains decreased as more carbohydrates were removed from the glycan. 19 Figure 1. Fc structure (A) indicating the locations of structural domains, two attached oligosaccharides colored cyan, and the position of Gly237 residues used to calculate the C H 2-C H 2 distance. For chain H, the C H 2 domain is red and the C H 3 domain is orange. For chain K, the C H 2 domain is green and the C H 3 domain is blue. The hinge, which was present in all glycoform variant simulations, has been left out of this representation for the sake of clarity. C H 2 domain schematic (B) with β-strands labeled and residue indices for strands in Kabat numbering scheme. Two phenylalanine residues 241 and 243, important for C H 2 domain stability, are located on β-strand A. The oligosaccharide attachment residue, Asn297, is located in loop C'E.  19 Desialylation and degalactosylation, however, did not affect C H 2 domain thermostability to the same degree (C H 2 Tm for G0 is 71.5 °C). 19 Because NAG9 does not form interactions with C H 2 domains, the loss of C H 2 domain thermostability upon the removal of NAG6 (see above for M3N2) can be attributed to key C H 2 domain stabilizing interactions formed between NAG6 and Phe243 (Kabat) of β-strand A. 34,35 No further destabilization of C H 2 domains was observed upon removal of MAN5 or MAN8 from M3N2 to produce MN2 (NAG1-NAG2-MAN4), (C H 2 Tm for MN2 is 67.3 °C). 19 However, full deglycosylation did lower the C H 2 melting temperature below that of MN2 (C H 2 Tm for DEGLCYO is 65.9 °C). 19 This provides supporting evidence for the importance of interactions between MAN4 and Phe241. 34,35 The significance of the MAN4-Phe241 and NAG6-Phe243 interactions in stabilizing C H 2 domains was also supported by the work of Voynov et al. who showed that mutations, F241S and F243S, significantly lowered C H 2 domain melting temperatures for fully glycosylated glycoforms (G0, G1, and G2). 28 Here, individual domain ensemble averages are reported to understand the stabilizing effect that full length glycans have on C H 2 and C H 3 thermostability. Simulation evaluations began by measuring the degree of structural change using root mean square deviation (RMSD), radius of gyration (Rg), and hydrophobic accessible surface area (HpASA) for glycoform variants during equilibrated simulations. For each of these global measures, the magnitudes of structural change were then compared with the thermodynamic melting temperature differences among glycoform variants or other available experimental data. To calculate domain ensemble averages for all glycoform conformation has several novel inter-glycan contacts that are not formed in the crystal structure (1HZH). Consistent with the above observations, structural differences among the glycoform variants during simulations suggest that glycan truncation/ removal can cause quaternary structural deformation as a result of the loss or disruption of a significant number of these inter-glycan contacts. Glycan truncation can also increase the structural deformation of C H 2 domains, demonstrating the importance of specific carbohydrates in stabilizing C H 2 domains. In elevated temperature simulations, glycan truncation appears to affect structure deformation differentially in locations (Helix-1 residues 246-251 and Helix-2 residues 310-315) that are far from the oligosaccharide attachment point (Asn 297). These helices form part of the FcRn binding site and their preservation is important in maintaining pharmacokinetic and pharmacodynamic properties of therapeutic mAb candidates. Observations from these glycan truncation simulations have improved our understanding on how glycan composition can impact pharmaceutical stability of protein therapeutics and demonstrated the usefulness of molecular dynamics simulations in biopharmaceutical development studies.

Results
Individual domain ensemble averages. The Fc structure is described in Figure 1, and Figure 2 presents the location of carbohydrates in glycan structure and glycan composition Figure 2. Structure, composition, and branches are labeled for oligosaccharides of G2FS2. The glycan has three terminal carbohydrates: FUC-2, and two N-acetyl neuraminic acids or sialylic acids, labeled S. Desialylating G2FS2 produces G2F. Alternate glycan compositions and structures are shown using a simplified notation. Note that MN2, not shown, can be formed by removing MAN5 and MAN8 from M3N2. ) across all glycoforms. This was expected because the average HpASA was calculated in the absence of the oligosaccharide and there are two hydrophobic residues (Phe241, Phe243) that become solvent-exposed when ignoring the glycan. In comparing the average HpASA for C H 2 domains across glycoform variants at 300 K, values were approximately the same given that the range of HpASAs for C H 2 domains among glycoform variant was ~125 Å 2 , which is a small percentage of the total average C H 2 domain HpASA (~3500 Å 2 ). HpASA standard deviations for C H 2 domains of glycoform variants were even smaller than the range ( Table 2). At 375 K, the average HpASA for C H 2 domain (chain K) of DEGLYCO is less than it was at 300 K (3509.5 Å 2 , 300 K; 3170.6 Å 2 , 375 K). This was due to the collapse of Fc quaternary structure at elevated temperature (Fc ensemble averages are described in the next section). All other domains of both chains (H and K) for glycoforms variants had similar average HpASA values between 300 K and 375 K simulations (± 100 Å 2 ). The lack of differences in average HpASA values among glycoform variants suggests that glycan truncation/removal does not impact the average HpASA of individual structural domains.
Overall, analyses of Fc individual structural domains during simulations indicate that tertiary structure changes in C H 2 and C H 3 domains were relatively small among different glycoform variants at both 300 K and 375 K ( Table 2). Furthermore, C H 3 domain tertiary structure changes upon glycan truncation/removal were smaller than those in C H 2 domains. This is expected because glycosylation sites are located within C H 2 domains. Moreover, the destabilization of C H 2 domains upon glycan truncation (G2F vs. M3N2) and deglycosylation (G2F variants, the last 90 ns of both trajectories at the same temperature were used. Salient observations are described below. At  For example, the C H 2-C H 2 distance in the crystal structure of MN2 is 4.7 Å shorter than it is in the structure for G2F. 35 Full deglycosylation shortened the C H 2-C H 2 distance even further. Based on these observations, it was concluded that progressive glycan truncation reduces Fc activity by shortening the distance between C H 2 domains, appearing to close the FcR (I-III) binding site and decrease FcR binding affinity. To investigate the effect glycan truncation has on overall Fc quaternary structure/activity during simulations, Fc average quaternary structure changes were determined at both temperatures. FcR binding site closure was monitored by computing the distance between Gly237 residues of β-strand A (see Materials and Methods). At 300 K, the C H 2-C H 2 distance did not change appreciably in G2F and M3N2 simulations (Fig. 3A). For DEGLYCO, C H 2 domains moved closer together such that the average C H 2-C H 2 distance was ~4 Å less than the average distance for G2F ( Table 2). The mutual approach of C H 2 domains in deglycosylated Fc simulations is consistent with observations from crystallographic studies (described above). At an elevated temperature (375 K), the C H 2-C H 2 distance actually increased for DEGLYCO over time. The increase was not due to domain separation as disulfide bonds in the hinge prevent this from happening. Instead, C H 2 domains collapsed onto C H 3 domains by twisting relative to one another, extending the distance between Gly237 residues (Fig. 3B). For DEGLCYO at 375 K, the extent of C H 2 domain collapse onto C H 3 domains can also be seen by comparing average HpASA values for both chains. At 300 K, the average HpASAs for C H 2 and C H 3 domains of chain K were 3509.5 Å 2 and 2367.2 Å 2 , respectively. Collapse of C H 2 domains onto C H 3 domains of chain K decreased the average HpASA for both domains in 375 K simulations (C H 2, 3170.6 Å 2 ; C H 3, 23041.1 Å 2 ). In M3N2 simulations at 375 K, C H 2 domains partially collapsed onto C H 3 domains, but the C H 2-C H 2 distance was restored to the initial simulation distance of 15 Å after 40 ns. C H 2 domains began collapsing again after 70 ns suggesting that the collapse of C H 2 domains in M3N2 simulations at 375 K could be a transient event. In G2F simulations at 375 K, the C H 2-C H 2 distance remained relatively constant and similar to the initial simulation distance (15 Å).
RMSD ensemble averages for the complete Fc fragment ( Table 2) were similar to the changes in C H 2-C H 2 distance during simulations (Fig. 3A and B). At 300 K, the structure of DEGLYCO deformed by 4.65 Å compared with the crystal structure, and M3N2 and G2F deformed by 4.25Å and 3.74 Å, respectively. The same increasing trend of average RMSD among glycoform variants was even more apparent at 375 K (DEGLYCO, 8.09 Å; M3N2, 6.97 Å; G2F, 5.60 Å). The Fc ensemble average RMSD for DEGLYCO was 2.5 Å greater than the average RMSD for G2F in elevated temperature simulations ( Table 2). The greater Fc quaternary structural deformations in elevated temperature simulations compared with 300 K simulations was likely due to the collapse of C H 2 domains onto C H 3 domains. By comparing Fc average RMSD differences among glycoform variants to their C H 2-C H 2 distance changes during simulations, it appears that these two descriptors of structural change are interrelated, i.e., the movements of the individual C H 2 domains with vs. DEGLYCO) was consistent with thermodynamic measurements summarized at the beginning of this section. Qualitative recapitulation of experiment results provides confidence in subsequent simulation analyses.
The Fc average radius of gyration (Rg) during DEGLYCO simulations was slightly larger than it was during M3N2 and G2F simulations at 300 K (DEGLYCO, 26.68 Å; M3N2, 26.48 Å; G2F, 26.30 Å), ( Table 2). This observation is unexpected because C H 2-C H 2 approach would seemingly decrease the overall volume of the quaternary structure (see above), but it is in agreement with experimental observations of a slight increase in hydrodynamic radius and susceptibility to papain cleavage upon deglycosylation. 20 In simulations at 375 K, the Fc average Rg for DEGLYCO was less than it was for G2F and M3N2 due to the collapse of C H 2 domains onto C H 3 domains (DEGLYCO, 24.12 Å; M3N2, 25.75 Å, G2F, 25.55 Å), ( Table 2). The collapse of C H 2 domains onto C H 3 domains in DEGLCYO simulations also affected the Fc average HpASA, which was calculated in the absence of the glycan for all glycoform variants ( Table 2). While the Fc average HpASA was approximately the same for all glycoform variants at 300 K (DEGLYCO, 14034.9 Å 2 ; M3N2, 14125.8 Å 2 ; G2F, 14015.6 Å 2 ), at 375 K, the average HpASA decreased due to the reduction in overall accessible surface area, resulting from C H 2 domain collapse (DELGYCO, 13185.5 Å 2 ; M3N2, 14081.0 Å 2 ; G2F, 14126.8 Å 2 , see Table 2). Overall, the observations from Fc and structural domain ensemble averages indicate that glycan truncation and deglycosylation can affect both C H 2 domain tertiary structure and Fc quaternary structure. The differences in quaternary structure deformation among glycoform variants were, however, significantly larger than differences in tertiary structure deformation. This suggests that individual C H 2 domains move with respect of each other and with respect to C H 3 domains. Similar observations have also been made in earlier simulations of a full length murine IgG2a mAb. 30 The consequence of these inter-domain movements is that the effector function activity can become impaired.
Local structure deformations upon glycan truncation/ removal. Glycan truncation/removal can also perturb residues in loops that constitute the FcR binding site such as BC, FG, and C'E loops (Fig. 1B). Crystal structures of G2F, G0F, and M3N2F suggest that removal of NAG6 perturbs the C'E loop. 35 Variants with an attached NAG6 carbohydrate (e.g., G0F) do not show differences in C'E loop perturbations when compared with the crystal structure of G2F. NMR studies have supported these findings and added that glycan truncation can also perturb the lower hinge region as well as the FG loop after full deglycosylation. 33 To monitor the local structure changes upon glycan truncation, average residue-wise RMSD values were computed. Since C H 3 domains were mostly unaffected by glycan truncation ( Table 2), this section focuses on the local structure changes in C H 2 domains only. The average residue-wise RMSD values are reported for β-strands and connecting loops by averaging between C H 2 domains of heavy chains, H and K (see Materials and Methods). The average residue-wise RMSD values were compared for similarities and differences at 300 K and 375 K among glycoform variants to identify local structure that is susceptible to glycan truncation and removal. temperature (300 K vs. 375 K) and the extent of glycan truncation. Helix-1 (residues 246-251) and Helix-2 (residues 310-315), which form part of the FcRn binding site, are located far from the glycan attachment site (Asn 297). Deformation of these helices could affect FcRn binding if these regions are unable to refold after temperature normalization. FcRn binding is an important determinant of pharmacokinetic and pharmacodynamic properties of therapeutic mAb candidates. 36 Glycan dynamics. Enzymatic sialylation of galactosylated glycan termini was shown to be anti-inflammatory in mice. [37][38][39] From a glycoengineering standpoint, enzymatic sialylation of normally galactosylated IgGs has the potential to improve therapy for patients suffering from inflammatory diseases such as rheumatoid arthritis. While a recent report suggests that sialylation produces an anti-inflammatory effect by biasing the FcR binding site toward the "closed" conformation, 40 it remains unclear how galactosylated termini are accessible for sialylation given their location in the cavity formed by C H 2 domains. Crystal structures of galactosylated Fc fragments show that glycan termini are normally buried between C H 2 domains. In early solution NMR measurements of Fc glycoforms, the α1-6 branch was completely immobilized to the surface of C H 2 domains while the α1-3 branch appeared to be more dynamic (Fig. 2). 41,42 Immobilization would seemingly prevent the α1-6 branch from enzymatic sialylation or interactions with receptors. A more recent NMR study indicates that the termini of both glycan branches are highly dynamic. 43,44 Here, we investigated the dynamics of both glycan branches at 300 K and 375 K using Root Mean Square Fluctuation (RMSF), ( Table 3). RMSF is a measure of how much a carbohydrate fluctuates about its average position. It is, therefore, an appropriate measure of glycan dynamic behavior. RMSF values for each carbohydrate were averaged between H and K glycans and across trajectories (see Methods).
At 300 K, the most mobile carbohydrates in G2F as measured by average RMSF were the branch termini GAL7 (1.0 Å), GAL10 (2.0 Å), and FUC2 (1.0 Å), (Fig. 2 and 5A). In comparing carbohydrates of the α1-6 and α1-3 branches, NAG9 and MAN8 (α1-3) fluctuated more than NAG6 and MAN5 (α1-6). Therefore, the α1-6 branch was less mobile than the α1-3 branch in 300 K simulations, agreeing with earlier NMR studies. Moreover, upon visual inspection of simulation trajectories, it was found that the α1-6 branch never lost all protein-glycan interactions as was proposed to explain α1-6 sialylation. 43,44 The retention of contacts between α1-6 carbohydrates and C H 2 domains can also be seen in Table 3. The NAG6-F243 contact and all GAL7 contacts were retained during simulations ~100% of the time (NAG6, GAL7 are located on the α1-6). It is possible, however, that 100 ns of simulation time may not be enough to see GAL7 lose its contacts, despite a reported α1-6 branch correlation time of 37 ns. 44 Average RMSF values between M3N2 and G2F were similar except at MAN5 (α1-6), (Fig. 5). The truncated glycan, M3N2, had similar carbohydrate mobilities to G2F in the more mobile α1-3, but the α1-6 branch was more mobile in M3N2 than in G2F due to the loss of glycan-protein interactions. Thus, GAL7 and NAG6 glycan-protein interactions at an elevated temperature (375 K), the perturbation differences in loop FG among variants did not exist (see below). Therefore, local structure deformation differences among glycoform variants were mostly localized to loop C'E loop. Furthermore, the localization of local structure perturbation differences to loop C'E suggests that differences in C H 2 domain ensemble average RMSDs among glycoform variants ( Table 2) could be a result of loop C'E deformation differences.
At 375 K, local structure perturbations in C H 2 domains were greater than at 300 K, and glycan truncation generally increased the extent of perturbation for local structure regions (Fig. 4A  and B). Comparing local structure perturbations between 300 K and 375 K, the local regions most affected by an increase in temperature were Helix-1 and Helix-2. At 375 K, both these helices deformed 2.0 Å more than they did at 300 K, indicating that these regions were highly susceptible to deformation at elevated temperatures. For DEGLCYO, the local regions most perturbed from the crystal structure were in Helix-1 (4.5 Å), loop BC (3.5 Å), β-strand C (3.0 Å), β-strand C' (3.0 Å), loop C'E (3.5 Å), Helix-2 (4.0 Å), and loop FG (2.5 Å). Comparing average residue-wise RMSDs among glycoform variants at 375 K, the regions of local structure with similar perturbations were in loops BC and FG loop, and all β-strands (Fig. 4A). Only β-strand C' differed among glycoform variants as it was heavily perturbed in one of the M3N2 simulations at 375 K. Local regions with differences in perturbation among glycoform variants were in Helix-1 (DEGLYCO, 4.0 Å; M3N2, 3.0 Å; G2F, 2.5 Å), loop AB (DEGLYCO, 4.5 Å; M3N2 and G2F, 3.0 Å), Helix-2 (DEGLYCO, 4.0 Å; M3N2, 3.0 Å; G2F, 2.5 Å), and loop C'E (DEGLYCO and M3N2, 3.5 Å; G2F, 2.5Å). Based on these differences, the extent of perturbation for these local structure regions generally increased due to glycan truncation/ removal. Taken together, observations from 300 K and 375 K simulations suggest that the amount of local structure perturbation within C H 2 domains depended both on the simulation NAG9, GAL10) facing in an outward direction. Overall, the inter-glycan contacts in the preferred conformation were shifted compared with those in the crystal (Fig. 6), and crystal contacts MAN4(H)-MAN8(K) and MAN5(H)-MAN8(K) were not well maintained ( Table 3). The critical stabilizing C H 2 domain contacts NAG6-F243, and MAN4-F241 were ~100% retained for both glycans. Comparatively, fucose-C'E loop contacts (FUC2-Y296 [H], FUC2-Q295 [K]) were partially lost during simulations, indicating much weaker interactions. Taken together, the preferred oligosaccharide conformation had significantly more inter-glycan contacts than observed in crystal structure (1HZH) and inter-glycan contacts were crucial to maintain Fc quaternary structure and/or activity. These observations may are important in restricting α1-6 branch carbohydrate dynamics. As in 300 K simulations, the most mobile carbohydrates at 375 K were the branch termini GAL7 (3.0 Å), GAL10 (3.0 Å), and FUC2 (2.5 Å). A comparison G2F carbohydrate mobilities in 300 K and 375 K simulations shows that there was an average 1.5 Å increase in average RMSF among carbohydrate moieties. At 375 K, MAN5 and NAG6 (α1-6) were more mobile than MAN8 and NAG9 (α1-3), (Fig. 5). Upon visual inspection of G2F simulation trajectories at 375 K, GAL7 did not leave the protein surface, but the greater mobility of the α1-6 branch than that of the α1-3 branch, at an elevated temperature, does suggest that this branch could expose its termini for sialylation over time.
Carbohydrate conformation/contacts. The crystal structure complex of FcRIII bound to the Fc fragment has revealed the existence of a single contact between oligosaccharides and FcRIII. 45 Because FcRIII activation is dependent on Fc glycan composition 19 (see section Fc quaternary structure ensemble averages), oligosaccharide moieties exert their influence indirectly through modification of the Fc conformation. As previously stated, glycan truncation can reduce FcR (I-III) binding by decreasing the distance between C H 2 domain, and thereby causing the FcR binding site to close. Therefore, contacts between oligosaccharides from glycans attached to different heavy chains (inter-glycan contacts) are crucial to maintain the Fc in an active state. Fc fragment crystal structures with intact glycans, however, show very few inter-glycan contacts with MAN8 mediating most of the inter-glycan contacts. 34,46,47 Moreover, Fc crystal structures obtained under different crystallization conditions have different inter-glycan interactions. 35 Thus, there is a need to investigate inter-glycan interactions and their dynamics to understand their effect on Fc conformation.
Here, the inter-glycan contacts from the crystal structure (1HZH) are compared with the preferred glycan-glycan conformation found during G2F simulations at 300 K ( Fig. 6A and B). The preferred glycan (chain H) -glycan (chain K) conformation, labeled 8240, was identified by clustering all G2F trajectory samples from both 100 ns simulations (see Methods). In total, 85% of all structures submitted for clustering had the preferred glycan-glycan conformation. The crystal structure (1HZH) with modeled GAL10(K) and FUC(H) showed only three inter-glycan contacts MAN4(H)-MAN8(K), MAN5(H)-MAN8(K), MAN8(H)-NAG3(K), (Fig. 6A; Table 3). The glycan-glycan conformation has both α1-3 branches facing down and away from the FcR binding site, and both α1-6 branch termini have GAL7 bound to the protein surface. 48 In the preferred glycanglycan conformation from simulations, both α1-3 branches were extended and faced outwards (Fig. 6B). The outward position of the α1-3 branch may allow sialylation. Compared with the crystal structure, several new inter-glycan contacts were observed in the preferred glycan-glycan structure. These are   57,58 Gray outlines show the surface of carbohydrates. Yellow dotted lines indicate inter-glycan contacts (identified at ≤ 4.5 Å for non-hydrogen atoms) between chains H and K in the crystal structure. Red dotted lines indicate inter-glycan contacts in the preferred glycanglycan conformation that were formed during simulations at 300 K and are not formed in the crystal structure for 1HZH.
www.landesbioscience.com mAbs 913 to mimic the thermal stress experienced by biopharmaceuticals as a result of occasional temperature incursions that may happen during manufacture, shipping and storage. These simulations are analogous to accelerated stability studies, which are routinely performed in pharmaceutical development. Importantly, the temperature values in molecular dynamics simulations should be carefully interpreted. In simulations, temperature is simply a measure of the amount of kinetic energy in the system. As such, elevating the temperature is a commonly used strategy to allow molecular systems to escape from deep energy minima. 49 However, elevating the simulation temperature is not a reliable strategy to reproduce experimental temperature-dependent phenomena such as changes in protein-protein binding affinity. On the other hand, the molecular origins that underpin protein conformational changes upon thermal stress can be reliably accelerated by arbitrarily increasing the temperature. [50][51][52] The molecular origins of these changes remain similar when simulations are performed at different elevated temperatures, but the time interval have important implications not only for basic understanding of Fc structure-function relationships, but also toward rational design of glycoengineered antibodies.

Discussion
In addition to simulations performed at 300 K, this study reports findings for molecular dynamics performed at an elevated temperature (375 K) used to enhance conformational sampling and perturb glycoform variants during simulations. In these elevated temperature simulations, the goal was not to unfold the IgG Fc molecular system. Protein folding/unfolding transitions occur over the millisecond to second timescales. One hundred nanoseconds of all-atom explicit solvent simulation time cannot achieve this for a molecular system as large as the Fc fragment, even when such studies are performed at an elevated temperature. Instead, the elevated temperature simulations in this study were performed to perturb the Fc molecular system as a means scheme. The inter-heavy chain disulfide bond formed by Cys229 residues, which is not formed in the 1HZH crystal structure, was modeled in using MOE. 53 Glycans for each glycoform variant were identical (M3N2 and G2F). GAL10 of chain K and FUC2 of chain H (Fig. 2), which are absent in the 1HZH crystal structure, were modeled for the G2F variant using MOE, with identical glycosidic dihedral angles to their corresponding carbohydrate, which is present in the opposite glycan. Histidines were doubly protonated, simulating a pH near or below 6, which is commonly used in formulation development. Each Fc system was solvated in a TIP3P water box with an initial box size of ~786 348 Å 3 for all glycoform variants. Chloride ions were added to neutralize all systems, resulting in average Clion concentration of ~16 mM per bath of water molecules. Energy minimizations and molecular dynamics simulations were performed using the NAMD package 54 with the AMBER99 all-atom forcefield 55 and the GLYCAM06 forcefield for carbohydrates. 56 Molecular systems were first subjected to 20 000 energy minimization steps with protein heavy atoms fixed, followed by 5000 steps of energy minimization without any restraints. Production runs were performed at 300 K and 375 K and each system was simulated for 100 ns. Trajectory samples were taken every 10 ps. For all glycoform variants, DEGLYCO, M3N2, G2F, two production simulations were performed at each temperature (Table 1). Overall, 12 production runs were performed during this study for a total simulation time of 1.2 μs.
Ensemble averages were calculated using the last 90 ns of each production run. The first 10 ns were used as an equilibration phase. System equilibration was assessed by monitoring energy, temperature, and pressure profiles. Domain ensemble averages were computed without the hinge or linker between C H 2-C H 3 domains. Reported ensemble averages are values averaged over simulation time (trajectory samples taken 10 ps apart) and between the two production runs at each temperature. RMSD and Rg values were calculated using all heavy atoms. HpASA was calculated using all carbon atoms and carbon bonded hydrogens in the absence of the glycan. The C H 2-C H 2 distance was computed between Gly237 (β-strand A) pairs, which is near the FcR (I-III) binding site (Fig. 1A).
For all trajectory samples, C H 2 structural domains were superimposed separately onto their corresponding crystal structure coordinates (chain H onto H, chain K onto K) and residue-wise RMSD values were computed using all atoms. Superimposing individual structural domains limits residue-wise RMSD values to the tertiary structure changes and excludes the effect of quaternary structural changes. Residue-wise RMSDs were then averaged across all trajectory structure samples (10 ps apart). Residue-wise RMSDs were also averaged between C H 2 domains (H and K) and across the two production runs at the same temperature (300 K and 375 K). RMSF values for glycans were computed by first superimposing glycans individually (chain H onto H, chain K onto K) onto the glycan coordinates found in the crystal structure (1HZH). RMSF values for glycan carbohydrates were also averaged between chains and production runs. necessary to observe these conformational changes is different due to varying the kinetic energy.
In the presented simulation analyses, glycoform variants were evaluated for differences in domain ensemble averages, FcR (I-III) binding site closure, local structure perturbations, and sugar dynamics/conformational preferences. Upon glycan truncation/removal, we find that simulations can recapitulate the thermodynamic trends for C H 2 and C H 3 domains. The closure of the FcR binding site was dependent on the degree of glycan truncation. During simulations at both 300 K and 375 K, most local structure deformation differences among glycoform variants were confined to the C'E loop, which contains the oligosaccharide attachment point, Asn297. In simulations at elevated temperature, Helix-1 and Helix-2 were the most perturbed local structures and were more perturbed for the deglycosylated glycoform variant than for the other variants. Within the glycan, the oligosaccharides of the α1-3 branch were more mobile than those of the α1-6 at room temperature (300K); however, at an elevated temperature (375K), the α1-6 branch was more mobile than the α1-3 branch. Furthermore, the preferred glycan-glycan conformation, identified by clustering glycan-glycan conformations in simulation trajectories, had significantly more inter-glycan contacts than those found in the crystal structure.
Observations of structural differences among glycoform variants suggest that glycan truncation and removal can cause quaternary structural deformation of the Fc region by C H 2 domain closure which results from the loss (DEGLYCO) or disruption (M3N2) of a significant number of inter-glycan contacts revealed by the preferred glycan-glycan conformation observed in simulations. These inter-glycan contacts are not formed in the crystal structure, PDB entry 1HZH. Glycan truncation also destabilized C H 2 domains by deforming the C'E loop at 300 K, suggesting that the absence of glycans can increase protein entropy, leading to domain destabilization. Glycan truncation also appears to affect structure deformation in locations (Helix-1 and Helix-2) that are far from the oligosaccharide attachment point in simulations at elevated temperature. Deformation of these local structures could affect FcRn binding if these regions are unable to refold after temperature normalization. Observations from these glycan truncation simulations have improved our understanding of how glycan composition can affect mAb stability and demonstrated the usefulness of molecular dynamics simulations in biopharmaceutical development studies.

Materials and Methods
Fc coordinates were taken from Protein Data Bank (PDB) entry 1HZH (heavy chains labeled H and K) that contains atomic coordinates for the full length human IgG1 b12 antibody. 48 The N and C termini were acetylated and amidated, respectively, to neutralize their charges. All inter-chain and intra-domain cysteine (Cys) pairs were disulfide bonded. Coordinates for the hinge (225-TCPPCPAPELLGG-237, Kabat numbering scheme 4 ) were included in all Fc glycoform variant simulations. All residues noted here are numbered according to the Kabat www.landesbioscience.com mAbs 915

Disclosure of Potential Conflicts of Interest
No potential conflict of interest was disclosed.