REVIEW article

Front. Mol. Biosci., 16 March 2021

Sec. Molecular Recognition

Volume 8 - 2021 | https://doi.org/10.3389/fmolb.2021.646808

NMR Characterization of the Interactions Between Glycosaminoglycans and Proteins

  • National Glycoengineering Research Center, Shandong Key Laboratory of Carbohydrate Chemistry and Glycobiology, Shandong University, Qingdao, China

Abstract

Glycosaminoglycans (GAGs) constitute a considerable fraction of the glycoconjugates found on cellular membranes and in the extracellular matrix of virtually all mammalian tissues. The essential role of GAG-protein interactions in the regulation of physiological processes has been recognized for decades. However, the underlying molecular basis of these interactions has only emerged since 1990s. The binding specificity of GAGs is encoded in their primary structures, but ultimately depends on how their functional groups are presented to a protein in the three-dimensional space. This review focuses on the application of NMR spectroscopy on the characterization of the GAG-protein interactions. Examples of interpretation of the complex mechanism and characterization of structural motifs involved in the GAG-protein interactions are given. Selected families of GAG-binding proteins investigated using NMR are also described.

Introduction

Glycosaminoglycans (GAGs) are linear acidic heteropolysaccharides that exist in all mammals and are formed by repeating disaccharide units composed of N-acetyl-hexosamine and hexuronic or hexose (Table 1; ). GAGs can have different sulfation patterns with different charge densities and heterogeneous monosaccharide compositions (). In addition to HA, GAGs are synthesized from the Golgi apparatus in the form of proteoglycans (). According to the disaccharide composition and sulfation pattern, GAGs can be divided into several groups, including heparin/heparan sulfate (HS), chondroitin sulfate (CS)/dermatan sulfate (DS), keratan sulfate (KS) and hyaluronic acid (HA) (). Heparin/HS is composed of repeating disaccharide units of glucosamine (GlcNAc) and glucuronic acid (GlcA) or iduronic acid (IdoA). The initial substrate is [→4)-β-D-GlcA-(1→4)- α-D-GlcNAc-(1→] n. GlcNAc can be substituted by sulfate groups at the amide, 3 or/and 6 hydroxyl groups, and the persulfation can be written as GlcNS3S6S. GlcA can be converted into IdoA by C5 epimerase, and both can be modified by 2-O-sulfation (written as IdoA2S or GlcA2S). CS consists of repeating disaccharide units of glucuronic acid (GlcA) and galactosamine (GalNAc). The initial substrate is [→4)-β-D-GlcA-(1→3)- β-D-GalNAc-(1→] n. CS can undergo sulfation modification similar to heparin except for N-sulfation. However, due to the difference in glycosidic linkage, 3-O-sulfation in heparin becomes 4-O-sulfation. DS is obtained by converting GlcA in CS by C5-epimerase into IdoA. KS consists of repeating disaccharide units of Gal and GlcNAc, both of which can be 6-O-sulfated (). HA is the only GAG that is not modified by sulfation and is not synthesized as proteoglycans. It is composed of repeating disaccharide units of GlcA and GlcNAc. According to the monosaccharide composition and sulfation pattern, GAG disaccharides can have 408 possible compositions ().

TABLE 1

GlycosaminoglycansDegree of sulfation per disaccharide unitMolecular weight rangeTissue distribution
Heparin/Heparan sulfate (HS)Heparin about 1.8∼2.4Heparin about 3∼30 kDaHeparin in liver, lungs and skin;
HS about 0.8∼1.8HS about 10∼100 kDaHS was widely distributed on the cell surface.
Chondroitin sulfate (CS)0.1–1.35∼50 kDacartilage, tendon, aorta, ligament
Dermatan sulfate (DS)< 115∼40 kDaskin, blood vessels, heart valves
Hyaluronic acid (HA)04∼12000 kDasynovial fluid, vitreous humour, ECM of loose connective tissue
Keratan sulfate (KS)< 15∼25 kDaKS I in cornea; KS II in cartilage aggregated; KS III in brain tissue

Structures and tissue distribution of glycosaminoglycans.

As an important component of the extracellular matrix (ECM), GAGs play important roles in the construction of biological systems and the transduction of biological signals (). Signal transduction occurs mainly through the interaction between GAGs and proteins, and these interactions are critical to the biological activity of these proteins. GAGs participate in a variety of physiological processes, including binding, activating and fixing a variety of protein ligands, such as growth factors, cytokines, chemokines, lipoproteins, proteases and their inhibitors, and other ECM components (; ; ). GAGs are also associated with many pathological processes, including degenerative neurological diseases (Alzheimer’s disease), cardiovascular diseases (thrombosis and atherosclerosis) and cancer (; ; ). In the invasion of viruses, GAGs also play roles that cannot be ignored (such as in herpes simplex virus and COVID-19) (). The interaction between GAGs and proteins occurs mainly through electrostatic forces. This puts forward requirements for amino acid sequences in proteins and meets some rules, such as the XBBXBX and XBBBXXBX heparin-binding sequences proposed by Cardin, where B is a basic amino acid and X is any amino acid (). However, long-term research has found that the interaction between GAGs and proteins is not simply determined by the primary structure sequence. A large number of studies have proven that hydrogen bonds and van der Waals forces sometimes even play roles far exceeding electrostatic forces in the interaction; a proper tertiary structure of the protein is also required (). This poses more serious and complex problems for studying the interactions between GAGs and proteins.

The interactions between GAGs and proteins are closely related to many factors, including saccharide unit composition, degree of sulfation, sulfation pattern, chain length, monosaccharide ring conformation and glycosidic linkage. The research methods used to characterize the interaction between GAGs and proteins mainly include gel electrophoresis (GE) (), affinity chromatography (AC) (), surface plasmon resonance (SPR) (), biological layer interferometry (BLI) (), isothermal titration (ITC) (), microarray methods (), crystal diffraction methods (X-ray) (), mass spectrometry (MS) (), and nuclear magnetic resonance spectroscopy (NMR) (). NMR is an insensitive technique compared with other analytical method for the study of interactions between biomolecules. The amount of sample needs to be in milligrams with high purity. In the study of proteins, NMR can characterize a protein with a molecular weight around 20 KD very well. However, proteins need to be isotope labeled by 15N and/or 13C when the molecular weight increases and can be studied up to 100 KD. The cross peaks will become broadening and overlapped severely for larger proteins. Even with the above limitations, NMR is still an irreplaceable technique in the characterization of the biomolecule interactions at the atomic level especially in the case of glycosamionoglycans. Both X-ray diffraction and NMR can provide more precise tertiary structure information, and they do not require sample derivatization and will not cause structural damage to the sample during the experiment. Due to the accuracy and refinement of the data, both types of data can be used for model construction. However, X-ray diffraction studies a crystal in solid state and provide only few conformations of the interaction. While, NMR studies a solution under physiological condition and records dynamic conformations during the whole interaction period. Glycosaminoglycans are very hard to obtain a crystal due to their high flexibilities and exchangeable conformations. The solution NMR can not only show the natural state of the complex, but also detect the change of the complex conformation on the ns-ms time scale (). Compared with the immobilization study of crystal diffraction, solution NMR can also be used for the dynamic study of interactions under physiological conditions.

