Biophysical Analyses for Probing Glycan-Protein Interactions.

Glycan-protein interactions occur at many physiological events, and the analyses are of considerable importance for understanding glycan-dependent mechanisms. Biophysical approaches including 3D structural analysis are essential for revealing glycan-protein interactions at the atomic level. The inherent diversity of glycans suits them to function as identification tags, e.g., distinguish self from the nonself components of pathogens. However, the complexity of glycans and poor affinities for interaction partners limit the usefulness of conventional analyses. To cope with such troublesome glycans, a logical sequence of biophysical analyses need to be developed. In this chapter, we introduce a workflow of glycan-protein interaction analysis consisting of six steps: preparation of lectin and glycan, screening of glycan ligand, determination of binding epitope, quantitative interaction analysis, 3D structural analysis, and molecular dynamics simulation. Our increasing knowledge and understanding of lectin-glycan interactions will hopefully lead to the design of glyco-based medicines and vaccines.

Introduction

pan class="Chemical">Glycan-protein interactions are crucial for numerous biological processes such as quality control of nascent proteins, cell-cell interaction, and pathogenic recognition. Of special interest is the breadth of n>n class="Chemical">glycan-mediated recognition events, which comes from the extremely high information-coding capacity of glycans. The diversity of oligosaccharides is inherent in their configurational variability, created by the broad range of monosaccharide species, anomeric configurations, multiple linkages, and branching. Glycan complexity is potentially suitable for coding biological information, and interest is growing to “decode” the sugar codes. To this end, investigating glycan-protein interaction at the atomic level is a necessity as only then will we understand the detailed mechanisms of these glycan-mediated biological activities. Recent techniques introduced include glycan microarray analysis and other high-throughput screening methods, and now we can begin to have a more systematic approach to analyzing glycan-protein interactions. We here introduce biophysical approaches to analyze the structure and function of glycan-protein interactions.

A Flowchart for Analyzing Glycan-Protein Interactions by Biophysical Methods

Biophysical investigation of pan class="Chemical">glycan-protein interactions is typically divided into six steps in our laboratory:

Large-scale preparation of protein and pan class="Chemical">glycan

Screening of pan class="Chemical">glycan ligand

Determination of binding epitope (pan class="Chemical">glycotope)

Quantitative pan class="Chemical">glycan-protein interaction analysis

3D structural study on proteipan class="Chemical">n-glycan interaction

Molecular dynamics simulation of pan class="Chemical">glycan-protein interaction

Crystallographic and nuclear magnetic resonance (NMR) analyses are direct methods for obtaining atomic details of pan class="Chemical">glycan-protein interactions if the binding epitope is known. In those cases where the binding epitopes are unknown, biophysical approaches such as listed above are required for structural studies. A computational approach combined with experiments is especially effective to investigate the flexible nature of glycans.

Large-Scale Preparation of Protein and Glycan

Biophysical studies require large amounts of protein and pan class="Chemical">glycan. In general, much effort needs to be paid for their prepan>ration. Examples of our achievements are shown below.

Preparation of Protein

pan class="Chemical">Carbohydrate recognition domain (lectin domain) or n>n class="Chemical">glyco-enzymes directly bind glycan, and it is convenient to use the isolated domain for biophysical interaction analysis. Normally, 5–10 mg of purified protein is required. Bacterial expression systems are frequently used to express lectin domains in recombinant form (Fig. 7.1a), but obtaining large amounts of pure lectin in a soluble form is often difficult using bacterial expression systems. For instance, C-type lectin and C-type lectin-like domains contain three conserved disulfide bridges, and the cysteine residues hamper soluble expression with correct folding. Refolding of protein from inclusion bodies may be needed to produce functional lectin domains. The protein is solubilized and refolded under controlled conditions. This refolding procedure has been effective for the purification of mammalian C-type lectins and C-type lectin-like domains (Nagae et al. 2013b, 2014a, b, 2016b).

Fig. 7.1

Large-scale preparation of recombinant proteins

(a) A flowchart of large-scale preparation of lectins for 3D structure and interaction studies

(b) Small-scale expression of differently fused MLS128-scFv in BL21 CodonPlus (left) and Origami B (right) strains analyzed by SDS-PAGE. M: Molecular weight of marker proteins; Lanes 1–8, odd number lanes are composition of total cell lysates, and even number lanes are composition of soluble cell lysates. Lanes 1 and 2 for scFv without fusion tag, Lanes 3 and 4 for thioredoxin (Trx)-fused scFv, Lanes 5 and 6 for DsbC-fused scFv, and Lanes 7 and 8 for PDI-fused scFv. The figure was adapted from a previous paper (Subedi et al. 2012) with permission

(c) Comparison of expression yield of Dectin-1 with and without GB1 fusion tag before and after concentration (to 1 mL) of purified sample. Results are presented as mean ± SD (n = 3) for pET28a murine Dectin-1, pCold murine Dectin-1, pCold-human Dectin-1 (codon-optimized), pCold murine Dectin-1 (optimized), and pCold GB1-fused murine Dectin-1 (codon-optimized). The figure was adapted from the paper (Dulal et al. 2016) with permission

(a) A flowcpan class="Chemical">hart of large-scale prepan>ration of lectins for 3D structure and interaction studies

(b) Small-scale expression of differently fused pan class="Chemical">MLS128-n>n class="Gene">scFv in BL21 CodonPlus (left) and Origami B (right) strains analyzed by SDS-PAGE. M: Molecular weight of marker proteins; Lanes 1–8, odd number lanes are composition of total cell lysates, and even number lanes are composition of soluble cell lysates. Lanes 1 and 2 for scFv without fusion tag, Lanes 3 and 4 for thioredoxin (Trx)-fused scFv, Lanes 5 and 6 for DsbC-fused scFv, and Lanes 7 and 8 for PDI-fused scFv. The figure was adapted from a previous paper (Subedi et al. 2012) with permission

(c) Comparison of expression yield of pan class="Gene">Dectin-1 with and without n>n class="Gene">GB1 fusion tag before and after concentration (to 1 mL) of purified sample. Results are presented as mean ± SD (n = 3) for pET28a murine Dectin-1, pCold murine Dectin-1, pCold-human Dectin-1 (codon-optimized), pCold murine Dectin-1 (optimized), and pCold GB1-fused murine Dectin-1 (codon-optimized). The figure was adapted from the paper (Dulal et al. 2016) with permission

Fusion of a tag protein can improve the solubility of target proteins. In our laboratory, several fusion tags pan class="Chemical">have been used, including n>n class="Gene">protein disulfide isomerase (PDI) and B1 domain of streptococcal protein G (GB1). Single chain Fv (scFv) originating from anti-carbohydrate MLS128 antibody was successfully expressed as a soluble form with PDI fusion (Subedi et al. 2012), while the other fusion tags such as thioredoxin or disulfide isomerase DsbC were not sufficient to produce soluble scFv (Fig. 7.1b). The expressed MLS128-scFv retains full binding activity toward synthetic carbohydrate antigens. The lectin domain of a C-type lectin-like receptor Dectin-1 was expressed as a GB1-fused protein (Dulal et al. 2016). The domain is correctly folded, is a monomer, and specifically binds a β-glucan ligand. The addition of the GB1 tag to the lectin domain greatly increased both yield (fivefold) and solubility (14-fold) (Fig. 7.1c).

