Direct Evidence for Hydrogen Bonding in Glycans: A Combined NMR and Molecular Dynamics Study

Measuring Hydroxyl Exchange Rates in Glycans Using a Synergistic Combination of Saturation Transfer and Relaxation Dispersion NMR

Molecular recognition events mediated by glycans play pivotal roles in controlling the fate of diverse biological processes such as cellular communication and the immune response. The affinity of glycans for their target receptors is governed primarily by the hydrogen bonds formed by hydroxyl groups decorating the glycan surface. Hydroxyl exchange rate constants are therefore vital parameters that report on glycan structure and dynamics. Here we present a strategy for characterizing hydroxyl hydrogen/deuterium (H/D) exchange in glycans that employs a synergistic combination of 13C chemical exchange saturation transfer (CEST) and Carr-Purcell-Meiboom-Gill relaxation dispersion (CPMG) NMR methods. We show that the combination of CEST and CPMG experiments facilitates the sensitive detection of the small (∼0.1 ppm) two-bond deuterium isotope shift on a 13C nucleus when the attached hydroxyl group fluctuates between protonated and deuterated states. This shift is leveraged for measuring site-specific kinetic H/D exchange rate constants as well as thermodynamic free energies of isotope fractionation. The CEST and CPMG modules are integrated with a selective J-cross-polarization scheme that provides the flexibility for rapid characterization of H/D exchange at a specific hydroxyl site. Moreover, our approach enables the precise isothermal measurement of hydroxyl exchange rate constants without the need for cumbersome isotope labeling. The H/D exchange rate constants of three different glycans assessed using this method highlight its potential for detecting transient intra- and intermolecular hydrogen bonds. In addition, the trends in H/D exchange rate constants establish site-specific steric accessibility as a key determinant of solvent exchange dynamics in glycans.

Primary Structure of Glycans by NMR Spectroscopy

Glycans, carbohydrate molecules in the realm of biology, are present as biomedically important glycoconjugates and a characteristic aspect is that their structures in many instances are branched. In determining the primary structure of a glycan, the sugar components including the absolute configuration and ring form, anomeric configuration, linkage(s), sequence, and substituents should be elucidated. Solution state NMR spectroscopy offers a unique opportunity to resolve all these aspects at atomic resolution. During the last two decades, advancement of both NMR experiments and spectrometer hardware have made it possible to unravel carbohydrate structure more efficiently. These developments applicable to glycans include, inter alia, NMR experiments that reduce spectral overlap, use selective excitations, record tilted projections of multidimensional spectra, acquire spectra by multiple receivers, utilize polarization by fast-pulsing techniques, concatenate pulse-sequence modules to acquire several spectra in a single measurement, acquire pure shift correlated spectra devoid of scalar couplings, employ stable isotope labeling to efficiently obtain homo- and/or heteronuclear correlations, as well as those that rely on dipolar cross-correlated interactions for sequential information. Refined computer programs for NMR spin simulation and chemical shift prediction aid the structural elucidation of glycans, which are notorious for their limited spectral dispersion. Hardware developments include cryogenically cold probes and dynamic nuclear polarization techniques, both resulting in enhanced sensitivity as well as ultrahigh field NMR spectrometers with a ¹H NMR resonance frequency higher than 1 GHz, thus improving resolution of resonances. Taken together, the developments have made and will in the future make it possible to elucidate carbohydrate structure in great detail, thereby forming the basis for understanding of how glycans interact with other molecules.

Epithelial Na+ channel and the glycocalyx: a sweet and salty relationship for arterial shear stress sensing

PURPOSE OF REVIEW

The ability of endothelial cells to sense mechanical force, and shear stress in particular, is crucial for normal vascular function. This relies on an intact endothelial glycocalyx that facilitates the production of nitric oxide (NO). An emerging arterial shear stress sensor is the epithelial Na+ channel (ENaC). This review highlights existing and new evidence for the interdependent activity of the glycocalyx and ENaC and its implications for vascular function.

RECENT FINDINGS

New evidence suggests that the glycocalyx and ENaC are physically connected and that this is important for shear stress sensing. The connection relies on N-glycans attached to glycosylated asparagines of α-ENaC. Removal of specific N-glycans reduced ENaC's shear stress response. Similar effects were observed following degradation of the glycocalyx. Endothelial specific viral transduction of α-ENaC increased blood pressure (∼40 mmHg). This increase was attenuated in animals transduced with an α-ENaC version lacking N-glycans.

SUMMARY