Nuclear magnetic resonance is widely used to study the conformation of GAGs alone or in complex with proteins (), but the information usually obtained indicates that there are multiple GAGs or complex structures in solution. According to NMR data, GAGs present different folds configurations in solution according to their type and environment (), such as the controversial 3-folds and 4-folds coexisting left-handed helix of HA (), which will directly affect the distribution of acidic groups in space. Generally speaking, the conformational changes of GAGs are mainly caused by two factors, one is the ring conformation of monosaccharides, and the second is the flexibility of the glycosidic linkages (). The conformation of the IdoA residue in heparin, HS and DS is different from that of the other three monosaccharides (GlcNAc, GalNAc, and GlcA). IdoA exist in the conformational equilibrium, with two chairs (1C4 and 4C1) and one shew-boat (2S0), instead of the fixed conformation 4C1 adopted in GalNAc, GlcNAc, or GlcA (). This gives these three different types of GAGs more flexible and various protein binding activities. This balance is affected by chain length, the degree of sulfation of adjacent monosaccharides, and its own 2-O sulfation (). When interacting with proteins, the conformational balance of IdoA will be tilted, such as binding to fibroblast growth factor-2 (FGF2), fibroblast growth factor-2 receptor (FGF2R), and eosinophil cationic protein (Ecp) (; ). In free state, when the conformational balance ratio is closer to the required binding state, the binding affinity is stronger (). Conversely, when the required conformation of the bound state cannot be achieved, the activity may be completely lost. But even if the protein has a clear tendency to a certain conformation of IdoA, there will generally be a conformational balance. The binding of AT III to heparin requires an absolute 2S0 conformation, but according to the NMR structure information, there is negligible 1C4 conformation in the whole binding process (). Even though IdoA brings more variable binding conformational selectivity, recent studies have shown that GlcA has a better effect on the overall conformation of GAGs (). In order to adapt to the ECM environment, the angle of the glycosidic linkages is allowed to change to a certain extent. The angle of the glycosidic linkages is affected by temperature, and the increase in temperature will result in a transition to the higher energy state (). When interacting with proteins, the glycosidic linkages can adopt proper orientations to meet the structural requirements during binding to proteins, and even cause the kinking of the GAGs polymer chain, thereby further enhancing the binding affinity (). Compared with the obvious conformational equilibrium of IdoA, sometimes GAGs have α/β isomeric equilibrium at the reducing end () and rapid intramolecular hydrogen bond exchange (). Due to the flexibility of GAGs, there may be multiple interaction modes at the same binding domain in the GAG-protein interaction process (). In the interaction between GAGs and proteins, the structure of the proteins is normally changed or stabilized. The weak interaction between GAGs and proteins undergoes on the ns-ms time scale, so the conformation of the protein in the system will change over time. Due to the structural heterogeneity and conformational flexibility of GAGs or the dynamic changes of the complex, it is also very difficult to construct a model of complexes in solution ().

Solution NMR can provide information about conformational changes and kinetic data during interactions between proteins and GAGs (). NMR can also reveal the effects of different temperatures, pH values, salt concentrations, and ligand concentrations on the binding activity. There are three main goals in using NMR to study GAG-protein interactions: the first is to detect the amino acids involved in binding from the perspective of proteins, the second is to analyze the saccharide and its groups involved in binding from the perspective of GAGs, and the third is to observe the conformational changes and kinetic information during binding from the perspective of the interaction. To achieve these three goals, three technologies, chemical shift perturbation (CSP), saturation transfer difference (STD), and exchange-transferred nuclear Overhauser effect (trNOE), are initially used (), while other technologies, such as saturation transfer double difference (STDD) (), paramagnetic relaxation enhancement (PRE) (), pseudocontact shifts (PCS) (), and exchange-transferred rotating-frame Overhauser effect (ROE), have been developed to compensate for the shortcomings of the former. The latest pulse sequences have been developed to provide a more detailed and accurate description of the binding process, such as the gradient spectroscopic observation of water ligands (waterLOGSY) () and heteronuclear in-phase single quantum coherence experiment (HISQC) (). In addition, solid-state NMR has also been applied to study interactions involving ligands with low solubility (; ). These techniques are based on four types of data: nuclear Overhauser effect (NOE), scalar coupling (J), residual dipole coupling (RDC) and chemical shift anisotropy (CSA). The purpose of this paper is to introduce some important findings of the application of NMR to the study of the interactions between GAGs and proteins (Table 2) and the review is classified according to the type of GAGs.

TABLE 2

