Structure and conformation of complex carbohydrates of glycoproteins, glycolipids, and bacterial polysaccharides

Taking the Scenic Route: Polyomaviruses Utilize Multiple Pathways to Reach the Same Destination

Members of the Polyomaviridae family differ in their host range, pathogenesis, and disease severity. To date, some of the most studied polyomaviruses include human JC, BK, and Merkel cell polyomavirus and non-human subspecies murine and simian virus 40 (SV40) polyomavirus. Although dichotomies in host range and pathogenesis exist, overlapping features of the infectious cycle illuminate the similarities within this virus family. Of particular interest to human health, JC, BK, and Merkel cell polyomavirus have all been linked to critical, often fatal, illnesses, emphasizing the importance of understanding the underlying viral infections that result in the onset of these diseases. As there are significant overlaps in the capacity of polyomaviruses to cause disease in their respective hosts, recent advancements in characterizing the infectious life cycle of non-human murine and SV40 polyomaviruses are key to understanding diseases caused by their human counterparts. This review focuses on the molecular mechanisms by which different polyomaviruses hijack cellular processes to attach to host cells, internalize, traffic within the cytoplasm, and disassemble within the endoplasmic reticulum (ER), prior to delivery to the nucleus for viral replication. Unraveling the fundamental processes that facilitate polyomavirus infection provides deeper insight into the conserved mechanisms of the infectious process shared within this virus family, while also highlighting critical unique viral features.

Interaction of poly-l-lysine coating and heparan sulfate proteoglycan on magnetic nanoparticle uptake by tumor cells

Background

Poly-l-lysine (PLL) enhances nanoparticle (NP) uptake, but the molecular mechanism remains unresolved. We asked whether PLL may interact with negatively charged glycoconjugates on the cell surface and facilitate uptake of magnetic NPs (MNPs) by tumor cells.

Methods

PLL-coated MNPs (PLL-MNPs) with positive and negative ζ-potential were prepared and characterized. Confocal and transmission electron microscopy was used to analyze cellular internalization of MNPs. A colorimetric iron assay was used to quantitate cell-associated MNPs (MNP_cell).

Results

Coadministration of PLL and dextran-coated MNPs in culture enhanced cellular internalization of MNPs, with increased vesicle size and numbers/cell. MNP_cell was increased by eight- to 12-fold in response to PLL in a concentration-dependent manner in human glioma and HeLa cells. However, the application of a magnetic field attenuated PLL-induced increase in MNP_cell. PLL-coating increased MNP_cell regardless of ζ-potential of PLL-MNPs, whereas magnetic force did not enhance MNP_cell. In contrast, epigallocatechin gallate and magnetic force synergistically enhanced PLL-MNP uptake. In addition, heparin, but not sialic acid, greatly reduced the enhancement effects of PLL; however, removal of heparan sulfate from heparan sulfate proteoglycans of the cell surface by heparinase III significantly reduced MNP_cell.

Conclusion

Our results suggest that PLL-heparan sulfate proteoglycan interaction may be the first step mediating PLL-MNP internalization by tumor cells. Given these results, PLL may facilitate NP interaction with tumor cells via a molecular mechanism shared by infection machinery of certain viruses.

The β-reducing end in α(2-8)-polysialic acid constitutes a unique structural motif

Over the years, structural characterizations of α(2-8)-polysialic acid (polySia) in solution have produced inconclusive results. Efforts for obtaining detailed information in this important antigen have focused primarily on the α-linked residues and not on the distinctive characteristics of the terminal ones. The thermodynamically preferred anomeric configuration for the reducing end of sialic acids is β, which has the [I]CO2- group equatorial and the OH ([I]OH2) axial, while for all other residues the CO2- group is axial. We show that this purportedly minor difference has distinct consequences for the structure of α(2-8)-polySia near the reducing end, as the β configuration places the [I]OH2 in a favorable position for the formation of a hydrogen bond with the carboxylate group of the following residue ([II]CO2-). Molecular dynamics (MD) simulations predicted the hydrogen bond, which we subsequently directly detected by NMR. The combination of MD and residual dipolar couplings shows that the net result for the structure of Sia2-βOH is a stable conformation with well-defined hydration and charge patterns, and consistent with experimental NOE-based hydroxyl and aliphatic inter-proton distances. Moreover, we provide evidence that this distinct conformation is preserved on Sia oligosaccharides, thus constituting a motif that determines the structure and dynamics of α(2-8)-polySia for at least the first two residues of the polymer. We suggest the hypothesis that this structural motif sheds light on a longtime puzzling observation for the requirement of 10 residues of α(2-8)-polySia in order to bind effectively to specific antibodies, about four units more than for analogous cases.

The sugar code: Why glycans are so important