These observations indicate that ENaC is connected to the glycocalyx and their activity is interdependent to facilitate arterial shear stress sensation. Future research focusing on how N-glycans mediate this interaction can provide new insights for the understanding of vascular function in health and disease.

Methods for Measuring Exchangeable Protons in Glycosaminoglycans

Recent NMR studies of the exchangeable protons of GAGs in aqueous solution, including those of the amide, sulfamate, and hydroxyl moieties, have demonstrated potential for the detection of intramolecular hydrogen bonds providing insights into secondary structure preferences. GAG amide protons are observable by NMR over wide pH and temperature ranges; however, specific solution conditions are required to reduce the exchange rate of the sulfamate and hydroxyl protons and allow their detection by NMR. Building on the vast body of knowledge on detection of hydrogen bonds in peptides and proteins, a variety of methods can be used to identify hydrogen bonds in GAGs including temperature coefficient measurements, evaluation of chemical shift differences between oligo- and monosaccharides, and relative exchange rates measured through line shape analysis and EXSY spectra. Emerging strategies to allow direct detection of hydrogen bonds through heteronuclear couplings offer promise for the future. Molecular dynamic simulations are important in this effort both to predict and confirm hydrogen bond donors and acceptors.

The Incorporation of Labile Protons into Multidimensional NMR Analyses: Glycan Structures Revisited

Glycan structures are often stabilized by a repertoire of hydrogen-bonded donor/acceptor groups, revealing longer-lived structures that could represent biologically relevant conformations. NMR provides unique data on these hydrogen-bonded networks from multidimensional experiments detecting cross-peaks resulting from through-bond (TOCSY) or through-space (NOESY) interactions. However, fast OH/H₂O exchange, and the spectral proximity among these NMR resonances, hamper the use of glycans’ labile protons in such analyses; consequently, studies are often restricted to aprotic solvents or supercooled aqueous solutions. These nonphysiological conditions may lead to unrepresentative structures or to probing a small subset of accessible conformations that may miss “active” glycan conformations. Looped, projected spectroscopy (L-PROSY) has been recently shown to substantially enhance protein NOESY and TOCSY cross-peaks, for ¹Hs that undergo fast exchange with water. This study shows that even larger enhancements can be obtained for rapidly exchanging OHs in saccharides, leading to the retrieval of previously undetectable 2D TOCSY/NOESY cross-peaks with nonlabile protons. After demonstrating ≥300% signal enhancements on model monosaccharides, these experiments were applied at 1 GHz to elucidate the structural network adopted by a sialic acid homotetramer, used as a model for α,2–8 linked polysaccharides. High-field L-PROSY NMR enabled these studies at higher temperatures and provided insight previously unavailable from lower-field NMR investigations on supercooled samples, involving mostly nonlabile nuclei. Using L-PROSY’s NOEs and other restraints, a revised structural model for the homotetramer was obtained combining rigid motifs and flexible segments, that is well represented by conformations derived from 40 μs molecular dynamics simulations.

NMR analysis and molecular dynamics conformation of α-1,6-linear and α-1,3-branched isomaltose oligomers as mimetics of α-1,6-linked dextran

The conformational preferences of several α-1,6-linear and α-1,3-branched isomalto-oligosaccharides were investigated by NMR and MD-simulations. Right-handed helical structure contributed to the solution geometry in isomaltotriose and isomaltotetraose with one nearly complete helix turn and stabilizing intramolecular hydrogen bonds in the latter by MD-simulation. Decreased helix contribution was observed in α-1,3-glucopyranosyl- and α-1,3-isomaltosyl-branched saccharide chains. Especially the latter modification was predicted to cause a more compact structure consistent with literature rheology measurements as well as with published dextranase-resistant α-1,3-branched oligosaccharides. The findings presented here are significant because they shed further light on the conformational preference of isomalto-oligosaccharides and provide possible help for the design of dextran-based drug delivery systems or for the targeted degradation of capsular polysaccharides by dextranases in multi-drug resistant bacteria.

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.

The use of MM/QM calculations of 13C chemical shifts in the conformational analysis of some monosaccharides and sucrose

The ¹³C NMR chemical shifts of some common pyranose monosaccharides in D₂O solution were predicted using a combined molecular mechanics (Pcmod 9.1/MMFF94) and ab initio (GIAO (B3LYP/DFT, 6-31G(d))) model.

A combined NMR, MD and DFT conformational analysis of 9-O-acetyl sialic acid-containing GM3 ganglioside glycan and its 9-N-acetyl mimic

