DFTMD studies of glucose and epimers: anomeric ratios, rotamer populations, and hydration energies

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.

Chitosan Polysaccharides from a Polarizable Multiscale Approach

We report simulations of chitosan polysaccharides in the aqueous phase, at infinite dilute conditions and zero ionic strength. Those simulations are performed by means of a polarizable multiscale modeling scheme that relies on a polarizable all atom force field to model solutes and on a polarizable solvent coarse grained approach. Force field parameters are assigned only from quantum chemistry ab initio data. We simulate chitosan monomer units, dimers and 50-long chains. Regarding the 50-long chains we simulate three sets of ten randomly built chain replica at three different pH conditions (corresponding to different chain protonation states, the chain degree of deacetylation is 85%). Our simulations show the persistence length of 50-long chitosan chains at strong acidic conditions (pH <5) to be 24 ± 2 nm (at weak/negligible ionic strength conditions), and to be 1 order of magnitude shorter at usual pH conditions. Our simulation data support the most recent simulation and experimental studies devoted to chitosan polysaccharides in the aqueous phase.

Between Two Chairs: Combination of Theory and Experiment for the Determination of the Conformational Dynamics of Xylosides

The conformational properties of monosaccharides constitute fundamental features of oligosaccharides. While the energy landscape of monosaccharides can be altered by a specific biochemical environment or by chemical modifications, the analysis of resulting dynamic conformational equilibria is not feasible by experimental means alone. In this work, a series of β-d-xylopyranosides is used to outline how a combination of experimental NMR parameters and computed molecular properties can be used to determine conformers and quantify the composition of conformational equilibria. We demonstrate that identifying the most stable conformers using energy calculations is challenging and computing of NMR shieldings is typically not sensitive enough. On the other hand, computed spin-spin coupling constants for the xyloside ring can be used to unambiguously assign experimental NMR data of dynamic conformational equilibria and quantify the ratio of different conformers in the mixture. As a proof of principle, this procedure allowed to analyze a hitherto unknown dynamic equilibrium of a diamino-xyloside as a precursor of a molecular switch.

Insight into the Mechanism of Reduced IgG/IgE Binding Capacity in Ovalbumin as Induced by Glycation with Monose Epimers through Liquid Chromatography and High-Resolution Mass Spectrometry

Ovalbumin (OVA) is one of the major food allergens in hen eggs. In this work, it was demonstrated that glycation with d-glucose and its epimers, including d-mannose, d-allose, d-galactose, and l-idose, could effectively attenuate the IgG/IgE binding of OVA, which was attributed to the covalent masking by sugars and to its structural changes. The glycation sites were determined, and their average degree of substitution was found using liquid chromatography coupled with high-resolution mass spectrometry. Fluctuations in OVA conformation were monitored by conventional spectrometry. Compared to those of OVA-Man and OVA-Glu, OVA-All, OVA-Gal, and OVA-Ido showed a higher glycation extent, and the alterations on their steric layouts were more drastic, suggesting that the configuration of hydroxyl groups at positions C-3, C-4, and C-5 in sugars might be important for the glycation reactivity; as such, their capabilities in binding with IgG/IgE decreased more significantly. Attempts were made to provide valuable information for in-depth understanding of the differences in biochemical functionality among epimeric sugars. These insights would be helpful for designing sweetened food products with a desirable level of safety.

Mutarotation of aldoses: Getting a deeper knowledge of a classic equilibrium enabled by computational analyses

The mutarotation equilibrium, by which reducing carbohydrates exist in solution as the α and β anomers of cyclic (furanoid and pyranoid) structures, along with open-chain (aldehyde and hydrate) forms, and whose ratios are depending on factors such as temperature, pH and solvent, portraits a phenomenon involved in numerous processes of chemical and biological importance. Herein, we have developed a DFT-based rationale that provides a broader landscape for anomerizations and ring-open chain interconversions, together with the pivotal role exerted not only by the aldehyde intermediate (essentially the only acyclic structure taken into account so far), but also the hydrate form (often more abundant at the equilibrium). These calculations reveal a more complex and richer scenario than was thought, and identify different mutarotation mechanisms that hinge on every monosaccharide. It is noteworthy that pyranose-furanose interconversion may actually occur without the intermediacy of open-chain forms. For the aldoses evaluated, namely d-glucose, d-ribose, and d-xylose, all structures involved in mutarotation undergo interconversion pathways, whose energy barriers calculated at the M06-2X/6-311++G(d,p) level, are in good agreement with previous experimental measurements.