The cell surface is the platform for presentation of biochemical signals that are required for intercellular communication. Their profile necessarily needs to be responsive to internal and external factors in a highly dynamic manner. The structural features of the signals must meet the criterion of high-density information coding in a minimum of space. Thus, only biomolecules that can generate many different oligomers ('words') from few building blocks ('letters') qualify to meet this challenge. Examining the respective properties of common biocompounds that form natural oligo- and polymers comparatively, starting with nucleotides and amino acids (the first and second alphabets of life), comes up with sugars as clear frontrunner. The enzymatic machinery for the biosynthesis of sugar chains can indeed link monosaccharides, the letters of the third alphabet of life, in a manner to reach an unsurpassed number of oligomers (complex carbohydrates or glycans). Fittingly, the resulting glycome of a cell can be likened to a fingerprint. Conjugates of glycans with proteins and sphingolipids (glycoproteins and glycolipids) are ubiquitous in Nature. This implies a broad (patho)physiologic significance. By looking at the signals, at the writers and the erasers of this information as well as its readers and ensuing consequences, this review intends to introduce a broad readership to the principles of the concept of the sugar code.

Structure of a protective epitope of group B Streptococcus type III capsular polysaccharide

Despite substantial progress in the prevention of group B Streptococcus (GBS) disease with the introduction of intrapartum antibiotic prophylaxis, this pathogen remains a leading cause of neonatal infection. Capsular polysaccharide conjugate vaccines have been tested in phase I/II clinical studies, showing promise for further development. Mapping of epitopes recognized by protective antibodies is crucial for understanding the mechanism of action of vaccines and for enabling antigen design. In this study, we report the structure of the epitope recognized by a monoclonal antibody with opsonophagocytic activity and representative of the protective response against type III GBS polysaccharide. The structure and the atomic-level interactions were determined by saturation transfer difference (STD)-NMR and X-ray crystallography using oligosaccharides obtained by synthetic and depolymerization procedures. The GBS PSIII epitope is made by six sugars. Four of them derive from two adjacent repeating units of the PSIII backbone and two of them from the branched galactose-sialic acid disaccharide contained in this sequence. The sialic acid residue establishes direct binding interactions with the functional antibody. The crystal structure provides insight into the molecular basis of antibody-carbohydrate interactions and confirms that the conformational epitope is not required for antigen recognition. Understanding the structural basis of immune recognition of capsular polysaccharide epitopes can aid in the design of novel glycoconjugate vaccines.

Vibrational Signatures of Isomeric Lithiated N-acetyl-D-hexosamines by Gas-Phase Infrared Multiple-Photon Dissociation (IRMPD) Spectroscopy

Three lithiated N-acetyl-D-hexosamine (HexNAc) isomers, N-acetyl-D-glucosamine (GlcNAc), N-acetyl-D-galactosamine (GalNAc), and N-acetyl-D-mannosamine (ManNAc) are investigated as model monosaccharide derivatives by gas-phase infrared multiple-photon dissociation (IRMPD) spectroscopy. The hydrogen stretching region, which is attributed to OH and NH stretching modes, reveals some distinguishing spectral features of the lithium-adducted complexes that are useful in terms of differentiating these isomers. In order to understand the effect of lithium coordination on saccharide structure, and therefore anomericity, chair configuration, and hydrogen bonding networks, the conformational preferences of lithiated GlcNAc, GalNAc, and ManNAc are studied by comparing the experimental measurements with density functional theory (DFT) calculations. The experimental results of lithiated GlcNAc and GalNAc show a good match to the theoretical spectra of low-energy structures adopting a ⁴ C ₁ chair conformation, consistent with this motif being the dominant conformation in condensed-phase monosaccharides. The epimerization effect upon going to lithiated ManNAc is significant, as in this case the ¹ C ₄ chair conformers give a more compelling match with the experimental results, consistent with their lower calculated energies. A contrasting computational study of these monosaccharides in their neutral form suggests that the lithium cation coordination with Lewis base oxygens can play a key role in favoring particular structural motifs (e.g., a ⁴ C ₁ versus ¹ C ₄ ) and disrupting hydrogen bond networks, thus exhibiting specific IR spectral features between these closely related lithium-chelated complexes. Graphical Abstract ᅟ.

Biomimetic mineralization of metal–organic frameworks around polysaccharides

Biomimetic mineralization exploits natural biomineralization processes for the design and fabrication of synthetic functional materials.

Ganglioside-Lipid and Ganglioside-Protein Interactions Revealed by Coarse-Grained and Atomistic Molecular Dynamics Simulations

Gangliosides are glycolipids in which an oligosaccharide headgroup containing one or more sialic acids is connected to a ceramide. Gangliosides reside in the outer leaflet of the plasma membrane and play a crucial role in various physiological processes such as cell signal transduction and neuronal differentiation by modulating structures and functions of membrane proteins. Because the detailed behavior of gangliosides and protein-ganglioside interactions are poorly known, we investigated the interactions between the gangliosides GM1 and GM3 and the proteins aquaporin (AQP1) and WALP23 using equilibrium molecular dynamics simulations and potential of mean force calculations at both coarse-grained (CG) and atomistic levels. In atomistic simulations, on the basis of the GROMOS force field, ganglioside aggregation appears to be a result of the balance between hydrogen bond interactions and steric hindrance of the headgroups. GM3 clusters are slightly larger and more ordered than GM1 clusters due to the smaller headgroup of GM3. The different structures of GM1 and GM3 clusters from atomistic simulations are not observed at the CG level based on the Martini model, implying a difference in driving forces for ganglioside interactions in atomistic and CG simulations. For protein-ganglioside interactions, in the atomistic simulations, GM1 lipids bind to specific sites on the AQP1 surface, whereas they are depleted from WALP23. In the CG simulations, the ganglioside binding sites on the AQP1 surface are similar, but ganglioside aggregation and protein-ganglioside interactions are more prevalent than in the atomistic simulations. Using the polarizable Martini water model, results were closer to the atomistic simulations. Although experimental data for validation is lacking, we proposed modified Martini parameters for gangliosides to more closely mimic the sizes and structures of ganglioside clusters observed at the atomistic level.