O-Acetylation of carbohydrates such as sialic acids is common in nature, but its role is not clearly understood due to the lability of O-acetyl groups. We demonstrated previously that 9-acetamido-9-deoxy-N-acetylneuraminic acid (Neu5Ac9NAc) is a chemically and biologically stable mimic of the 9-O-acetyl-N-acetylneuraminic acid (Neu5,9Ac2) of the corresponding sialoglycans. Here, a systematic nuclear magnetic resonance (NMR) spectroscopic and molecular dynamics (MD) simulation study was undertaken for Neu5,9Ac2-containing GM3 ganglioside glycan (GM3-glycan) and its Neu5Ac9NAc analog. GM3-glycan with Neu5Ac as the non-O-acetyl form of Neu5,9Ac2 was used as a control. Complete 1H and 13C NMR chemical shift assignments, three-bond 1H-13C trans-glycosidic coupling constants (3JCH), accurate 1H-1H coupling constants (3JHH), nuclear Overhauser effects and hydrogen bonding detection were carried out. Results show that structural modification (O- or N-acetylation) on the C-9 of Neu5Ac in GM3 glycan does not cause significant conformational changes on either its glycosidic dihedral angles or its secondary structure. All structural differences are confined to the Neu5Ac glycerol chain, and minor temperature-dependent changes are seen in the aglycone portion. We also used Density Functional Theory (DFT) quantum mechanical calculations to improve currently used 3JHH Karplus relations. Furthermore, OH chemical shifts were assigned at -10°C and no evidence of an intramolecular hydrogen bond was observed. The results provide additional evidence regarding structural similarities between sialosides containing 9-N-acetylated and 9-O-acetylated Neu5Ac and support the opportunity of using 9-N-acetylated Neu5Ac as a stable mimic to study the biochemical role of 9-O-acetylated Neu5Ac.

Conformational Studies of Oligosaccharides

The conformation of a molecule strongly affects its function, as demonstrated for peptides and nucleic acids. This correlation is much less established for carbohydrates, the most abundant organic materials in nature. Recent advances in synthetic and analytical techniques have enabled the study of carbohydrates at the molecular level. Recurrent structural features were identified as responsible for particular biological activities or material properties. In this Minireview, recent achievements in the structural characterization of carbohydrates, enabled by systematic studies of chemically defined oligosaccharides, are discussed. These findings can guide the development of more potent glycomimetics. Synthetic carbohydrate materials by design can be envisioned.

In Situ Differentiation of Multiplex Noncovalent Interactions Using SERS and Chemometrics

Probing changes of noncovalent interactions is crucial to study the binding efficiencies and strengths of (bio)molecular complexes. While surface-enhanced Raman scattering (SERS) offers unique molecular fingerprints to examine such interactions in situ, current platforms are only able to recognize hydrogen bonds because of their reliance on manual spectral identification. Here, we differentiate multiple intermolecular interactions between two interacting species by synergizing plasmonic liquid marble-based SERS platforms, chemometrics, and density functional theory. We demonstrate that characteristic 3-mercaptobenzoic acid (probe) Raman signals have distinct peak shifts upon hydrogen bonding and ionic interactions with tert-butylamine, a model interacting species. Notably, we further quantify the contributions from each noncovalent interaction coexisting in different proportions. As a proof-of-concept, we detect and categorize biologically important nucleotide bases based on molecule-specific interactions. This will potentially be useful to study how subtle changes in biomolecular interactions affect their structural and binding properties.

Conformational properties of l-fucose and the tetrasaccharide building block of the sulfated l-fucan from Lytechinus variegatus

Although the 3D structure of carbohydrates is known to contribute to their biological roles, conformational studies of sugars are challenging because their chains are flexible in solution and consequently the number of 3D structural restraints is limited. Here, we investigate the conformational properties of the tetrasaccharide building block of the Lytechinus variegatus sulfated fucan composed of the following structure [l-Fucp4(SO₃^-)-α(1-3)-l-Fucp2,4(SO₃^-)-α(1-3)-l-Fucp2(SO₃^-)-α(1-3)-l-Fucp2(SO₃^-)] and the composing monosaccharide unit Fucp, primarily by nuclear magnetic resonance (NMR) experiments performed at very low temperatures and using H₂O as the solvent for the sugars rather than using the conventional deuterium oxide. By slowing down the fast chemical exchange rates and forcing the protonation of labile sites, we increased the number of through-space ¹H-¹H distances that could be measured by NMR spectroscopy. Following this strategy, additional conformational details of the tetrasaccharide and l-Fucp in solution were obtained. Computational molecular dynamics was performed to complement and validate the NMR-based measurements. A model of the NMR-restrained 3D structure is offered for the tetrasaccharide.