Performance of Density-Functional Tight-Binding in Comparison to Ab Initio and First-Principles Methods for Isomer Geometries and Energies of Glucose Epimers in Vacuo and Solution

Density functional theory (DFT) is a widely used methodology for the computation of molecular and electronic structure, and we confirm that B3LYP and the high-level ab initio G3B3 method are in excellent agreement for the lowest-energy isomers of the 16 glucose epimers. Density-functional tight-binding (DFTB) is an approximate version of DFT with typically comparable accuracy that is 2 to 3 orders of magnitude faster, therefore generally very suitable for processing large numbers of complex structures. Conformational isomerism in sugars is well known to give rise to a large number of isomer structures. On the basis of a comprehensive study of glucose epimers in vacuo and aqueous solution, we found that the performance of DFTB is on par to B3LYP in terms of geometrical parameters excluding hydrogen bonds and isomer energies. However, DFTB underestimates both hydrogen bonding interactions as well as torsional barriers associated with rotations of the hydroxy groups, resulting in a counterintuitive overemphasis of hydrogen bonding in both gas phase as well as in water. Although the associated root mean squared deviation from B3LYP within epimer isomer groups is only on the order of 1 kcal/mol, this deviation affects the correct assignment of major isomer ordering, which span less than 10 kcal/mol. Both second- as well as third-order DFTB methods are exhibiting similar deviations from B3LYP. Even after the inclusion of empirical dispersion corrections in vacuum, these deviations remain for a large majority of isomer energies and geometries when compared to dispersion-corrected B3LYP.

Probing deoxysugar conformational preference: A comprehensive computational study investigating the effects of deoxygenation

Deoxysugars are intrinsic components in a number of antibiotics, antimicrobials, and therapeutic agents that often dictate receptor binding, improve efficacy, and provide a diverse toolbox in modifying glycoconjugate function due to an extensive number of unique isomers and inherent conformational flexibility. Hence, this work provides a comprehensive examination of the conformational effects associated with deoxygenation of the pyranose ring. Both the location and degree of deoxygenation were evaluated by interrogating the energetic landscape for a number of mono- and dideoxyhexopyranose derivatives using DFT methods (M05-2X/cc-pVTZ(-f)). Both anomeric forms and in some cases, the alternate chair form, have been investigated in the gas phase. As was documented in a preceding study, variation of the C-6 oxidation state has been shown to affect the anomeric preference of select glucose stereoisomers. Similar results were also observed for several deoxysugar isomers in this work, wherein the alternate anomer was favored upon reduction to the 6-deoxyhexose derivative or oxidation to the hexonic acid. Additionally, comparison of relative Gibbs free energies revealed C-3 deoxygenation imparts greater instability compared to C-2 or C-4 deoxygenation, as indicated by an increase in free energy for 3-deoxysugars. A polarizable continuum solvation model was also applied to empirically validate theoretical results for several deoxysugars, wherein good agreement with both carbon (σ = 1.6 ppm) and proton (σ = 0.20 ppm) NMR shifts was observed for the majority of isomers. Solvated and gas phase anomeric ratios were also calculated and compared favorably to reported literature values, although some discrepancies are noted.

Quantification of anomeric structural changes of glucose solutions using near-infrared spectra

Glucose is the most abundant carbohydrate found in living organisms. It exists as two anomers: α-D-glucose and β-D-glucose, which differ in how the hydroxyl group on the C1 carbon is directed. In solutions, the ratio between α- and β-D-glucose is typically 4:6 but can vary depending on the surrounding ions or temperature. In this study, we obtained near-infrared (NIR) spectra of the glucose anomers based on concentration, and analyzed the spectral difference between each anomer by spectra subtraction and principal component analysis, respectively. Moreover, by simultaneously measuring the optical rotation and NIR spectra from dissolution to equilibration, we showed that NIR spectra quantitatively estimated the specific rotations of glucose solutions using partial least-squares regression in the 1100-1800 nm wavelength range. All the analytical results indicated that the absorption at 1742 nm possess the potential to distinguish each glucose anomer quantitatively. Therefore, we addressed the prediction of the specific rotation by the absorption at 1742 nm, and demonstrated that the absorption normalized by line subtraction showed the high correlation with measured specific rotation. The absorption at 1742 nm reflects structural changes of the glucose anomers in solution. Our spectroscopy study not only provides spectral information about glucose anomers, which are the most fundamental chemical compounds in organisms, but also shows the possibility to detect the anomer ratio in vivo for the fields of agriculture and medicine by taking advantage of NIR.

