Neutron diffraction and simulation studies of the exocyclic hydroxymethyl conformation of glucose

Outstanding Properties of the Hydration Shell around β-d-Glucose: A Computational Study

Ab initio molecular dynamics (AIMD) simulations have been performed on aqueous solutions of four simple sugars, α-d-glucose, β-d-glucose, α-d-mannose, and α-d-galactose. Hydrogen-bonding (HB) properties, such as the number of donor- and acceptor-type HB-s, and the lengths and strengths of hydrogen bonds between sugar and water molecules, have been determined. Related electronic properties, such as the dipole moments of water molecules and partial charges of the sugar O atoms, have also been calculated. The hydrophilic and hydrophobic shells were characterized by means of spatial distribution functions. β-d-Glucose was found to form the highest number of hydrophilic and the smallest number of hydrophobic connections to neighboring water molecules. The average sugar-water H-bond length was the shortest for β-d-glucose, which suggests that these are the strongest such H-bonds. Furthermore, β-d-glucose appears to stand out in terms of the symmetry properties of both its hydrophilic and hydrophobic hydration shells. In summary, in all aspects considered here, there seems to be a correlation between the distinct characteristics of β-d-glucose reported here and its outstanding solubility in water. Admittedly, our findings represent only some of the important factors that influence the solubility.

Do non-thermal effects exist in microwave heating of glucose aqueous solutions? Evidence from molecular dynamics simulations

The existence of microwave non-thermal effects in food processing is still debated. In this study, molecular dynamics (MD) simulations were performed to investigate the conformation, electrostatic profiles, and intramolecular hydrogen bonds (intra-HB) of glucose in aqueous solution under alternating electric fields of microwaves ranging from 0 to 10⁹ V/m at 2.45 GHz. The results showed a field-induced threshold of 10⁹ V/m. At the threshold, alternating microwaves reoriented the flexible moieties and thus enhanced the intra-HB. The conformational transition among gg, gt, and tg conformers at 10⁹ V/m possibly resulted from the uneven electrostatic potential and the increased intra-HB. In practice, the maximum electric field of microwaves is several times weaker than the threshold, verifying the absence of microwave non-thermal effects for glucose molecules in food processing. This study provides a novel strategy to evaluate the potential non-thermal effects of microwaves in food processing and the related underlying food safety issues.

Structural Comparison between Sucrose and Trehalose in Aqueous Solution

The two sugar molecules sucrose and trehalose are both considered as stabilizing molecules for the purpose of preserving biological materials during, for example, lyophilization or cryo-preservation. Although these molecules share a similar molecular structure, there are several important differences in their properties when they interact with water, such as differences in solubility, viscosity, and glass transition temperature. In general, trehalose has been shown to be more efficient than other sugar molecules in preserving different biological molecules against stress, and thus by investigating how these two disaccharides differ in their water interaction, it is possible to further understand what makes trehalose special in its stabilizing properties. For this purpose, the structure of aqueous solutions of these disaccharides was studied by using neutron and X-ray diffraction in combination with empirical potential structure refinement (EPSR) modeling. The results show that there are surprisingly few differences in the overall structure of the solutions, although there are indications for that trehalose perturbs the water structure slightly more than sucrose.

Effects of varying the 6-position oxidation state of hexopyranoses: a systematic comparative computational analysis of 48 monosaccharide stereoisomers