Inter-residual Hydrogen Bonding in Carbohydrates Unraveled by NMR Spectroscopy and Molecular Dynamics Simulations

Carbohydrates, also known as glycans in biological systems, are omnipresent in nature where they as glycoconjugates occur as oligo- and polysaccharides linked to lipids and proteins. Their three-dimensional structure is defined by two or three torsion angles at each glycosidic linkage. In addition, transglycosidic hydrogen bonding between sugar residues may be important. Herein we investigate the presence of these inter-residue interactions by NMR spectroscopy in D₂ O/[D₆ ]DMSO (70:30) or D₂ O and by molecular dynamics (MD) simulations with explicit water as solvent for disaccharides with structural elements α-d-Manp-(1→2)-d-Manp, β-d-GlcpNAc-(1→2)-d-Manp, and α-d-Glcp-(1→4)-β-d-Glcp, all of which have been suggested to exhibit inter-residue hydrogen bonding. For the disaccharide β-d-GlcpNAc-(1→2)-β-d-Manp-OMe, the large extent of O5'⋅⋅⋅HO3 hydrogen bonding as seen from the MD simulation is implicitly supported by the ¹ H NMR chemical shift and ³ J_HO3,H3 value of the hydroxy proton. In the case of α-d-Glcp-(1→4)-β-d-Glcp-OMe, the existence of a transglycosidic hydrogen bond O2'⋅⋅⋅HO3 was proven by the presence of a cross-peak in ¹ H,¹³ C HSQC-TOCSY experiments as a result of direct TOCSY transfer between HO3 of the reducing end residue and H2' (detected at C2') of the terminal residue. The occurrence of inter-residue hydrogen bonding, albeit transient, is judged important for the stabilization of three-dimensional structures, which may be essential in maintaining a conformational state for carbohydrate-protein interactions of glycans to take place in biologically important environments.

Novel NMR Avenues to Explore the Conformation and Interactions of Glycans

This perspective article is focused on the presentation of the latest advances in NMR methods and applications that are behind the exciting achievements in the understanding of glycan receptors in molecular recognition events. Different NMR-based methodologies are discussed along with their applications to scrutinize the conformation and dynamics of glycans as well as their interactions with protein receptors.

Investigation of the Amide Proton Solvent Exchange Properties of Glycosaminoglycan Oligosaccharides

One-dimensional ¹H NMR experiments were conducted for aqueous solutions of glycosaminoglycan oligosaccharides to measure the amide proton temperature coefficients and activation energy barriers for solvent exchange and evaluate the effect of pH on the solvent exchange properties. A library of mono- and oligosaccharides was prepared by enzymatic depolymerization of amide-containing polysaccharides and by chemical modification of heparin and heparan sulfate saccharides including members that contain a 3- O-sulfated glucosamine residue. The systematic evaluation of this saccharide library facilitated assessment of the effects of structural characteristics, such as size, sulfation number and site, and glycosidic linkage, on amide proton solvent exchange rates. Charge repulsion by neighboring negatively charged sulfate and carboxylate groups was found to have a significant impact on the catalysis of amide proton solvent exchange by hydroxide. This observation leads to the conclusion that solvent exchange rates must be interpreted within the context of a given chemical environment. On their own, slow exchange rates do not conclusively establish the involvement of a labile proton in a hydrogen bond, and additional supporting experimental evidence such as reduced temperature coefficients is required.

Synthesis and O-Glycosidic Linkage Conformational Analysis of 13C-Labeled Oligosaccharide Fragments of an Antifreeze Glycolipid