Type of proteinName of proteinType of GAGparticipating binding residues and Secondary structureAffinity (Kd)References
ChemokineCCL5Heparin40S loop (R44KNR47), α helix (K55, K56)18 μM
CS40S loop(R44KNR47), N loop (R17, L19, I15)0.25 μM
CXCL1Heparin/HSN terminus(R8), N-loop (H19, K21),40S turn (K45, R48), β3-strand (R49), C-helix (K60, K61, K65)50 μM; ;
CS/DSN terminus(R8), N-loop (H19, K21), 40S turn (R48),4 μM
CXCL2HeparinN-loop (R17, K21), 40S turn (K45), C-helix (K61, K65, K69)25 μM
CXCL5Heparin/HSN-loop (H23, K25), 40S turn (K49), β3 strand (K52), C-helix (Lys64, Lys65, Lys69, Lys76)30 μM;
CS/DSN-loop (H23, K25), 40S turn (K49), β3 strand (K52)3 μM
CXCL7HeparinN-loop (H15, K17), β3-strand(R44), C-helix (R54, K57, K61)
CXCL8HeparinN-loop (K15, H18, K20, K23), C -helix (R60, K64, R68), β3-strand (R47), 50S loop (K54)μM
CXCL11HeparinC- helix (K57SKQAR62)
CXCL12HeparinC-helix (R12, A40), 20S loop (K24), 40S loop (N46)μM
CXCL13HSC-helix (K60, R64, R67, H68), C-loop (K84, R85, R86)19 nM
CXCL14Heparin10S loop(I12), β2-strand (I36, T37), 40S loop(K54), C -helix (R72),
CS/DS10S loop(I12), 40S loop(K54), C -helix (R72)
Growth factorFGF1Heparinβ1–β2 loop (N18), β8–β9 loop (N92), β10–β11 loop (K113), β11 strand (K118), β11–β12 loop (Q127, K128)nM
FGF2Heparinβ1 strand (K27), β1–β2 loop (N28), β8–β9 loop (N102), β10–β11 loop (R121), β11 strand (K126), β11–β12 loop (Q135, K136)nM
FGF7Heparinβ3 strand (R18), 40s loop(N92), β10(N114), 110s loop(Q115), 120s loop (V120, K124, Q129, K130, T131)
SerpinAT IIIHeparinN-terminal end (K11, R13), A helix (R46, R47), D helix (K114, K125, R129, R132, K133, K136)20 nM
Type II cytokinesIL-10Heparin/CS/DSD helix (K99, R102, R104, R106), 110S loop (R107, K117, K119)0.41 mM
IFNγHeparinC-terminal end (D1:K125TGKRKR131, D2:R137GRR140)1.63 nM
Roundabout 1HS80s loop (K81), 130s loop (V133, H134, G135, R136, K137), βA strand (I167, R169)
CytokinePleiotrophinCSC-terminal TSR domain β-sheet (K60, K61, K69, K91, K92, K84, K86, K107)90 μM
Link proteinCD44HAβ1 strand (K38), 40s loop (R41, Y42), 70s loop (R78, Y79), 90s loop (N100, N101), 150s loop (R150), β9 strand (R154), 160s loop(R162)μM
TSG-6HA10s loop (K11, Y12), 40s loop (H45, C47), β3 strand (A49), β3 strand (Y59), 60s loop (V62, K63), 80s loop (Y78, R81)μM
Viral pathogenviral CCL2Heparin10s loop(R18), 40s loop (K45, R46, R48)113 mM
DefensinsHuman β-defensin 2Heparin/DS20s loop(R22RYK25), β3 strand (K39), 40s loop(K40)5 mM
RNase AEosinophil cationic proteinHeparinα1 helix (R7, Q14, H15), β1 strand (Q40), loop4(H64), β6 strand(H128)15 μM

GAG binding proteins.

Heparin/Heparan Sulfate

Heparin is the most negatively charged polymer found in nature, and it is also the most studied in the GAG family (). One way to distinguish between heparin and HS is based on whether the mature body is still connected to the core protein. HS will be secreted out of the cell in the form of glycoproteins, most of which are fixed on the cell membrane to mediate many intercellular signaling pathways. Heparin is cleaved by β-endoglucuronidase and is combined with alkaline protease in the form of oligosaccharide chains to be stored in secretory granules (). The binding of heparin to protein mostly relies on its own high electronegativity and the positively charged domains in the protein. Hydrogen bonds and van der Waals forces also play important roles in the binding process. Moreover, the binding of heparin and protein is sometimes ion-dependent. For example, the binding of Langerin and heparin is mainly Ca2 + -dependent, although there are additional non-Ca2 + -dependent binding sites (; ; ). HS can be divided into a high-sulfation domain (NS domain) and a low-sulfation domain (NA domain). Heparin essentially contains all possible sulfation modification structures of the NS domain due to the degree of high sulfation. Most of the biological functions of HS are concentrated in the NS domain, although the NA domain is more flexible and more suitable for bending. Due to the early large-scale clinical application of heparin, it was relatively easy to obtain. Early research mainly used heparin as a substitute for HS to carry out functional and structural studies. In approximately the past thirty years, the study of the interaction between heparin and various proteins has become a hot spot, and the gradual maturity of chemical enzyme synthesis has given this field new vitality. Heparin can induce the oligomerization or heteromerization of proteins, which can prevent proteins from being hydrolyzed by protein-degrading enzymes and increase or decrease the possibility of their binding to receptors.

Antithrombin III (AT III) is an absolutely conserved serine protease with two different glycosylation forms (α, β), consisting of three β-sheets (A-C) and nine α-helices (A-I) (). Heparin is a cofactor of the antithrombin-mediated coagulation cascade, and the interaction between them directly affects the activities of factors IXa, Xa and IIa (). Choay, J used chemical enzymatic synthesis of various heparin-related oligosaccharides to determine that the minimum specific sequence required for binding to AT III was the pentasaccharide A1GA2IA3 (Figure 1), which is also the only specific recognition sequence for heparin and protein binding found thus far (; ). Although the specific pentasaccharide can meet the requirement of binding to AT III, it can only inhibit the activity of Xa. Inhibiting thrombin activity requires a heparin chain containing more than 16 saccharides, which can form a ternary complex with antithrombin and thrombin (). The interaction between heparin and AT III was described as a three-state, two-step kinetic process (Figure 2; ), which assumed that AT III was in a balance of ‘native unactivated,’ ‘intermediate-activated’ and ‘fully activated’ states under physiological conditions (). First, A1GA2 was driven by K125 and K114 to combine with the C- terminus of helix D in “native unactivated” AT III, and the reducing end faced the N-terminus (). Then, accompanied by conformational changes in AT III (helix D extension, reactive center loop exposure, and closure of sheet A) and heparin (IdoA from equilibrium conformation between 1C4 and 2S0 to complete 2S0), each unit in the pentasaccharide was further combined with AT III (). The combined complex can interact with the target protease or enzymatically decompose, and heparin is dissociated accordingly. In the electrostatic binding of heparin and AT III, several sulfate groups of heparin-specific pentasaccharide (N-SO3 for A2 and A3, 6-O-SO3 for A1, and 3-O-SO3 for A2) and carboxyl groups were irreplaceable ().

FIGURE 1

FIGURE 2

). The models of the three states are derived from X-ray. The reactive center loop (RCL) (red), sheet A (green), and helix D (gray blue) and the helix extension (dark blue) are highlighted in each state.