Effects of temperature, concentration, and isomer on the hydration structure in monosaccharide solutions

Water-monosaccharide coupled interactions are essential for the function, stability, and dynamics of all glycans. Using molecular dynamics simulations, we investigated the effects of temperature, concentration, and monosaccharide isomer on the hydration structure and water dynamics in the hydration shell of monosaccharides in solution. We found that perturbations of the hydrogen-bond (H-bond) network in the first hydration shell around each monosaccharide molecule can be separated into two regions: one rich in water molecules with donor H-bonds (in the 2.4-2.8 Å region) and the other rich in water molecules with abundant acceptor H-bonds (in the 2.8-3.3 Å region). Moreover, we investigated the dependencies of clustering and conversion of the conformers of the monosaccharides on temperature and concentration. Increasing the concentration enhances monosaccharide clustering in all the monosaccharide solutions, while cluster formation does not depend on temperature. In the clusters, some water molecules in the hydration shell are replaced with monosaccharide oxygen atoms, which contributes to the shrinkage of the hydration shell with increasing monosaccharide concentration. The monosaccharides basically adopt one of two conformers, the stable chair or the unstable boat conformer. We revealed that the hydration structures of the boat and chair conformers were dramatically different. As the temperature increases, the content of the chair conformer decreases. Thus, the conversion of conformers strongly affects the hydration structure around the monosaccharide. These results are critical to understand the important roles of the hydration structure in glycan solutions.

Statistical Mechanics-Based Theoretical Investigation of Solvation Effects on Glucose Anomer Preferences

The importance of solvation effects on the stability of glucose anomers has been studied by the combination of quantum mechanics and statistical mechanics, namely, the reference interaction site model self-consistent field spatial electron density distribution. The preferences of α- and β-glucose in H₂O are well reproduced with the obtained ratio of 35:65 for α- and β-glucose, respectively. Indirect interactions and bulk effects, described by the Onsager model, are relatively small compared to the direct solute-solvent interactions, especially in [DMIM]Cl and dimethyl sulfoxide. From the decomposition of solvation free energy and solvation structures, it can be seen that the interactions with the solvent molecules greatly contribute to the anomer preferences.

Molecular dynamics simulations of hexopyranose ring distortion in different force fields

The four classical, biomolecular force fields designed to study hexopyranose-based carbohydrates (GROMOS 56a6CARBO/56a6CARBO_R, GROMOS 53a6GLYC, CHARMM and GLYCAM06) have been tested in the context of ring-inversion properties. These properties were evaluated for both unfunctionalized monomers of all hexopyranoses of the d series and for residues in a chain composed of uniform units connected by α(1→4) and β(1→4) glycosidic linkages. The results indicate that the tested force fields differ in their predictions of the ring-inversion properties of both monomers and residues in a chain. The comparison with the available experimental data and with the semi-empirical Angyal scheme reveals that, at the level of monomers, GROMOS 56a6CARBO, GROMOS 53a6GLYC and CHARMM correctly reproduce the ring-inversion free energies. However, due to the lack of analogous reference data we cannot state which force field is more or less accurate in the context of ring distortion of residues in a chain. Therefore, the use of ab initio potentials is recommended in the prospective, quantitative studies on the related subject.

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.

Potential Energy Surface-Based Automatic Deduction of Conformational Transition Networks and Its Application on Quantum Mechanical Landscapes of d-Glucose Conformers

This paper describes our approach that is built upon the potential energy surface (PES)-based conformational analysis. This approach automatically deduces a conformational transition network, called a conformational reaction route map (r-map), by using the Scaled Hypersphere Search of the Anharmonic Downward Distortion Following method (SHS-ADDF). The PES-based conformational search has been achieved by using large ADDF, which makes it possible to trace only low transition state (TS) barriers while restraining bond lengths and structures with high free energy. It automatically performs sampling the minima and TS structures by simply taking into account the mathematical feature of PES without requiring any a priori specification of variable internal coordinates. An obtained r-map is composed of equilibrium (EQ) conformers connected by reaction routes via TS conformers, where all of the reaction routes are already confirmed during the process of the deduction using the intrinsic reaction coordinate (IRC) method. The postcalculation analysis of the deduced r-map is interactively carried out using the RMapViewer software we have developed. This paper presents computational details of the PES-based conformational analysis and its application to d-glucose. The calculations have been performed for an isolated glucose molecule in the gas phase at the RHF/6-31G level. The obtained conformational r-map for α-d-glucose is composed of 201 EQ and 435 TS conformers and that for β-d-glucose is composed of 202 EQ and 371 TS conformers. For the postcalculation analysis of the conformational r-maps by using the RMapViewer software program we have found multiple minimum energy paths (MEPs) between global minima of ¹C₄ and ⁴C₁ chair conformations. The analysis using RMapViewer allows us to confirm the thermodynamic and kinetic predominance of ⁴C₁ conformations; that is, the potential energy of the global minimum of ⁴C₁ is lower than that of ¹C₄ (thermodynamic predominance) and that the highest energy of those of all the TS structures along a route from ⁴C₁ to ¹C₄ is lower than that of ¹C₄ to ⁴C₁ (kinetic predominance).

Revision of the GROMOS 56A6(CARBO) force field: Improving the description of ring-conformational equilibria in hexopyranose-based carbohydrates chains

This article describes a revised version 56A6(CARBO_R) of the GROMOS 56A6(CARBO) force field for hexopyranose-based carbohydrates. The simulated properties of unfunctionalized hexopyranoses are unaltered with respect to 56A6CARBO . In the context of both O1 -alkylated hexopyranoses and oligosaccharides, the revision stabilizes the regular (4) C1 chair for α-anomers, with the opposite effect for β-anomers. As a result, spurious ring inversions observed in α(1→4)-linked chains when using the original 56A6(CARBO) force field are alleviated. The (4) C1 chair is now the most stable conformation for all d-hexopyranose residues, irrespective of the linkage type and anomery, and of the position of the residue along the chain. The methylation of a d-hexopyranose leads to a systematic shift in the ring-inversion free energy ((4) C1 to (1) C4 ) by 7-8 kJ mol(-1), positive for the α-anomers and negative for the β-anomers, which is qualitatively compatible with the expected enhancement of the anomeric effect upon methylation at O1. The ring-inversion free energies for residues within chains are typically smaller in magnitude compared to those of the monomers, and correlate rather poorly with the latter. This suggests that the crowding of ring substituents upon chain formation alters the ring flexibility in a nonsystematic fashion. In general, the description of carbohydrate chains afforded by 56A6(CARBO_R) suggests a significant extent of ring flexibility, i.e., small but often non-negligible equilibrium populations of inverted chairs, and challenges the "textbook" picture of conformationally locked carbohydrate rings.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

Polarizable empirical force field for hexopyranose monosaccharides based on the classical Drude oscillator

A polarizable empirical force field based on the classical Drude oscillator is presented for the hexopyranose form of selected monosaccharides. Parameter optimization targeted quantum mechanical (QM) dipole moments, solute-water interaction energies, vibrational frequencies, and conformational energies. Validation of the model was based on experimental data on crystals, densities of aqueous-sugar solutions, diffusion constants of glucose, and rotational preferences of the exocylic hydroxymethyl of d-glucose and d-galactose in aqueous solution as well as additional QM data. Notably, the final model involves a single electrostatic model for all sixteen diastereomers of the monosaccharides, indicating the transferability of the polarizable model. The presented parameters are anticipated to lay the foundation for a comprehensive polarizable force field for saccharides that will be compatible with the polarizable Drude parameters for lipids and proteins, allowing for simulations of glycolipids and glycoproteins.

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.

Vacuum-ultraviolet electronic circular dichroism study of methyl α-D-glucopyranoside in aqueous solution by time-dependent density functional theory

The vacuum-ultraviolet (VUV) electronic circular dichroism (ECD) spectrum of methyl α-D-glucopyranoside (methyl α-D-Glc) was measured down to 163 nm in aqueous solution using a synchrotron-radiation VUV-ECD spectrophotometer. The spectrum exhibited two characteristic ECD peaks around 170 nm, which depend on the trans (T) and gauche (G) configurations of the hydroxymethyl group at C-5. To elucidate the influences of the T and G configurations on the spectrum, the ECD spectra of three rotamers (α-GT, α-GG, and α-TG) of methyl α-D-Glc were calculated using time-dependent density functional theory (TDDFT) combined with molecular dynamics simulation. A linear combination of the ECD spectra of these three rotamers, which differ markedly from each other, produced a methyl α-D-Glc spectrum similar to that observed experimentally. The spectrum was assignable to the n-σ* transitions of the ring oxygen and methoxy oxygen with minor contributions from the hydroxyl oxygen. The differences in α-GT, α-GG, and α-TG spectra were attributed to fluctuations of the configurations of the hydroxymethyl group at C-5 and the hydroxyl group at C-4, which strongly affected the orientations of intramolecular hydrogen bonds around the ring oxygen. These findings demonstrate that combining VUV-ECD and TDDFT is useful for structural characterization of saccharides in aqueous solution.

Dependence of pyranose ring puckering on anomeric configuration: methyl idopyranosides

In the aldohexopyranose idose, the unique presence of three axial ring hydroxyl groups causes considerable conformational flexibility, rendering it challenging to study experimentally and an excellent model for rationalizing the relationship between puckering and anomeric configuration. Puckering in methyl α- and β-L-idopyranosides was predicted from kinetically rigorous 10 μs simulations using GLYCAM11 and three explicit water models (TIP3P, TIP4P, and TIP4P-EW). In each case, computed pyranose ring three-bond (vicinal) (1)H-(1)H spin couplings ((3)J(H,H)) trended with NMR measurements. These values, calculated puckering exchange rates and free energies, were independent of the water model. The α- and β-anomers were (1)C(4) chairs for 85 and >99% of their respective trajectories and underwent (1)C(4)→(4)C(1) exchange at rates of 20 μs(-1) and 1 μs(-1). Computed α-anomer (1)C(4)↔(4)C(1) puckering rates depended on the exocyclic C6 substituent, comparing hydroxymethyl with carboxyl from previous work. The slower kinetics and restricted pseudorotational profile of the β-anomer were caused by water occupying a cavity bounded by the anomeric 1-O-methyl and the C6 hydroxymethyl groups. This finding rationalizes the different methyl α- and β-L-idopyranoside (3)J(H,H) values. Identifying a relationship between idopyranose anomeric configuration, microsecond puckering, and water structure facilitates engineering of biologically and commercially important derivatives and underpins deciphering presently elusive structure-function relationships in the glycome.

Specific rotation as a property to validate monosaccharide conformations

Specific rotation ([α](D)) values were calculated for the 15 conformations of xylopyranose that are the most stable in the gas phase and in aqueous solution. The effects of different theoretical descriptions and the medium on the geometry of the conformers and the [α](D) values are evaluated. Differences in [α](D) values found for the same conformer in all descriptions used were smaller than those found between any two different conformers in the same description. The difference between [α](D) values is prominent, even for two conformations that are distinguished from each other only by the orientation of one secondary hydroxyl group. This finding suggests that [α](D) values may potentially be used in the validation of monosaccharide conformations that are theoretically sampled.

Surface-confined conformers and coassembly-induced conformer resolution

Stereoisomerism is a fundamental chemistry issue and has been intensively investigated because of its importance in organic chemistry, biology, and pharmacology. Molecules with freely rotatable single bonds have many interconvertable conformers. Herein, we report the surface-adsorption-induced conformer resolution by employing azobenzene-3,3-dicarboxylic acid (ADA-33) as a model compound. Two linear assembly phases composed of trans conformers on a highly oriented pyrolytic graphite (HOPG) surface are observed by scanning tunneling microscopy. With the codeposition with 1-octanoic acid (OA), only one trans conformer of ADA-33 can be recognized by OA to form a two-component assembly with alternately arranged ADA-33 and OA stripes, which can be attributed to the epitaxial assembly of ADA-33 and OA on the HOPG surface, and weak hydrogen bonding exists between conformer I and OA molecules. The results are of significance with respect to the discrimination and resolution of conformers on a solid surface and provide molecular insights into the coadsorption assembly on the surface.

DFTMD studies of β-cellobiose: conformational preference using implicit solvent

Previous DFT in vacuo studies on the conformational preferences for cellobiose showed that upon optimization the φ(H)-anti conformations were of lower energy than the syn forms. Upon optimization using an implicit solvation method, COSMO, the syn or observed form was still not predicted to be of lower energy than the φ(H)-anti form, even though optimization after addition of several explicit water molecules did show a relative energy difference favoring the syn form. In order to examine the predictive ability of COSMO on this carbohydrate, constant energy dynamics, DFTMD, simulations were carried out on low energy syn and φ(H)-anti conformations with and without COSMO included during the dynamics. The resulting analysis confirmed that when COSMO is included in the dynamics, the syn conformations become energetically favored over the φ(H)-anti forms suggesting that both solvent and entropy play roles in dictating the solution conformation of cellobiose. Analysis of the dynamic runs includes distributions of selected dihedral angles versus time, conformational transitions, and populations of some quasi-planar, boat, skew forms during the simulations.