Perspectives on Anti-Glycan Antibodies Gleaned from Development of a Community Resource Database

Antibodies are used extensively for a wide range of basic research and clinical applications. While an abundant and diverse collection of antibodies to protein antigens have been developed, good monoclonal antibodies to carbohydrates are much less common. Moreover, it can be difficult to determine if a particular antibody has the appropriate specificity, which antibody is best suited for a given application, and where to obtain that antibody. Herein, we provide an overview of the current state of the field, discuss challenges for selecting and using antiglycan antibodies, and summarize deficiencies in the existing repertoire of antiglycan antibodies. This perspective was enabled by collecting information from publications, databases, and commercial entities and assembling it into a single database, referred to as the Database of Anti-Glycan Reagents (DAGR). DAGR is a publicly available, comprehensive resource for anticarbohydrate antibodies, their applications, availability, and quality.

Multivalent Carbohydrate-Lectin Interactions: How Synthetic Chemistry Enables Insights into Nanometric Recognition

Glycan recognition by sugar receptors (lectins) is intimately involved in many aspects of cell physiology. However, the factors explaining the exquisite selectivity of their functional pairing are not yet fully understood. Studies toward this aim will also help appraise the potential for lectin-directed drug design. With the network of adhesion/growth-regulatory galectins as therapeutic targets, the strategy to recruit synthetic chemistry to systematically elucidate structure-activity relationships is outlined, from monovalent compounds to glyco-clusters and glycodendrimers to biomimetic surfaces. The versatility of the synthetic procedures enables to take examining structural and spatial parameters, alone and in combination, to its limits, for example with the aim to produce inhibitors for distinct galectin(s) that exhibit minimal reactivity to other members of this group. Shaping spatial architectures similar to glycoconjugate aggregates, microdomains or vesicles provides attractive tools to disclose the often still hidden significance of nanometric aspects of the different modes of lectin design (sequence divergence at the lectin site, differences of spatial type of lectin-site presentation). Of note, testing the effectors alone or in combination simulating (patho)physiological conditions, is sure to bring about new insights into the cooperation between lectins and the regulation of their activity.

Linkage and anomeric differentiation in trisaccharides by sequential fragmentation and variable-wavelength infrared photodissociation

Carbohydrates and their derivatives play important roles in biological systems, but their isomeric heterogeneity also presents a considerable challenge for analytical techniques. Here, a stepwise approach using infrared multiple-photon dissociation (IRMPD) via a tunable CO2 laser (9.2-10.7 μm) was employed to characterize isomeric variants of glucose-based trisaccharides. After the deprotonated trisaccharides were trapped and fragmented to disaccharide C2 fragments in a Fourier transform ion cyclotron resonance (FTICR) cell, a further variable-wavelength infrared irradiation of the C2 ion produced wavelength-dependent dissociation patterns that are represented as heat maps. The photodissociation patterns of these C2 fragments are shown to be strikingly similar to the photodissociation patterns of disaccharides with identical glycosidic bonds. Conversely, the photodissociation patterns of different glycosidic linkages exhibit considerable differences. On the basis of these results, the linkage position and anomericity of glycosidic bonds of disaccharide units in trisaccharides can be systematically differentiated and identified, providing a promising approach to characterize the structures of isomeric oligosaccharides.

A tensor-free method for the structural and dynamical refinement of proteins using residual dipolar couplings

Residual dipolar couplings (RDCs) are parameters measured in nuclear magnetic resonance spectroscopy that can provide exquisitely detailed information about the structure and dynamics of biological macromolecules. We describe here a method of using RDCs for the structural and dynamical refinement of proteins that is based on the observation that the RDC between two atomic nuclei depends directly on the angle ϑ between the internuclear vector and the external magnetic field. For every pair of nuclei for which an RDC is available experimentally, we introduce a structural restraint to minimize the deviation from the value of the angle ϑ derived from the measured RDC and that calculated in the refinement protocol. As each restraint involves only the calculation of the angle ϑ of the corresponding internuclear vector, the method does not require the definition of an overall alignment tensor to describe the preferred orientation of the protein with respect to the alignment medium. Application to the case of ubiquitin demonstrates that this method enables an accurate refinement of the structure and dynamics of this protein to be obtained.

Glycans - the third revolution in evolution