Bacterial expression systems are the first pan class="CellLine">choice for protein expression but are not always successful for n>n class="Chemical">glycosylated proteins. In these cases, mammalian expression systems may be necessary to express correctly folded proteins with posttranslational modifications. HEK293 cells, CHO cells, and their variant cells are frequently used as hosts. For example, protein O-mannosyl kinase (POMK), which phosphorylates the mannose residue on α-dystroglycan, has three N-glycosylation sites within the catalytic domain, and two sites out of three are partially N-glycosylated. To abolish the heterogeneity originating from partial N-glycan occupancy, two asparagine residues were mutated to glutamine and the mutated POMK expressed in HEK293 variant cells lacking N-acetylglucosaminyl transferase I (GnT-I), to produce uniformly glycosylated protein (Man5GlcNAc2-glycan attached), which was successfully crystallized and the atomic structure elucidated (Nagae et al. 2017b).

Preparation of Glycan

It is another cpan class="Chemical">hallenge to obtain sufficient amounts of glycan and glycoconjugates suitable for biophysical assays. Prior to interaction analysis, it is necessary to design and prepared a suitable glycan ligand. Ligand design is an extremely important step. The design may arise from information in the literature, or else the binding specificity needs to be ascertained experimentally (Steps 2 and 3). Several methods can be used to prepare glycan ligands: (1) purification of a free natural glycan or liberation of a specific glycan from natural glycoconjugates ( glycoprotein, glycolipid, etc.), (2) chemical synthesis, and (3) chemoenzymatic synthesis. Liberation of glycan from natural glycoprotein can yield large amounts of glycan provided that there is a plentiful supply of the natural glycoprotein. N-glycans attached onto glycoproteins can be liberated by chemical methods (e.g., hydrazinolysis) or by enzymatic methods (e.g., PNGase F) and then labeled with a fluorescent tag 2-aminopyridine for detection. Pyridylamino derivatives of oligosaccharides can be purified by HPLC, and their chemical structures confirmed by solution NMR spectroscopy and mass spectrometry. For example, a biantennary N-glycan bearing bisecting GlcNAc was successfully purified with a yield of 120 μg from 200 mg of IgG and used for the analysis of a plant lectin E4-PHA (Nagae et al. 2014c, 2016b).

pan class="Chemical">Glycan ligands “designed” through synthetic approaches are of great importance in structural studies. n>n class="Chemical">N-glycan core units in particular cannot be prepared easily by enzymatic approaches, and chemical synthesis plays a significant role (Hanashima et al. 2014b). Synthetic bisected hexasaccharide has been used as ligand for the C-type lectin receptor murine DCIR2, the legume lectin PHA-E, and the plant lectin Calsepa. The use of designed synthetic ligands has been invaluable for revealing the ligand recognition mode of these lectins by X-ray crystallography and NMR.

Recombinant and/or chemoenzymatic methods pan class="Chemical">have been developed to obtain sufficient amounts of n>n class="Chemical">glycopeptides for structural analysis. It is known that several lectins recognize not only the glycan part but also the aglycon part (Nagae and Yamaguchi 2015), and here it is best to use the appropriate glycoconjugate (e.g., glycopeptide, glycolipid, etc.) as the ligand molecule. A C-type lectin receptor CLEC-2 binds to O-glycosylated podoplanin, and binding requires both disialyl core 1 glycan (NeuAcα2-3Galβ1–3[NeuAcα2–6] GalNAc) and neighboring amino acid residues. The O-glycosylated podoplanin glycopeptide was prepared by overexpression using engineered yeast cells and in vitro sialylation (Kato et al. 2008). The substrate of POMK was chemoenzymatically synthetized, composed of core M3 trisaccharide unit ( GalNAcβ1–3GlcNAcβ1-4Man) and α-dystroglycan peptide. Mono-mannosylated peptide was chemically synthesized in advance, and then additional disaccharide units were successively introduced by two glycosyltransferases (Nagae et al. 2017b).

Screening of Glycan Ligand

It is very important to understand how pan class="Chemical">glycans function in recognition and signaling systems within an organism and the responses with microbes and pathogens (Cummings and Pierce 2014). Heterogeneous n>n class="Chemical">glycans are potentially a vast source of information, and glycan microarrays embedding such diverse glycans are now becoming one of the key technologies to provide information on recognition determinants (glycotopes) of carbohydrate-binding proteins (Feizi and Chai 2004; Geissner and Seeberger 2016).

Glycan Microarray

When doing the pan class="Chemical">glycan microarray analysis, it is essential to immobilize a series of n>n class="Chemical">glycans onto the plate. The immobilization is one of the critical issues in the development of a glycan microarray. There are several ways to immobilize carbohydrate probes onto a suitable surface with retention of function. The chemically modified or unmodified glycans can be non-covalently or covalently attached onto the surface (Fig. 7.2a). Neoglycolipid (NGL) is an artificial glycolipid used to construct glycan microarrays via non-covalent interaction. In NGL technology, lipid-conjugated oligosaccharides prepared from natural sources or chemically synthesized are immobilized onto nitrocellulose-coated microarray slides (Palma et al. 2014). In covalent immobilization, various chemical reactions have been developed (Hyun et al. 2017). Among them, the conjugation reaction using amine-linked glycans is most widely used, and the Consortium for Functional Glycomics uses this chemistry for a glycan microarray platform (Fig. 7.2b). Several different functional groups which react with amines are utilized such as cyanuric chloride, N-hydroxysuccinimide ester (NHS), and epoxide. Measuring lectin binding requires labeling of the lectin with a fluorescent probe or the use of an antibody. Even so carbohydrate-lectin interaction is often too weak to be detected. In such cases, oligomerization of the lectin is a strategy to enhance inherent weak binding. Interactions may be detected by introducing dimeric glutathione S-transferase (GST) into the lectin as a fusion tag, as in the case of a mammalian lectin ZG16p-binding immobilized glycans (Kanagawa et al. 2014; Hanashima et al. 2015). Alternatively, it is possible to enhance binding by increasing the amount of surface glycans (Hanashima et al. 2015). However, it should be noted that the density of immobilized glycans may affect the specificity of lectins (Horan et al. 1999). Instead of fluorescence detection, label-free glycan biosensors have been developed to push the limit of detection (Hushegyi and Tkac 2014).

Fig. 7.2

Glycan microarray analysis for the lectin specificity

(a) Non-covalent and covalent glycan immobilization for glycan microarray analysis

(b) On-array inhibition assays for analyzing the binding of GST-fused human ZG16p to selected ligands which are non-covalently immobilized on the plate as neoglycolipids. Man-α-OMe, Gal-β-OMe (upper panel), oligosaccharide fractions (>20-mer) of heparin (Hep), and hyaluronic acid (HA) (lower panel) were used as inhibitors. This figure was adapted from a previous paper (Kanagawa et al. 2014) with permission

pan class="Chemical">Glycan microarray analysis for the lectin specificity

(a) Nopan class="Chemical">n-covalent and covalent n>n class="Chemical">glycan immobilization for glycan microarray analysis

(b) Opan class="Chemical">n-array inhibition assays for analyzing the binding of n>n class="Gene">GST-fused human ZG16p to selected ligands which are non-covalently immobilized on the plate as neoglycolipids. Man-α-OMe, Gal-β-OMe (upper panel), oligosaccharide fractions (>20-mer) of heparin (Hep), and hyaluronic acid (HA) (lower panel) were used as inhibitors. This figure was adapted from a previous paper (Kanagawa et al. 2014) with permission

