Solvent-Modulated Influence of Intramolecular Hydrogen-Bonding on the Conformational Properties of the Hydroxymethyl Group in Glucose and Galactose: A Molecular Dynamics Simulation Study

Conformational Preferences at the Glycosidic Linkage of Saccharides in Solution as Deduced from NMR Experiments and MD Simulations: Comparison to Crystal Structures

Glycans are central to information content and regulation in biological systems. These carbohydrate molecules are active either as oligo- or polysaccharides, often in the form of glycoconjugates. The monosaccharide entities are joined by glycosidic linkages and stereochemical arrangements are of utmost importance in determining conformation and flexibility of saccharides. The conformational preferences and population distributions at the glycosidic torsion angles φ and ψ have been investigated for O-methyl glycosides of three disaccharides where the substitution takes place at a secondary alcohol, viz., in α-l-Fucp-(1→3)-β-d-Glcp-OMe, α-l-Fucp-(1→3)-α-d-Galp-OMe and α-d-Glcp-(1→4)-α-d-Galp-OMe, corresponding to disaccharide structural elements present in bacterial polysaccharides. Stereochemical differences at or adjacent to the glycosidic linkage were explored by solution state NMR spectroscopy using one-dimensional 1 H,1 H-NOESY NMR experiments to obtain transglycosidic proton-proton distances and one- and two-dimensional heteronuclear NMR experiments to obtain 3 JCH transglycosidic coupling constants related to torsion angles φ and ψ. Computed effective proton-proton distances from molecular dynamics (MD) simulations showed excellent agreement to experimentally derived distances for the α-(1→3)-linked disaccharides and revealed that for the bimodal distribution at the ψ torsion angle for the α-(1→4)-linked disaccharide experiment and simulation were at variance with each other, calling for further force field developments. The MD simulations disclosed a highly intricate inter-residue hydrogen bonding pattern for the α-(1→4)-linked disaccharide, including a nonconventional hydrogen bond between H5' in the glucosyl residue and O3 in the galactosyl residue, supported by a large downfield 1 H NMR chemical shift displacement compared to α-d-Glcp-OMe. Comparison of population distributions of the glycosidic torsion angles φ and ψ in the disaccharide entities to those of corresponding crystal structures highlighted the potential importance of solvation on the preferred conformation.

Multiscale structure and precipitation mechanism of debranched starch precipitated by different alcohols

Alcohol solution is a cheap, simple, and effective precipitating solvent frequently used for separating debranched starch (DBS), yet little is known about the precipitation mechanism of DBS by different alcohols. This study precipitated DBS from pullulanase-hydrolyzed starch using ethanol, n-butanol, and isopentanol. The multiscale structures of DBS were characterized, including chain length, single/double helix, and crystalline. The chain conformation and precipitation mechanism of DBS in different alcohols was investigated using molecular dynamics (MD) simulation. DBS precipitated by n-butanol contained the largest proportion of short chain (DP6-24, 83.2 %), the highest V-type crystallinity (21.1 %), and the largest single-helix content (24.7 %). A single helix conformation of DBS chain was determined in alcohols, where alcohol molecules entered the helix cavity. Intra/inter-molecular hydrogen bonds stabilized the helix, with a large number of hydrogen bonds leading to strong molecular interaction and stable helical structure. The solvent accessible surface area of DBS chain decreased by 7.88-19.32 % in alcohols, and the radial distribution function revealed that the first solvent layer of DBS chain at 0.29 nm was closely related to hydrogen bonding. This study provides a basis for the choice of precipitation solvent for preparing DBS with different chain lengths and physicochemical properties.

Formation and molecular dynamics simulation of inclusion complex of large-ring cyclodextrin and 4-terpineol

In the present study, the formation and structure of the inclusion compound of large-ring cyclodextrin and 4-terpineol were obtained through different experiments and molecular dynamics (MD) simulation. The analysis of FTIR, ¹ H-NMR, and thermodynamic results confirmed the formation of clathrates. Analysis of molecular structure (root-mean-square deviation and radius of gyration), solubility, and interaction energy (Coul, H bond) based on MD simulations further clarified the nature of the clathrate and the conformational changes caused by guest molecules as well as inclusion complexes process trends. The inclusion complex reportedly has a new crystal structure with improved thermal stability. PRACTICAL APPLICATION: This is the first work to demonstrate the complex formation between 4-terpineol and large-ring cyclodextrin by molecular dynamics simulation. Molecular dynamics simulation confirmed the formation of inclusion complexes theoretically. Conformational changes of the molecules and the formation of complexes with improved thermal stability were observed. Complexing with large-ring cyclodextrin can be used as an effective means to encapsulate the aroma/flavor compounds.

Reinvestigation of N,N-Diacetylimido-Protected 2-Aminothioglycosides in O-Glycosylation: Intermolecular Hydrogen Bonds Contributing to 1,2-Orthoamide Formation

N,N-Diacetylimido protection of 2-aminoglycosides is an elegant strategy but has had limited applications due to unexpected side reactions in glycosylation. We found that high acid concentrations could diminish the side reactions. We observed intermolecular hydrogen bonding among alcohols and acids could disrupt. Assuming that intermolecular hydrogen bonding accelerates the formation of 1,2-orthoamides and disrupting intermolecular hydrogen bonds could turn to the desired glycosylation, we successfully employed sulfenyl triflate pre-activation in the glycosylation of a broad scope of alcohol acceptors, as well as in a one-pot synthesis of a protected human milk oligosaccharide, lacto-N-neotetraose.

Driving torsion scans with wavefront propagation

The parameterization of torsional/dihedral angle potential energy terms is a crucial part of developing molecular mechanics force fields. Quantum mechanical (QM) methods are often used to provide samples of the potential energy surface (PES) for fitting the empirical parameters in these force field terms. To ensure that the sampled molecular configurations are thermodynamically feasible, constrained QM geometry optimizations are typically carried out, which relax the orthogonal degrees of freedom while fixing the target torsion angle(s) on a grid of values. However, the quality of results and computational cost are affected by various factors on a non-trivial PES, such as dependence on the chosen scan direction and the lack of efficient approaches to integrate results started from multiple initial guesses. In this paper, we propose a systematic and versatile workflow called TorsionDrive to generate energy-minimized structures on a grid of torsion constraints by means of a recursive wavefront propagation algorithm, which resolves the deficiencies of conventional scanning approaches and generates higher quality QM data for force field development. The capabilities of our method are presented for multi-dimensional scans and multiple initial guess structures, and an integration with the MolSSI QCArchive distributed computing ecosystem is described. The method is implemented in an open-source software package that is compatible with many QM software packages and energy minimization codes.

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.

Extension of the GROMOS 56a6CARBO/CARBO_R Force Field for Charged, Protonated, and Esterified Uronates

An extension of the GROMOS 56a6_{CARBO/CARBO_R} force field for hexopyranose-based carbohydrates is presented. The additional parameters describe the conformational properties of uronate residues. The three distinct chemical states of the carboxyl group are considered: deprotonated (negatively charged), protonated (neutral), and esterified (neutral). The parametrization procedure was based on quantum-chemical calculations, and the resulting parameters were tested in the context of (i) flexibility of the pyranose rings under different pH conditions, (ii) conformation of the glycosidic linkage of the (1 → 4)-type for uronates with different chemical states of carboxyl moieties, (iii) conformation of the exocyclic (i.e., carboxylate and lactol) moieties, and (iv) structure of the Ca²⁺-linked chain-chain complexes of uronates. The presently proposed parameters in combination with the 56a6_{CARBO/CARBO_R} set can be used to describe the naturally occurring polyuronates, composed either of homogeneous (e.g., glucuronans) or heterogeneous (e.g., pectins, alginates) segments. The results of simulations relying on the new set of parameters indicate that the conformation of glycosidic linkage is nearly unaffected by the chemical state of the carboxyl group, in contrary to the ring conformational equilibria. The calculations for the poly(α-d-galacturonate)-Ca²⁺ and poly(α-l-guluronate)-Ca²⁺ complexes show that both parallel and anitiparallel arrangements of uronate chains are possible but differ in several structural aspects.

Validation and Comparison of Force Fields for Native Cyclodextrins in Aqueous Solution

Molecular dynamics simulations of native α-, β-, and γ-cyclodextrin in aqueous solution have been conducted with the goal to investigate the performance of the CHARMM36 force field, the AMBER-compatible q4md-CD force field, and five variants of the GROMOS force field. The properties analyzed are structural parameters derived from X-ray diffraction and NMR experiments as well as hydrogen bonds and hydration patterns, including hydration free enthalpies. Recent revisions of the torsional-angle parameters for carbohydrate systems within the GROMOS family of force fields lead to a significant improvement of the agreement between simulated and experimental NMR data. Therefore, we recommend using the variant 53A6_GLYC instead of 53A6 and 56A6_{CARBO_R} or 2016H66 instead of 56A6_CARBO to simulate cyclodextrins in solution. The CHARMM36 and q4md-CD force fields show a similar performance as the three recommended GROMOS parameter sets. A significant difference is the more flexible nature of the cyclodextrins modeled with the CHARMM36 and q4md-CD force fields compared to the three recommended GROMOS parameter sets.