NMR studies of two ¹³C-labeled disaccharides and a tetrasaccharide were undertaken that comprise the backbone of a novel thermal hysteresis glycolipid containing a linear glycan sequence of alternating [βXyl p-(1→4)-βMan p-(1→4)] _n dimers. Experimental trans-glycoside NMR J-couplings, parameterized equations obtained from density functional theory (DFT) calculations, and an in-house circular statistics package ( MA'AT) were used to derive conformational models of linkage torsion angles ϕ and ψ in solution, which were compared to those obtained from molecular dynamics simulations. Modeling using different probability distribution functions showed that MA'AT models of ϕ in βMan(1→4)βXyl and βXyl(1→4)βMan linkages are very similar in the disaccharide building blocks, whereas MA'AT models of ψ differ. This pattern is conserved in the tetrasaccharide, showing that linkage context does not influence linkage geometry in this linear system. Good agreement was observed between the MA'AT and MD models of ψ with respect to mean values and circular standard deviations. Significant differences were observed for ϕ, indicating that revision of the force-field employed by GLYCAM is probably needed. Incorporation of the experimental models of ϕ and ψ into the backbone of an octasaccharide fragment leads to a helical amphipathic topography that may affect the thermal hysteresis properties of the glycolipid.

O-Benzoyl side-chain conformations in 2,3,4,6-tetra-O-benzoyl-β-D-galactopyranosyl-(1→4)-1,2,6-tri-O-benzoyl-β-D-glucopyranose (ethyl acetate solvate) and 1,2,4,6-tetra-O-benzoyl-β-D-glucopyranose (acetone solvate)

The crystal structures of 2,3,4,6-tetra-O-benzoyl-β-D-galactopyranosyl-(1→4)-1,2,6-tri-O-benzoyl-β-D-glucopyranose ethyl acetate hemisolvate, C61H50O18·0.5C4H8O2, and 1,2,4,6-tetra-O-benzoyl-β-D-glucopyranose acetone monosolvate, C34H28O10·C3H6O, were determined and compared to those of methyl β-D-galactopyranosyl-(1→4)-β-D-glucopyranoside (methyl β-lactoside) and methyl β-D-glucopyranoside hemihydrate, C7H14O6·0.5H2O, to evaluate the effects of O-benzoylation on bond lengths, bond angles and torsion angles. In general, O-benzoylation exerts little effect on exo- and endocyclic C—C and endocyclic C—O bond lengths, but exocyclic C—O bonds involved in O-benzoylation are lengthened by 0.02–0.04 Å depending on the site of substitution. The conformation of the O-benzoyl side-chains is highly conserved, with the carbonyl O atom either eclipsing the H atom attached to a 2°-alcoholic C atom or bisecting the H—C—H bond angle of an 1°-alcoholic C atom. Of the three bonds that determine the side-chain geometry, the C—O bond involving the alcoholic C atom exhibits greater rotational variability than the remaining C—O and C—C bonds involving the carbonyl C atom. These findings are in good agreement with recent solution NMR studies of the O-acetyl side-chain conformation in saccharides.

A Solution NMR Approach To Determine the Chemical Structures of Carbohydrates Using the Hydroxyl Groups as Starting Points

An efficient NMR approach is described for determining the chemical structures of the monosaccharide glucose and four disaccharides, namely, nigerose, gentiobiose, leucrose and isomaltulose. This approach uses the 1H resonances of the -OH groups, which are observable in the NMR spectrum of a supercooled aqueous solution, as the starting point for further analysis. The 2D-NMR technique, HSQC-TOCSY, is then applied to fully define the covalent structure (i.e., the topological relationship between C-C, C-H, and O-H bonds) that must be established for a novel carbohydrate before proceeding to further conformational studies. This process also leads to complete assignment of all 1H and 13C resonances. The approach is exemplified by analyzing the monosaccharide glucose, which is treated as if it were an "unknown", and also by fully assigning all the NMR resonances for the four disaccharides that contain glucose. It is proposed that this technique should be equally applicable to the determination of chemical structures for larger carbohydrates of unknown composition, including those that are only available in limited quantities from biological studies. The advantages of commencing the structure elucidation of a carbohydrate at the -OH groups are discussed with reference to the now well-established 2D-/3D-NMR strategy for investigation of peptides/proteins, which employs the -NH resonances as the starting point.

Direct Observation of Carbohydrate Hydroxyl Protons in Hydrogen Bonds with a Protein

Hydroxyl proton resonances of uniformly ¹³C-labeled Manα(1-2)Manα(1-2)ManαOMe (Man₃) bound to cyanovirin-N (CV-N) were detected at ambient temperature in aqueous solution by NMR spectroscopy. The directions of the hydroxyl groups were determined on the basis of NOEs, and a previously unknown hydrogen-bonding network between Man₃ and CV-N was discovered. This is the first report on detecting hydroxyl protons of a protein-bound carbohydrate in aqueous solution by NMR. Approaches such as those presented here may open the door for accurately determining intermolecular hydrogen bonds in carbohydrate-protein complexes.

Glycan OH Exchange Rate Determination in Aqueous Solution: Seeking Evidence for Transient Hydrogen Bonds

Hydrogen bonds (Hbonds) are important stabilizing forces in biomolecules. However, for glycans in aqueous solution, direct NMR detection of Hbonds is elusive because of their transient nature. Here, we present Isotope-based Natural-abundance TOtal correlation eXchange SpectroscopY (INTOXSY), a new ¹H-¹³C heteronuclear single quantum coherence-total correlation spectroscopy based method, to extract OH groups' exchange rate constants (k_ex) for molecules in natural ¹³C abundance and show that OH Hbonds can be inferred from "slower" H/D k_ex. We evaluate k_ex measured with INTOXSY in light of those extracted with line-shape analysis. Subsequently, we use a set of common glycans to establish a k_ex reference basis set and to infer the existence of transient Hbonds involving OH donor groups. Then, we report k_ex values for a series of mono- and disaccharides, as well as for oligosaccharides sialyl Lewis X and β-cyclodextrin, and compare the results with those from the reference set to extract Hbond information. Finally, we utilize NMR experimental data in conjunction with molecular dynamics simulations to establish donor and acceptor Hbond pairs. Our exchange rate measurements indicate that OH/OD exchange rates, k_HD, values <10 s^-1 are consistent with transient Hbond OH groups and potential acceptor groups can be uncovered through MD simulations.

Flexibility at a glycosidic linkage revealed by molecular dynamics, stochastic modeling, and (13)C NMR spin relaxation: conformational preferences of α-L-Rhap-α-(1 → 2)-α-L-Rhap-OMe in water and dimethyl sulfoxide solutions

The monosaccharide L-rhamnose is common in bacterial polysaccharides and the disaccharide α-L-Rhap-α-(1 → 2)-α-L-Rhap-OMe represents a structural model for a part of Shigella flexneri O-antigen polysaccharides. Utilization of [1'-(13)C]-site-specific labeling in the anomeric position at the glycosidic linkage between the two sugar residues facilitated the determination of transglycosidic NMR (3)JCH and (3)JCC coupling constants. Based on these spin-spin couplings the major state and the conformational distribution could be determined with respect to the ψ torsion angle, which changed between water and dimethyl sulfoxide (DMSO) as solvents, a finding mirrored by molecular dynamics (MD) simulations with explicit solvent molecules. The (13)C NMR spin relaxation parameters T1, T2, and heteronuclear NOE of the probe were measured for the disaccharide in DMSO-d6 at two magnetic field strengths, with standard deviations ≤1%. The combination of MD simulation and a stochastic description based on the diffusive chain model resulted in excellent agreement between calculated and experimentally observed (13)C relaxation parameters, with an average error of <2%. The coupling between the global reorientation of the molecule and the local motion of the spin probe is deemed essential if reproduction of NMR relaxation parameters should succeed, since decoupling of the two modes of motion results in significantly worse agreement. Calculation of (13)C relaxation parameters based on the correlation functions obtained directly from the MD simulation of the solute molecule in DMSO as solvent showed satisfactory agreement with errors on the order of 10% or less.

Direct NOE simulation from long MD trajectories

A software package, MD2NOE, is presented which calculates Nuclear Overhauser Effect (NOE) build-up curves directly from molecular dynamics (MD) trajectories. It differs from traditional approaches in that it calculates correlation functions directly from the trajectory instead of extracting inverse sixth power distance terms as an intermediate step in calculating NOEs. This is particularly important for molecules that sample conformational states on a timescale similar to molecular reorientation. The package is tested on sucrose and results are shown to differ in small but significant ways from those calculated using an inverse sixth power assumption. Results are also compared to experiment and found to be in reasonable agreement despite an expected underestimation of water viscosity by the water model selected.

Uncovering Nonconventional and Conventional Hydrogen Bonds in Oligosaccharides through NMR Experiments and Molecular Modeling: Application to Sialyl Lewis-X

We describe the direct NMR detection of a C-H···O nonconventional hydrogen bond (Hbond) and provide experimental and theoretical evidence for conventional Hbonds in the pentasaccharide sialyl Lewis-X (sLe(X)-5) between 5 and 37 °C in water. Extensive NMR structural studies together with molecular dynamics simulations offer strong evidence for significant local dynamics in the Le(X) core and for previously undetected conventional Hbonds in rapid equilibrium that modulate structure. These NMR studies also showed temperature-dependent (1)H and (13)C line broadening. The resulting model emerging from this study is more complex than a simple rigid core description of Le(X)-like molecules and improves our understanding of stabilizing interactions in glycans.

Molecular dynamics simulations and NMR spectroscopy studies of trehalose-lipid bilayer systems

The disaccharide trehalose (TRH) strongly affects the physical properties of lipid bilayers. We investigate interactions between lipid membranes formed by 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and TRH using NMR spectroscopy and molecular dynamics (MD) computer simulations. We compare dipolar couplings derived from DMPC/TRH trajectories with those determined (i) experimentally in TRH using conventional high-resolution NMR in a weakly ordered solvent (bicelles), and (ii) by solid-state NMR in multilamellar vesicles (MLV) formed by DMPC. Analysis of the experimental and MD-derived couplings in DMPC indicated that the force field used in the simulations reasonably well describes the experimental results with the exception for the glycerol fragment that exhibits significant deviations. The signs of dipolar couplings, not available from the experiments on highly ordered systems, were determined from the trajectory analysis. The crucial step in the analysis of residual dipolar couplings (RDCs) in TRH determined in a bicelle-environment was access to the conformational distributions derived from the MD trajectory. Furthermore, the conformational behavior of TRH, investigated by J-couplings, in the ordered and isotropic phases is essentially identical, indicating that the general assumptions in the analyses of RDCs are well founded.

Synthesis and Physicochemical Characterization of D-Tagatose-1-Phosphate: The Substrate of the Tagatose-1-Phosphate Kinase in the Phosphotransferase System-Mediated D-Tagatose Catabolic Pathway of Bacillus licheniformis

We report the first enzymatic synthesis of D-tagatose-1-phosphate (Tag-1P) by the multicomponent phosphoenolpyruvate:sugar phosphotransferase system (PEP-PTS) present in tagatose-grown cells of Klebsiella pneumoniae. Physicochemical characterization by (31)P and (1)H nuclear magnetic resonance spectroscopy reveals that, in solution, this derivative is primarily in the pyranose form. Tag-1P was used to characterize the putative tagatose-1-phosphate kinase (TagK) of the Bacillus licheniformis PTS-mediated D-tagatose catabolic pathway (Bli-TagP). For this purpose, a soluble protein fusion was obtained with the 6 His-tagged trigger factor (TF(His6)) of Escherichia coli. The active fusion enzyme was named TagK-TF(His6). Tag-1P and D-fructose-1-phosphate are substrates for the TagK-TF(His6) enzyme, whereas the isomeric derivatives D-tagatose-6-phosphate and D-fructose-6-phosphate are inhibitors. Studies of catalytic efficiency (kcat/Km) reveal that the enzyme specificity is markedly in favor of Tag-1P as the substrate. Importantly, we show in vivo that the transfer of the phosphate moiety from PEP to the B. licheniformis tagatose-specific Enzyme II in E. coli is inefficient. The capability of the PTS general cytoplasmic components of B. subtilis, HPr and Enzyme I to restore the phosphate transfer is demonstrated.

Methods for measuring exchangeable protons in glycosaminoglycans

Recent NMR studies of the exchangeable protons of GAGs in aqueous solution, including those of the amide, sulfamate, and hydroxyl moieties, have demonstrated potential for the detection of intramolecular hydrogen bonds, providing insights into secondary structure preferences. GAG amide protons are observable by NMR over wide pH and temperature ranges; however, specific solution conditions are required to reduce the exchange rate of the sulfamate and hydroxyl protons and allow their detection by NMR. Building on the vast body of knowledge on detection of hydrogen bonds in peptides and proteins, a variety of methods can be used to identify hydrogen bonds in GAGs including temperature coefficient measurements, evaluation of chemical shift differences between oligo- and monosaccharides, and relative exchange rates measured through line shape analysis and EXSY spectra. Emerging strategies to allow direct detection of hydrogen bonds through heteronuclear couplings offer promise for the future. Molecular dynamic simulations are important in this effort both to predict and confirm hydrogen bond donors and acceptors.

NMR structure analysis of uniformly 13C-labeled carbohydrates

In this study, a set of nuclear magnetic resonance experiments, some of them commonly used in the study of (13)C-labeled proteins and/or nucleic acids, is applied for the structure determination of uniformly (13)C-enriched carbohydrates. Two model substances were employed: one compound of low molecular weight [(UL-(13)C)-sucrose, 342 Da] and one compound of medium molecular weight ((13)C-enriched O-antigenic polysaccharide isolated from Escherichia coli O142, ~10 kDa). The first step in this approach involves the assignment of the carbon resonances in each monosaccharide spin system using the anomeric carbon signal as the starting point. The (13)C resonances are traced using (13)C-(13)C correlations from homonuclear experiments, such as (H)CC-CT-COSY, (H)CC-NOESY, CC-CT-TOCSY and/or virtually decoupled (H)CC-TOCSY. Based on the assignment of the (13)C resonances, the (1)H chemical shifts are derived in a straightforward manner using one-bond (1)H-(13)C correlations from heteronuclear experiments (HC-CT-HSQC). In order to avoid the (1) J CC splitting of the (13)C resonances and to improve the resolution, either constant-time (CT) in the indirect dimension or virtual decoupling in the direct dimension were used. The monosaccharide sequence and linkage positions in oligosaccharides were determined using either (13)C or (1)H detected experiments, namely CC-CT-COSY, band-selective (H)CC-TOCSY, HC-CT-HSQC-NOESY or long-range HC-CT-HSQC. However, due to the short T2 relaxation time associated with larger polysaccharides, the sequential information in the O-antigen polysaccharide from E. coli O142 could only be elucidated using the (1)H-detected experiments. Exchanging protons of hydroxyl groups and N-acetyl amides in the (13)C-enriched polysaccharide were assigned by using HC-H2BC spectra. The assignment of the N-acetyl groups with (15)N at natural abundance was completed by using HN-SOFAST-HMQC, HNCA, HNCO and (13)C-detected (H)CACO spectra.

NMR of glycans: shedding new light on old problems

The diversity in molecular arrangements and dynamics displayed by glycans renders traditional NMR strategies, employed for proteins and nucleic acids, insufficient. Because of the unique properties of glycans, structural studies often require the adoption of a different repertoire of tailor-made experiments and protocols. We present an account of recent developments in NMR techniques that will deepen our understanding of structure-function relations in glycans. We open with a survey and comparison of methods utilized to determine the structure of proteins, nucleic acids and carbohydrates. Next, we discuss the structural information obtained from traditional NMR techniques like chemical shifts, NOEs/ROEs, and coupling-constants, along with the limitations imposed by the unique intrinsic characteristics of glycan structure on these approaches: flexibility, range of conformers, signal overlap, and non-first-order scalar (strong) coupling. Novel experiments taking advantage of isotopic labeling are presented as an option for overcoming spectral overlap and raising sensitivity. Computational tools used to explore conformational averaging in conjunction with NMR parameters are described. In addition, recent developments in hydroxyl detection and hydrogen bond detection in protonated solvents, in contrast to traditional sample preparations in D2O for carbohydrates, further increase the tools available for both structure information and chemical shift assignments. We also include previously unpublished data in this context. Accurate determination of couplings in carbohydrates has been historically challenging due to the common presence of strong-couplings. We present new strategies proposed for dealing with their influence on NMR signals. We close with a discussion of residual dipolar couplings (RDCs) and the advantages of using (13)C isotope labeling that allows gathering one-bond (13)C-(13)C couplings with a recently improved constant-time COSY technique, in addition to the commonly measured (1)H-(13)C RDCs.

Hydroxyl-proton hydrogen bonding in the heparin oligosaccharide Arixtra in aqueous solution

Heparin is best known for its anticoagulant activity, which is mediated by the binding of a specific pentasaccharide sequence to the protease inhibitor antithrombin-III (AT-III). Although heparin oligosaccharides are thought to be flexible in aqueous solution, the recent discovery of a hydrogen bond between the sulfamate (NHSO3(-)) proton and the adjacent 3-O-sulfo group of the 3,6-O-sulfated N-sulfoglucosamine residue of the Arixtra (fondaparinux sodium) pentasaccharide demonstrates that definable elements of local structure are accessed. Molecular dynamics simulations of Arixtra suggest the presence of additional hydrogen bonds involving the C3-OH groups of the glucuronic acid and 2-O-sulfo-iduronic acid residues. NMR measurements of temperature coefficients, chemical shift differences, and solvent exchange rate constants provide experimental confirmation of these hydrogen bonds. We note that the extraction of rate constants from cross-peak buildup curves in 2D exchange spectroscopy is complicated by the presence of radiation damping in aqueous solution. A straightforward model is presented that explicitly takes into account the effects of radiation damping on the water proton relaxation and is sufficiently robust to provide an accurate measure of the proton exchange rate between the analyte hydroxyl protons and water.