Knowledge of multi-dimensional carbohydrate structure is essential when delineating structure-function relationships in the development of analytical techniques such as ion mobility-mass spectrometry and of carbohydrate-based therapeutics, as well as in rationally modifying the chemical and physical properties of drugs and materials based on sugars. Although monosaccharides are conventionally presumed to adopt the canonical ⁴C₁ chair conformation, it is not well known how altering the substituent identity around the pyranose ring affects the favored conformational state. This work provides a comprehensive and systematic computational comparison of all eight aldohexose isomers in the gas phase with reduction and oxidation at the C-6 position using density functional theory (M05-2X/cc-pVTZ(-f)//B3LYP/6-31G**) to determine the conformational and anomeric preference for each sugar in the gas phase. All 6-deoxyhexose and aldohexose isomers favored the ⁴C₁ chair conformation, while oxidation at C-6 showed a shift in equilibrium to favor the ¹C₄ chair for β-alluronic acid, β-guluronic acid, and β-iduronic acid. The anomeric preference was found to be significantly affected by a remote change in oxidation state, with the alternate anomer favored for several isomers. These findings provide a fundamental platform to empirically test steric and electronic effects of pyranose substituents, with the goal of formulating straightforward rules that govern carbohydrate reactivity and drive quicker, more efficient syntheses. Graphical abstract A systematic comparative conformational analysis of all eight aldohexose isomers using DFT methods (M05-2X/cc-pVTZ(-f)) reveals changes in anomeric and ring conformational preference upon reduction or oxidation at the C-6 position for several sugars.

Structure-activity relationships in carbohydrates revealed by their hydration

One of the more intriguing aspects of carbohydrate chemistry is that despite having very similar molecular structures, sugars have very different properties. For instance, there is a sensible difference in sweet taste between glucose and trehalose, even though trehalose is a disaccharide that comprised two glucose units, suggesting a different ability of these two carbohydrates to bind to sweet receptors. Here we have looked at the hydration of specific sites and at the three-dimensional configuration of water molecules around three carbohydrates (glucose, cellobiose, and trehalose), combining neutron diffraction data with computer modelling. Results indicate that identical chemical groups can have radically different hydration patterns depending on their location on a given molecule. These differences can be linked with the specific activity of glucose, cellobiose, and trehalose as a sweet substance, as building block of cellulose fiber, and as a bioprotective agent, respectively. This article is part of a Special Issue entitled "Recent Advances in Bionanomaterials" Guest Editors: Dr. Marie-Louise Saboungi and Dr. Samuel D. Bader.

Are all polar molecules hydrophilic? Hydration numbers of nitro compounds and nitriles in aqueous solution

Typical highly polar nitro compounds and nitriles should be classified ashydroneutralcompounds because of their hydration numbers of zero.

Sodium ion interactions with aqueous glucose: insights from quantum mechanics, molecular dynamics, and experiment

In the last several decades, significant efforts have been conducted to understand the fundamental reactivity of glucose derived from plant biomass in various chemical environments for conversion to renewable fuels and chemicals. For reactions of glucose in water, it is known that inorganic salts naturally present in biomass alter the product distribution in various deconstruction processes. However, the molecular-level interactions of alkali metal ions and glucose are unknown. These interactions are of physiological interest as well, for example, as they relate to cation-glucose cotransport. Here, we employ quantum mechanics (QM) to understand the interaction of a prevalent alkali metal, sodium, with glucose from a structural and thermodynamic perspective. The effect on β-glucose is subtle: a sodium ion perturbs bond lengths and atomic partial charges less than rotating a hydroxymethyl group. In contrast, the presence of a sodium ion significantly perturbs the partial charges of α-glucose anomeric and ring oxygens. Molecular dynamics (MD) simulations provide dynamic sampling in explicit water, and both the QM and the MD results show that sodium ions associate at many positions with respect to glucose with reasonably equivalent propensity. This promiscuous binding nature of Na(+) suggests that computational studies of glucose reactions in the presence of inorganic salts need to ensure thorough sampling of the cation positions, in addition to sampling glucose rotamers. The effect of NaCl on the relative populations of the anomers is experimentally quantified with light polarimetry. These results support the computational findings that Na(+) interacts similarly with α- and β-glucose.

The conformation of a ribose derivative in aqueous solution: a neutron-scattering and molecular dynamics study

The structure of aqueous solutions of methyl β-D-ribofuranoside was investigated by coupling molecular dynamics (MD) simulations and neutron scattering measurements with isotopic substitution. Using a sample of the sugar isotopically-labeled at a single unique position, neutron scattering structure factors and radial distribution functions can be compared with MD simulations constrained to different conformations to determine which conformer best fits the experimental results. Three different simulations were performed with the methyl ether group of the sugar unconstrained and constrained in each of its staggered orientations. The results of the unconstrained simulation showed that the methyl ester group occupied predominantly the 300° position, which is in agreement with the diffraction experimental results. This result suggests that the molecular mechanics force field used in the simulation adequately describes the conformation of the 1-methyl ether group in the methyl β-D-ribofuranoside.

Are all polar molecules hydrophilic? Hydration numbers of ketones and esters in aqueous solution

Hydration numbers of typical polar compounds like ketones and esters in aqueous solution were precisely determined using high-frequency dielectric relaxation techniques up to a frequency of 50 GHz at 25 °C. Because the hydration number is one of the most quantitative parameters to demonstrate how much are molecules hydrophilic, it is a critical parameter to determine the hydrophilicity of compounds. Hydration numbers of some ketones bearing carbonyl groups were determined to be ca. 0 irrespective of the species of molecules. Moreover, hydration numbers of some esters were also evaluated to be ca. 0 as well as the ketones. These findings suggested that there is no hydrogen bond formation between the ester group and water molecules, nor is there the hydrogen bond formation between the carbonyl group and water molecules. Consequently, esters and ketones bearing typical polar groups are not classified into hydrophilic compounds, but into "hydroneutral" compounds positioned between hydrophilic and hydrophobic ones. Molecular motions of the examined polar molecules in aqueous solution were well described with single Debye-type rotational relaxation modes without strong interaction between solute and water molecules, and also between solute molecules because of the obtained Kirkwood factor close to unity. This independent rotational mode for the polar compounds results from the hydroneutral characteristics caused by the relationship n(H) = 0.

Simulation and Neutron Diffraction Studies of Small Biomolecules in Water

Modern biophysics has benefited greatly from the use of X-ray and neutron diffraction from ordered single crystals of proteins and other macromolecules to give highly detailed pictures of these molecules in the solid state. However, the most biologically relevant environments for these molecules are liquid solutions, and their liquid state properties are sensitive to details of the liquid structuring. The best experimental method for studying such structuring is also neutron diffraction, but of course, the inherent disorder of the liquid state means that these experiments cannot hope to achieve the level of informational detail available from single crystal diffraction. Nonetheless, recent advances in neutron beam intensity, beam stability, and detector sensitivity mean that it should be possible, at least in principle, to use such measurements to extract information about structuring in much more complex systems than have previously been studied. We describe a series of neutron diffraction studies of isotopically labeled molecules in aqueous solution which, when combined with results from computer simulations, can be used to extract conformational information of the hydration of the molecules themselves, essentially opening up new avenues of investigation in structural biology.

Observation of pyridine aggregation in aqueous solution using neutron scattering experiments and MD simulations

Neutron diffraction with isotopic substitution (NDIS) experiments have been used to examine the structuring of aqueous solutions of pyridine. A new method is described for extracting the structure factors relating to intermolecular correlations from neutron scattering experiments on liquid solutions of complex molecular species. This approach performs experiments at different concentrations and exploits the intramolecular coordination number concentration invariance (ICNCI) to separate the internal and intermolecular contributions to the total intensities. The ability of this method to deconvolute molecular and intermolecular correlations is tested and demonstrated using simulated neutron scattering results predicted from molecular dynamics simulations of aqueous solutions of the polyatomic solute pyridine in which the inter- and intramolecular terms are known. The method is then implemented using neutron scattering measurements on solutions of pyridine. The results confirm that pyridine shows a significant propensity to aggregate in solution and demonstrate the prospects for the future application of the ICNCI approach to the study of large polyatomic solutes using next-generation neutron sources and detectors.

Pathways and populations: stereoelectronic insights into the exocyclic torsion of 5-(hydroxymethyl)tetrahydropyran

High level ab initio computations in vacuum and with the IEFPCM implicit solvent model are carried out on 5-(hydroxymethyl)tetrahydropyran to investigate the effects of water on the exocyclic torsional surface. Rotamer populations evaluated from the omega(C-C-C-O), theta(C-C-C-O) solvent surface agree almost quantitatively with experimental values for the closely related methyl 4-deoxy-alpha-D-xylohexopyranoside. Potentials of mean force obtained from the two surfaces show substantial solvent stabilization of the TG (omega = 180 +/- 60 degrees) rotamer and the barriers at omega= 120 and 240 degrees but solvent destabilization at the cis barrier (omega = 0 degrees). Natural bond orbital analyses indicate that energetics of these effects are largely explained by overstabilization of the vacuum GT (omega= 60 +/- 60 degrees) and GG (omega = 300 +/- 60 degrees) rotamers. Solvent stabilization of theta conformations provides entropic stabilization.