Further research using NMR focused on the specific role of each monosaccharide in the binding of heparin to AT III and the effect of extended pentasaccharide on the binding. The ratio of the 2S0 conformation in IdoA in the A1GA2IA3 sequence was 20% higher than that in the general heparin sequence (). In the three different chemically synthesized heparin pentasaccharides, the pentasaccharide had anticoagulant activity only when IdoA was in 2S0 (). Therefore, the proportion of 2S0 of IdoA in the heparin pentasaccharide sequence was one of the factors affecting the binding rate, which was affected by the degree of sulfation of glucosamine on both sides and its own 2-O-SO3 (). Although the absence of 2-O-SO3 in IdoA had no significant effect on the binding conformation, it resulted in a decrease in the proportion of free state 2S0 and a two-fold decrease in affinity (). At the same time, the flexibility of IdoA provided unlimited possibilities for the binding of heparin to protein. A recent study used IdoA2S instead of GlcA in the AT III binding sequence (). The results showed that IdoA2S, which replaced GlcA, was in a pure 1C4 conformation when bound, and the affinity was tripled, which provided a basis for the application of bovine heparin. The unique structure of bovine heparin also provided unique ideas for the study of the specific mechanism of anticoagulation between heparin and AT III (). The 3-O-SO3 and 6-O-SO3 also had significant effects on the conformational balance of IdoA (; ). The contribution of A2’s 3-O-SO3 to binding was in not only the conformation of heparin but also the formation of ‘intermediate-activated’ AT III (; ). Octa-7 (Figure 1), an octasaccharide with extended reducing end, showed that adding an extra 3-O-SO3 to the A3 would increase the ratio of 2S0 in I by approximately 15%. The additional 3-O-SO3 formed new ionic bonds with R46 and R47. The extended disaccharide also had a certain contribution to the binding (by interacting with E113 and R24), and the binding force of otcasccharide and AT III was 40% higher than that of the specific pentasaccharide sequence and AT III. In the binding state, I and extended nonreducing end IdoA2S was completely in 2S0. In a similar structure (OCTA-1), due to the lack of 3-O-SO3 in the reducing end of A3, the extended IdoA2S was completely in 1C4 when bound, resulting in a substantial decrease in affinity (). When extended reducing end IdoA2S’s 2-O-SO3 was removed (OCTA-2), the affinity increased slightly (). In addition, there was little interaction between the reducing end extended disaccharide and AT III. In the other two octasaccharides with GlcA or IdoA as the extended nonreducing end (OCTA-3, OCTA-4), there was a significant polarization of affinity. The affinity of octasaccharide with GlcA as the nonreducing end was one order of magnitude higher than that with IdoA, which was in pure 2S0. In recent years, the appearance of low-molecular-weight heparin has become a research hotspot due to its unique fragments produced by cleavage or hydrolysis on anticoagulation. In Guerrini’s study, the affinity of two octasaccharides (OCTA-5, OCTA-6) containing specific pentasaccharide sequences derived from enoxaparin in binding with AT III decreased by 60-fold compared with the hexasaccharide with a complete pentasaccharide sequence. Because of the special pentasaccharide unit, the binding of the reducing end became weaker (). The interaction difference of the octasaccharides with AT III showed that the substitution of different groups on heparin not only affected the binding strength with AT III but also changed the conformation during binding.

Heparin plays a key role in the regional aggregation and oligomerization of fibroblast growth factor (FGF), protecting it from denaturation and degradation and inducing its binding to the receptor (FGFR) (). FGF is a growth factor family with 23 members, and its structure is highly related (12 β strands form the classic β-trefoil structure) (). The receptor proteins of FGF include four categories (FGFR1-4), which are composed of three immunoglobulin (Ig)-like domains, which can be subdivided into seven categories according to the difference in Ig3 (). FGFR Ig2 is a key site for the binding of FGF and FGFR mediated by heparin (). In the study of the effect of FGF and heparin, acidic fibroblast growth factor (aFGF, FGF1) and basic fibroblast growth factor (bFGF, FGF2) were the most classic models (). Studies have shown that the binding of heparin to FGF does not change the FGF conformation, and the binding domain is mainly located at the β1-2 and β10-11 strands (). Although there is clear evidence in the study of Crystallography, in the free state, 116-120 (131-136) of FGF1 (FGF2) constitute βXI structure (). However, Moy’s NMR study on the structure of FGF2 in solution showed that there was no evidence to prove the existence of βXI (). It is speculated that this is the structural change caused by the combination with HSPG, and this change is very important for the combination. This was confirmed in the subsequent NMR structural study of FGF1, Ogura pointed out that in the binding state, the 116-120 sequence has an obvious tendency of β-chain structure (). In addition, K125 in FGF2 and K118 in FGF1 had high affinity in binding with heparin. Therefore, the β11 chain was considered to be the key structure for the binding of FGF to heparin. In the combination of FGF2 and heparin, 2-O-SO3 and N-SO3 were necessary (), and additional 6-O-SO3 was required for FGF1 (). However, in the study using 48 kinds of heparin disaccharides to bind FGF1, 3-O-SO3 provided a stronger binding ability, and further C6 sulfation seemed to have a negative effect on the binding (). In the study of the binding of heparin to FGF, 1C4 might have been the more favorable conformation (; ). Interestingly, a recent study showed that specific AT-binding sequences can bind to FGFR2 Ig2 as a high-affinity complex, and IdoA remained in a high proportion of 2S0 (). Some experiments have shown that the combination of FGF and heparin seem to require a certain regular sequence of monosaccharide units or a special sulfation pattern (). The mirror image of the carbohydrate structure also caused a significant reduction or loss of activity (). For FGF1, only a single 6-sulfated tetrasaccharide was needed to induce its dimerization (). However, for FGF2 to be fully activated, heparin fragments of approximately decasaccharide might be required (), although there was also evidence that tetrasaccharides could induce FGF2 dimerization (). Heparin can induce FGF dimerization, but whether it is a critical step is controversial. Some NMR data showed that heparin, which formed a high-affinity complex with FGF, did not induce the dimerization of FGF but still had high activity ().

