Internal motions and hydration of sucrose in a diluted water solution

Investigation of the Hydration Behavior of Different Sugars by Time Domain-NMR

The hydration behavior of sugars varies from each other and examining the underlying mechanism is challenging. In this study, the hydration behavior of glucose, fructose, allulose (aka rare sugar), and sucrose have been explored using different Time Domain Nuclear Magnetic Resonance (TD-NMR) approaches (relaxation times, self-diffusion, and Magic Sandwich Echo (MSE)). For that purpose, the effects of different sugar concentrations (2.5%, 5%, 10%, 15%, 20%, 30%, and 40%) (w/v) and hydration at different times for 1 day were investigated by T₂ relaxation times and self-diffusion coefficients. Crystallinity values of the solid and hydrated sugars were also determined with MSE. Change in T₂ relaxation times with concentration showed that the fastest binding with water (parallel with the shortest T₂ values) was observed for sucrose for all concentrations followed by glucose, fructose, and allulose. Furthermore, dependency of T₂ relaxation times with hydration time showed that sucrose was the fastest in binding with water followed by glucose, fructose, and allulose. The study showed that allulose, one of the most famous rare sugars that is known to be a natural low-calorie sugar alternative, had the lowest interaction with water than the other sugars. TD-NMR was suggested as a practical, quick, and accurate technique to determine the hydration behavior of sugars.

Translational Dynamics of Lipidated Ras Proteins in the Presence of Crowding Agents and Compatible Osmolytes

Ras proteins are small GTPases and are involved in transmitting signals that control cell growth, differentiation, and proliferation. Since the cell cytoplasm is crowded with different macromolecules, understanding the translational dynamics of Ras proteins in crowded environments is crucial to yielding deeper insight into their reactivity and function. Herein, the translational dynamics of lipidated N-Ras and K-Ras4B is studied in the bulk and in the presence of a macromolecular crowder (Ficoll) and the compatible osmolyte and microcrowder sucrose by fluorescence correlation spectroscopy. The results reveal that N-Ras forms dimers due to the presence of its lipid moiety in the hypervariable region, whereas K-Ras4B remains in its monomeric form in the bulk. Addition of a macromolecular crowding agent gradually favors clustering of the Ras proteins. In 20 wt % Ficoll N-Ras forms trimers and K-Ras4B dimers. Concentrations of sucrose up to 10 wt % foster formation of N-Ras trimers and K-Ras dimers as well. The results can be rationalized in terms of the excluded-volume effect, which enhances the association of the proteins, and, for the higher concentrations, by limited-hydration conditions. The results of this study shed new light on the association state of these proteins in a crowded environment. This is of particular interest for the Ras proteins, because their solution state-monomeric or clustered-influences their membrane-partitioning behavior and their interplay with cytosolic interaction partners.

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.

Evaluation of the regioselective delactonization of tri-sialic acid lactone by in-solution molecular dynamics simulation

An approximate model for the delactonization of tri-sialic acid lactone is presented with two water-layers that led to neutral hydrolysis of δ-lactone. The hydrolytic reactivity was studied with a 10-ns in-solution molecular dynamics simulation. The initial step of this hydrolysis involves a reactant water nucleophile complex via a proton transfer with another water molecule. Therefore, the probability of water molecules localized at the hydrolytic center correlates to the hydrolysis of δ-lactone. The stepwise delactonization of α2,8-(NeuAc)(3) lactone results/resulted from water concentration discrepancy near the carbonyl carbon of lactones in two water oxygen···carbonyl carbon shells, and the distances of OC···O(water) layers were 2.8 Å and 5.1 Å. Based on in-solution molecular dynamics study, the motion of water molecules over the re-face of the carbonyl groups was used for the quantitative description of the residence probability, p, whose value is 0.11 for lactone I and 0.33 for lactone II. The geometric criteria used to determine the residence statistics are (1) the distance of water-oxygen···carbonyl carbon in less than 5.1 Å and (2) the cone angle, θ, of carbonyl OC···O(water) in the range of 85-115°. As expected, a higher residence probability at lactone II led to its faster hydrolysis. Both the radial g(r) and angular p(θ) pair distribution functions of water oxygen and carbonyl groups of lactones ensure a better surrounding hydration encounter for lactone II. In contrast, water molecules around lactone I are deduced due to a steric hindrance by the turn structure of α2,8-(NeuAc)(3) lactone.

Use of sugar on the healing of diabetic ulcers: a review

With the advent of several innovative wound care management tools, the choice of products and treatment modalities available to clinicians continues to expand. High costs associated with wound care, especially diabetic foot wounds, make it important for clinician scientists to research alternative therapies and optimally incorporate them into wound care protocols appropriately. This article reviews using sugar as a treatment option in diabetic foot care and provides a guide to its appropriate use in healing foot ulcers. In addition to a clinical case study, the physiological significance and advantages of sugar are discussed.

Solvation of Glucose, Trehalose, and Sucrose by the Soft Sticky Dipole-Quadrupole-Octupole Water Model