The development and maintenance of a complex organism composed of trillions of cells is an extremely complex task. At the molecular level every process requires a specific molecular structures to perform it, thus it is difficult to imagine how less than tenfold increase in the number of genes between simple bacteria and higher eukaryotes enabled this quantum leap in complexity. In this perspective article we present the hypothesis that the invention of glycans was the third revolution in evolution (the appearance of nucleic acids and proteins being the first two), which enabled the creation of novel molecular entities that do not require a direct genetic template. Contrary to proteins and nucleic acids, which are made from a direct DNA template, glycans are product of a complex biosynthetic pathway affected by hundreds of genetic and environmental factors. Therefore glycans enable adaptive response to environmental changes and, unlike other epiproteomic modifications, which act as off/on switches, glycosylation significantly contributes to protein structure and enables novel functions. The importance of glycosylation is evident from the fact that nearly all proteins invented after the appearance of multicellular life are composed of both polypeptide and glycan parts.

The power of coarse graining in biomolecular simulations

Computational modeling of biological systems is challenging because of the multitude of spatial and temporal scales involved. Replacing atomistic detail with lower resolution, coarse grained (CG), beads has opened the way to simulate large-scale biomolecular processes on time scales inaccessible to all-atom models. We provide an overview of some of the more popular CG models used in biomolecular applications to date, focusing on models that retain chemical specificity. A few state-of-the-art examples of protein folding, membrane protein gating and self-assembly, DNA hybridization, and modeling of carbohydrate fibers are used to illustrate the power and diversity of current CG modeling.

A guide into glycosciences: How chemistry, biochemistry and biology cooperate to crack the sugar code

BACKGROUND

The most demanding challenge in research on molecular aspects within the flow of biological information is posed by the complex carbohydrates (glycan part of cellular glycoconjugates). How the 'message' encoded in carbohydrate 'letters' is 'read' and 'translated' can only be unraveled by interdisciplinary efforts.

SCOPE OF REVIEW

This review provides a didactic step-by-step survey of the concept of the sugar code and the way strategic combination of experimental approaches characterizes structure-function relationships, with resources for teaching.

MAJOR CONCLUSIONS

The unsurpassed coding capacity of glycans is an ideal platform for generating a broad range of molecular 'messages'. Structural and functional analyses of complex carbohydrates have been made possible by advances in chemical synthesis, rendering production of oligosaccharides, glycoclusters and neoglycoconjugates possible. This availability facilitates to test the glycans as ligands for natural sugar receptors (lectins). Their interaction is a means to turn sugar-encoded information into cellular effects. Glycan/lectin structures and their spatial modes of presentation underlie the exquisite specificity of the endogenous lectins in counterreceptor selection, that is, to home in on certain cellular glycoproteins or glycolipids.

GENERAL SIGNIFICANCE

Understanding how sugar-encoded 'messages' are 'read' and 'translated' by lectins provides insights into fundamental mechanisms of life, with potential for medical applications.

Complex carbohydrates as a possible source of high energy to formulate functional feeds

Carbohydrates (CHOs) are the most abundant organic compounds found in living organisms and are a great source of metabolic energy, both for plants and animals. Besides of CHOs great potential to solve animal's energy requirements and diminishing high feed cost, we first must to understand its digestibility and assimilation to avoid several inconvenients. Today, CHOs feed animal inclusions are of great concern about cost-benefits, animal's health status, and environmental pollution. In this chapter, we make a brief description about sugars (DP1-2), oligosaccharides (DP3-9), polysaccharides (DP ≥10), and their essential characteristics to understand the role of marine and terrestrial CHOs in animal nutrition. Subsequently, we talk about basic concepts, CHOs functional benefits, suggestions about their application and successful cases. This information will contribute to produce a new generation of high-quality and energetic functional feed formulations for livestock and aquaculture farms; which must be of low cost, healthy, and environmentally friendly, with the inclusion of prebiotics and probiotics.

BSA-boronic acid conjugate as lectin mimetics

We report bovine serum albumin (BSA)-boronic acid (BA) conjugates as lectin mimetics and their glyco-capturing capacity. The BSA-BA conjugates were synthesized by amidation of carboxylic acid groups in BSA with aminophenyl boronic acid in the presence of EDC, and were characterized by Alizarin Red S (ARS) assay and SDS-PAGE gel. The BSA-BA conjugates were immobilized onto maleimide-functionalized silica beads and their sugar capturing capacity and specificity were confirmed by ARS displacement assay. Further, surface plasmon resonance (SPR) analysis of the glyco-capturing activity of the BSA-BA conjugates was conducted by immobilizing BSA-BA onto SPR gold chip. Overall, we demonstrated a BSA-BA-based lectin mimetics for glyco-capturing applications. These lectin mimetics are expected to provide an important tool for glycomics and biosensor research and applications.

Conformation and dynamics at a flexible glycosidic linkage revealed by NMR spectroscopy and molecular dynamics simulations: analysis of β-L-Fucp-(1→6)-α-D-Glcp-OMe in water solution

