Solution NMR conformation of glycosaminoglycans

Structural Insights into the Cofactor Role of Heparin/Heparan Sulfate in Binding between the SARS-CoV-2 Spike Protein and Host Angiotensin-Converting Enzyme II

The viral entry process of the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) requires heparin and heparan sulfates from the cell surface, functioning as a cofactor for human angiotensin-converting enzyme 2 (ACE2) for recognizing the receptor-binding domain (RBD) of the spike (S) protein on the surface of the virion. In the present study, the binding poses of an oligosaccharide with four repeating units of GlcNS6S-IdoA2S (octa) predicted by Vina-Carb in the RBD binding site were employed in molecular dynamics (MD) simulations to provide atomic details for studying the cofactor mechanism. The molecular model in the MD simulations reproduced the length- and sequence-dependent behavior observed from the microarray experiments and revealed an important planar U-turn shape for HP/HS binding to RBD. The model for octa with this shape in the ACE2-RBD complex enhanced the interactions in the binding interface. The comparisons with the ACE2-RBD complex suggested that the presence of octa in the RBD binding site blocked the movements in a loop region at the distal end of the RBD binding interface and promoted the contacts of this loop region with the ACE2 N-terminus helix. This study shed light on the atomic and dynamic details for HP/HS interacting with RBD and provided insights into their cofactor role in the ACE2-RBD interactions.

Investigation of the structure of regulatory proteins interacting with glycosaminoglycans by combining NMR spectroscopy and molecular modeling - the beginning of a wonderful friendship

The interaction of regulatory proteins with extracellular matrix or cell surface-anchored glycosaminoglycans (GAGs) plays important roles in molecular recognition, wound healing, growth, inflammation and many other processes. In spite of their high biological relevance, protein-GAG complexes are significantly underrepresented in structural databases because standard tools for structure determination experience difficulties in studying these complexes. Co-crystallization with subsequent X-ray analysis is hampered by the high flexibility of GAGs. NMR spectroscopy experiences difficulties related to the periodic nature of the GAGs and the sparse proton network between protein and GAG with distances that typically exceed the detection limit of nuclear Overhauser enhancement spectroscopy. In contrast, computer modeling tools have advanced over the last years delivering specific protein-GAG docking approaches successfully complemented with molecular dynamics (MD)-based analysis. Especially the combination of NMR spectroscopy in solution providing sparse structural constraints with molecular docking and MD simulations represents a useful synergy of forces to describe the structure of protein-GAG complexes. Here we review recent methodological progress in this field and bring up examples where the combination of new NMR methods along with cutting-edge modeling has yielded detailed structural information on complexes of highly relevant cytokines with GAGs.

Unveiling the Molecular Structure of Antimalarial Drugs Chloroquine and Hydroxychloroquine in Solution through Analysis of 1H NMR Chemical Shifts

Chloroquine (CQ) and hydroxychloroquine (HCQ) have been standard antimalarial drugs since the early 1950s, and very recently, the possibility of their use for the treatment of COVID-19 patients has been considered. To understand the drug mode of action at the submicroscopic level (atoms and molecules), molecular modeling studies with the aid of computational chemistry methods have been of great help. A fundamental step in such theoretical investigations is the knowledge of the predominant drug molecular structure in solution, which is the real environment for the interaction with biological targets. Our strategy to access this valuable information is to perform density functional theory (DFT) calculations of ¹H NMR chemical shifts for several plausible molecular conformers and then find the best match with experimental NMR profile in solution (since it is extremely sensitive to conformational changes). Through this procedure, after optimizing 30 trial distinct molecular structures (ωB97x-D/6-31G(d,p)-PCM level of calculation), which may be considered representative conformations, we concluded that the global minimum (named M24), stabilized by an intramolecular N-H hydrogen bond, is not likely to be observed in water, chloroform, and dimethyl sulfoxide (DMSO) solution. Among fully optimized conformations (named M1 to M30, and MD1 and MD2), we found M12 (having no intramolecular H-bond) as the most probable structure of CQ and HCQ in water solution, which is a good approximate starting geometry in drug-receptor interaction simulations. On the other hand, the preferred CQ and HCQ structure in chloroform (and CQ in DMSO-d₆) solution was assigned as M8, showing the solvent effects on conformational preferences. We believe that the analysis of ¹H NMR data in solution can establish the connection between the macro level (experimental) and the sub-micro level (theoretical), which is not so apparent to us and appears to be more appropriate than the thermodynamic stability criterion in conformational analysis studies.

NMR Characterization of the Interactions Between Glycosaminoglycans and Proteins

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.

A review of NMR analysis in polysaccharide structure and conformation: Progress, challenge and perspective

