Abstract:Glycosaminoglycans (GAGs) are negatively charged polysaccharides found on cell surfaces, where they regulate transport pathways of foreign molecules toward the cell. The structural and functional diversity of GAGs is largely attributed to varied sulfation patterns along the polymer chains, which makes understanding their molecular recognition mechanisms crucial. Molecular dynamics (MD) simulations, thanks to their unmatched microscopic resolution, have the potential to be a reference tool for exploring the patterns responsible for biologically relevant interactions. However, the capability of molecular dynamics force fields used in biosimulations to accurately capture sulfation-specific interactions is not well established, partly due to the intrinsic properties of GAGs that pose challenges for most experimental techniques. In this work, we evaluate the performance of molecular dynamics force fields for sulfated GAGs by studying ion pairing of Ca2+ to sulfated moieties─N-methylsulfamate and methylsulfate─that resemble N- and O-sulfation found in GAGs, respectively. We tested available nonpolarizable (CHARMM36 and GLYCAM06) and explicitly polarizable (Drude and AMOEBA) force fields, and derived new implicitly polarizable models through charge scaling (prosECCo75 and GLYCAM-ECC75) that are consistent with our developed "charge-scaling" framework. The calcium-sulfamate/sulfate interaction free energy profiles obtained with the tested force fields were compared against reference ab initio molecular dynamics (AIMD) simulations, which serve as a robust alternative to experiments. AIMD simulations indicate that the preferential Ca2+ binding mode to sulfated GAG groups is solvent-shared pairing. Only our scaled-charge models agree satisfactorily with the AIMD data, while all other force fields exhibit poorer agreement, sometimes even qualitatively. Surprisingly, even explicitly polarizable force fields display a notable disagreement with the AIMD data, likely attributed to difficulties in their optimization and possible inherent limitations in depicting high-charge-density ion interactions accurately. Finally, the underperforming force fields lead to unrealistic aggregation of sulfated saccharides, which qualitatively disagrees with our understanding of the soft glycocalyx environment. Our results highlight the importance of accurately treating electronic polarization in MD simulations of sulfated GAGs and caution against over-reliance on currently available models without thorough validation and optimization.

Molecular Dynamics Simulation of the Cooperative Effect by Different Force Fields in Monosodium Glutamate Aqueous Solution

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Simulation Study of Ion Pairing in Concentrated Aqueous Salt Solutions with a Polarizable Force Field

Protein Simulation Using Coarse-Grained Two-Bead Multipole Force Field with Polarizable Water Models.

Interaction of divalent cations and amino acids in bulk water: Molecular simulations with neural network potentials

Systematic comparison of force fields for microscopic simulations of NaCl in aqueous solutions: Diffusion, free energy of hydration and structural properties

Spatiotemporal characterization of water diffusion anomalies in saline solutions using machine learning force field

Comparative Study of the Water Response to External Force at Nanoscale and Mesoscale

Combined Fragment-Based Machine Learning Force Field with Classical Force Field and Its Application in the NMR Calculations of Macromolecules in Solutions.

Molecular Dynamics Simulation for Dmf Aqueous Solution

Developing and Benchmarking Sulfate and Sulfamate Force Field Parameters via Ab Initio Molecular Dynamics Simulations To Accurately Model Glycosaminoglycan Electrostatic Interactions

Tmff-A Two-Bead Multipole Force Field for Coarse-Grained Molecular Dynamics Simulation of Protein

Changes in the Local Conformational States Caused by Simple Na+ and K+ Ions in Polyelectrolyte Simulations: Comparison of Seven Force Fields with and without NBFIX and ECC Corrections

GCMe: Efficient implementation of the Gaussian core model with smeared electrostatic interactions for molecular dynamics simulations of soft matter systems

Examination of the Quality of Various Force Fields and Solvation Models for the Equilibrium Simulations of GA88 and GB88

Spatio-Temporal Characterization of Water Diffusion Anomalies in Saline Solutions Using Machine Learning Force Field

Developing and Benchmarking Sulfate and Sulfamate Force Field Parameters for Glycosaminoglycans via Ab Initio Molecular Dynamics Simulations

Multiscale Theory in the Molecular Simulation of Electrolyte Solutions

Molecular Dynamics in Physiological Solutions: Force Fields, Alkali Metal Ions, and Ionic Strength

A Coarse-Grained Force Field Parameterized for MgCl2 and CaCl2 Aqueous Solutions.

A Coarse-Grained Force Field Parameterized for MgCl2 and CaCl2 Aqueous Solutions