The intrinsic flexibility of carbohydrates facilitates different 3D structures in response to altered environments. At glycosidic (1→6)-linkages, three torsion angles are variable, and herein the conformation and dynamics of β-L-Fucp-(1→6)-α-D-Glcp-OMe are investigated using a combination of NMR spectroscopy and molecular dynamics (MD) simulations. The disaccharide shows evidence of conformational averaging for the ψ and ω torsion angles, best explained by a four-state conformational distribution. Notably, there is a significant population of conformations having ψ = 85° (clinal) in addition to those having ψ = 180° (antiperiplanar). Moderate differences in (13)C R1 relaxation rates are found to be best explained by axially symmetric tumbling in combination with minor differences in librational motion for the two residues, whereas the isomerization motions are occurring too slowly to be contributing significantly to the observed relaxation rates. The MD simulation was found to give a reasonably good agreement with experiment, especially with respect to diffusive properties, among which the rotational anisotropy, D∥/D⊥, is found to be 2.35. The force field employed showed too narrow ω torsion angles in the gauche-trans and gauche-gauche states as well as overestimating the population of the gauche-trans conformer. This information can subsequently be used in directing parameter developments and emphasizes the need for refinement of force fields for (1→6)-linked carbohydrates.

Shaping up for structural glycomics: a predictive protocol for oligosaccharide conformational analysis applied to N-linked glycans

•

Aqueous 10 μs simulations of N-linked mannosyl cores and sialyl Lewis (sLe) antennae are validated.

•

Sequence dependent glycosidic linkage and pyranose ring μs motions are implicated in bioactivity.

•

Stacked pyranoses in sLe^a and sLe^x are predicted to be atypically rigid on μs timescales.

•

In a 25 μs simulation of sLe^x, all known conformers were sampled within the initial 10 μs of dynamics.

•

Unbiased 10 μs simulations are proposed as a route to systematic and accurate glycomic 3D-analysis.

The human glycome comprises a vast untapped repository of 3D-structural information that holds the key to glycan recognition and a new era of rationally designed mimetic chemical probes, drugs, and biomaterials. Toward routine prediction of oligosaccharide conformational populations and exchange rates at thermodynamic equilibrium, we apply hardware-accelerated aqueous molecular dynamics to model μs motions in N-glycans that underpin inflammation and immunity. In 10 μs simulations, conformational equilibria of mannosyl cores, sialyl Lewis (sLe) antennae, and constituent sub-sequences agreed with prior refinements (X-ray and NMR). Glycosidic linkage and pyranose ring flexing were affected by branching, linkage position, and secondary structure, implicating sequence dependent motions in glycomic functional diversity. Linkage and ring conformational transitions that have eluded precise quantification by experiment and conventional (ns) simulations were predicted to occur on μs timescales. All rings populated non-chair shapes and the stacked galactose and fucose pyranoses of sLe^a and sLe^x were rigidified, suggesting an exploitable 3D-signature of cell adhesion protein binding. Analyses of sLe^x dynamics over 25 μs revealed that only 10 μs were sufficient to explore all aqueous conformers. This simulation protocol, which yields conformational ensembles that are independent of initial 3D-structure, is proposed as a route to understanding oligosaccharide recognition and structure–activity relationships, toward development of carbohydrate-based novel chemical entities.

Understanding the specificity of human Galectin-8C domain interactions with its glycan ligands based on molecular dynamics simulations

Human Galectin-8 (Gal-8) is a member of the galectin family which shares an affinity for β-galactosides. The tandem-repeat Gal-8 consists of a N- and a C-terminal carbohydrate recognition domain (N- and C-CRD) joined by a linker peptide of various length. Despite their structural similarity both CRDs recognize different oligosaccharides. While the molecular requirements of the N-CRD for high binding affinity to sulfated and sialylated glycans have recently been elucidated by crystallographic studies of complexes with several oligosaccharides, the binding specificities of the C-CRD for a different set of oligosaccharides, as derived from experimental data, has only been explained in terms of the three-dimensional structure for the complex C-CRD with lactose. In this study we performed molecular dynamics (MD) simulations using the recently released crystal structure of the Gal-8C-CRD to analyse the three-dimensional conditions for its specific binding to a variety of oligosaccharides as previously defined by glycan-microarray analysis. The terminal β-galactose of disaccharides (LacNAc, lacto-N-biose and lactose) and the internal β-galactose moiety of blood group antigens A and B (BGA, BGB) as well as of longer linear oligosaccharide chains (di-LacNAc and lacto-N-neotetraose) are interacting favorably with conserved amino acids (H53, R57, N66, W73, E76). Lacto-N-neotetraose and di-LacNAc as well as BGA and BGB are well accommodated. BGA and BGB showed higher affinity than LacNAc and lactose due to generally stronger hydrogen bond interactions and water mediated hydrogen bonds with α1-2 fucose respectively. Our results derived from molecular dynamics simulations are able to explain the glycan binding specificities of the Gal-8C-CRD in comparison to those of the Gal-8N -CRD.

Separation of isomeric disaccharides by traveling wave ion mobility mass spectrometry using CO2 as drift gas

The use of CO(2) as a massive and polarizable drift gas is shown to greatly improve peak-to-peak resolution (R(p-p) ), as compared with N(2) , for the separation of disaccharides in a Synapt G2 traveling wave ion mobility cell. Near or baseline R(p-p) was achieved for three pairs of sodiated molecules of disaccharide isomers, that is, cellobiose and sucrose (R(p-p)  = 0.76), maltose and sucrose (R(p-p)  = 1.04), and maltose and lactose (R(p-p)  = 0.74). Ion mobility mass spectrometry using CO(2) as the drift gas offers therefore an attractive alternative for fast and efficient separation of isomeric disaccharides.

