Protein-Ligand Dissociation Simulated by Parallel Cascade Selection Molecular Dynamics.

We investigated the dissociation process of tri-N-acetyl-d-glucosamine from hen egg white lysozyme using parallel cascade selection molecular dynamics (PaCS-MD), which comprises cycles of multiple unbiased MD simulations using a selection of MD snapshots as the initial structures for the next cycle. Dissociation was significantly accelerated by PaCS-MD, in which the probability of rare event occurrence toward dissociation was enhanced by the selection and rerandomization of the initial velocities. Although this complex was stable during 1 µs of conventional MD, PaCS-MD easily induced dissociation within 100-101 ns. We found that velocity rerandomization enhances the dissociation of triNAG from the bound state, whereas diffusion plays a more important role in the unbound state. We calculated the dissociation free energy by analyzing all PaCS-MD trajectories using the Markov state model (MSM), compared the results to those obtained by combinations of PaCS-MD and umbrella sampling (US), steered MD (SMD) and US, and SMD and the Jarzynski equality, and experimentally determined binding free energy. PaCS-MD/MSM yielded results most comparable to the experimentally determined binding free energy, independent of simulation parameter variations, and also gave the lowest standard errors.

Analytical characterization of human milk oligosaccharides - potential applications in pharmaceutical analysis.