Water structure around sugars modeled by partial charges is compared for soft-sticky dipole-quadrupole-octupole (SSDQO), a fast single-site multipole model, and commonly used multi-site models in Monte Carlo simulations. Radial distribution functions and coordination numbers of all the models indicate similar hydration by hydrogen-bond donor and acceptor waters. However, the new optimized SSDQO1 parameters as well as TIP4P-Ew and TIP5P predict a "lone-pair" orientation for the water accepting the sugar hydroxyl hydrogen bond that is more consistent with the limited experimental data than the "dipole" orientation in SPC/E, which has important implications for studies of the cryoprotectant properties of sugars.

The glass transition temperatures of amorphous trehalose–water mixtures and the mobility of water: an experimental and in silico study

Isothermal-isobaric molecular dynamics simulations are used to calculate the specific volume of models of trehalose and three amorphous trehalose-water mixtures (2.9%, 4.5% and 5.3% (w/w) water, respectively) as a function of temperature. Plots of specific volume versus temperature exhibit a characteristic change in slope when the amorphous systems change from the glassy to the rubbery state and the intersection of the two regression lines provides an estimate of the glass transition temperature T(g). A comparison of the calculated and experimental T(g) values, as obtained from differential scanning calorimetry, shows that despite the predicted values being systematically higher (about 21-26K), the trend and the incremental differences between the T(g) values have been computed correctly: T(g)(5.3%(w/w))<T(g)(4.5%(w/w))<T(g)(2.9%(w/w))<T(g)(0.0%(w/w)). The mobility of water has been investigated over temperature ranges covering the rubbery and the glassy phases of the trehalose-water mixtures by calculating the diffusion coefficients of water. The temperature dependence of the diffusion coefficient changes in the region of the glass transition and can be used as well to estimate T(g) values. The activation energies for water diffusion were found to be independent of the amount of water in amorphous trehalose.

Methods for peak assignment in low-resolution multidimensional NMR cross-correlation relaxometry

Several NMR protocols are presented for assigning peaks in complex T1-T2 spectra, including the effects of varying the spectrometer frequency and the CPMG pulsing rate. Extensions into a third dimension based on chemical-shift; diffusion- and field-cycled weighted T1-T2 cross-correlation methods are also explored as a means of peak assignment. We illustrate the power of these novel techniques with reference to simple aqueous sucrose solutions, but the methodology should be generally applicable.

Water T2 relaxation in sugar solutions

1H spin-spin relaxation times of water were measured with the CPMG sequence in dilute aqueous solutions of glucitol, mannitol, glycerol, glycol, the methyl D-pyranosides of alpha-glucose, beta-glucose, alpha-galactose, beta-galactose, alpha-xylose, beta-xylose, beta-arabinose and sucrose, alpha,alpha-trehalose, beta-maltose, maltotriose and maltoheptaose. The relaxation-time dispersion was measured by varying the CPMG pulse spacing, tau. These data were interpreted by means of the Carver-Richards model in which exchange between water protons and labile solute hydroxyl protons provides a significant contribution to the relaxation. From the dependences on temperature and tau, parameters characteristic of the pool of hydroxyls belonging to a given solute were extracted by nonlinear regression, including: the fraction of exchangeable protons, P, the chemical-shift difference between water protons and hydroxyl protons, deltaomega, the intrinsic spin-spin relaxation time, T2, and the chemical exchange rate, k. These solute-specific parameters are related, respectively, to the concentration, identity, mobility and exchange life-time of the hydroxyl site. At 298 K, values of deltaomega, T2 and k were found to be of the order of 1 ppm, 100 ms and 1000 s(-1), respectively. Effects of molecular size, conformation and solute concentration were investigated. The exchange mechanism was characterised by Eyring activation enthalpies and entropies with values in the ranges 50-70 kJ mol(-1) and -10 to 60 J K(-1)mol(-1), respectively.

A hydration study of (1-->4) and (1-->6) linked alpha-glucans by comparative 10 ns molecular dynamics simulations and 500-MHz NMR

The hydration behavior of two model disaccharides, methyl-alpha-D-maltoside (1) and methyl-alpha-D-isomaltoside (2), has been investigated by a comparative 10 ns molecular dynamics study. The detailed hydration of the two disaccharides was described using three force fields especially developed for modeling of carbohydrates in explicit solvent. To validate the theoretical results the two compounds were synthesized and subjected to 500 MHz NMR spectroscopy, including pulsed field gradient diffusion measurements (1: 4.0. 10(-6) cm(2). s(-1); 2: 4.2. 10(-6) cm(2). s(-1)). In short, the older CHARMM-based force field exhibited a more structured carbohydrate-water interaction leading to better agreement with the diffusional properties of the two compounds, whereas especially the alpha-(1-->6) linkage and the primary hydroxyl groups were inaccurately modeled. In contrast, the new generation of the CHARMM-based force field (CSFF) and the most recent version of the AMBER-based force field (GLYCAM-2000a) exhibited less structured carbohydrate-water interactions with the result that the diffusional properties of the two disaccharides were underestimated, whereas the simulations of the alpha-(1-->6) linkage and the primary hydroxyl groups were significantly improved and in excellent agreement with homo- and heteronuclear coupling constants. The difference between the two classes of force field (more structured and less structured carbohydrate-water interaction) was underlined by calculation of the isotropic hydration as calculated by radial pair distributions. At one extreme, the radial O em leader O pair distribution function yielded a peak density of 2.3 times the bulk density in the first hydration shell when using the older CHARMM force field, whereas the maximum density observed in the GLYCAM force field was calculated to be 1.0, at the other extreme.