Carbohydrate force fields

Carbohydrates present a special set of challenges to the generation of force fields. First, the tertiary structures of monosaccharides are complex merely by virtue of their exceptionally high number of chiral centers. In addition, their electronic characteristics lead to molecular geometries and electrostatic landscapes that can be challenging to predict and model. The monosaccharide units can also interconnect in many ways, resulting in a large number of possible oligosaccharides and polysaccharides, both linear and branched. These larger structures contain a number of rotatable bonds, meaning they potentially sample an enormous conformational space. This article briefly reviews the history of carbohydrate force fields, examining and comparing their challenges, forms, philosophies, and development strategies. Then it presents a survey of recent uses of these force fields, noting trends, strengths, deficiencies, and possible directions for future expansion.

Differentiation of the stereochemistry and anomeric configuration for 1-3 linked disaccharides via tandem mass spectrometry and 18O-labeling

Collision-induced dissociation (CID) of deprotonated hexose-containing disaccharides (m/z 341) with 1-2, 1-4, and 1-6 linkages yields product ions at m/z 221, which have been identified as glycosyl-glycolaldehyde anions. From disaccharides with these linkages, CID of m/z 221 ions produces distinct fragmentation patterns that enable the stereochemistries and anomeric configurations of the non-reducing sugar units to be determined. However, only trace quantities of m/z 221 ions can be generated for 1-3 linkages in Paul or linear ion traps, preventing further CID analysis. Here we demonstrate that high intensities of m/z 221 ions can be built up in the linear ion trap (Q3) from beam-type CID of a series of 1-3 linked disaccharides conducted on a triple quadrupole/linear ion trap mass spectrometer. (18)O-labeling at the carbonyl position of the reducing sugar allowed mass-discrimination of the "sidedness" of dissociation events to either side of the glycosidic linkage. Under relatively low energy beam-type CID and ion trap CID, an m/z 223 product ion containing (18)O predominated. It was a structural isomer that fragmented quite differently than the glycosyl-glycolaldehydes and did not provide structural information about the non-reducing sugar. Under higher collision energy beam-type CID conditions, the formation of m/z 221 ions, which have the glycosyl-glycolaldehyde structures, were favored. Characteristic fragmentation patterns were observed for each m/z 221 ion from higher energy beam-type CID of 1-3 linked disaccharides and the stereochemistry of the non-reducing sugar, together with the anomeric configuration, were successfully identified both with and without (18)O-labeling of the reducing sugar carbonyl group.

Solvation properties of N-acetyl-β-glucosamine: molecular dynamics study incorporating electrostatic polarization

N-Acetyl-β-glucosamine (NAG) is an important moiety of glycoproteins and is involved in many biological functions. However, conformational and dynamical properties of NAG molecules in aqueous solution, the most common biological environment, remain ambiguous due to limitations of experimental methods. Increasing efforts are made to probe structural properties of NAG and NAG-containing macromolecules, like peptidoglycans and polymeric chitin, at the atomic level using molecular dynamics simulations. In this work, we develop a polarizable carbohydrate force field for NAG and contrast simulation results of various properties using this novel force field and an analogous nonpolarizable (fixed charge) model. Aqueous solutions of NAG and its oligomers are investigated; we explore conformational properties (rotatable bond geometry), electrostatic properties (dipole moment distribution), dynamical properties (self-diffusion coefficient), hydrogen bonding (water bridge structure and dynamics), and free energy of hydration. The fixed-charge carbohydrate force field exhibits deviations from the gas phase relative rotation energy of exocyclic hydroxymethyl side chain and of chair/boat ring distortion. The polarizable force field predicts conformational properties in agreement with corresponding first-principles results. NAG-water hydrogen bonding pattern is studied through radial distribution functions (RDFs) and correlation functions. Intermolecular hydrogen bonding between solute and solvent is found to stabilize NAG solution structures while intramolecular hydrogen bonds define glycosidic linkage geometry of NAG oligomers. The electrostatic component of hydration free energy is highly dependent on force field atomic partial charges, influencing a more favorable free energy of hydration in the fixed-charge model compared to the polarizable model.

Structure and binding analysis of Polyporus squamosus lectin in complex with the Neu5Ac{alpha}2-6Gal{beta}1-4GlcNAc human-type influenza receptor

Glycan chains that terminate in sialic acid (Neu5Ac) are frequently the receptors targeted by pathogens for initial adhesion. Carbohydrate-binding proteins (lectins) with specificity for Neu5Ac are particularly useful in the detection and isolation of sialylated glycoconjugates, such as those associated with pathogen adhesion as well as those characteristic of several diseases including cancer. Structural studies of lectins are essential in order to understand the origin of their specificity, which is particularly important when employing such reagents as diagnostic tools. Here, we report a crystallographic and molecular dynamics (MD) analysis of a lectin from Polyporus squamosus (PSL) that is specific for glycans terminating with the sequence Neu5Acα2-6Galβ. Because of its importance as a histological reagent, the PSL structure was solved (to 1.7 Å) in complex with a trisaccharide, whose sequence (Neu5Acα2-6Galβ1-4GlcNAc) is exploited by influenza A hemagglutinin for viral adhesion to human tissue. The structural data illuminate the origin of the high specificity of PSL for the Neu5Acα2-6Gal sequence. Theoretical binding free energies derived from the MD data confirm the key interactions identified crystallographically and provide additional insight into the relative contributions from each amino acid, as well as estimates of the importance of entropic and enthalpic contributions to binding.