In the study of the FGF-FGFR-heparin binding model (Figure 3), the crystal study gave two hypotheses: a 2:2:1 trans-binding model and a 2:2:2 cis-binding model (). NMR research in recent years has explained the formation process of the 2:2:2 model. Nieto used FGF1 and FGFR2 Ig2 and two heparin oligosaccharides to study the mechanism (). In the activity experiment, FGF1 and FGF2 had different requirements for heparin. In deheparinized cells, FGF2 activity was completely lost. However, after pretreatment of the cells with heparin, the activity recovered. FGF1 requires the presence of an additional heparin-like stabilizer myo-inositol hexasulfate (MIHS). It is speculated that the role of heparin in FGF1 was not limited to mediating the binding of FGF and FGFR. There was a second binding site in the FGF-FGFR complex, which was a clear cis-dimer binding model mark. Subsequent speculation suggested that the signaling pathway should be regarded as follows: FGFR dimerization was initially induced by GAGs, and then FGF and the ternary complex formed a higher-order aggregate and activated the subsequent enzyme cascade. Schieborr investigated the interactions among FGF1/FGF2, FGFR4 Ig2, and three different heparin polysaccharides (). The experimental results showed that the hexasaccharide could meet all the binding site requirements for inducing FGF dimerization, but the stability of the resulting complex was extremely poor. STD experiments showed that the combination of octasaccharide and FGF2 had a positive synergistic effect, but due to the lack of heparin structure data, the exact mechanism needs further experimental verification. Heparin was proven to have an extremely low dimerization ability for inducing FGFR4 Ig2, which was clear proof of the trans-dimer model in the description by . However, the NMR data suggested there was a secondary binding site in the FGF-FGF Ig2 complex, which was again a clear cis-dimer binding model. Schieborr proposed that hexasaccharides and octasaccharide could mediate FGF2 signaling pathways under different mechanisms, and the positive synergistic effect of octasaccharide was due to the different residues involved in the binding. However, while there should theoretically be an FGF/FGFR/heparin 4:2:2 complex in the pathway, there were no data to support its existence. The existence of the FGF/FGFR/heparin 2:2:1 model was clearly supported by Brown’s ITC data, but no NMR evidence was obtained ().

FIGURE 3

CXCL12 has six different splicing variants (CXCL12α-φ) in humans and is the only CXC chemokine with differential gene splicing (). The complex of CXCL12 and the receptor CXCR4 mediates many physiological functions, including physiological processes such as hematopoiesis, embryonic development, vascular repair, and inflammation (). CD26, a leukocyte-activating antigen, can be cleaved CXCL12 between the N-terminal P2 and V3 residues (). The cleaved product has a reduced affinity for CXCR4 and cannot activate it any more. Research on the binding domain of CXCL12 and heparin/HS can be traced back to 1999. The K24HLK27 base sequence in the β1-strand of the β-sheet, conforming to the BBXB rule, was verified in a mutation experiment (). Sadir believed that R41 and R43 in the β2 strand were additional binding sites, in addition to K1 at the N-terminus as a potential binding site (). The binding between heparin/HS and K1 in CXCL12 was believed to protect CXCL12 from being cleaved by CD26 (). Murphy first used X-ray crystallography to study the interaction between CXCL12 and heparin/HS and proposed two binding domains in CXCL12: one at the interface of the dimer and the other in the N-loop region and the N-terminal helix similar to the binding domain in CXCL8 (). Using 13C-labeled octasaccharides in the NMR experiment, Laguri determined that the heparin-binding sequence was related to the GlcN-3, GlcA-4, and GlcN-5 units of the octasaccharides (). N-sulfation and 6-O-sulfation are essential for binding. The nonreducing end monosaccharide and reducing end disaccharide of the octasaccharide formed additional contact with the N-terminus of CXCL12 (R8 and R12 are the most prominent), and a consistent molecular binding model was constructed. However, Ziarek proposed a controversial molecular model (). He believed that heparin and two CXCL12 molecules should drive the formation of the polymer in an almost orthogonal conformation, instead of the previously proposed interface of two CXCL12 molecules (composed of a β1 strand and the N-terminus). The data indicated that the binding site in CXCL12 should be on the six-strand of the β-sheet, while the N-terminus was not involved. The main residues involved in binding included K20, K24, K27, K41, K43 and R47, while A8 and A12 provided additional binding. It was proposed that the reason why heparin protected CXCL12 from CD26 cleavage was not the preemptive combination but the coverage of K1 caused by dimerization. Panitz’s study proved that the interaction affinity between heparin and CXCL12 was much higher than that of other GAGs, and the degree of sulfation was not the only factor influencing the binding (). The binding sites in CXCL12 with other GAGs were similar to heparin, with the exception of a second binding site for CS compared to heparin (R20, A21, N30, K64).

Type II cytokines have six secondary structure elements (A-F) to form an α-helical structure, of which A, C, D, and F adopt the classic four-helix topology, while B and E exist as the connecting structure (). Interleukin-10 (IL-10), interferon-γ (IFNγ) and interleukin-26 (IL-26) are the three proteins in this family that exist in the form of dimers. Although IL-10 and IFN-γ had the same protein folding mode, their binding with heparin split into two completely different manners. STD data indicated that when IL-10 bound to heparin, the degree of sulfation rather than the site had a greater impact on the binding (), although the effect of 6-O-SO3 on affinity was 2-3 times greater than the effects of N-SO3 and 2-O-SO3. Data showed that there was a hydrogen bond or strong van der Waals force between IL-10 and the methyl group in the N-acetyl residue of the saccharides. As the heparin chain length increases, the affinity increases. When the chain length reached eight sugars, the affinity suddenly increased. It was calculated using STD data that when IL-10 bound to a heparin oligosaccharide with more than eight sugars, the Hill coefficient was approximately 2. This indicated that heparin and each monomer of the IL-10 dimer were bound, and the binding was synergistically positive. It was speculated that the binding site in IL-10 was located at the C-terminus of the D helix and the basic amino acid cluster L101RLRLRRCHRF111 of the adjacent DE loop. This heparin-binding domain existed in both monomers, which also supported the positive synergistic combination of octasaccharide and IL-10. NOE data showed that the conformation of a tetrasaccharide in the binding center did not change much. Further PCS data confirmed that the binding domain of IL-10 with heparin was in the 101-111 basic amino acid cluster (). This domain is absolutely conserved in IL-10 from various sources, and it is also located in the binding domain of IL-10R2 and IL-10. The reason why GAG had an inhibitory effect on IL-10 might be due to the low-affinity IL-10R2 competing with heparin for binding.

Unlike IL-10, the binding domain of IFN-γ with heparin was located at the C-terminus. IFN-γ had four clusters of enriched basic amino acids, but only two C-terminal domains, K125-R131 (D1) and R137-R140 (D2), interacted with heparin (). NOE data showed that the interaction between the protein and heparin had no effect on the conformation of the protein, and only the electrostatic force contributed to the binding without any other interaction force. The increase in sugar chain length increased not only the affinity between heparin and IFNγ but also the bending degree of the whole sugar chain. The binding of IFNγ to heparin protected the D1 domain from protease hydrolysis, and D1 acts as the main binding domain to heparin. ITC experiments have shown that D2 is not necessary for the binding of IFNγ to heparin, but removing D2 will increase the binding of IFNγ to heparin (). Further studies have shown that the combination of D1 with heparin was mainly a thermodynamic process, while the combination of D2 with heparin was a kinetic process (). The main function of D2 was to strengthen the binding of IFNγ with heparin. The binding of the C-terminus of IFNγ to heparin is a two-step process. First, D1 bound to heparin, and the binding site was oriented. Then, D2 combined with heparin to strengthen the binding. The binding of IFNγ to its receptor includes two domains, one of which is the C-terminus. Therefore, HSPG on the cell surface competed with the IFNγ receptor for binding; and the addition of exogenous heparin could also reduce the IFNγ concentration on the cell surface. The inhibitory effect of heparin on the activity of certain proteins might be due to its competition with the protein receptor for binding, which led to the decreased or even disappearance of the binding affinity between the receptor and the protein. IL-10 inhibits the activity of IFNγ, so its mechanism might be more complicated. Studying the interaction between GAGs and proteins of a specific sequence may help to develop a more thorough understanding of the mechanism.