Nuclear magnetic resonance (NMR) has been widely used as an analytical chemistry technique to investigate the molecular structure and conformation of polysaccharides. Combined with 1D spectra, chemical shifts and coupling constants in both homo- and heteronuclear 2D NMR spectra are able to infer the linkage and sequence of sugar residues. Besides, NMR has also been applied in conformation, quantitative analysis, cell wall in situ, degradation, polysaccharide mixture interaction analysis, as well as carbohydrates impurities profiling. This review summarizes the principle and development of NMR in polysaccharides analysis, and provides NMR spectra data collections of some common polysaccharides. It will help to promote the application of NMR in complex polysaccharides of biochemical interest, and provide valuable information on commercial polysaccharide products.

Exploiting the Amazing Diversity of Natural Source-Derived Polysaccharides: Modern Procedures of Isolation, Engineering, and Optimization of Antiviral Activities

Naturally occurring polysaccharide sulfates are highly diverse, owning variations in the backbone structure, linkage pattern and stereochemistry, branching diversity, sulfate content and positions of sulfate group(s). These structural characteristics bring about diverse sulfated polymers with dissimilar negative charge densities and structure-activity relationships. Herein, we start with a short discussion of techniques needed for extraction, purification, chemical sulfation, and structural characterization of polysaccharides. Processes of isolation and sulfation of plant-derived polysaccharides are challenging and usually involve two steps. In this context, we describe an integrated extraction-sulfation procedure that produces polysaccharide sulfates from natural products in one step, thereby generating additional pharmacological activities. Finally, we provide examples of the spectrum of natural source-derived polysaccharides possessing specific features of bioactivity, in particular focusing on current aspects of antiviral drug development and drug-target interaction. Thus, the review presents a detailed view on chemically engineered polysaccharides, especially sulfated derivatives, and underlines their promising biomedical perspectives.

The conformation of the idopyranose ring revisited: How subtle O-substituent induced changes can be deduced from vicinal 1H-NMR coupling constants

The idopyranose ring plays a pivotal role in the conformational, dynamical, and intermolecular binding aspects of glycosaminoglycans like heparin and dermatan sulfate and it was early on assigned a role in the Sugar Code governing biological recognition processes. There is consensus that next to the two canonical ¹C₄ and ⁴C₁ chair conformations, the conformational space accessible to the idopyranose ring entails a ²S_O skew-boat conformation, but the equilibrium between these three ring puckers has evaded satisfactory quantification. In this study a meta-analysis of X-ray solid-state data and vicinal NMR coupling constants is presented, based on the Truncated Fourier Puckering (TFP) formalism and the generalized Karplus (CAGPLUS) equation. This approach yields a model-free, granular and consistent reckoning of 159 idopyranose solution puckering equilibria studied by NMR and allows us to reproduce the involved 636 NMR vicinal couplings with an overall residual RMS(J_obs-J_calc) of 0.184 Hz. Our analyses show that for all ring systems examined, the idopyranosyl chair conformations take up the same ring pucker irrespective of the ring substituent pattern or a vast variety in experimental conditions. Instead, it is the (skew-)boat conformation that adapts to the substitution pattern of the idopyranose ring or a specific sulfation pattern of neighboring saccharides. All idopyranose rings are involved in conformational equilibria that subsume the aforementioned conformers which turn out to differ only a few kJ/mole in conformational energy. Thus, the plasticity and flexibility of idopyranose remains intact under practically all circumstances and, as the glycosidic linkages in heparin are considered to be relatively stiff, the iduronic moiety functions as the linchpin of heparin flexibility thereby being rather a "space(r)" than a "letter" in the alleged Sugar Code alphabet.

Characterization and Quantification of Highly Sulfated Glycosaminoglycan Isomers by Gated-Trapped Ion Mobility Spectrometry Negative Electron Transfer Dissociation MS/MS

Glycosaminoglycans (GAGs) play vital roles in many biological processes and are naturally present as complex mixtures of polysaccharides with tremendous structural heterogeneity, including many structural isomers. Mass spectrometric analysis of GAG isomers, in particular highly sulfated heparin (Hep) and heparan sulfate (HS), is challenging because of their structural similarity and facile sulfo losses during analysis. Herein, we show that highly sulfated Hep/HS isomers may be resolved by gated-trapped ion mobility spectrometry (gated-TIMS) with negligible sulfo losses. Subsequent negative electron transfer dissociation (NETD) tandem mass spectrometry (MS/MS) analysis of TIMS-separated Hep/HS isomers generated extensive glycosidic and cross-ring fragments for confident isomer differentiation and structure elucidation. The high mobility resolution and preservation of labile sulfo modifications afforded by gated-TIMS MS analysis also allowed relative quantification of highly sulfated heparin isomers. These results show that the gated-TIMS-NETD MS/MS approach is useful for both qualitative and quantitative analysis of highly sulfated Hep/HS compounds in a manner not possible with other techniques.