Tuning the cavity of cyclodextrins: altered sugar adaptors in protein pores

Cyclodextrins (CDs) have been widely used in host-guest molecular recognition because of their chiral and hydrophobic cavities. For example, β-cyclodextrin (βCD) lodged as a molecular adaptor in protein pores such as α-hemolysin (αHL) is used for stochastic sensing. Here, we have tuned the cavity and overall size of βCD by replacing a single oxygen atom in its ring skeleton by a disulfide unit in two different configurations to both expand our ability to detect analytes and understand the interactions of βCD with protein pores. The three-dimensional structures of the two stereoisomeric CDs have been determined by the combined application of NMR spectroscopy and molecular simulation and show distorted conformations as compared to natural βCD. The interactions of these synthetic βCD analogues with mutant αHL protein pores and guest molecules were studied by single-channel electrical recording. The dissociation rate constants for both disulfide CDs from the mutant pores show ∼1000-fold increase as compared to those of unaltered βCD, but are ∼10-fold lower than the dissociation rate constants for βCD from wild-type αHL. Both of the skeleton-modified CDs show altered selectivity toward guest molecules. Our approach expands the breadth in sensitivity and diversity of sensing with protein pores and suggests structural parameters useful for CD design, particularly in the creation of asymmetric cavities.

A molecular dynamics study of the inclusion complexes of C60 with some cyclodextrins

The strongly hydrophobic C(60) fullerene is a carbon allotrope of huge interest in materials science and in pharmaceutical chemistry that can be solubilized in water either by extensive chemical functionalization or by inclusion in appropriate carriers such as the cyclodextrins with formation of host-guest complexes. Here we report a molecular dynamics study of the complexes formed in solution by C(60) with gamma- and delta-cyclodextrins. The most stable host-guest complex stoichiometry is determined to be 2:1 through simulations in vacuo and in explicit water by the stepwise addition of the cyclodextrins to C(60). No a priori assumption about the inclusion stoichiometry and geometry is made. The equilibrium fluctuations of the complexes that can affect the system stability are also investigated within the molecular dynamics runs.

Conformational properties of glucose-based disaccharides investigated using molecular dynamics simulations with local elevation umbrella sampling

Explicit-solvent molecular dynamics (MD) simulations of the 11 glucose-based disaccharides in water at 300K and 1bar are reported. The simulations were carried out with the GROMOS 45A4 force-field and the sampling along the glycosidic dihedral angles phi and psi was artificially enhanced using the local elevation umbrella sampling (LEUS) method. The trajectories are analyzed in terms of free-energy maps, stable and metastable conformational states (relative free energies and estimated transition timescales), intramolecular H-bonds, single molecule configurational entropies, and agreement with experimental data. All disaccharides considered are found to be characterized either by a single stable (overwhelmingly populated) state ((1-->n)-linked disaccharides with n=1, 2, 3, or 4) or by two stable (comparably populated and differing in the third glycosidic dihedral angle omega ; gg or gt) states with a low interconversion barrier ((1-->6)-linked disaccharides). Metastable (anti-phi or anti-psi) states are also identified with relative free energies in the range of 8-22 kJ mol(-1). The 11 compounds can be classified into four families: (i) the alpha(1-->1)alpha-linked disaccharide trehalose (axial-axial linkage) presents no metastable state, the lowest configurational entropy, and no intramolecular H-bonds; (ii) the four alpha(1-->n)-linked disaccharides (n=1, 2, 3, or 4; axial-equatorial linkage) present one metastable (anti-psi) state, an intermediate configurational entropy, and two alternative intramolecular H-bonds; (iii) the four beta(1-->n)-linked disaccharides (n=1, 2, 3, or 4; equatorial-equatorial linkage) present two metastable (anti-phi and anti-psi) states, an intermediate configurational entropy, and one intramolecular H-bond; (iv) the two (1-->6)-linked disaccharides (additional glycosidic dihedral angle) present no (isomaltose) or a pair of (gentiobiose) metastable (anti-phi) states, the highest configurational entropy, and no intramolecular H-bonds. The observed conformational preferences appear to be dictated by four main driving forces (ring conformational preferences, exo-anomeric effect, steric constraints, and possible presence of a third glycosidic dihedral angle), leaving a secondary role to intramolecular H-bonding and specific solvation effects. In spite of the weak conformational driving force attributed to solvent-exposed H-bonds in water (highly polar protic solvent), intramolecular H-bonds may still have a significant influence on the physico-chemical properties of the disaccharide by decreasing its hydrophilicity. Along with previous work, the results also complete the suggestion of a spectrum of approximate transition timescales for carbohydrates up to the disaccharide level, namely: approximately 30 ps (hydroxyl groups), approximately 1 ns (free lactol group, free hydroxymethyl groups, glycosidic dihedral angleomega in (1-->6)-linked disaccharides), approximately 10 ns to 2 micros (ring conformation, glycosidic dihedral angles phi and psi). The calculated average values of the glycosidic torsional angles agree well with the available experimental data, providing validation for the force-field and simulation methodology employed.