Specific empirical free energy function for automated docking of carbohydrates to proteins

We present an automated docking protocol specifically optimized to predict the structure and affinity of a protein-carbohydrate complex. A scoring function was developed based on a training set of 30 protein-carbohydrate complexes of known structure and affinity. Combinations of several models for hydrogen bonding, torsional entropy loss, and solvation were tested for their ability to fit the training set data, and the best model was used with AutoDock. The electrostatic empirical coefficient is larger than in a previously obtained model using a training set comprised of various types of protein-ligand complexes, indicating that electrostatic interactions play a more important role in determining the affinity between a carbohydrate and a protein. The differences in the relative weighting of the empirical coefficients in the model yields predicted free energies for the training set with a standard error of 1.403 kcal/mol. The new scoring function was tested on 17 Aspergillus niger glucoamylase inhibitors for which binding energies had been determined experimentally. Free energies of complex formation were predicted with a residual standard error of 1.101 kcal/mol. The new scoring function therefore provides a robust method for predicting free energies of formation and optimal conformations of carbohydrate-protein complexes.

Crystal structure of Pterocarpus angolensis lectin in complex with glucose, sucrose, and turanose

The crystal structure of the Man/Glc-specific seed lectin from Pterocarpus angolensis was determined in complex with methyl-alpha-d-glucose, sucrose, and turanose. The carbohydrate binding site contains a classic Man/Glc type specificity loop. Its metal binding loop on the other hand is of the long type, different from what is observed in other Man/Glc-specific legume lectins. Glucose binding in the primary binding site is reminiscent of the glucose complexes of concanavalin A and lentil lectin. Sucrose is found to be bound in a conformation similar as seen in the binding site of lentil lectin. A direct hydrogen bond between Ser-137(OG) to Fru(O2) in Pterocarpus angolensis lectin replaces a water-mediated interaction in the equivalent complex of lentil lectin. In the turanose complex, the binding site of the first molecule in the asymmetric unit contains the alphaGlc1-3betaFruf form of furanose while the second molecule contains the alphaGlc1-3betaFrup form in its binding site.

Sugar for wounds

Sugar in its pure form, or incorporated into a paste containing an adhesive hydropolymer (gum), is a non-toxic treatment for a variety of wounds. Not only does it provide a suitable clean environment for angiogenesis to take place, but it will debride the wound surface and reduce odour. The presence of an adhesive hydropolymer seems to prevent hypergranulation, scarring and contraction.

Molecular dynamics study on lipid A from Escherichia coli: insights into its mechanism of biological action

Structural properties of the Escherichia coli lipid A moiety were analysed by means of molecular mechanics and molecular dynamics simulations and compared to synthetic monophospho and dephospho analogues with different biological activities in the Limulus assay. The conformation of glucosamine disaccharide headgroup, order and packing of fatty acid chains, solvation of phosphate groups, coordination by water molecules, sodium counterions and models of cationic amino acid side chains were described in terms of mean values, mean residence times, radial distribution functions, coordination numbers, solvation and interaction energies. Solvation and polar interactions of the phosphate groups were correlated to known biological activities the lipid A variants. The observed relationship between the biological effect and the number and position of the phosphate groups were explained with the help of simple mechanistic models of lipid A action. The possible mechanism of action involving specific binding of lipid A disaccharide headgroup to cationic residues of a receptor model was compared with an alternative mechanism, which assumes a relationship between the ability to adopt non-lamellar supramolecular structures and the biological activity. Conclusions are drawn about the probable mode of lipid A action. Implications for rational drug design of endotoxin-neutralising agents are discussed.

The mean hydration of carbohydrates as studied by normalized two-dimensional radial pair distributions

The hydration of carbohydrates plays a key role in many biological processes. Molecular dynamics simulations provide an effective tool for investigating the hydration of complex solutes such as carbohydrates. In this article we devise an algorithm for the calculation of two-dimensional radial pair distributions describing the probability of finding a water molecule in a site defined by two reference atoms. The normalized 2D radial pair distribution is proposed as an effective tool for investigating and comparing localized or ordered water sites around flexible molecules such as carbohydrates when analyzing molecular dynamics simulations and the utility of 2D radial pair distributions is demonstrated using sucrose as an example. In this relatively simple structure, 2D radial pair distributions were able to characterize and quantify the importance of two unique interresidue hydration sites in which a water molecule is forming a bridge between the glycopyranosyl and fructofuranosyl residues. The approach is proposed to be a valuable tool for comparing and understanding the hydration of flexible biomolecules such as carbohydrates.