Abstract:The conformational features of the glycosidic linkage are the most important variable to consider when studying di-, oligo-, and polysaccharide molecules using molecular dynamics (MD) simulations. The accuracy of the theoretical model describing this degree of freedom influences the quality of the results obtained from MD calculations based on this model. This article focuses on the following two issues related to the conformation of the glycosidic linkage. First, we describe the results of a comparative analysis of the predictions of three carbohydrate-dedicated classical force fields for MD simulations, namely, CHARMM, GLYCAM, and GROMOS, in the context of different parameters of structural and energetic nature related to the conformation of selected types of glycosidic linkages, α(1 → 4), β(1 → 3), and β(1 → 4), connecting glucopyranose units. This analysis revealed several differences, mainly concerning the energy levels of the secondary and tertiary conformers and the linkage flexibility within the dominant exo-syn conformation for α(1 → 4) and β(1 → 3) linkages. Some aspects of the comparative analysis also included the newly developed, carbohydrate-dedicated Martini 3 coarse-grained force field. Second, to overcome the time-scale problem associated with sampling slow degrees of freedom in polysaccharide chains during MD simulations, we developed a coarse-grained (CG) model based on the data from MD simulations and designed for Monte Carlo modeling. This model (CG MC) is based on information from simulations of short saccharide chains, effectively sampled in atomistic MD simulations, and is capable of extrapolating local conformational properties to the case of polysaccharides of arbitrary length. The CG MC model has the potential to estimate the conformations of very long polysaccharide chains, taking into account the influence of secondary and tertiary conformations of glycosidic linkages. With respect to the comparative analysis of force fields, the application of CG MC modeling showed that relatively small differences in the predictions of individual force fields with respect to a single glycosidic linkage accumulate when considering their effect on the structure of longer chains, leading to drastically different predictions with respect to parameters describing the polymer conformation, such as the persistence length.

A first-principles exploration of the conformational space of neutral and sodiated di-saccharides assisted by semi-empirical methods and neural network potentials

Searching low-energy conformers of neutral and protonated di-, tri-, and tetra-glycine using first-principles accuracy assisted by the use of neural network potentials

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

Exploring conformational energy landscape of glassy disaccharides by CPMAS 13C NMR and DFT/GIAO simulations. I. Methodological aspects

Conformational analysis and flexibility of carbohydrates using the CICADA approach with MM3

The Conformation of Glycosidic Linkages According to Various Force Fields: Monte Carlo Modeling of Polysaccharides Based on Extrapolation of Short-Chain Properties

First-Principles Molecular Structure Search with a Genetic Algorithm

Optimizing Nonbonded Interactions of the OPLS Force Field for Aqueous Solutions of Carbohydrates: How to Capture Both Thermodynamics and Dynamics

Accelerated lignocellulosic molecule adsorption structure determination

Precise AIE‐Based Ternary Co‐Assembly for Saccharide Recognition and Classification

A vibrational molecular force field of model compounds with biological interest. IV. Parameters for the different glycosidic linkages of oligosaccharides

Exploring conformational energy landscape of glassy disaccharides by CPMAS 13C NMR and DFT/GIAO simulations. II. Enhanced molecular flexibility in amorphous trehalose

Exploring the Molecular Conformation Space by Soft Molecule–Surface Collision

Modern Semiempirical Electronic Structure Methods and Machine Learning Potentials for Drug Discovery: Conformers, Tautomers, and Protonation States

A Scalable Graph Neural Network Method For Developing An Accurate Force Field Of Large Flexible Organic Molecules

Single-Molecule Discrimination of Saccharides Using Carbon Nitride Nanopores.

Revealing elusive conformations of sucrose from hydrogen bond J-coupling in H2O: A combined NMR and quantum mechanics study

Hydrogen bonding and dominant conformations of hydrated sugar analogue complexes using tetrahydrofurfuryl alcohol as the model sugar molecule.

Polysaccharide conformations measured by solution state X-ray scattering

Direct imaging of carbohydrate stereochemistry