Effects of sulfation and the environment on the structure of chondroitin sulfate studied via Raman optical activity

Raman optical activity reflects differences in the secondary structure of chondroitin caused by its sulfation.

Solution Conformation of Heparin Tetrasaccharide. DFT Analysis of Structure and Spin⁻Spin Coupling Constants

Density functional theory (DFT) has provided detailed information on the molecular structure and spin⁻spin coupling constants of heparin tetrasaccharide (GlcNS,6S-IdoA2S-GlcNS,6S-IdoA2S-OMe) representing the predominant heparin repeating-sequence. The fully optimised molecular structures of two tetrasaccharide conformations (differing from each other in the conformational form of the sulphated iduronic acid residue⁻one ¹C₄ and the other ²S₀) were obtained using the B3LYP/6-311+G(d,p) level of theory and applying explicit water molecules to simulate the presence of a solvent. The theoretical data provided insight into variations of the bond lengths, bond angles and torsion angles, formations of intra- and intermolecular hydrogen bonds and ionic interactions. Optimised molecular structures indicated the formation of a complex hydrogen bond network, including interresidue and intraresidue bonds. The ionic interactions strongly influence the first hydration shell and, together with hydrogen bonds, play an important role in shaping the 3D tetrasaccharide structure. DFT-derived indirect three⁻bond proton⁻proton coupling constants (³J_H-C-C-H) showed that the best agreement with experiment was obtained with a weighted average of 67:33 (¹C₄:²S₀) of the IdoA2S forms. Detailed analysis of Fermi-contact contributions to ³J_H-C-C-H showed that important contributions arise from the oxygen lone pairs of neighbouring oxygen atoms. The analysis also showed that the magnitude of diamagnetic spin⁻orbit contributions are sufficiently large to determine the magnitude of some proton⁻proton coupling constants. The data highlight the need to use appropriate quantum-chemical calculations for a detailed understanding of the solution properties of heparin oligosaccharides.

Differential calixarene receptors create patterns that discriminate glycosaminoglycans

A well-designed fluorescence displacement sensing array based on calixarene receptors realizes the discrimination of glycosaminoglycans.

Computational Chemistry Tools in Glycobiology: Modelling of Carbohydrate–Protein Interactions

Molecular modelling provides a major impact in the field of glycosciences, helping in the characterisation of the molecular basis of the recognition between lectins from pathogens and human glycoconjugates, and in the design of glycocompounds with anti-infectious properties. The conformational properties of oligosaccharides are complex, and therefore, the simulation of these properties is a challenging task. Indeed, the development of suitable force fields is required for the proper simulation of important problems in glycobiology, such as the interatomic interactions responsible for oligosaccharide and glycoprotein dynamics, including O-linkages in oligo- and polysaccharides, and N- and O-linkages in glycoproteins. The computational description of representative examples is discussed, herein, related to biologically active oligosaccharides and their interaction with lectins and other proteins, and the new routes open for the design of glycocompounds with promising biological activities.

Solution Structure of Heparin Pentasaccharide: NMR and DFT Analysis

High-resolution NMR and density functional theory (DFT) calculations have been applied to analysis of heparin pentasaccharide 3D structure in aqueous solution. The fully optimized molecular geometry of two pentasaccharide conformations (differing from each other in the form, one (1)C4 and the other (2)S0, of the sulfated iduronic acid residue) were obtained using the B3LYP/6-311+G(d,p) level of theory in the presence of solvent, the latter included as explicit water molecules. The presented approach enabled insight into variations of the bond lengths, bond angles, and torsion angles, formations of intra- and intermolecular hydrogen bonds, and ionic interactions in the two pentasaccharide conformations. A rather complex hydrogen bond network is formed, including inter-residue and intraresidue bonds between the NH group in the GlcN,3,6S with oxygens linked to C-2 at the IdoA2S residue and the glycosidic O-1 and the neighboring OSO3(-) group linked to C-3 in the same residue. On the other hand, because the first hydration shell is strongly influenced by strong ion-ion and ion-dipole interactions between sodium ions, sulfates, carboxylates, and -OH groups, ionic interactions play an important role in the stabilization of the 3D structure. The DFT-computed three-bond proton-proton coupling constants also showed that best agreement with experiment was obtained with a weighted average of 15:85 ((1)C4/(2)S0) of the sulfated iduronic acid forms indicating that the ratio is even more shifted toward the (2)S0 form than previously supposed. The DFT-computed pentasaccharide conformation differs from the previously published data, with the main changes at the glycosidic linkages, namely, the ψ1 torsion angles and the ϕ3 angle. The comparison of the glycosidic linkage torsion angle values in solution with the antithrombin-pentasaccharide complex also indicates that the pentasaccharide conformation changes upon binding to antithrombin III. The data supports the assumption that the protein selects the more populated (2)S0 conformer of heparin pentasaccharide and, consequently, the binding process of heparin pentasaccharide with antithrombin III is energetically more favorable than formerly expected.