Chondroitin Sulfate

According to the type of uronic acid and sulfation, common CS can be divided into five categories: nonsulfated chondroitin sulfate (CS-O), 4-O-sulfated chondroitin sulfate (CS-A), 6-O-sulfated chondroitin sulfate (CS-C), 2, 4-O-disulfated chondroitin sulfate (CS-D), and 4,6-O-disulfated chondroitin sulfate (CS-E) (). CS-B (DS) has all of the sulfation modification types of the above five types of CS, but its uronic acid is epimerized into IdoA. Oversulfated chondroitin sulfate (OSCS) was sulfated at all sites that could be sulfated, and it was one of the culprits that triggered the “heparin crisis” in 2008 (). There is a special kind of 3-O-sulfated chondroitin sulfate (CS-K) in marine organisms that has a high affinity for growth factors ().

In the interaction with chemokines, the main function of GAG was to locally aggregate chemokines to increase their binding to G-coupled protein receptors and to form a concentration gradient required for the migration of leukocytes, among which HS was dominant (). However, CS also played an important role in the interaction with certain chemokines, such as the chemokine CCL5 (regulated upon activation of normal T cell expressed and secreted factor, RANTES). CS plays an important role in a variety of biological pathways mediated by CCL5, such as inducing T cell apoptosis and monocyte blockade. Deshauer studied the interaction between two CS hexasaccharides and CCL5 and used TEMPO to label CS for PRE experiments to study the binding sites in depth (). In the titration of CCL5 with CS444 (GlcA-GalNAc4S- GlcA-GalNAc4S- GlcA-GalNAc4S), there were obvious chemical shift changes in the 40S loop, the N-terminus and the N loop (Figure 4). At a ratio of 1:1, the chemical shift had no significant change. When CS644 (GlcA-GalNAc6S- GlcA-GalNAc4S- GlcA-GalNAc4S) is used for titration, there are only small chemical shift disturbances at these three binding sites. However, when the ratio of CS644:CCL5 was more than 1:1, R17 and L19 in the N loop showed obvious chemical shift disturbances. In the PRE experiment, CS444 data showed that its reducing end was close to the 40S loop BBXB sequences. However, CS644 had additional chemical shift changes at Y3, A16, and R21, indicating that CS644 was also close to the 20S loop, N-loop and N-terminus, which suggested that the combination of CS644 and CCL5 was more heterogeneous. It can be seen that the type of GAG, the degree of sulfation and the ring conformation had a huge influence on the binding conformation between GAG and protein, which was also reflected in Pichert’s CXCL8 and CS hexasaccharide interaction study ().

FIGURE 4

Midkine (MK) and pleiotropic protein (PTN) form the MK/PTN cytokine family, which is a heparin-binding nerve growth factor. They are highly similar in structure and share more than 50% of the amino acid sequence (). They consist of two TSR domains with a hinge connection. Each domain consists of three antiparallel β-strands to form β-sheets. The C-terminal domain (CTD) of PTN is the main CS-binding domain, which has an affinity far greater than that of the N-terminal domain (NTD) (). CTD has two basic residue clusters (cluster 1: K69, K91, K92 and cluster 2: K84, K86, K107). The electrostatic potential diagram showed that the two sides of the β-sheet can be coplanar. According to the PRE data, CS-A preferred cluster 2, while CS-E preferred cluster 1. The data showed that K54 in NTD was close to the paramagnetic center, but NTD had only a few residues with side chains and HN atom transfer perturbation. The hydrophobic hinge can arrange two lysines (K60 and K61) near CTD cluster 1 to participate in the binding of CS. Although there was no clear reason to prove the effect of the C-terminus on the binding of CS to PTN, the affinity of CS-A, but not CS-E, to the C-terminal truncated PTN was greatly reduced. CS-E had a greater affinity than CS-A, which might be the reason why the PTN/MK family was associated with many tumorous inflammations (). Unlike PTN, according to STD data, CS-E can simultaneously bind to the two domains of the midkine ().

Tumor necrosis factor-stimulated gene-6 (TSG-6) is a classic HA-binding protein that shows different binding modes with CS compared to HA (). The combination of CS and Link- TSG-6 had at least two binding sites, and 4-O-sulfation was preferred. The slow exchange site was similar to the HA-binding site, but there were still some differences due to the sulfation pattern of CS. STD data indicated that there was a second group of rapid exchange binding sites, which were close to the heparin-binding site according to the model based on PRE data. The change in the relaxation rate ratio R2/R1 indicated that the initial combination of CS and Link- TSG-6 can induce dimerization. The dimerization interface and the CS binding site were located on opposite sides, so CS plays a neutralizing role rather than functioning as a bridge in inducing dimerization.

Dermatan Sulfate

Although DS was similar in structure to CS, the existence of IdoA gave it unparalleled structural flexibility. For example, in combination with hepatocyte growth factor/scattering factor (HGF/SF), the presence or absence of IdoA was the key to the combination of GAG with HGF/SF (). The binding mode of DS and NK1 (HGF/SF heparin-binding domain) was similar to that of heparin, although the affinity was slightly lower. The binding was concentrated in the N domain. Although crystallographic data proved that the K1 domain was involved in binding, this binding was based on the premise of dimerization. However, the NMR data showed that in solution, the low-molecular-weight GAGs would not induce its dimerization.