Human breast milk is the gold standard for infant feeding and the best possible nourishment a new-born could have. Breastfeeding is the natural way to provide optimal nutritional, immunological and emotional nurturing for the healthy growth and development of infants. Human milk is a complex and dynamic biofluid comprised of many hundreds to thousands of distinct bioactive structures, among which one of the most abundant substances are the non-conjugated complex carbohydrates referred to as human milk oligosaccharides (HMOs). Due to their structural diversity and abundance, HMOs possess many beneficial biological functions. In order to understand human milk composition and HMO functions, state-of-the-art glycomic methods are inevitable. The industrial, large scale chemoenzymatic production of the most abundant HMOs became a reality in the last years and it evokes the need for straightforward and genuine analytical procedures to monitor the synthetic process and the quality of the products. It is obvious, that HMOs represent the next breakthrough in infant nutrition, as the addition of HMOs (such as 2'-fucosyllactose or lacto-N-neotetraose) to infant- and follow-on formulas, processed cereal-based food and baby foods for infants and young children etc. will revolutionize this field. This review highlights the potential applications of HMOs in the (bio)pharmaceutical industry, also summarizes the analytical methods available for the characterization of HMOs. An overview of the structure and function of HMOs along with their determination methods in complex matrices are provided. Various separation methods including liquid- and gas chromatography and capillary electrophoresis for the characterization and novel approaches for the quantitation of HMOs are discussed.

The carbohydrate-binding promiscuity of Euonymus europaeus lectin is predicted to involve a single binding site.

Euonymus europaeus lectin (EEL) is a carbohydrate-binding protein derived from the fruit of the European spindle tree. EEL was first identified for its erythrocyte agglutinating properties and specificity for B and H blood groups. However, a detailed molecular picture of the structural basis of carbohydrate recognition by EEL remains to be developed. In this study, we performed fluorescence titrations of a range of carbohydrates against EEL. Binding of EEL to a wide range of carbohydrates was observed, including a series of blood group-related carbohydrates, mannosides, chitotriose and sialic acid. Affinity was strongest for carbohydrates with H-related structures and the B trisaccharide. A homology model of EEL was produced from templates identified using the HHPred server, which employs hidden Markov models (HMMs) to identify templates. The HMM approach identified that the best templates for EEL were proteins featuring a ricin B-like (R-type) fold. Separate templates were used to model the core and binding site regions of the lectin. Through the use of constrained docking and spatial comparison with a template ligand, binding modes for the carbohydrate ligands were predicted. A relationship between the experimental binding energies and the computed binding energies of the selected docked poses was determined and optimized. Collectively, our results suggest that EEL utilizes a single site for recognition of carbohydrates terminating in a variety of monosaccharides.

Integrating ReSET with glycosyl iodide glycosylation in step-economy syntheses of tumor-associated carbohydrate antigens and immunogenic glycolipids.

Carbohydrates mediate a wide range of biological processes, and understanding these events and how they might be influenced is a complex undertaking that requires access to pure glycoconjugates. The isolation of sufficient quantities of carbohydrates and glycolipids from biological samples remains a significant challenge that has redirected efforts toward chemical synthesis. However, progress toward complex glycoconjugate total synthesis has been slowed by the need for multiple protection and deprotection steps owing to the large number of similarly reactive hydroxyls in carbohydrates. Two methodologies, regioselective silyl exchange technology (ReSET) and glycosyl iodide glycosylation have now been integrated to streamline the synthesis of the globo series trisaccharides (globotriaose and isoglobotriaose) and α-lactosylceramide (α-LacCer). These glycoconjugates include tumor-associated carbohydrate antigens (TACAs) and immunostimulatory glycolipids that hold promise as immunotherapeutics. Beyond the utility of the step-economy syntheses afforded by this synthetic platform, the studies also reveal a unique electronic interplay between acetate and silyl ether protecting groups. Incorporation of acetates proximal to silyl ethers attenuates their reactivity while reducing undesirable side reactions. This phenomenon can be used to fine-tune the reactivity of silylated/acetylated sugar building blocks.

Synthesis and assessment of globotriose-chitosan conjugate, a novel inhibitor of shiga toxins produced by Escherichia coli.

Shiga toxin (Stx)-producing Escherichia coli (STEC) causes diarrhea and colitis in humans that can develop into a life-threatening hemolytic uremic syndrome (HUS). Developing efficient means of controlling STEC diseases, for which no drugs or vaccines are currently available, remains a high priority. We report here the construction and development of chitosan conjugates bearing the Stx ligand trisaccharide globotriose to demonstrate their potential as STEC disease treatment agents. The synthesis was accomplished by grafting a globotriose derivative containing an aldehyde-functionalized aglycone to chitosan amino groups. The obtained globotriose-chitosan conjugate bound with high affinity to Stx and efficiently neutralized its toxicity on Vero cells. Moreover, Stx levels in the gut of infected mice receiving oral doses of the conjugate were greatly diminished, enabling the mice to resist a fatal STEC challenge. The conjugate appears to function as a Stx adsorbent in the gut, preventing toxin entry into the bloodstream and consequent development of HUS. As such, the conjugate could act as a novel agent against STEC disease.

Synthesis and characterization of a potential prebiotic trisaccharide from cheese whey permeate and sucrose by Leuconostoc mesenteroides dextransucrase.

The production of new bioactive oligosaccharides is currently garnering much attention for their potential use as functional ingredients. This work addresses the enzymatic synthesis and NMR structural characterization of 2-α-D-glucopyranosyl-lactose derived from sucrose:lactose and sucrose:cheese whey permeate mixtures by using a Leuconostoc mesenteroides B-512F dextransucrase. The effect of synthesis conditions, including concentration of substrates, molar ratio of donor/acceptor, enzyme concentration, reaction time, and temperature, on the formation of transfer products is evaluated. Results indicated that cheese whey permeate is a suitable material for the synthesis of 2-α-D-glucopyranosyl-lactose, giving rise to yields around 50% (in weight respect to the initial amount of lactose) under the optimum reaction conditions. According to its structure, this trisaccharide is an excellent candidate for a new prebiotic ingredient, due to the reported high resistance of α-(1→2) linkages to the digestive enzymes in humans and animals, as well as to its potential selective stimulation of beneficial bacteria in the large intestine mainly attributed to the two linked glucose units located at the reducing end that reflects the disaccharide kojibiose (2-α-D-glucopyranosyl-D-glucose). These findings could contribute to broadening the use of important agricultural raw materials, such as sucrose or cheese whey permeates, as renewable substrates for enzymatic synthesis of oligosaccharides of nutritional interest.

Conformational analysis and dynamics of mannobiosides and mannotriosides using Monte Carlo/stochastic dynamics simulations.

The conformation and dynamics of alpha-(1-->2)-mannobioside, alpha-(1-->6)-mannobioside, and of the trisaccharide alpha-Man-(1-->2)-alpha-Man-(1-->6)-alpha- Man-OMe were studied using Monte Carlo/stochastic dynamics (MC/SD) simulations, the AMBER* force field, and the GB/SA implicit water solvation model. The results are in agreement with available experimental data.

Conformational studies of Lewis X and Lewis A trisaccharides using NMR residual dipolar couplings.

The conformations of the histo-blood group carbohydrate antigens Lewis X (Le(x)) and Lewis A (Le(a)) were studied by NMR measurements of one-bond C-H residual dipolar couplings in partially oriented liquid crystal solutions. A strategy for rapid calculation of the difference between theoretical and experimental dipolar couplings of a large number of model structures generated by computer simulations was developed, resulting in an accurate model structure for the compounds. Monte Carlo simulations were used to generate models for the trisaccharides, and orientations of each model were sought that could reproduce the experimental residual dipolar coupling values. For both, Le(a) and Le(x), single low energy models giving excellent agreement with experiment were found, implying a compact rigidly folded conformation for both trisaccharides. The new approach was also applied to the pentasaccharides lacto-N-fucopentaose 2 (LNF-2) and lacto-N-fucopentaose 3 (LNF-3) proving its consistency and robustness. For describing the conformation of tightly folded oligosaccharides, a definition for characterization of ring planes in pyranoside chairs is proposed and applied to the analysis of the relation between the fucose and galactose residues in the epitopes, revealing the structural similarity between them.

A conformational study of the trisaccharide beta-D-Glcp-(1-->2)[beta-D-Glcp-(1-->3)]alpha-D-Glcp-OMe by NMR NOESY and TROESY experiments, computer simulations, and X-ray crystal structure analysis.

Proton-proton cross-relaxation rates have been measured for the trisaccharide beta-D-Glcp-(l --> 2)[beta-D-Glcp-(1 --> 3)]alpha-D-Glcp-OMe in D2O as well as in D2O/[D6]DMSO 7:3 solution at 30 degrees C by means of one-dimensional NMR pulsed field gradient 1H,1H NOESY and TROESY experiments. Interatomic distances for the trisaccharide in D2O were calculated from the cross-relaxation rates for two intraresidue and three interglycosidic proton pairs, using the isolated spin-pair approximation. In the solvent mixture one intraresidue and three interglycosidic distances were derived without the use of a specific molecular model. In this case the distances were calculated from the cross-relaxation rates in combination with "model-free" motional parameters previously derived from 13C relaxation measurements. The proton-proton distances for interglycosidic pairs were compared with those averaged from Metropolis Monte Carlo and Langevin Dynamics simulations with the HSEA, PARM22, and CHEAT95 force fields. The crystal structure of the trisaccharide was solved by analysis of X-ray data. Interresidue proton pairs from the crystal structure and those observed by NMR experiments were similar. However, the corresponding proton-proton distances generated by computer simulations were longer. For the (1 --> 2) linkage the glycosidic torsion angles of the crystal structure were found in a region of conformational space populated by all three force fields, whereas for the (1 --> 3) linkage they occupied a region of low population density, as seen from the simulations.

Molecular dynamics simulation and nuclear magnetic resonance studies of the terminal glucotriose unit found in the oligosaccharide of glycoprotein precursors.

The trisaccharide alpha-d-Glcp-(1 --> 2)-alpha-d-Glcp-(1 --> 3)-alpha-d-Glcp-OMe, a model for the terminal glucotriose in Glc(3)Man(9)GlcNAc(2) in glycoprotein precursors, has been investigated by computer simulations and NMR spectroscopy. Molecular dynamics simulations were performed for 1 ns in aqueous solution and 20 ns in vacuo using the CHARMM-based force fields PARM22 and CHEAT95. An additional Monte Carlo simulation with the HSEA force field was also carried out. Experimental NMR data in water solution was obtained from measurement of long-range (1)H,(13)C heteronuclear trans-glycosidic coupling constants, (3)J(H,C), using one-dimensional Hadamard spectroscopy. Calculation of the (3)J(H,C) values from the simulations showed a varying degree of agreement to experimental data. It could be shown from simulation that the φ torsion angles differed, which was corroborated by the NMR measurements. Analyses were done of radial distribution functions and of hydrogen bonds. It was suggested that intermolecular hydrogen bonds were present, but in contrast to simulation the results from NMR spectroscopy did not support any major contribution. Hence, their influence on the conformation of the trisaccharide is rather small. Comparison of (1)H NMR chemical shifts for the trisaccharide and the glucotriose in Glc(3)Man(8)GlcNAc revealed high similarity. However, the derived conformation of the model substance in this work differed at one glycosidic torsion angle compared to the glucotriose on a large oligosaccharide.

Conformational analysis of an alpha-galactosyl trisaccharide epitope involved in hyperacute rejection upon xenotransplantation.

alpha-Galactosyl epitopes are carbohydrate structures bearing an alpha-Gal-(1-->3)-Gal terminus (alpha-Gal epitopes). The interaction of these epitopes on the surface of animal cells with anti alpha-Gal antibodies in human serum is believed to be the main cause in antibody-mediated hyperacute rejection in xenotransplantation. In this paper, conformational analysis of an N-linked alpha-D-Galp-(1-->3)-beta-D-Galp-(1-->4)-beta-D-Glcp trisaccharide epitope was conducted in terms of each monosaccharide residue conformation, primary hydroxymethyl group configuration, and interglycosidic conformations. Selective 2D J-delta INEPT experiments have been carried out at three different temperatures to evaluate three-bond, long-range 13C-1H coupling constants for the crucial alpha-(1-->3) linkage. The NMR experimental data were complemented by theoretical calculations. The flexibility and dynamics of the trisaccharide have been studied by Metropolis Monte Carlo simulations. Ensemble-averaged three-bond, long-range 13C-1H coupling constants and nuclear Overhauser effects were in good agreement with the experimental data. The alpha-(1-->3) glycosidic linkage has shown a restricted flexibility as indicated by NMR spectroscopy and molecular modeling.

A conformational study of the vicinally branched trisaccharide beta-D-glcp-(1 --> 2)[beta-D-glcp-(1 --> 3)]alpha-D-Manp-OMe by nuclear Overhauser effect spectroscopy (NOESY) and transverse rotating-frame Overhauser effect spectroscopy (TROESY) experiments: comparison to Monte Carlo and Langevin dynamics simulations.

The trisaccharide beta-D-Glcp-(1 --> 2)[beta-D-Glcp-(1 --> 3)]alpha-D-Manp-OMe, a model for branching regions in oligosaccharides, has been investigated by one-dimensional DPFGSE (1)H, (1)H nuclear Overhauser effect spectroscopy (NOESY) and transverse rotating-frame Overhauser effect spectroscopy (TROESY) experiments at 30 degrees C in water and in the solvent mixture water : dimethyl sulfoxide (7 : 3). Cross-relaxation rates were obtained from the nmr experiments and interpreted as proton-proton distances. From Metropolis Monte Carlo and Langevin dynamics simulations, distances were calculated and compared to those obtained from experiment. Using the previously determined dynamics from carbon-13 nmr relaxation measurements of the trisaccharide in the solvent mixture, intraresidue proton distances could be obtained that were in excellent to reasonable agreement with those calculated from simulations. In water, the isolated spin-pair approximation was used for comparison of interproton distances. The experimentally derived distances in both solvents showed that the trans-glycosidic distances were shorter between the anomeric proton of the glucosyl group and the proton at the linkage position, respectively, than to the proton on the adjacent carbon on the mannosyl residue. The interresidue distances calculated from the computer simulations, performed with three different force fields, namely HSEA, PARM22, and CHEAT95, resulted in the reverse order in all cases but one.

MMC and LD simulations of alpha-D-Glcp-(1-->2)-alpha-D-Glcp-(1-->3)-alpha-D-Glcp-OMe. A model for the terminal trisaccharide in glycoprotein precursors.

The conformational flexibility and the dynamics of alpha-D-Glcp-(1-->2)-alpha-D-Glcp(1-->3)-alpha-D-Glcp-OMe (I) has been investigated by Metropolis-Monte Carlo with the HSEA (Hard Sphere Exo-Anomeric) force field and Langevin dynamics simulations employing two different CHARMm (Chemistry at HARvard Molecular Mechanics) force fields, CHEAT95 and PARM22. The conformational space spanned by the molecule is similar for the two former force fields but differ significantly for the latter. Hydrogen bonding between O2" and O4 of the title compound is analysed in comparison to NMR and preliminary results from X-ray powder diffraction studies.

Synthesis and structural studies of branched 2-linked trisaccharides related to blood group determinants.

A series of trisaccharide glycosides, Fuc-(1 reversible 2)-beta-Gal-(1 reversible 3)-beta-X-OMe (X = GlcNAc, Glc, 2-deoxy-Glc) related to the blood group determinant Le(d) have been synthesised both as their alpha- and beta-Fuc anomers together with the component disaccharide starting compounds. The conformational properties of the six trisaccharides together with their parent disaccharides have been investigated by NMR spectroscopy (proton and carbon chemical shifts and proton NOEs) in combination with computer modeling using the Monte Carlo approach and the HSEA force field using the GEGOP programme. The interaction between the terminal fucose unit and the reducing unit was probed by substitution of bulky NAc group with hydroxyl and deoxy substituents, respectively. Compounds with severe steric interactions were identified by the non-additivity of their carbon chemical shifts. This was subsequently confirmed by the detailed conformational assessment by NOE spectroscopy and computer modeling. The most severe contacts arose in the alpha-L-Fuc derivatives, whereas the beta-linked L-Fuc derivatives only in one case exhibit severe steric interaction as probed by the NMR parameters.

Assessing glycosidic linkage flexibility: conformational analysis of the repeating trisaccharide unit of Aeromonas salmonicida.

A detailed conformational analysis was performed for the synthetic branched trisaccharide beta-D-ManNAc-(1-->4)-[alpha-D-Glc-(1-->3)]-L-Rha 1 which represents the repeating unit of the O-antigenic polysaccharide of Aeromonas salmonicida. The study was based on 26 experimental NOE curves from 1D transient NOE experiments, employing Gaussian-shaped inversion pulses at 600 MHz. Eight of the NOE curves were interglycosidic and thus useful for an analysis of glycosidic linkage orientations. Metropolis Monte Carlo (MMC) simulations and minimum-energy calculations with the program GEGOP were used to obtain theoretical NOE curves which were compared to the experimental ones. MMC simulations with different temperature parameters of 310, 600, 900 and 2000 K allowed identification of NOEs which are sensitive towards different conformation distributions--not only different conformations--at both glycosidic linkages in 1. A comparison of trisaccharide 1 with the constituent disaccharides beta-D-ManNAc-(1-->4)-L-Rha 2 and alpha-D-Glc-(1-->3)-L-Rha 3 revealed effects of branching on glycosidic linkage flexibility. A quantitative evaluation was facilitated by the introduction of entropy-related flexibility parameters. Our study indicates a notable restriction of flexibility, especially at the (1-->3) linkage in 1. Although overall flexibility in 1 is reduced as compared to the constituent disaccharides 2 and 3, it cannot be neglected altogether. In summary, combined transient NOE experiments and MMC simulations provide a simple approach to analyse glycosidic linkage flexibility.