Sulfation and cation effects on the conformational properties of the glycan backbone of chondroitin sulfate disaccharides

Chondroitin sulfate (CS) is one of several glycosaminoglycans that are major components of proteoglycans. A linear polymer consisting of repeats of the disaccharide -4GlcAβ1-3GalNAcβ1-, CS undergoes differential sulfation resulting in five unique sulfation patterns. Because of the dimer repeat, the CS glycosidic "backbone" has two distinct sets of conformational degrees of freedom defined by pairs of dihedral angles: (ϕ1, ψ1) about the β1-3 glycosidic linkage and (ϕ2, ψ2) about the β1-4 glycosidic linkage. Differential sulfation and the possibility of cation binding, combined with the conformational flexibility and biological diversity of CS, complicate experimental efforts to understand CS three-dimensional structures at atomic resolution. Therefore, all-atom explicit-solvent molecular dynamics simulations with Adaptive Biasing Force sampling of the CS backbone were applied to obtain high-resolution, high-precision free energies of CS disaccharides as a function of all possible backbone geometries. All 10 disaccharides (β1-3 vs β1-4 linkage × five different sulfation patterns) were studied; additionally, ion effects were investigated by considering each disaccharide in the presence of either neutralizing sodium or calcium cations. GlcAβ1-3GalNAc disaccharides have a single, broad, thermodynamically important free-energy minimum, whereas GalNAcβ1-4GlcA disaccharides have two such minima. Calcium cations but not sodium cations bind to the disaccharides, and binding is primarily to the GlcA -COO(-) moiety as opposed to sulfate groups. This binding alters the glycan backbone thermodynamics in instances where a calcium cation bound to -COO(-) can act to bridge and stabilize an interaction with an adjacent sulfate group, whereas, in the absence of this cation, the proximity of a sulfate group to -COO(-) results in two like charges being both desolvated and placed adjacent to each other and is found to be destabilizing. In addition to providing information on sulfation and cation effects, the present results can be applied to building models of CS polymers and as a point of comparison in studies of CS polymer backbone dynamics and thermodynamics.

Impact of sulfation pattern on the conformation and dynamics of sulfated fucan oligosaccharides as revealed by NMR and MD

Sulfated fucans from sea urchin egg jelly express well-defined chemical structures that vary with species. This species specificity regulates the sperm acrosome reaction, a critical step to assure intra-specific fertilization. In addition, these polysaccharides are involved in other biological activities such as anticoagulation. Although sulfation patterns are relevant to the levels of response in both activities, conformation and dynamics of these glycans are also contributing factors. However, data about these features of sulfated fucans are very rare. To address this, we have employed nuclear magnetic resonance experiments combined with molecular dynamics on structurally defined oligosaccharides derived from two sulfated fucans. The results have indicated that the oligosaccharides are flexible in solution. Ring conformation of their composing units displays just the (1)C4 chair configuration. In a particular octasaccharide, composed of two tetrasaccharide sequences, inter-residual hydrogen bonds play a role to decrease dynamics in these repeating units. Conversely, the linking disaccharide [-3)-α-L-Fucp-2(OSO3(-))-(1-3)-α-L-Fucp-4(OCO3(-))-(1-] located right between the two tetrasaccharide units has amplified motions suggested to be promoted by electrostatic repulsion of sulfates on opposite sides of the central glycosidic bond. This conjunction of information about conformation and dynamics of sulfated fucan oligosaccharides provides new insights to explain how these glycans behave free in solution and influenced by sulfation patterns. It may also serve for future studies concerning structure-function relationship of sulfated fucans, especially those involving sea urchin fertilization and anticoagulation.

Keratan sulfate: an up-to-date review

Keratan sulfate (KS) is a glycosaminoglycan (GAG) type consisted of a sulfated poly-N-acetyl lactosamine chain. Besides acting as a constitutive molecule of the extracellular matrices, this GAG also plays a role as a hydrating and signaling agent in cornea and cartilage tissues. Inasmuch, KS is widely explored in the pharmaceutical industry. This review will cover the major achievements described in the literature of 2010-2014 concerning this GAG. Discussion about KS' roles in physiopathological conditions, as target or therapeutic molecule in diseases, methods of analysis and detection as well as KS-related enzymes, metabolism and developmental biology is properly provided.

NMR-based dynamics of free glycosaminoglycans in solution

Dynamical behaviors of glycosaminoglycans, as here illustrated with a hyaluronan oligosaccharide, are key regulators of biological functions.