Sepuru used medium-length GAG to study the interaction with CXCL1 or CXCL5 in the presence of monomers and dimers through CSP experiments (). The two binding sites in CXCL1 with HS were on the opposite sides of the protein, the α-domain (H19, K21, K45, K60, K61, K65) and the β-domain (R8, K29, R48, K49). The results showed that CXCL1 and HS were combined in a ratio of 1:2, and ITC experiments verified this result. The binding sites of CXCL1 with CS and DS are located in the γ-domain (R8, H19, K21, K45, K49). The binding domain of CXCL5 with GAG was similar to that of CXCL1, but there was no obvious specificity for GAG species. Neither CXCL1 nor CXCL5 bound to GAG involved helices, which was different from the previous proposal that helices are an important binding site for the interaction of chemokines that activate CXCR2 with GAG. In the HADDOCK model, the interaction between DS and CXCL1 involved two sulfate groups, two carboxyl groups and two N-acetyl groups, and the interaction model with CXCL5 involved two sulfate groups, one N-acetyl and one hydroxyl group. The molecular docking models of CS and DS with different structures were quite different. They involved different residue-binding groups and positions. This was consistent with the differences in the interaction morphology of GAG with different structures proposed previously. This was also reflected in the combination of CXCL14 and DS (). The binding of DS and heparin with CXCL14 occurred in the C-terminal helix, part of the N-terminus and the transition between the second and third β-sheets (Y44-Q47). However, the maximum perturbation in the combination of DS and CXCL14 was associated with R72, while I36 and T37 were more affected in terms of heparin. DS and CS also had significant differences in N-terminal disturbances. The interaction between DS and protein was also dependent on chain length and sulfation pattern. In the study of the interaction between tau protein and DS, tau was favored for 6-O-sulfation (). Disulfated DS had a higher affinity than monosulfated DS, although the affinity of both was less than that of heparin.

Decorin binding protein B (DBPB) bound to DS in a different binding mode than DBPA, mainly through the linker between helices 1 and 2, the C-terminal tail, and the alkaline patch (). In the PRE experiment, there were no clear data indicating that the C-terminal tail was involved in binding. It was speculated that this was because the binding occurs at the nonreducing end of DS, while the TEMPO label was at the reducing end of DS. The mutation data showed that the three sites all had a promoting effect on binding, and the C-terminus played a key role in binding. The most obvious difference between DBPB and DBPA was only the C-terminal disulfide bond, which again emphasized the influence of protein structure on binding. Due to the lack of disulfide bonds, the C-terminus could exist in multiple conformations when combined with DS, which was also thermodynamically favorable. Although the BXBB sequence in DBPA remained highly dynamic in DBPB, it did not contribute much to the binding due to the exposure of the C-terminus and the position of the linker in DBPB.

Hyaluronic Acid

Hyaluronic acid has a different synthesis site (plasma membrane) and a different synthesis form (non-glycoprotein) compared to other GAGs. HA will not undergo further modification; thus, the interaction between it and the protein seems to be structurally specific. The hydrogen bonds and intramolecular hydrogen bonds with water molecules gave it a complex β-sheet structure (). In the double helix structure of HA, every two monosccharide flip 180°. HA, as a structural scaffold, widely exists in the epithelial tissue, connective tissue and nerve tissue of vertebrates and regulates the physical and chemical processes of tissue hydration and penetration. The interaction between HA and HA-binding protein (hyaluroadhesin) mediates various physiological activities, such as cell signal transduction, wound repair, tissue regeneration, leukocyte rolling adhesion and inflammation (). Most HA-binding proteins belong to the link protein superfamily. Some other proteins (such as receptor for hyaluronan-mediated motility, RHAMM) and peptides (thymosin α1, Tα1) bound to HA are independent of the link module ().

The 14 human link proteins can be divided into three categories (A, B, C) according to their structural composition (). TSG-6 was the most typical type A Link protein, and its HA-binding domain (HABD) was the only Link module (Figure 5; ). The link module was composed of 100 amino acids and structured by two β-sheets and two α-helices, which were stabilized by two extremely conserved disulfide bonds. The two β-sheets were composed of four and two β-strands. Type B Link protein used CD44 as a template. It extended the β-sheet at the C- and N-termini on the basis of type A (adding four β strands), and the HABD of type B was redefined (). The type C link protein was composed of two links in series, both of which participate in binding with HA. This subcategory included aggrecan, versican and HAPLN1-4, but detailed research on its structure is lacking. The binding of HA and protein had very strict requirements on the tertiary structure of the protein. This was most obvious in the type C Link protein, which did not interact with GAGs other than HA. In one study, three link modules were connected in series, but the binding activity with HA was completely lost ().

FIGURE 5

Kahmann proposed that the binding of Link-TSG-6 and HA was concentrated in the β4/β5 loop. The association was accompanied by the rearrangement of C47 and C68 disulfide bonds (). In the previously proposed B(X)7B rule motif (R5EARSGKYK13), R5 and K13 had no obvious evidence of involvement in binding, but K11 was the main binding residue. In Blundell’s subsequent research, it was shown that the folding of the link module remains unchanged during the combination (). The largest structural change was found in β4/β5. K11 also changed its orientation and became more oriented. For Y59 and Y58, the benzene rings did not rotate due to ring stacking. Due to the derived polarity of the binding, the two ends of the binding were located at K11 and R81. Higman proposed that in the free state, the β4/β5 loop of TSG-6 was highly dynamic. In this state, there was a conformation that exposes aromatic residues and captured HA by stacking interactions and then rearranged structural elements, such as the β4/β5 loop (). There were two structural elements that were obviously solidified, one of which was G10 located at the corner of α1/β1, and the other was K54 of β3/β4. K54 was far from the HA-binding site but played an important role in the binding of heparin to TSG-6. Its solidification explained the problem that HA and heparin could not bind to TSG-6 at the same time, although they have different binding sites.

In the 2014 study, HA and hybrid HA of different lengths were used to study the interaction with Link-TSG-6 (). Although the heptasaccharide with the reducing end of GlcA (HA7AA) had a complete binding structure, the entropy was unfavorable. Therefore, the octasaccharide with the reducing end of GlcNAc (HA8AN) was defined as the minimum unit required for binding. HSQC data clearly showed that HA8NA and HA7AA had two binding modes, with the reducing end GlcA bound to K63/H45 as the dominant one. The affinity of HA8NA was twice that of HA8AN, while the affinity of the two heptasaccharides had no such difference. The reason for the difference in specific affinity is unknown. In the binding model of HA8AN and TSG-6, H45 and K63 appear to be new binding residues. They bound to the reducing terminal disaccharide of the octasaccharide to make the binding tighter. The binding of HA and Link-TSG-6 was mainly through ionic interactions, ring-stacking interactions, hydrogen bonding, van der Waals forces and hydrophobic repulsion. Since the binding occurred on two interfaces, this imposed an inevitable requirement for the distortion of the two glycosidic bonds between the fifth and seventh residues. For heptasaccharides, the significant reduction in the affinity of hexasaccharides might be due to the lack of multiple groups of binding, resulting in instability of the distortion of glycosidic bonds. The CS part of hybrid HA will also be distorted during binding, but due to the lack of structural elements and the lack of hydrogen bonds during binding, the affinity was far lower than that of HA. However, due to the existence of binding, this provided a certain explanation for the chondroprotective function of TSG-6. CS, Heparin and HA had different binding modes with TSG-6, giving TSG-6 complex biological functions.

