Pedro M. Nieto- Glycosystems Laboratory, Instituto de Investigaciones Científicas, cicCartuja, CSIC/USE, Seville, Spain
Carbohydrates are biologically ubiquitous and are essential to the existence of all known living organisms. Although they are better known for their role as energy sources (glucose/glycogen or starch) or structural elements (chitin or cellulose), carbohydrates also participate in the recognition events of molecular recognition processes. Such interactions with other biomolecules (nucleic acids, proteins, and lipids) are fundamental to life and disease. This review focuses on the application of NMR methods to understand at the atomic level the mechanisms by which sugar molecules can be recognized by proteins to form complexes, creating new entities with different properties to those of the individual component molecules. These processes have recently gained attention as new techniques have been developed, while at the same time old techniques have been reinvented and adapted to address newer emerging problems.
Introduction
Glycobiology, which can be defined as the study of the structure, chemistry, biosynthesis, and biological functions of glycans and their derivatives, is fundamental in many critical biological processes (Varki, 1993, 2017). A challenge of structural glycobiology is to reconcile the large variety of 3D shapes that carbohydrates can assume with the high degree of selectivity found between closely related glycans (DeMarco and Woods, 2008). In addition, the binding constants between individual carbohydrates and proteins are generally low. Under these conditions, most of the interactions are in the fast exchange region of the NMR chemical shift scale and a single averaged set of signals is detected. This behavior can be exploited for structural elucidation. The rationale behind transient or transfer experiments is that if the equilibrium is fast enough and the property is dependent on the correlation time, it is possible to observe properties characteristic of the bound state in the averaged ligand signals, due to the larger correlation time of the bound form (Hyde et al., 1980; Feeney, 2000; Meyer and Peters, 2003).
In order to understand the evolution of NMR methods for the study of glycans, the scarcity of simple methods for obtaining isotopically-labeled glycans needs to be kept in mind. When such methods are available, the use of isotope-filtered/edited NMR experiments should be considered. When there is no such alternative, however, transfer techniques need to be employed. These can be applied to systems in fast equilibrium in the relaxation time scale and to properties weighted by the correlation time. Under these circumstances, the average of the NMR property obtained is biased toward the minor population of the bound carbohydrate due to its larger correlation time, even in the presence of an excess of free ligand. This is the basis of the transfer NOE technique and its analogs (Hyde et al., 1980; Ni, 1994). The advent of methodological improvements, monodimensional analogs with better signal to noise ratio, has effectively relaunched the transfer techniques (Stott et al., 1997). Other parameters that can be used to study molecular complexes in a transient state are T1sel (selective longitudinal relaxation time), or T2 (transversal relaxation time) and ligand-detected 1H relaxation dispersion because they depend on the correlation time, which is a function of the size of the molecule. Theoretical descriptions of T1 and T2 and of their application to the analysis of binding constants have been reviewed by Stockman and Dalvit (2002) and Peng et al. (Lepre et al., 2004).
A further experiment that benefits from the difference in the correlation times is the STD (saturation transfer difference) experiment (Mayer and Meyer, 1999; Meyer and Peters, 2003). This experiment uses the faster transfer of saturation from the receptor, caused by its long correlation time, that is further transferred to the ligand within the complex and that finally is observed in signals corresponding to the free-state. Since its formulation, the STD technique has experienced a rapid growth, both in its range of applications and in the number of labs using it. Finally, WaterLOGSY uses protonated water molecules to distinguish those that are in fast exchange from those that are buried in the interface with the receptor and therefore reflects the interaction surfaces between receptor and ligand.
NMR Characteristics of Free Carbohydrates
In general, the principal sources of glycans are via chemical synthesis or natural product isolation. In each case, they are almost completely restricted to non-isotopically-labeled compounds. As, from an NMR viewpoint, there are few experimental restrictions for unlabeled-carbohydrates, the analysis conditions must nevertheless be carefully documented and quantified.
Coupling Constants
Three-bond H – H coupling constants depend on the dihedral angles, and consequently they are key to determining the ring conformation. In general, most hexoses are monoconformational and are generally in a 4C1 or 1C4 chair conformation. Due to the cyclic nature of the sugars, many redundant interprotonic three-bond coupling constants are available and it is relatively simple to define the conformation from the Cremer–Pople polar coordinates (Cremer and Pople, 1975). Residual dipolar coupling constants provide also data about the relative orientation of two vectors in partially oriented media (Martin-Pastor and Bush, 2001). In carbohydrates, few examples of conformational equilibrium are known; however, one such example is the iduronate residue, which can exist in at least three conformations, 1C4, 4C1, and 2SO, easily distinguishable by their coupling constants. Generally, internal iduronate residues are found in fast equilibrium between 1C4 and 2SO (Figure 1A,B) (Mulloy et al., 1993). This equilibrium is fundamental in the molecular recognition events in which heparin is involved (Canales et al., 2005) and also has been studied by residual dipolar coupling (Jin et al., 2009)
Figure 1. Structure of a heparin trisaccharide with the central iduronate in: (A) 1C4 and (B) 2SO conformation; (C) NOESY 500 MHz and 600 ms tmix; (D) transfer NOESY (800 MHz, 200 ms tmix)—the labels correspond to the assignment of cross-peaks—notice the peak I2 H2—H5, a NOESY peak exclusive to the 2SO conformation; and (E) STD—affinity factor growing curves for the extraction of the initial rates for quantitative analysis (Muñoz-García et al., 2015).
NOE, Distances: The ISPA (Isolated Spin Pair Approach)
NOE (and other related parameters such as ROE, transverse NOE, and off-resonance ROE)-based distance constraints are often used to study the glycosidic angle behavior and can thus be used for the determination of the 3D shapes of saccharides. The scarcity of experimental restrictions in the analysis of carbohydrates implies a more accurate quantification of the NOE. Fortunately, for small molecules, distances are easily quantifiable from the NOE by means of an analysis of the growing NOE curve at several mixing times (Macura et al., 1986). Then, assuming equal motional behaviors for an undetermined distance and for a known reference, the initial NOE growth rate (cross-relaxation rate) relationship is inversely proportional to the sixth power of the interprotonic distance (Neuhaus and Williamson, 2000).
Applications of this methodology have increased since the arrival of modern selective 1D NOESY methods based on the double pulsed field gradients spin echo technique (Stott et al., 1997). This scheme has advantages over the bidimensional analogs; it is faster and has better resolution. Moreover, as the protons are relaxing by T1sel instead of T1, the growing curves have longer linearity, and thus better-fitting and more precise results can be obtained (Hu and Krishnamurthy, 2006; Munoz-Garcia et al., 2013). The technique has also been applied to off-resonance ROESY, so as to obtain reference-independent distances. Using this approach, it has been possible to obtain the interglycosidic distances for a strongly anisotropic heparin hexasaccharide in which the interprotonic correlation time depends upon the orientation of the vector relative to the axis of the molecule (Munoz-Garcia et al., 2013).
NMR Characteristics of Bound Small Glycans
The ability of NMR to observe and analyze specific signals of individual atoms, focusing analysis on a particular aspect of a complex without the need to solve its entire structure, is the most important consideration.
A bound ligand has NMR properties that are governed by its geometry in the bound state and these can therefore be different to those of the ligand in its free state. In a situation of fast equilibrium in the NMR timescale, the magnitude of the parameter that is measured corresponds to the weighted average of the NMR properties of the free and the bound states. For the special case of correlation time-dependent properties, the weight of the properties derived from the complexed ligand is so large that it is predominant in the averaged values, even at low molar ratio (Ni, 1994; Neuhaus and Williamson, 2000).
Transfer NOE (TR-NOE)
Most conformational analyses of sugars have been based on NOE data. For the determination of the structures of ligands bound to receptor proteins, transfer NOE (TR-NOE)-based experiments are an excellent tool (Ni, 1994). In suitable conditions, a fast equilibrium in the NMR timescale is established between the free ligand (a small to medium-sized fast-tumbling molecule with a short correlation time) and the ligand bound to the high-molecular-weight receptor (a large and slow-tumbling molecule with a long correlation time) and therefore behaving as the receptor, with a large and fast-growing negative NOE. This averaged situation is reflected in the observation of a large-molecule negative NOE in the signals of the ligand (see Figures 1C,D). In general, modifications of the magnetic field strength or changes in the sample conditions such as temperature or viscosity can be used for fine adjustment if the free ligand is not clearly in the fast-tumbling regime.
In such a case, it is possible to obtain the bound ligand NOE by the analysis of the ligand signals (Ni, 1994; Neuhaus and Williamson, 2000). As a result, the qualitative analysis considering the changes between the free ligand conformation and that of the bound ligand is direct. On the other hand, the quantitative interpretation of the spectrum suffers from the potential effects of spin diffusion. In this situation, complete relaxation matrix calculations (CORCEMA or Mardigras (Borgias and James, 1990; Jayalakshmi and Krishna, 2002) should be used.
In principle, any experiment that measures interprotonic relaxation rates can be used for the analysis of the structure of the bound ligand. Thus, ROESY has been used to improve the experiment, taking advantage of its lower sensitivity to spin diffusion (Poveda et al., 1997). As the ROESY experiment suffers from some spurious effects, transverse-ROESY (Hwang and Shaka, 1992), is better suited as it is less prone to spin-lock artifacts. In addition, the quiet-NOESY experiment has been used. In this experiment, in the middle of the mixing time a biselective signal inversion is applied and a selective band inversion pulse is included, which avoids diffusion through the use of inverted spins (Zwahlen et al., 1994).
Special mention should be made of the monodimensional analogs (Stott et al., 1997). By selecting a single signal, not only is a better resolution achieved due to the change from 2D to 1D but a new dependence on T1sel is created and the NOE exhibits a longer linearity, a feature that is particularly useful when precise distances are required (Hu and Krishnamurthy, 2006).
STD-NMR
The STD-NMR (Mayer and Meyer, 1999) experiment consists of the difference between two experiments. It is undertaken with low-power irradiation during the relaxation delay on a sample in equilibrium that comprises a large excess of ligand(s) (from 10:1 up to 1,000:1) relative to the receptor (a large molecule), which is present at low concentration (nM to μM). In one experiment, recorded as the reference, dummy irradiation far from the signals is performed, while in the other experiment some signals from the receptor are selectively irradiated. If binding occurs, magnetization from the receptor is transferred to the ligand through close contacts with the receptor in the complex, and the effect will emerge in the difference experiment (Mayer and Meyer, 1999; Meyer and Peters, 2003).
Since its formulation, several applications have been explored for STD-NMR. The first application proposed for the technique was ligand screening; this consists of the deconvolution of a library of potential binding molecules by using several compounds at once in each experiment (Henrichsen et al., 1999). The molecules that bind better will show STD signals. Epitope mapping was the next application proposed. The experiment can be used to determine the relative importance of diverse regions of the ligand in the interaction with the receptor (Mayer and Meyer, 2001). The technique can also be used for affinity constant evaluation. The first application in which calculation of a binding constant was described used the Cheng and Prusoff equation, since the STD values are biased by the relaxation properties of the ligand proton being considered (Meyer and Peters, 2003). Latterly, the initial growth rate of the STD affinity factor has been used, in order to avoid this dependence upon the relaxation of the protons, see Figure 1E (Angulo et al., 2010). One of the potential complications in the STD is peak overlap that can make difficult the interpretation of the results. One of the solutions relies on adding another transfer step, then the peaks will be spread in two dimensions: STD-TOCSY (Mayer and Meyer, 1999) and STD-HMQC (Vogtherr and Peters, 2000). Similarly, spectral editing can also be performed using other nuclei 15N or 19F (Kövér et al., 2007; Diercks et al., 2009).
Using methods based on the relaxation matrix, it is possible to calculate the theoretical STD of a given complex. Thus, a quantitative analysis of the ligand within the complex can be carried out using a 3D structure. Based on a 3D model of the complex, using some NMR parameters and considering the saturation time and the irradiation frequency among others, CORCEMA-ST can estimate the different values of saturation transferred from the protein to a particular proton of the ligand (Jayalakshmi and Krishna, 2004, 2005). An example is given in Figure 2. This methodology has also been extended to the case of multiple binding modes within the same binding site (Angulo et al., 2008) and to the iterative refinement of a complex structure (Jayalakshmi and Krishna, 2005).
Figure 2. (A) Detail of the binding-site structure of a mannose pseudodisaccharide bound to DC-SIGN ECD (ExtraCellular Domain), showing the interactions between the Ca2+ ion and the terminal mannose through hydroxyls 3 and 4, and STD growing rates for selected protons; (B) experimental; and (C) CORCEMA-calculated using the structure shown in (A) (Thepaut et al., 2013).
Selective T1 (T1sel) and T2
T1sel and T2 also depend on correlation times and in principle they can be used to study association processes for fast equilibrium between ligand and receptor, observing the ligand signals at low concentration of the receptor. These values can be used for the analysis of binding constants by performing a titration (García-Jiménez et al., 2017). A series of T1sel /T2 values are obtained at different ratios of ligand to receptor, and by fitting these to the variation in the ligand/receptor ratio, KD or IC50 can be extracted. Recently the application of 1H relaxation dispersion measurements has been described for weak binding systems using both T2 and T1ρ (Moschen et al., 2015; Trigo-Mouriño et al., 2017).
WaterLOGSY
The WaterLOGSY experiment (Dalvit et al., 2000) takes into account explicitly the water molecules surrounding the macromolecular complex. It is easy to implement and consists essentially of a selective excitation of the water molecules in a sample prepared in H2O, followed by a NOESY mixing period and a water suppression scheme for detection (Dalvit et al., 2001).
This experiment traces the mechanisms that follow the magnetization from excited water to ligand protons. It can be a direct transfer of magnetization from bulk water to ligand, through exchangeable protons from the receptor or via water molecules buried into the receptor interface. These molecules have tumbling times closer to the receptor than to the bulk water or ligand. Due to the slow correlation times of the bound ligand or the buried water molecules the peaks through the two later mechanisms will have the opposite sign as the direct transfer of magnetization of the free compounds (Stockman and Dalvit, 2002).
The WaterLOGSY competition experiment is used for ligand library deconvolution, in which several relay pathways are used constructively, transferring bulk water magnetization to the ligand in a selective manner (Dalvit et al., 2002). In this experiment, the resonances of non-binding compounds appear with opposite sign, corresponding to NOE in the narrowing limit, and they tend to be weaker than those of the interacting ligands.
New Developments
A recent application of the WaterLOGSY is the LOGSY titration (Geist et al., 2017). From the slopes of the WaterLOGSY vs the protein concentration curves, LOGSY-titration factors can be obtained, that can be used for epitope mapping alternative to STD-AF0. In this case, a great advantage is that isolated water molecules are always in the narrowing limit of the NOE curve.
Complementary information can be also obtained using paramagnetic tags attached to carbohydrates. At least two methods can be exploited: PCS, where changes in chemical shifts are produced and PRE, where additional relaxation is caused. In both cases the effects are proportional to the distance to the paramagnetic center. The paramagnetic center can be either a Lanthanide (Gao et al., 2016; Canales et al., 2017) or nitroxide carrying spin labels like TEMPO.
Integrating the Results. From NMR Data to 3D Structures
The process of deriving a 3D structure from NMR data is not straightforward and it requires the application of force field modeling, with experimental constraints taken from the experimental data. For carbohydrates, accurate distances have to be used and complemented by high-quality theoretical data; this is because an accurate description of the flexibility is required. This often implies the use of explicit solvent MD (molecular dynamics) calculations for the free ligand.
Several force field programs have been developed for carbohydrates, the most widely used being Amber, making use of the parameters for carbohydrates developed by Woods (GLYCAM) (Kirschner et al., 2008); CHARMM (Brooks et al., 2009), optimized by Mackerel; and GROMACS (Hess et al., 2008). These force fields are compatible with other force fields developed for proteins and peptides and can be used in mixed systems calculations.
Using the experimental NMR data from a complex, we could obtain a 3D structure of the ligand within the complex using transfer NOESY evaluation, or at least qualitative data relative to the potential conformational changes upon binding (Ni, 1994). Finally, CORCEMA-ST software is capable of calculating the STD values for individual protons (Krishna and Jayalakshmi, 2008), and of verifying the bound ligand conformation. STD0-AF have been used for the definition of the structure, position, and conformation of the ligand within the complex (Enríquez-Navas et al., 2011). This approach can also be used to tackle more complex situations, such as multiple binding modes in one or more than one site (Angulo et al., 2008). Finally, it has been used to refine NMR structures, using STD-NMR intensity-restrained optimization (Jayalakshmi and Krishna, 2005).
Future Perspectives
NMR provides tools for the analysis and study of low-resolution 3D structures of protein–carbohydrate complexes. A particular advantage is that there are no upper limits to the size of the receptor and the analyses can be performed even with complete cells. In order to improve the reliability of the results, however, efforts should be made to achieve an integrated semi-quantitative method combining all the transfer NMR techniques that can deliver spatial resolution, namely transfer NOESY, STD-NMR, and WaterLOGSY. Recently, promising new methods based on PCS (Pseudo Contact Shifts) have emerged. In the case of the lanthanide paramagnetic tags using PCS they can be applied to protein–carbohydrate complexes with the advantage of the larger dispersion of the signals induced by the lanthanide linked to the carbohydrate (Canales et al., 2017) or the protein (Gao et al., 2016).
Author Contributions
The author confirms being the sole contributor of this work and approved it for publication.
Funding
Research described in this paper has been performed with financial support from MINECO grants (CTQ2012-32605, CTQ2015-70134-P), Junta de Andalucía (FQM-1303), co-financed by European Regional Development Funds (ERDF). The author acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).
Conflict of Interest Statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The author is grateful to his colleagues (J. L. de Paz, J. Rojo, S. Gil-Caballero, M.J. Garcia, and M. Torres) for their careful reading of the manuscript and their participation in the studies cited.
References
Angulo, J., Diaz, I., Reina, J. J., Tabarani, G., Fieschi, F., Rojo, J., et al. (2008). Saturation Transfer Difference (STD) NMR spectroscopy characterization of dual binding mode of a mannose disaccharide to DC-SIGN. ChemBioChem 9, 2225–2227. doi: 10.1002/cbic.200800361
Angulo, J., Enriquez-Navas, P. M., and Nieto, P. M. (2010). Ligand-Receptor Binding Affinities from Saturation Transfer Difference (STD) NMR spectroscopy: the binding isotherm of STD initial growth rates. Chemistry 16, 7803–7812. doi: 10.1002/chem.200903528
Borgias, B. A., and James, T. L. (1990). Mardigras - A procedure for matrix analysis of relaxation for discerming geometry of an aqueous structure. J. Magn. Reson. 87, 475–487.
Brooks, B. R., Brooks, C. L. III, Mackerell, A. D. Jr., Nilsson, L., Petrella, R. J., Roux, B., et al. (2009). CHARMM: the biomolecular simulation program. J. Comput. Chem. 30, 1545–1614. doi: 10.1002/jcc.21287
Canales, A., Angulo, J., Ojeda, R., Bruix, M., Fayos, R., Lozano, R., et al. (2005). Conformational flexibility of a synthetic glycosylaminoglycan bound to a fibroblast growth factor. FGF-1 recognizes both the (1)C(4) and (2)S(O) conformations of a bioactive heparin-like hexasaccharide. J. Am. Chem. Soc. 127, 5778–5779. doi: 10.1021/ja043363y
Canales, A., Boos, I., Perkams, L., Karst, L., Luber, T., Karagiannis, T., et al. (2017). Breaking the limits in analyzing carbohydrate recognition by NMR spectroscopy: resolving branch-selective interaction of a tetra-antennary N-glycan with lectins. Angew. Chem. Int. 56, 14987–14991. doi: 10.1002/anie.201709130
Cremer, D., and Pople, J. A. (1975). General definition of ring puckering coordinates. J. Am. Chem. Soc. 97, 1354–1358. doi: 10.1021/ja00839a011
Dalvit, C., Flocco, M., Knapp, S., Mostardini, M., Perego, R., Stockman, B. J., et al. (2002). High-throughput NMR-based screening with competition binding experiments. J. Am. Chem. Soc. 124, 7702–7709. doi: 10.1021/ja020174b
Dalvit, C., Fogliatto, G., Stewart, A., Veronesi, M., and Stockman, B. (2001). WaterLOGSY as a method for primary NMR screening: practical aspects and range of applicability. J. Biomol. NMR 21, 349–359. doi: 10.1023/A:1013302231549
Dalvit, C., Pevarello, P., Tato, M., Veronesi, M., Vulpetti, A., and Sundstrom, M. (2000). Identification of compounds with binding affinity to proteins via magnetization transfer from bulk water. J. Biomol. NMR 18, 65–68. doi: 10.1023/A:1008354229396
DeMarco, M. L., and Woods, R. J. (2008). Structural Glycobiology: a game of snakes and ladders. Glycobiology 18, 426–440. doi: 10.1093/glycob/cwn026
Diercks, T., Ribeiro, J. P., Cañada, F. J., André, S., Jiménez-Barbero, J., and Gabius, H.-J. (2009). Fluorinated carbohydrates as lectin ligands: versatile Sensors in 19F-detected saturation transfer difference NMR spectroscopy. Chemistry 15, 5666–5668. doi: 10.1002/chem.200900168
Enríquez-Navas, P. M., Marradi, M., Padro, D., Angulo, J., and Penadés, S. (2011). A solution NMR study of the interactions of oligomannosides and the Anti-HIV-1 2G12 antibody reveals distinct binding modes for Branched Ligands*. Chemistry 17, 1547–1560. doi: 10.1002/chem.201002519
Feeney, J. (2000). NMR studies of ligand binding to dihydrofolate reductase. Angew. Chem. 39, 290–312. doi: 10.1002/(SICI)1521-3773(20000117)39:2<290::AID-ANIE290>3.0.CO;2-1
Gao, Q., Chen, C. Y., Zong, C., Wang, S., Ramiah, A., Prabhakar, P., et al. (2016). Structural Aspects of Heparan Sulfate Binding to Robo1-Ig1-2. ACS Chem. Biol. 11, 3106–3113. doi: 10.1021/acschembio.6b00692
García-Jiménez, M. J., Gil-Caballero, S., Canales, Á., Jiménez-Barbero, J., de Paz, J. L., and Nieto, P. M. (2017). Interactions between a heparin trisaccharide library and FGF-1 analyzed by NMR methods. Int. J. Mol. Sci. 18:1293. doi: 10.3390/ijms18061293
Geist, L., Mayer, M., Cockcroft, X.-L., Wolkerstorfer, B., Kessler, D., Engelhardt, H., et al. (2017). Direct NMR probing of hydration shells of protein ligand interfaces and its application to drug design. J. Med. Chem. 60, 8708–8715. doi: 10.1021/acs.jmedchem.7b00845
Henrichsen, D., Ernst, B., Magnani, J. L., Wang, W. T., Meyer, B., and Peters, T. (1999). Bioaffinity NMR spectroscopy: identification of an E-selectin antagonist in a substance mixture by transfer NOE. Angew. Chem. Int. 38, 98–102.
Hess, B., Kutzner, C., van der Spoel, D., and Lindahl, E. (2008). GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theor. Comput. 4, 435–447. doi: 10.1021/ct700301q
Hu, H., and Krishnamurthy, K. (2006). Revisiting the initial rate approximation in kinetic NOE measurements. J. Magn. Reson. 182, 173–177. doi: 10.1016/j.jmr.2006.06.009
Hwang, T. L., and Shaka, A. J. (1992). Cross relaxation without TOCSY - transverse rotating.frame overhauser spectroscopy. J. Am. Chem. Soc. 114, 3157–3159. doi: 10.1021/ja00034a083
Hyde, E. I., Birdsall, B., Roberts, G. C. K., Feeney, J., and Burgen, A. S. V. (1980). Proton nuclear magnetic-resonance saturation transfer studies of coenzyme binding to lactobacillus-casei dihydrofolate-reductase. Biochemistry 19, 3738–3746. doi: 10.1021/bi00557a015
Jayalakshmi, V., and Krishna, N. R. (2002). Complete relaxation and conformational exchange matrix (CORCEMA) analysis of intermolecular saturation transfer effects in reversibly forming ligand-receptor complexes. J. Magn. Reson. 155, 106–118. doi: 10.1006/jmre.2001.2499
Jayalakshmi, V., and Krishna, N. R. (2004). CORCEMA refinement of the bound ligand conformation within the protein binding pocket in reversibly forming weak complexes using STD-NMR intensities. J. Magn. Reson. 168, 36–45. doi: 10.1016/j.jmr.2004.01.017
Jayalakshmi, V., and Krishna, N. R. (2005). Determination of the conformation of trimethoprim in the binding pocket of bovine dihydrofolate reductase from a STD-NMR intensity-restrained CORCEMA-ST optimization. J. Am. Chem. Soc. 127, 14080–14084. doi: 10.1021/ja054192f
Jin, L., Hricovini, M., Deakin, J. A., Lyon, M., and Uhrin, D. (2009). Residual dipolar coupling investigation of a heparin tetrasaccharide confirms the limited effect of flexibility of the iduronic acid on the molecular shape of heparin. Glycobiology 19, 1185–1196. doi: 10.1093/glycob/cwp105
Kirschner, K. N., Yongye, A. B., Tschampel, S. M., Gonzalez-Outeirino, J., Daniels, C. R., Foley, B. L., et al. (2008). GLYCAM06: a generalizable Biomolecular force field. Carbohydrates. J. Comput. Chem. 29, 622–655. doi: 10.1002/jcc.20820
Kövér, K. E., Groves, P., Jiménez-Barbero, J., and Batta, G. (2007). Molecular recognition and screening using a 15N group selective STD NMR method. J. Am. Chem. Soc. 129, 11579–11582. doi: 10.1021/ja073291l
Krishna, N. R., and Jayalakshmi, V. (2008). Quantitative analysis of STD-NMR spectra of reversibly forming ligand-receptor complexes. Top. Curr. Chem. 15–54. doi: 10.1007/128_2007_144
Lepre, C. A., Moore, J. M., and Peng, J. W. (2004). Theory and applications of NMR-based screening in pharmaceutical research. Chem. Rev. 104, 3641–3675. doi: 10.1021/cr030409h
Macura, S., Farmer, B. T., and Brown, L. R. (1986). An Improved method for the determination of cross-relaxation rates from NOE data. J. Magn. Reson. 70, 493–499. doi: 10.1016/0022-2364(86)90143-5
Martin-Pastor, M., and Bush, C. A. (2001). Refined structure of a flexible heptasaccharide using 1H-13C and 1H-1H NMR residual dipolar couplings in concert with NOE and long range scalar coupling constants. J. Biomol. NMR 19, 125–139. doi: 10.1023/A:1008327926009
Mayer, M., and Meyer, B. (1999). Characterization of ligand binding by saturation transfer difference NMR spectroscopy. Angew. Chem. Int. 38, 1784–1788.
Mayer, M., and Meyer, B. (2001). Group epitope mapping by saturation transfer difference NMR to identify segments of a ligand in direct contact with a protein receptor. J. Am. Chem. Soc. 123, 6108–6117. doi: 10.1021/ja0100120
Meyer, B., and Peters, T. (2003). NMR Spectroscopy techniques for screening and identifying ligand binding to protein receptors. Angew. Chem. Int. 42, 864–890. doi: 10.1002/anie.200390233
Moschen, T., Wunderlich, C. H., Spitzer, R., Levic, J., Micura, R., Tollinger, M., et al. (2015). Ligand-detected relaxation dispersion NMR spectroscopy: dynamics of preQ1–RNA Binding. Angew. Chem. Int. 54, 560–563. doi: 10.1002/ange.201409779
Mulloy, B., Forster, M. J., Jones, C., and Davies, D. B. (1993). Nmr and molecular-modeling studies of the solution conformation of heparin. Biochem. J. 293, 849–858. doi: 10.1042/bj2930849
Muñoz-García, J. C., Chabrol, E., Vivès, R. R., Thomas, A., de Paz, J. L., Rojo, J., et al. (2015). Langerin–heparin interaction: two binding sites for small and large ligands as revealed by a combination of NMR spectroscopy and cross-linking mapping experiments. J. Am. Chem. Soc. 137, 4100–4110. doi: 10.1021/ja511529x
Munoz-Garcia, J. C., Solera, C., Carrero, P., de Paz, J. L., Angulo, J., and Nieto, P. M. (2013). 3D structure of a heparin mimetic analogue of a FGF-1 activator. A NMR and molecular modelling study. Organ. Biomol. Chem. 11, 8269–8275. doi: 10.1039/c3ob41789a
Neuhaus, D., and Williamson, M. P. (2000). The Nuclear Overhauser Effect in Structural and Conformational Analysis. New York, NY: Wiley-VCH.
Ni, F. (1994). Recent developments in transferred NOE methods. Prog. Nuclear Magn. Reson. Spectrosc. 26(Pt 6), 517–606. doi: 10.1016/0079-6565(94)90000-0
Poveda, A., Asensio, J. L., MartinPastor, M., and JimenezBarbero, J. (1997). Solution conformation and dynamics of a tetrasaccharide related to the Lewis(x) antigen deduced by NMR relaxation measurements. J. Biomol. NMR 10, 29–43. doi: 10.1023/A:1018395627017
Stockman, B. J., and Dalvit, C. (2002). NMR screening techniques in drug discovery and drug design. Prog. Nuclear Magn. Reson. Spectrosc. 41, 187–231. doi: 10.1016/S0079-6565(02)00049-3
Stott, K., Keeler, J., Van, Q. N., and Shaka, A. J. (1997). One-dimensional NOE experiments using pulsed field gradients. J. Magn. Reson. 125, 302–324. doi: 10.1006/jmre.1997.1110
Thepaut, M., Guzzi, C., Sutkeviciute, I., Sattin, S., Ribeiro-Viana, R., Varga, N., et al. (2013). Structure of a glycomimetic ligand in the carbohydrate recognition domain of C-type Lectin DC-SIGN. Structural requirements for selectivity and ligand design. J. Am. Chem. Soc. 135, 2518–2529. doi: 10.1021/ja3053305
Trigo-Mouriño, P., Griesinger, C., and Lee, D. (2017). Label-free NMR-based dissociation kinetics determination. J. Biomol. NMR 69, 229–235. doi: 10.1007/s10858-017-0150-5
Varki, A. (1993). Biological roles of oligosaccharides - all of the theories are correct. Glycobiology 3, 97–130. doi: 10.1093/glycob/3.2.97
Vogtherr, M., and Peters, T. (2000). Application of NMR based binding assays to identify key hydroxy groups for intermolecular recognition. J. Am. Chem. Soc. 122, 6093–6099. doi: 10.1021/ja0001916
Keywords: NMR, protein–carbohydrate interaction, STD-NMR, transfer-NOESY, transient interactions
Citation: Nieto PM (2018) The Use of NMR to Study Transient Carbohydrate—Protein Interactions. Front. Mol. Biosci. 5:33. doi: 10.3389/fmolb.2018.00033
Received: 08 January 2018; Accepted: 26 March 2018;
Published: 11 April 2018.
Edited by:
Rivka Isaacson, King's College London, United KingdomReviewed by:
Gyula Batta, University of Debrecen, HungaryPeter J. Simpson, Imperial College London, United Kingdom
Copyright © 2018 Nieto. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pedro M. Nieto, pedro.nieto@iiq.csic.es