pan class="Chemical">Glycan microarrays have various applications. Microarrays usually simply display relative binding affinity, but quantitative parameters (e.g., dissociation constant, K D) can be obtained by applying various concentrations (Liang et al. 2007). The estimated K D values are comparable to those obtained by other methods such as surface plasmon resonance (SPR). Moreover, effects of heterogeneous glycan mixing can be evaluated by immobilizing a mixture of glycans onto one spot (Liang et al. 2011). On-array inhibition assays can show if two different ligands bind competitively to the lectin (sharing a common binding site) or independently (each using a different binding site) (Fig. 7.2b) (Kanagawa et al. 2014). Glycan microarrays can also evaluate direct binding of mammalian cells and viruses (Nimrichter et al. 2004; Song et al. 2011), as well as biologically active glycans, if cell signaling can be linked to a fluorescent probe (Pai et al. 2016). This could result in rapid progress toward elucidating pathway-specific glycans.

Other Screening Methods to Evaluate Lectin-Glycan Interaction

pan class="Chemical">Glycan microarray analysis often may overestimate the strength of weak binding, and weak binding is often seen in n>n class="Chemical">carbohydrate-lectin interactions. In addition, non-specific binding is usually present, which contributes a background signal. Therefore it is preferable to include other screening methods to validate the results. Several methods have already been developed to screen glycans which bind target lectins. The thermal shift assay is widely used for ligand screening. In this approach, melting temperature is measured in the absence and the presence of ligands and any shifts assessed (Vedadi et al. 2006). A temperature shift is indicative of increased stability acquired upon ligand binding. The purified protein is mixed with ligand in the presence of a fluorescent probe, such as SYPRO Orange, that binds to the denatured state. This method is simple and can be used for rapid ligand screening. Frontal affinity chromatography (FAC) can quantitatively evaluate weak interactions between lectins and glycans, as explained in the following section. The FAC method can also be used to screen for ligand binding. NMR is an excellent technique to determine the binding epitope of a lectin-glycan interaction. When the binding epitope is determined, a suitable ligand of minimum size is designed for use in X-ray crystallography. NMR can also be applied to high-throughput ligand screening (Meyer and Peters 2003) and be advantageously used with a mixture of binder and non-binders.

Determination of Binding Epitope

The binding epitope of a pan class="Chemical">glycan (n>n class="Chemical">glycotope) can be deduced from the results of glycan microarray and FAC analysis using a series of glycans. More directly, NMR spectroscopy can provide the glycotope at atomic level. NMR spectroscopy has demonstrated its potential with a wide range of applications for studies of structure and function of biological molecules. The unique potential of NMR spectroscopy is its ability to generate high-resolution data not only on structure but also on dynamics and binding interactions. There are many types of NMR experiments that can be applied for very different experimental conditions both for the sample and the NMR parameters. NMR experiments provide various parameters such as chemical shifts, spin-spin couplings (J couplings), nuclear Overhauser effect (NOEs), and relaxation times. These parameters provide the atomic information on chemical environment, dihedral angles, distance, and mobility, respectively. Here two NMR techniques, saturation transfer difference NMR (STD-NMR) and deuterium-induced isotope shift (DIS), are described below.

Saturation Transfer Difference NMR in Lectin–Glycan Interactions

Saturation transfer difference NMR (STD-NMR) spectroscopy is a method for studying protein class="Chemical">n-ligand interactions by observing signals emanating from the ligand (Mayer and Meyer 2001; n>n class="Chemical">Haselhorst et al. 2009). When the target protein is saturated with selective pulses, the saturation is transferred to a bound ligand through intermolecular 1H-1H cross relaxation (Fig. 7.3a). Protons which are spatially close to the protein surface receive a higher degree of saturation compared with protons further away, and those on a non-binding ligand do not receive any saturation from the protein and are therefore not attenuated in the on-resonance spectrum. The subtraction of two spectra, one with protein resonance saturation (on-resonance) and the other without saturation of the protein resonance (off-resonance), results in a final STD-NMR spectrum showing only signals from bound ligands. For binding epitope analyses, the signal intensities of each peak are evaluated from an amplification factor (AF) = (I off – I on)/I off, where I on is the on-resonance signal intensity and I off is off-resonance signal intensity. Association constants can be estimated from a plot of AF against ligand concentration or initial growth rates of STD AF (Hanashima et al. 2010; Angulo et al. 2010). STD-NMR requires only a small amount of native protein, and no expensive protein labeling is necessary. STD-NMR spectroscopy can be utilized to detect binding of ligands with K D values in the 10−2–10−8 M range. Often screening for optimum STD-NMR conditions is required, such as probe temperature and optimal ligand-to-protein ratio (Fig. 7.3b). D2O solvent is preferred over H2O, to avoid direct and water-mediated saturation artifacts. High-affinity ligands typically reside longer within the protein binding site, undergo slow chemical exchange, and are thus not detectable by STD-NMR spectroscopy.

Fig. 7.3

Saturation transfer difference NMR (STD-NMR) analysis for group epitope mapping

(a) Basic principle of saturation transfer difference NMR (STD-NMR). The ligand protons which receive saturation from the protein are highlighted in red. From the degree of ligand saturation, one can map the group epitope which is recognized by the protein

(b) STD-NMR spectra of (Neu5Ac)6 with an anti-polysialic acid antibody 12E3 at 10 °C (ii), 20 °C (iii), and 30 °C (iv) with the control 1H-NMR spectrum (off-resonance) at 10 °C (i). STD-NMR spectra were collected by irradiating iteratively at 7 ppm (on resonance)/ 40 ppm (off resonance). x; the signals from low-molecular-weight impurities, which are nearly null in STD-NMR spectra. This figure was adapted from a previous paper (Hanashima et al. 2013) with permission

(c) STD-NMR spectrum (lower) and reference 1H-NMR spectrum (upper) of the di-galactosylated biantennary glycan bearing bisecting GlcNAc and core fucose (top panel) in the presence of PHA-E (glycan: subunit molar ratio is 1:2). Normalized amplification factor (GlcNAc-2 signal is set to 100%) is shown in parenthesis. The signal indicated with an asterisk originates from PHA-E protein. This figure was adapted from a previous paper (Nagae et al. 2014c) with permission. Crystal structure of PHA-E in complex with bisected N-glycan (PDB code: 5AVA) is shown in right panel

(b) STD-NMR spectra of (pan class="Chemical">Neu5Ac)6 with an anti-n>n class="Chemical">polysialic acid antibody 12E3 at 10 °C (ii), 20 °C (iii), and 30 °C (iv) with the control 1H-NMR spectrum (off-resonance) at 10 °C (i). STD-NMR spectra were collected by irradiating iteratively at 7 ppm (on resonance)/ 40 ppm (off resonance). x; the signals from low-molecular-weight impurities, which are nearly null in STD-NMR spectra. This figure was adapted from a previous paper (Hanashima et al. 2013) with permission

(c) STD-NMR spectrum (lower) and reference 1H-NMR spectrum (upper) of the di-galactosylated biantennary pan class="Chemical">glycan bearing bisecting n>n class="Chemical">GlcNAc and core fucose (top panel) in the presence of PHA-E (glycan: subunit molar ratio is 1:2). Normalized amplification factor (GlcNAc-2 signal is set to 100%) is shown in parenthesis. The signal indicated with an asterisk originates from PHA-E protein. This figure was adapted from a previous paper (Nagae et al. 2014c) with permission. Crystal structure of PHA-E in complex with bisected N-glycan (PDB code: 5AVA) is shown in right panel

1D STD-NMR analysis is often pan class="Chemical">hampered by severe overlapping of non>n class="Chemical">n-anomeric signals. For instance, the ring proton signals from the two branches (α1-3 and α1-6 branches) of bisected biantennary N-glycan largely overlap in 1D STD-NMR analysis. In contrast, the methyl proton signals from GlcNAc and core fucose are sharp and well separated for certain glycans. Hence these signals can be used as a reliable probe for branch-specific epitope mapping as exemplified by an analysis of the interaction between biantennary bisected glycan and a legume lectin PHA-E (Nagae et al. 2014c). The methyl proton signals of GlcNAc from outer branches and inner chitobiose unit show almost comparable intensities, indicating that the entire glycan contributes equally to binding with PHA-E (Fig. 7.3c). Consistent with this result, the crystal structure of bisected glycan-PHA-E complex showed that most sugar residues of the glycan ligand are recognized by PHA-E (Fig. 7.3c, right panel). Another approach is the use of 2D 1H-13C STD-heteronuclear single quantum coherence ( HSQC) to avoid the severe signal overlap in 1D NMR spectra. 2D STD- spectra have been collected to investigate the interaction between oligosialic acid and anti-oligosialic acid antibodies, 12E3 and A2B5 (Hanashima et al. 2013). STD-NMR is also applicable to analyze multiple binding modes of a unique ligand in a single binding site (Angulo and Nieto 2011). The duplicate sugar-binding modes of C-type lectin, DC-SIGN (Angulo et al. 2008), and anti-HIV antibody, 2G12 (Enríquez-Navas et al. 2012), have been analyzed in this way.

Application of 13C-NMR Using Deuterium-Induced Isotope Shift

In solution NMR experiments, it is difficult to directly observe the hydroxyl protons of pan class="Chemical">glycans because of their inherently fast excn>n class="Chemical">hange with water protons. Deuterium secondary isotope shifts on 13C-chemical shifts (deuterium-induced 13C isotope shifts, DIS) have been observed at the geminal 13C-signal of exchangeable hydroxyl protons (Pfeffer et al. 1979). In the presence of H2O and D2O, the line shape of the 13C-OH/D signal is highly dependent on the H/D exchange rate (Fig. 7.4a). In a very slow H/D exchange environment, a deuterium-induced 13C isotope shift (DIS) is observed in corresponding isotopomers, 13C-OH and 13C-OD, with a 0.09–0.15 ppm difference, whereas in the case of fast exchange, the equilibrium makes each isotopomer indistinguishable and gives a singlet signal at an averaged chemical shift.

Fig. 7.4

Hydrogen exchange analysis using deuterium-induced isotope shift

(a) Deuterium-induced 13C isotope shifts for analyzing the proton exchange rate of sugar hydroxyl groups. k ex, exchange rates of hydroxyl protons with water. The chemical shift difference (isotope shift) is ~0.15 ppm, which is not dependent on the magnetic field

(b) 13C-NMR spectra of 40 mM LewisX tetrasaccharide with expanded secondary carbon area. LewisX at 5 °C (top), and LewisX with 1.0 M calcium chloride at 5 °C (bottom). The sample was dissolved in 10 mM sodium acetate buffer (pH 6.0) composed of H2O:D2O = 1:1

(c) 13C-NMR spectra of short β1–3-glucans DP7 (red) and laminarin (black) at 5, 10, 15, 20, and 25 °C. 13C-NMR signal focuses on the C4 position. All ligands were dissolved in the same buffer shown in (B). The figures were adapted from previous papers (Hanashima et al. 2011, 2014a) with modifications

pan class="Chemical">Hydrogen excn>n class="Chemical">hange analysis using deuterium-induced isotope shift

(a) pan class="Chemical">Deuterium-induced n>n class="Chemical">13C isotope shifts for analyzing the proton exchange rate of sugar hydroxyl groups. k ex, exchange rates of hydroxyl protons with water. The chemical shift difference (isotope shift) is ~0.15 ppm, which is not dependent on the magnetic field

(b) pan class="Chemical">13C-NMR spectra of 40 mM n>n class="Gene">LewisX tetrasaccharide with expanded secondary carbon area. LewisX at 5 °C (top), and LewisX with 1.0 M calcium chloride at 5 °C (bottom). The sample was dissolved in 10 mM sodium acetate buffer (pH 6.0) composed of H2O:D2O = 1:1

(c) pan class="Chemical">13C-NMR spectra of short β1–3-glucans DP7 (red) and laminarin (black) at 5, 10, 15, 20, and 25 °C. n>n class="Chemical">13C-NMR signal focuses on the C4 position. All ligands were dissolved in the same buffer shown in (B). The figures were adapted from previous papers (Hanashima et al. 2011, 2014a) with modifications

This NMR technique pan class="Chemical">has been applied in the weak n>n class="Chemical">carbohydrate-carbohydrate interaction of LewisX (Hanashima et al. 2011). LewisX, a trisaccharide of galactose, fucose and GlcNAc, is often found at the termini of glycolipids and glycoproteins and initiates cell-cell interactions by forming a LewisX-Ca2+ complex (Eggens et al. 1989). Calcium-dependent intermolecular interactions of the LewisX-Ca2+ complex are observed, via the proton exchange rates of each hydroxyl group in an equimolar H2O/D2O solution (Fig. 7.4b). The conformation of the glycan ligand can also be analyzed by this method, including that of β(1,3)-glucan dependent on the degree of polymerization (DP) (Hanashima et al. 2014a). With β-glucan, the exchange of 13C-OH/D is dramatically retarded as DP extends, suggesting changes in quaternary structure (Fig. 7.4c). Recently a new two-dimensional 1H-13C NMR method, named INTOXSY, was developed to estimate the H/D exchange rate constants of OH groups based on DIS (Battistel et al. 2017).

Quantitative Lectin-Glycan Interaction Analysis

Like other proteipan class="Chemical">n-ligand interactions, lectin>n class="Chemical">n-glycan interactions are characterized by several biophysical properties, such as the dissociation constant K d, thermodynamics parameters such as enthalpy ΔH and entropy ΔS, and kinetics (e.g., k on and k off rate constants). Many quantitative biophysical methods exist to analyze lectin-glycan interactions such as frontal affinity chromatography (FAC) analysis, isothermal titration calorimetry (ITC), and surface plasmon resonance (SPR). All these methods can provide K d; however, a suitable method needs to be selected for each purpose. Among them, ITC can be performed without labeling or immobilization of lectin or glycan. A lot of sample (typically 10 mmol for ligand and protein) is usually needed. In FAC analysis, immobilization of lectin is required, and the glycan must be labeled for detection. SPR also requires immobilization of lectin or glycan, but amount of sample is typically small. Combining complementary methods provides more insights than a single method.

Frontal Affinity Chromatography (FAC)

Frontal affinity chromatography (FAC) is a powerful method to assess lectipan class="Chemical">n-glycan interactions as it is applicable to weak interactions (Tateno et al. 2007). In FAC analysis, n>n class="Chemical">glycans are applied to a lectin-immobilized column and the elution profile monitored. Dissociation constants for each glycan are calculated from the equation:where V and V 0 are the elution volumes of the column and negative control, respectively. B t is the effective lectin content expressed in mole, and [A]0 is the initial glycan molar concentration. The amount of analyte needed depends on the affinity of analyte and ligand: higher affinities require smaller amounts.

An example of FAC analysis is tpan class="Chemical">hat of the binding of various n>n class="Chemical">N-glycans to a mannose-binding Jacalin-related lectin Calsepa. The lectin was immobilized on the column, and a panel of pyridylaminated (PA) or p-nitrophenyl (pNP) glycans were successively passed through the column. In this way, Calsepa was found to specifically bind bisected N-glycans (Nagae et al. 2017a) (Fig. 7.5a). The specificity is explained by the crystal structure of bisected glycan-Calsepa complex (Fig. 7.5a, right). The ligand is in a flipped-back conformation, which bisected glycan prefers, and the unique ligand conformation enables the 2:1 sandwich binding mode.

Fig. 7.5

Quantitative analyses of lectin-glycan interaction

(a) FAC analysis of a mannose-binding Jacalin-related lectin Calsepa toward N-glycans. Association constant (K A) of Calsepa lectin for each glycan is shown. Monosaccharide symbols follow the SNFG (symbol nomenclature for glycans) system (Varki et al. 2015). The FAC data is adapted from a paper (Nagae et al. 2017a) with permission. Crystal structure of Calsepa in complex with bisected glycan (corresponding to #104) is shown in right panel (PDB code: 5AV7)

(b) Isothermal titration calorimetry (ITC) analysis of scFv735 using polysialic acid (PolySia, DP 80–130) (left). Bar graph and table of thermodynamic parameters for DP4, DP5, DP6, and PolySia are shown in middle panel. Crystal structure of scFv735 in complex with octa- sialic acid (DP8) is shown in the right panel (PDB code: 3WBD)

(c) Surface plasmon resonance (SPR) analyses of the protein-protein interaction between p24β1 and p24δ1 GOLD domains. (i) Sensorgram showing the interaction between immobilized p24δ1 GOLD domain and analyte p24β1 GOLD domain. p24β1 GOLD domain aliquots (0, 5, 25, 50, 100, 200, and 400 μM) were injected for 60 s at a flow rate of 30 μl/min. Subtracted sensorgrams with blank lanes are shown. (ii) Sensorgram showing the interaction between immobilized p24β1 GOLD domain and analyte p24δ1 GOLD domain. The dilution series of p24δ1 GOLD domain are the same as those of p24β1 GOLD domain shown in (i). The SPR data was adapted from a paper (Nagae et al. 2016a) with permission

Quantitative analyses of lectipan class="Chemical">n-glycan interaction

(a) FAC analysis of a pan class="Chemical">mannose-binding Jacalin-related lectin Calsepa toward N-glycans. Association constant (K A) of Calsepa lectin for each glycan is shown. Monosaccharide symbols follow the SNFG (symbol nomenclature for glycans) system (Varki et al. 2015). The FAC data is adapted from a paper (Nagae et al. 2017a) with permission. Crystal structure of Calsepa in complex with bisected glycan (corresponding to #104) is shown in right panel (PDB code: 5AV7)

(b) Isothermal titration calorimetry (ITC) analysis of pan class="Gene">scFv735 using n>n class="Chemical">polysialic acid (PolySia, DP 80–130) (left). Bar graph and table of thermodynamic parameters for DP4, DP5, DP6, and PolySia are shown in middle panel. Crystal structure of scFv735 in complex with octa- sialic acid (DP8) is shown in the right panel (PDB code: 3WBD)

(c) Surface plasmon resonance (SPR) analyses of the proteipan class="Chemical">n-protein interaction between n>n class="Gene">p24β1 and p24δ1 GOLD domains. (i) Sensorgram showing the interaction between immobilized p24δ1 GOLD domain and analyte p24β1 GOLD domain. p24β1 GOLD domain aliquots (0, 5, 25, 50, 100, 200, and 400 μM) were injected for 60 s at a flow rate of 30 μl/min. Subtracted sensorgrams with blank lanes are shown. (ii) Sensorgram showing the interaction between immobilized p24β1 GOLD domain and analyte p24δ1 GOLD domain. The dilution series of p24δ1 GOLD domain are the same as those of p24β1 GOLD domain shown in (i). The SPR data was adapted from a paper (Nagae et al. 2016a) with permission

Isothermal Titration Calorimetry (ITC)

Isothermal titration calorimetry (ITC) is a technique tpan class="Chemical">hat directly determines thermodynamic parameters of protein-ligand interactions (Dam et al. 2016). In a single experiment, ITC measures n (number of binding sites on the protein), ΔH ( enthalpy of binding), and K A ( association constant). From the K A value, the free energy of binding (ΔG) and the entropy of binding (ΔS) can be calculated from the following equation:where ΔG, ΔH, and ΔS are the changes in free energy, enthalpy, and entropy of binding, respectively. T is the absolute temperature and R is the gas constant (8.314 J/mol/K).

An example is the ITC analysis of an anti-pan class="Chemical">polysialic acid antibody 735. This antibody shows high affinity toward α2–8-linked n>n class="Chemical">polysialic acid as the degree of polymerization (DP) increases. To investigate the DP-dependent affinity enhancement mechanism, a series of oligo (DP4–6) and polysialic acids (DP 80–130) and single chain Fv originating from antibody 735 (scFv735) were evaluated by ITC (Nagae et al. 2013a) (Fig. 7.5b). As DP increases so does the affinity of scFv735. Thermodynamic parameters revealed that a large entropy loss upon interaction (TΔS = −72.7 kcal/mol for polysialic acid) is compensated for by a larger enthalpy gain (ΔH = −80.1 kcal/mol). Importantly, the crystal structure of scFv735-ligand complex indicates that antibody 735 can recognize internal residues of polysialic acids (Fig. 7.5b, right). Therefore, increasing DP may provide a rapid cycle of dissociation-association and thereby enhance the affinity (Nagae and Yamaguchi 2014).

Typical lectin class="Chemical">n-glycan interactions are exothermic, but occasionally an endothermic reaction is observed, e.g., interaction between n>n class="Chemical">heparan sulfate and heparin-binding hemagglutinin (Huang et al. 2017). Kinetic parameters can also be obtained from ITC experiments (Burnouf et al. 2012; Vander Meulen et al. 2016).

Surface Plasmon Resonance (SPR)

Surface plasmon resonance (SPR) is one of the most widely used analytical methods for detecting biomolecular interactions and kinetics. It monitors interactions between an immobilized ligand on a thin pan class="Chemical">metal layer and a soluble-free analyte by detecting cn>n class="Chemical">hanges in the resonance angle due to an increase in concentration of analyte at the surface. SPR can monitor real-time interaction with biosensors by measuring the resonance units against time.

Sensorgrams can be utilized for the determination of kinetic parameters, association/pan class="Chemical">dissociation constants (K A/K D), on/off rates (k on/k off), as well as thermodynamics parameters. There are two major methods to calculate kinetic parameters, first is by directly fitting the curve sensorgram to ideal binding models (more applicable to slow association/dissociation binding) and the second is steady-state analysis (for fast association/dissociation interactions). Maximum resonance values of a series of analyte concentrations (R eq) are plotted against analyte concentration. In simple 1:1 steady-state binding model, the plot is fitted to the following equation:where C is concentration of analyte, R max is maximum binding response, and K D is the equilibrium dissociation constant. For example, the heterophilic interaction of p24 Golgi dynamics (GOLD) domains is shown between p24β1 and p24δ1. These GOLD domains interact with each other to form a hetero p24 family protein complex which is involved in the transport of glycosylphosphatidylinositol (GPI)-anchored proteins. The sensorgram shows fast association/dissociation binding and was analyzed by a 1:1 steady-state model (Nagae et al. 2016a) (Fig. 7.5c). To avoid immobilization artifacts, two experiments were performed: first p24β1 was immobilized and p24δ1 injected and then the vice versa experiment performed. The estimated K D values of the two runs are comparable (1.2 × 10−4 and 1.5 × 10−4 M). SPR is widely used in lectin-glycan interactions, either alone or combined with other methods such as the enzyme-linked immunosorbent assay (ELISA), ITC, and thermal shift assays (Wesener et al. 2015; Ribeiro et al. 2016; Machon et al. 2017). When the lectin is immobilized and glycan injected, the SPR detection limit needs to be confirmed in advance. To overcome the sensitivity problem, biotinylated carbohydrates are immobilized onto streptavidin or neutravidin-coated sensor chips. Diluted lectins are injected and affinities calculated .

3D Structural Studies on Lectin-Glycan Interactions

Biophysical data of lectipan class="Chemical">n-glycan interaction is strengthened when combined with structural data obtained from X-ray crystallography and NMR spectroscopy. X-ray crystallography is the favored method to determine the 3D structures of protein-ligand complexes at atomic resolution. However, glycans are highly flexible in solution, and the electron density is sometimes missing even in the presence of lectin. In contrast, solution NMR analyses such as titration assays and transferred nuclear Overhauser effects (trNOE) can give a more realistic picture under physiological conditions. Titration assays provide affinities, while trNOEs give intramolecular distance constraints of the bound glycan. This chapter focuses on NMR analysis of lectin-glycan interaction.

NMR Titration Study

For any binding scheme , a ligand excpan class="Chemical">hanges between free and bound states with rate constant k (expressed per second). In NMR spectra, resonance frequencies of bound and free ligand are different, with Δν representing the chemical shift difference (expressed in Hz). In a fast excn>n class="Chemical">hange regime, k >> Δν, the signal will appear in a population-weighted position between free and bound resonances. This allows one to monitor binding using the chemical shift change. In 1H-15N HSQC experiments, the dissociation constant K D can be obtained by plotting the weighted chemical shift change, Δδ, as a function of the carbohydrate/protein molar ratio. Δδ of each concentration is calculated with the following equation:where Δδ H and Δδ N are the observed chemical shift changes (ppm) of 1H and 15N and α is a scaling factor (typically 0.1~0.2). Then the dissociation constants is estimated by fitting the curve:where Δδ max is the maximum chemical shift change, [L] the total ligand concentration, [P] the total protein concentration, and K D the dissociation constant for a 1:1 binding equilibrium.

In a slow excpan class="Chemical">hange regime, k < < Δν, two distinct resonances are observed, corresponding to the free and bound states. Increasing concentrations of carbohydrate increases the signal of the bound form and decreases that of the free. When the exchange rate is comparable to the frequency difference between the two states (intermediate exchange regime), k ~ Δν, extensive line broadening occurs, and sometimes the signals are beyond detection. Occasionally, it may be possible to escape this undesirable intermediate exchange situation by changing experimental conditions, aiming to either accelerate or slow the exchange.

A first example is the interaction of pan class="Species">human soluble lectin n>n class="Gene">ZG16p and phosphatidylinositol mannosides (PIMs) (Hanashima et al. 2015). PIM2 was titrated into a 15N-labeled ZG16p solution, and signal perturbations in the 1H-15N HSQC spectra were recorded (Fig. 7.6a). The binding is a fast exchange process, since each set of specific signals featured a gradual chemical shift change. The binding interface of PIM2 can be inferred from a map of the hot spots on the 3D structure of ZG16p, and the dissociation constant for PIM2 was 3.0 mM from a fit of the titration curve.

Fig. 7.6

NMR titration studies and trNOE-based conformational analyses

(a) NMR titration study of the interaction between ZG16p and PIM2. (Left panel) 1H-15N HSQC spectra of uniformly 15N-labeled ZG16p in titration with PIM2 glycan (Black, 0 equiv.; red, 5 equiv.; and green, 20 equiv. of PIM2). Blue arrows indicate directions of the chemical shift changes. The figure was adapted from a paper (Hanashima et al. 2015) with permission.Fig. 7.6 (continued) (Right panel) Mapping surface residues in the PIM2 glycan interaction on the crystal structure of human ZG16p (PDB ID; 3APA). The signal from I149 (green) was broadened upon PIM2 binding

(b) NMR titration study of the interaction between a C-type lectin receptor mDCIR2 CRD and glycans. The pyridylaminated bisected and non-bisected N-glycans (glycans a and b) used in this study are shown in middle panel. A part of 1H-NMR spectra of 20 μM mDCIR2 CRD in the absence of glycan (top), in the presence of 40 μM glycan (second), and 40–120 μM glycan (third to the last). The figure was adapted and modified from a previous paper (Nagae et al. 2013b) with permission. Crystal structure of mDCIR2 CRD in complex with bisected glycan corresponding to glycan is shown in right panel (PDB code: 3VYK)

(c) Basic principle of transferred NOE (trNOE). If the glycan ligand is bound to protein for a sufficiently short time (in fast exchange), one can observed the bound NOEs (reflecting the interproton distances in the bound state) by using the free ligand signals

(d) NOE build up curves of Man-3 H1 and Man-3 O-methyl proton signals in the presence (left) and absence (right) of Calsepa lectin. Man-3 H1 signal was selectively inverted using a 180° rectangular pulse with 40-ms duration. This figure was adapted from a paper (Nagae et al. 2016c) with permission

(e) 1D selective NOESY spectra of bisected glycan in the presence (upper) and absence (lower) of Calsepa. Strong intra-residue trNOE was observed from Man-3 H1 to Man-3 H2, and long-range Fig. 7.6 (continued) trNOEs were also observed from Man-3 H1 to GlcNAc-5′ (α1–6 branch) H1 and to Man-4’ H2 signals (Right panel). Proton-proton distances were indicated between Man-3 H1 and GlcNAc-5’ H1 and between Man-3 H1 and Man-4’ H2. Structures of back-fold conformations are derived from Calsepa complex (PDB code: 5AV7). Hydrogen atoms are generated with PyMOL. This figure was adapted and modified from a paper (Nagae et al. 2016c) with permission

(f) 2D 1H-1H NOESY spectrum of PIM1 glycan (4.6-fold excess) in the presence of wild-type ZG16p collected with a mixing time of 250 ms at 10 °C (upper panel). Key inter-residual correlation was provided between Man-H1 and Ino-H2 and intra-residual correlations (Man-H1 and Man-H2) as negative NOEs. The atomic distance of inter-residue Man-H1-Ino-H2 was determined as 2.2 Å from the relative intensity of the signal. No correlation was observed when using the inactive ZG16p mutant (D151N) (lower panel). This figure is reproduced from the paper (Hanashima et al. 2015) with permission

(a) NMR titration study of the interaction between pan class="Gene">ZG16p and n>n class="Gene">PIM2. (Left panel) 1H-15N HSQC spectra of uniformly 15N-labeled ZG16p in titration with PIM2 glycan (Black, 0 equiv.; red, 5 equiv.; and green, 20 equiv. of PIM2). Blue arrows indicate directions of the chemical shift changes. The figure was adapted from a paper (Hanashima et al. 2015) with permission.Fig. 7.6 (continued) (Right panel) Mapping surface residues in the PIM2 glycan interaction on the crystal structure of human ZG16p (PDB ID; 3APA). The signal from I149 (green) was broadened upon PIM2 binding

(b) NMR titration study of the interaction between a C-type lectin receptor pan class="Gene">mDCIR2 CRD and n>n class="Chemical">glycans. The pyridylaminated bisected and non-bisected N-glycans (glycans a and b) used in this study are shown in middle panel. A part of 1H-NMR spectra of 20 μM mDCIR2 CRD in the absence of glycan (top), in the presence of 40 μM glycan (second), and 40–120 μM glycan (third to the last). The figure was adapted and modified from a previous paper (Nagae et al. 2013b) with permission. Crystal structure of mDCIR2 CRD in complex with bisected glycan corresponding to glycan is shown in right panel (PDB code: 3VYK)

(c) Basic principle of transferred NOE (trNOE). If the pan class="Chemical">glycan ligand is bound to protein for a sufficiently short time (in fast excn>n class="Chemical">hange), one can observed the bound NOEs (reflecting the interproton distances in the bound state) by using the free ligand signals

(d) NOE build up curves of Mapan class="Chemical">n-3 H1 and Man-3 O-methyl proton signals in the presence (left) and absence (right) of Calsepa lectin. Man-3 H1 signal was selectively inverted using a 180° rectangular pulse with 40-ms duration. This figure was adapted from a paper (Nagae et al. 2016c) with permission

(e) 1D selective NOESY spectra of bisected pan class="Chemical">glycan in the presence (upper) and absence (lower) of n>n class="Chemical">Calsepa. Strong intra-residue trNOE was observed from Man-3 H1 to Man-3 H2, and long-range Fig. 7.6 (continued) trNOEs were also observed from Man-3 H1 to GlcNAc-5′ (α1–6 branch) H1 and to Man-4’ H2 signals (Right panel). Proton-proton distances were indicated between Man-3 H1 and GlcNAc-5’ H1 and between Man-3 H1 and Man-4’ H2. Structures of back-fold conformations are derived from Calsepa complex (PDB code: 5AV7). Hydrogen atoms are generated with PyMOL. This figure was adapted and modified from a paper (Nagae et al. 2016c) with permission

(f) 2D 1H-1H NOESY spectrum of pan class="Gene">PIM1 n>n class="Chemical">glycan (4.6-fold excess) in the presence of wild-type ZG16p collected with a mixing time of 250 ms at 10 °C (upper panel). Key inter-residual correlation was provided between Man-H1 and Ino-H2 and intra-residual correlations (Man-H1 and Man-H2) as negative NOEs. The atomic distance of inter-residue Man-H1-Ino-H2 was determined as 2.2 Å from the relative intensity of the signal. No correlation was observed when using the inactive ZG16p mutant (D151N) (lower panel). This figure is reproduced from the paper (Hanashima et al. 2015) with permission

The second example is the interaction of a C-type lectin receptor pan class="Gene">mDCIR2 and bisected n>n class="Chemical">N-glycan (Nagae et al. 2013b). Bisected and non-bisected N-glycans (glycans and ) were analyzed. Concentration-dependent spectral shifts in the 1H-NMR spectra were only observed in the presence of bisected glycan (Fig. 7.6b). The interaction involved slow exchange in terms of chemical shift, and the dissociation constant was estimated to be 3.0 × 10−5 M, a ligand concentration at which there is an equal amount of ligand-free and ligand-bound lectin. The specificity is well explained from the crystal structure of mDCIR2-bisected glycan complex, in which the bisected GlcNAc residue is directly interacting with the protein (Fig. 7.6b, right).

Transferred Nuclear Overhauser Effects (trNOE)

The nuclear Overpan class="Chemical">hauser effect (NOE) describes the phenomenon of magnetization transfer from one n>n class="Gene">spin to another. This process proceeds through cross relaxation originating from dipolar interaction between two spins that are spatially close (less than 5 Å apart). Roughly, the NOE is related to the distance between the two protons and shows a 1/r 6 dependence.

If the pan class="Chemical">glycan ligand binds to lectin with slow excn>n class="Chemical">hange, both intermolecular lectin-glycan NOEs and intramolecular NOEs of the glycan are important in elucidating the bound conformation of the ligand. However, in most lectin-glycan interactions, the exchange is fast because the affinity is weak. In such cases, transferred- NOESY (tr-NOESY) is the suitable measurement for conformational analysis of bound glycan ligand. trNOE analysis is generally valid for ligands that bind in the μM-mM range (Fig. 7.6c).

TrNOE experiments were applied to the conformational analysis of bisected n class="Chemical">N-n>n class="Chemical">glycan bound to a lectin, Calsepa (Nagae et al. 2016c). The NOE build-up curve shows linearity up to 200 ms for the glycan-Calsepa mixture and to 500 ms for glycan alone (Fig. 7.6d). The mixing time was then chosen as 200 ms for the trNOE experiments. In the presence of Calsepa lectin, long-range inter-residual trNOEs were observed between core β-mannose (Man-3) and GlcNAc-5′ on the outer branch (Fig. 7.6e). The presence of this long-range trNOE indicates that the glycan assumes a compact back-fold conformation when bound to Calsepa lectin. In a second example, the interaction between ZG16p and phosphatidylinositol mannoside glycans (PIM1 and PIM2) was investigated by trNOE (Hanashima et al. 2015). In the presence of ZG16p, several inter- and intra-residual NOE correlations were observed in PIM1 (Fig. 7.6f). There were none with the inactive ZG16p mutant D151N. From these correlations, the conformation of PIM1 in ZG16p-bound state was elucidated.

The structural information obtained from trNOE is usually limited by low proton density and by the inherently fpan class="Gene">lexible nature of n>n class="Chemical">carbohydrate chains. Thus, obtained distances are typically further evaluated by molecular dynamics simulations and/or X-ray crystallographic analysis.

Other NMR Techniques

Residual dipolar coupling (RDC) measured in liquid crystalline alignment medium offer a viable alternative to traditional NOE-based approaches for structural analysis. RDC provides a way to constrain the relative orientation of two molecules in compn class="Gene">lex with each other by aligning their independently determined order tensors (Jain 2009). Proton>n class="Chemical">n-proton and proton-carbon RDC of heparin tetrasaccharide has been measured (Jin et al. 2009), which, when combined with molecular dynamics simulations, clarified ring conformer and global shape of the glycan.

Paramagnetic probes attached to the reducing ends of pan class="Chemical">oligosaccharides cause paramagnetic relaxation enn>n class="Chemical">hancement (PRE) and/or pseudocontact shift (PCS) which can resolve the peak overlap problem (Kato and Yamaguchi 2015). These spatial perturbations can be sources of long-range atomic distance information (~10 Å), which complements the local conformational information derived from J-couplings and NOEs (~5 Å).

Molecular Dynamics Simulation Studies of Lectin-Glycan Interactions

Theoretical approaches such as docking models and molecular dynamics (MD) simulation complement experimental information and pan class="Chemical">have great possibilities for visualizing the highly mobile n>n class="Chemical">glycans. Docking simulations of lectin-glycan complexes provide a theoretical foundation and reduce ambiguities in experimental results. We validated the accuracy of docking simulation by applying it to a crystal structure of POMK-ligand complex, in which the electron density of ligand was partially observed (Fig. 7.7a). The position of the trisaccharide ligand in the docking model coincides well with the fragmented electron density in the crystal structure, underscoring the usefulness of docking and MD simulations in determining bound structures (Nagae et al. 2017b).

Fig. 7.7

Molecular dynamics simulations of protein-glycan complexes

(a) Top five docking poses of the core M3 trisaccharide ( GalNAcβ1–3GlcNAcβ1-4Man, shown in yellow) in POMK binding site (i-v) by AutoDock Vina. The mannose residue assigned in the crystal structure (PDB code: 5GZ9) is overlaid in cyan. The table contains docking results for each model (i)-(v). Calculated binding affinity (in kcal/mol) together with root mean square deviation (r.m.s.d. in Å) from top scoring docking pose (i) is shown for all five models. This figure is adapted from the paper (Nagae et al. 2017b) with permission

(b) Three major conformations (Conformers 1, 2, and 3) of the α1–6 branch of glycan A bound to Orysata lectin in MD simulations. The simulation was performed using the crystal structure of Orysata in complex with glycan A (PDB code: 5XFH). This figure is adapted from the paper (Nagae et al. 2017a) with permission

(c) Bar graph and table of calculated average enthalpy (∆H), entropy (T∆S), and total binding free energy (∆G) for the biantennary glycans A–D. Experimental dissociation constant (K D = 1/K A) determined by FAC analysis is also indicated. N.D.: K D was not determined due to weak or no binding

Molecular dynamics simulations of proteipan class="Chemical">n-glycan comppan class="Gene">lexes

(a) Top five docking poses of the core M3 pan class="Chemical">trisaccharide ( GalNAcβ1–3n>n class="Chemical">GlcNAcβ1-4Man, shown in yellow) in POMK binding site (i-v) by AutoDock Vina. The mannose residue assigned in the crystal structure (PDB code: 5GZ9) is overlaid in cyan. The table contains docking results for each model (i)-(v). Calculated binding affinity (in kcal/mol) together with root mean square deviation (r.m.s.d. in Å) from top scoring docking pose (i) is shown for all five models. This figure is adapted from the paper (Nagae et al. 2017b) with permission

(b) Three major conformations (Conformers 1, 2, and 3) of the α1–6 branch of pan class="Chemical">glycan A bound to Orysata lectin in MD simulations. The simulation was performed using the crystal structure of Orysata in compn>n class="Gene">lex with glycan A (PDB code: 5XFH). This figure is adapted from the paper (Nagae et al. 2017a) with permission

(c) Bar graph and table of calculated average entpan class="Chemical">halpy (∆H), entropy (T∆S), and total binding free energy (∆G) for the biantennary n>n class="Chemical">glycans A–D. Experimental dissociation constant (K D = 1/K A) determined by FAC analysis is also indicated. N.D.: K D was not determined due to weak or no binding

MD simulation is also useful for validation of possible binding modes with the estimation of binding free energy. Orysata lectin, one of the pan class="Chemical">mannose-binding Jacalin-related lectins, shows a preference for complex-type glycans with a α1-6 branch extension rather than a α1-3 branch. It was not easy to understand the structural basis for the preference for the α1-6 branch from the “snap-shot” crystal structure of Orysata in complex with biantennary glycan. MD simulations yield the dynamic behavior of the glycan attached to Orysata lectin (Fig. 7.7b). Several docking models of Orysata in complex with a series of N-glycans were constructed, and the calculated binding free energies suggest that multiple binding modes exist in the solution, and the binding preference seems somehow to originate from an average of multiple binding modes that favors α1-6 branch binding (Fig. 7.7c). MD simulations are especially useful for analyzing such multiple ligand-binding modes .

Future Perspective: Application of Cryo-electron Microscopy Analysis to Investigate Lectin-Glycan Interaction

Recent improvement in cryo- electron microscopy (cryo-EM) is a great breakthrough for structural biology (Cheng et al. 2015). EM analysis pan class="Chemical">has an inherent advantage for investigating the global conformations of macromolecules, rather tn>n class="Chemical">han atomic details of local interfaces. The structural changes of glycans on glycoproteins affect the physiological functions of carrier glycoproteins. This has been clearly exemplified by EM analysis. For example, EM has shown that mutation or elimination of N-glycans affects the conformational equilibrium of α5β1 integrin (Li et al. 2017). Likewise, the inter-domain angle of LDL-receptor-related protein 6 (LRP6) is regulated by the evolutionarily conserved N-glycan (Matoba et al. 2017).

In X-ray crystallography, the trimming of surface pan class="Chemical">glycans on glycoproteins by glycosidase is usually necessary to produce diffraction quality crystals (Chang et al. 2007). Mutations at N-glycosylation sites can also produce highly uniform proteins (Nagae et al. 2017b). However, such treatments prevent visualization of the atomic structure of glycans, although some native glycoproteins provide good glycan electron density due to extensive intramolecular interactions between glycan and carrier glycoproteins (Nagae and Yamaguchi 2012). In contrast, there is no need to eliminate the surface glycans of glycoproteins in single particle analysis of cryo-EM. Thus, cryo-EM has great potential to provide the atomic structure of native glycans. Actually, cryo-EM structures of glycoproteins, such as γ-secretase and human coronavirus, often include the coordinates of attached glycans (Bai et al. 2015; Walls et al. 2016) (Fig. 7.8a and b). Even tetra-antennary N-glycan structures are visible in the EM map of the HIV-1 envelope glycoprotein trimer (Env) (Lee et al. 2016) (Fig. 7.8c).

Fig. 7.8

Electron microscopic images of carbohydrates covalently attached onto proteins

(a) N-glycan core attached onto N55 of human nicastrin, a component of γ-secretase complex ((Bai et al. 2015), PDB code: 5A63, EMDB: 3061). The density map is depicted at 5.0 σ level cutoff. Schematic representation of the observed glycan is shown in right panel

(b) High-mannose-type N-glycan attached onto N426 of human coronavirus NL63 (HCoV-NL63) ((Walls et al. 2016), PDB code, 5SZS; EMDB, 8331). The density map is contoured at 7.0 σ level

(c) Tetra-antennary glycan attached onto N637 of HIV-1 Env trimer ((Lee et al. 2016), PDB code, 5FUU; EMDB, 3308). The density map is depicted in cyan mesh contoured at 5.5 σ level

Electron microscopic images of pan class="Chemical">carbohydrates covalently attached onto proteins

(a) pan class="Chemical">N-glycan core attached onto N55 of n>n class="Species">human nicastrin, a component of γ-secretase complex ((Bai et al. 2015), PDB code: 5A63, EMDB: 3061). The density map is depicted at 5.0 σ level cutoff. Schematic representation of the observed glycan is shown in right panel

(b) High-pan class="Chemical">mannose-type n>n class="Chemical">N-glycan attached onto N426 of human coronavirus NL63 (HCoV-NL63) ((Walls et al. 2016), PDB code, 5SZS; EMDB, 8331). The density map is contoured at 7.0 σ level

(c) Tetra-antennary pan class="Chemical">glycan attached onto N637 of n>n class="Disease">HIV-1 Env trimer ((Lee et al. 2016), PDB code, 5FUU; EMDB, 3308). The density map is depicted in cyan mesh contoured at 5.5 σ level

Due to limitation of molecular weight, it is normally difficult to directly visualize lectipan class="Chemical">n-glycan interactions using single particle analysis by cryo-EM. However, population analysis of many particles is the great advantage of EM analysis. Although the flexibility of the glycan part remains a serious problem for visualization (averaging problem), we can expect that cryo-EM will become the method of choice for analyzing the structure of native glycan on glycoproteins. The quaternary structural analyses of elaborate complex systems have recently been reported, such as human MHC-I peptide-loading complex (PLC), which includes TAP1, TAP2, tapasin, MHC-I heavy chain, β2m, calreticulin, and ERp57 in the complex (Blees et al. 2017). It is likely that multicomponent glycoprotein complexes and lectin receptor-glycoprotein complexes will be analyzed by cryo-EM at atomic resolution, coupled with other biophysical methods in the very near future.