The HABD in CD44 was mainly located in the link module, C-terminal extension and α1-helix. Two N-linked glycosylation sites (N25 and N100) were also located in the HABD (). Teriete pointed out that octasaccharide might be the smallest unit that satisfies all binding requirements (). All binding sites were located on the same plane, but due to the scattered distribution, there might be two incompatible binding modes. One used N100/N101 to R150/R154, similar to the combination of TSG-6 and HA. The other used K38/R162 as the terminal binding, and the binding was farther away from the charged area. The data showed that the binding is accompanied by a structural rearrangement. Takeda proposed that the parallel sheets of β8 and β0 involved rearrangement, which might be related to the special structure of β8 (). More thorough structural changes were located at the C-terminal extensions of α3 and β9, and their structure changed from a regular to a randomized structure after the combination. This result was in conflict with crystal studies, which showed that binding did not involve changes in C-terminal extension (). But unlike other studies, the protein used by Banerji is of mouse origin. And in the model established in this study, the complex is in two conformational equilibrium (type A and B, Figure 6). The difference between the two conformations is the orientation of R45 (human CD44 R41). Ogino also proposed that CD44 was in the balance of two conformations in the unbound or bound state (). In the unbound state, it had a regular structure and low HA affinity, which was conducive to cell rolling. In the combined state, it was mainly a random structure with high HA affinity, which was conducive to cell adhesion. The balance of these two states was conducive to the physiological activity of CD44-mediated cell rolling.

FIGURE 6

In terms of RHAMM, two amino acid clusters were mainly involved in binding with HA: the first was the proposed BX7B structure (K531-K541), and the second was K553-K562 (). Studies have shown that the second binding site plays a major role in binding. Studies on Tα1 indicated that the binding is mainly related to its terminal L16KEKK20 (). The combination of HA and these two substances occurred mainly through electrostatic forces, which was different from the role of HA with TSG-6 and CD44. The combination of HA and CD44 was mainly through hydrogen bonding and van der Waals forces, while the combination with TSG-6 was mainly through electrostatic forces and aromatic accumulation.

Kertan Sulfate

Kertan sulfate is the only GAG without any acidic uronic acid residue, and its interaction with proteins mainly depends on structural characteristics and sulfation modification. KS is mainly distributed in the cornea and cartilage tissue and is divided into three categories (I-III) according to the distribution and connection with glycoproteins (). KS plays an important role in brain development, neurodevelopment and regeneration, implantation and fertilization and maintains the balance of tissue hydration properties (; ; ). KS has many protein partners, including tyrosine protein kinases, inflammatory cytokines, growth factors, chemokines, cytoskeletal cells, and lectins. Only a few studies of the interaction between KS and protein have been investigated using NMR ().

Galectin 3 (Gal-3) seems to be one of KS’s most tacit partners, and its distribution is extremely close to that of KS. The interaction between full-length Gal-3 and KS has been studied using HSQC; the disturbance was found to be in the β1, β3, β4, β5, β6, and β10 strands, and the β10 strand was the most important strand. The binding domain can be on the S- and F-faces in Gal-3. When the N-terminal tail of Gal-3 was truncated, KS interaction on the S-face became more obvious. The presence of other negatively charged regions did not affect the binding between KS and the Gal-3 S-face according to MD data. In the binding state, the conformations of the F-face and the N-terminal tail were changed. The binding was mainly concentrated on the left side of the S-face, which facilitated its combination with other proteins or heteropolymerization with other galectins. However, the pulse field gradient NMR data showed that KS did not induce oligomerization of Gal-3. Desulfated KS had far less affinity than KS, and the chemical shift disturbances on the F-face and N-terminal tail were greatly reduced.

Conclusion

Glycosaminoglycans, as common glycoproteins in biological systems, are involved in many physiological and pathological processes. The study of their structure and interaction with proteins has received extensive attention, but the study of molecular perspectives is only the tip of the iceberg. This not only is due to the delay of carbohydrate research but is also related to the limitations of technology. The information produced by NMR is incomparable to all other technologies. For example, it can provide information about the binding affinity constant, on/off chemical exchange rate, binding site and atomic information, but high-precision research is more demanding for technology. In particular, regarding the special existence of GAG, its highly complex structure not only endows it with rich biological functions but also brings incomparable difficulties for research. The study of the interactions between GAG and proteins using NMR is based on complete structural characterizations of GAG and/or proteins, which face huge obstacles. Biosynthesis carriers of GAG are difficult to find, while chemical and enzymatic syntheses are limited to a few scientists. This in turn makes it difficult to obtain isotope-labeled GAG. Because the binding of GAG and protein has obvious multibinding characteristics, it will cause oligomerization and even precipitation. The application of NMR technology is mainly limited by several factors, including the length of the oligosaccharides, the molecular weight of the proteins, and the concentration range and stability of the complex. However, with the renewal and iteration of technology, the rise of high magnetic flux nuclear magnetic spectrometry and enzymatic chemical synthesis has injected a steady stream of vitality into interaction research. The study of the interaction between GAG and proteins is helpful for understanding various physiological and pathological mechanisms and has a huge impetus for drug development.

Statements

Author contributions

CB and LJ participated in preparation, creation, initial draft writing and review of this article. Both authors contributed to the article and approved the submitted version.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Summary

Keywords

glycosaminoglycans, proteins, interaction, NMR, conformation

Citation

Bu C and Jin L (2021) NMR Characterization of the Interactions Between Glycosaminoglycans and Proteins. Front. Mol. Biosci. 8:646808. doi: 10.3389/fmolb.2021.646808

Received

28 December 2020

Accepted

24 February 2021

Published

16 March 2021

Volume

8 - 2021

Edited by

Fuming Zhang, Rensselaer Polytechnic Institute, United States

Reviewed by

Xiaojun Sun, University of Jinan, China; Yasuteru Shigeta, University of Tsukuba, Japan

Updates

Copyright

*Correspondence: Lan Jin,

This article was submitted to Molecular Recognition, a section of the journal Frontiers in Molecular Biosciences

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics