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.

Polarizable AMOEBA force field predicts thin and dense hydration layer around monosaccharides for accurate simulations of carbohydrates in solution

Drude Polarizable Force Field Parametrization of Carboxylate and N-acetyl Amine Carbohydrate Derivatives.

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

Polarizable Mean-Field Model of Water for Biological Simulations with Amber and Charmm force fields

Overcoming the limitations of the MARTINI force field in Molecular Dynamics simulations of polysaccharides

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

Parametrizing A Polarizable Force Field from Ab Initio Data. I. the Fluctuating Point Charge Model

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

Polarizable Multipole-Based Force Field for Dimethyl and Trimethyl Phosphate

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

Accounting for electronic polarization in nonpolarizable force fields

Modeling polyelectrolyte hydration from a multi scale polarizable pseudo particle solvent coarse grained approach

A Polarizable Atomic Multipole-Based Force Field for Molecular Dynamics Simulations of Anionic Lipids

Constant-pH Simulations with the Polarizable Atomic Multipole AMOEBA Force Field

Validation of Polarizable Force Field Parameters for Nucleic Acids by Inter-Molecular Interactions

Balanced polarizable Drude force field parameters for molecular anions: phosphates, sulfates, sulfamates, and oxides

Accurate transferable polarization model derived from the monomer electron density

Refinement of the Drude Polarizable Force Field for Hexose Monosaccharides: Capturing Ring Conformational Dynamics with Enhanced Accuracy

A New Parameterization of an All-Atom Force Field for Cellulose

Multiscale Polarizable Coarse-Graining Water Models on Cluster-Level Electrostatic Dipoles.

The multi-scale polarizable pseudo-particle solvent coarse-grained approach: From NaCl salt solutions to polyelectrolyte hydration