Small molecule beta-amyloid inhibitors that stabilize protofibrillar structures in vitro improve cognition and pathology in a mouse model of Alzheimer's disease

Beta-amyloid (Abeta) peptides are thought to play a major role in the pathogenesis of Alzheimer's disease. Compounds that disrupt the kinetic pathways of Abeta aggregation may be useful in elucidating the role of oligomeric, protofibrillar and fibrillar Abeta in the etiology of the disease. We have previously reported that scyllo-inositol inhibits Abeta(42) fibril formation but the mechanism(s) by which this occurs has not been investigated in detail. Using a series of scyllo-inositol derivatives in which one or two hydroxyl groups were replaced with hydrogen, chlorine or methoxy substituents, we examined the role of hydrogen bonding and hydrophobicity in the structure-function relationship of scyllo-inositol-Abeta binding. We report here that all scyllo-inositol derivatives demonstrated reduced effectiveness in preventing Abeta(42) fibrillization compared with scyllo-inositol, suggesting that scyllo-inositol interacts with Abeta(42) via key hydrogen bonds that are formed by all hydroxyl groups. Increasing the hydrophobicity of scyllo-inositol by the addition of two methoxy groups (1,4-di-O-methyl-scyllo-inositol) produced a derivative that stabilized Abeta(42) protofibrils in vitro. Prophylactic administration of 1,4-di-O-methyl-scyllo-inositol to TgCRND8 mice attenuated spatial memory impairments and significantly decreased cerebral amyloid pathology. These results suggest that Abeta aggregation can be targeted at multiple points along the kinetic pathway for the improvement of Alzheimer's disease-like pathology.

A 3D-structural model of unsulfated chondroitin from high-field NMR: 4-sulfation has little effect on backbone conformation

The glycosaminoglycan chondroitin sulfate is essential in human health and disease but exactly how sulfation dictates its 3D-structure at the atomic level is unclear. To address this, we have purified homogenous oligosaccharides of unsulfated chondroitin (with and without (15)N-enrichment) and analysed them by high-field NMR to make a comparison published chondroitin sulfate and hyaluronan 3D-structures. The result is the first full assignment of the tetrasaccharide and an experimental 3D-model of the hexasaccharide (PDB code 2KQO). In common with hyaluronan, we confirm that the amide proton is not involved in strong, persistent inter-residue hydrogen bonds. However, in contrast to hyaluronan, a hydrogen bond is not inferred between the hexosamine OH-4 and the glucuronic acid O5 atoms across the beta(1-->3) glycosidic linkage. The unsulfated chondroitin bond geometry differs slightly from hyaluronan by rotation about the beta(1-->3) psi dihedral (as previously predicted by simulation), while the beta(1-->4) linkage is unaffected. Furthermore, comparison shows that this glycosidic linkage geometry is similar in chondroitin-4-sulfate. We therefore hypothesise that both hexosamine OH-4 and OH-6 atoms are solvent exposed in chondroitin, explaining why it is amenable to sulfation and hyaluronan is not, and also that 4-sulfation has little effect on backbone conformation. Our conclusions exemplify the value of the 3D-model presented here and progress our understanding of glycosaminoglycan molecular properties.

Validating a strategy for molecular dynamics simulations of cyclodextrin inclusion complexes through single-crystal X-ray and NMR experimental data: a case study

A theoretical and experimental study about the formation and structure of the inclusion complex (-)-menthyl-O-beta-D-glucopyranoside 1 with beta-cyclodextrin (beta-CD) 2 is presented as paradigmatic case study to test the results of molecular dynamics (MD) simulations. The customary methodological approach-the use of experimental geometrical parameters as restraints for MD runs-is logically reversed and the calculated structures are a posteriori compared with those obtained from NMR spectroscopy in D(2)O solution and single crystal X-ray diffraction so as to validate the simulation procedure. The guest molecule 1 allows for a broad repertoire of intermolecular interactions (dipolar, hydrophobic, hydrogen bonds) concurring to stabilize the host-guest complex, thus providing the general applicability of the simulation procedure to cyclodextrin physical chemistry. Many starting geometries of the host-guest association were chosen, not assuming any a priori inclusion. The simulation protocol, involving energy minimization and MD runs in explicit water, yielded four possible inclusion geometries, ruling out higher-energy outer adducts. By analysis of the average energy at room temperature, the most stable geometry in solution was eventually obtained, while the kinetics of formation showed that it is also kinetically favored. The reliability of such geometry was thoroughly checked against the NOE distances via the pair distribution functions, that is, the statistical distribution of intermolecular distances among selected diagnostic atoms calculated from the MD trajectories at room temperature. An analogous procedure was adopted both with implicit solvent and in vacuo. The most stable geometry matched that found with explicit solvent but major differences were observed in the relative stability of the metastable complexes as a consequence of the lack of hydration on the polar moiety of the guest. Finally, a control set of geometrical parameters of the thermodynamically favored complex matched the corresponding one obtained from the X-ray structure, while local conformational differences were indicative of packing effects.