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, with their unmatched microscopic perspective, have the potential to be a reference tool for exploring the patterns responsible for biologically relevant interactions. However, the capability of molecular dynamics models (i.e., force fields) used in biosimulations to accurately capture sulfation-specific interactions is not well established. 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 nonpolarizable (CHARMM36 and GLYCAM06), explicitly polarizable (Drude and AMOEBA), and implicitly polarizable through charge scaling (prosECCo75 and GLYCAM-ECC75) force fields. The Ca--sulfamate/sulfate interaction free energy profiles obtained with the tested force fields were compared against reference ab initio molecular dynamics (AIMD) simulations. AIMD reveals that the preferential Ca2+ binding mode to sulfated GAG groups is solvent-shared pairing, and only the charge-scaled models agree satisfactorily with the AIMD data. All other force fields exhibit poorer performance, sometimes even qualitatively. Surprisingly, even explicitly polarizable force fields display a notable shortfall in their performance, 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, qualitatively distorting 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.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to evaluate the performance of existing molecular dynamics force - field models in simulating sulfated glycosaminoglycans (GAGs), especially the ability of these models to capture specific interactions of sulfate groups. Specifically, the researchers evaluated the performance of different types of force - field models (non - polarized, explicitly polarized, and implicitly polarized) by studying the coordination patterns of calcium ions (Ca²⁺) with sulfated and sulfonated groups (simulating N - and O - sulfation respectively), and compared them with the results of high - precision ab initio molecular dynamics (AIMD) simulations. ### Research Background GAGs are a class of negatively charged polysaccharides that exist on the cell surface and regulate the transport pathways of foreign molecules into the cell interior. The structural and functional diversity of GAGs is mainly attributed to the sulfation patterns on their polymer chains, which makes it crucial to understand their molecular recognition mechanisms. Although molecular dynamics (MD) simulations provide a microscopic perspective and are helpful for exploring the patterns of biologically relevant interactions, the ability of current molecular dynamics models (i.e., force fields) used for biological simulations to accurately capture sulfation - specific interactions has not been fully verified. ### Research Methods The researchers selected N - methylsulfamate (NMSM) and methylsulfate (MSO₄) as model molecules to simulate N - and O - sulfated groups in GAGs respectively. They tested non - polarized (CHARMM36 and GLYCAM06), explicitly polarized (Drude and AMOEBA), and force - field models with implicit polarization achieved through charge scaling (prosECCo75 and GLYCAM - ECC75). By using the umbrella sampling and accelerated weight histogram (AWH) methods, the researchers calculated the binding free energy of Ca²⁺ with NMSM and MSO₄ and compared it with the results of AIMD simulations. ### Main Findings 1. **Binding Pattern**: AIMD simulations showed that the preferred binding pattern of Ca²⁺ with sulfated GAG groups is solvent - shared pairing, and only the charge - scaling models consistently predicted this pattern with the AIMD data. 2. **Model Performance**: Other force - field models showed poor performance and sometimes even qualitative errors. Even the explicitly polarized force fields also showed obvious deficiencies, which may be related to their optimization difficulties and the inherent limitations in depicting high - charge - density ion interactions. 3. **Aggregation Phenomenon**: Poor - performing force fields led to unrealistic aggregation of sulfated sugar molecules, distorting the understanding of the glycocalyx environment. ### Conclusion The study emphasizes the importance of accurately handling electronic polarization in MD simulations of GAGs and warns against over - relying on existing models that have not been thoroughly verified and optimized. To ensure the accuracy of future large - scale simulations, it is necessary to use charge - scaling force fields and further develop or improve existing GAGs models.

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

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

Extension and validation of the GLYCAM force field parameters for modeling glycosaminoglycans

Benchmarking Water Models in Molecular Dynamics of Protein–Glycosaminoglycan Complexes

Coarse-Grained Model of Glycosaminoglycans for Biomolecular Simulations

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

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

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

AUTOMATED FORCE FIELD PARAMETERIZATION FOR NON-POLARIZABLE AND POLARIZABLE ATOMIC MODELS BASED ON AB INITIO TARGET DATA.

Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties

Simulating Bacterial Membrane Models at the Atomistic Level: A Force Field Comparison

Comparative Assessment of Water Models in Protein-Glycan Interaction: Insights from Alchemical Free Energy Calculations and Molecular Dynamics Simulations

Benchmarking polarizable and non-polarizable force fields for Ca2+-peptides against a comprehensive QM dataset

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

Sulfation Effects of Chondroitin Sulfate to Bind a Chemokine in Aqueous Medium: Conformational Heterogeneity and Dynamics from Molecular Simulation

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

An Automated, High-Throughput Method for Interpreting the Tandem Mass Spectra of Glycosaminoglycans

Development of a fused-sphere SAFT-γ Mie force field for poly(vinyl alcohol) and poly(ethylene)

Ab Initio Molecular Dynamics Simulations of the S N 1/S N 2 Mechanistic Continuum in Glycosylation Reactions

Development of the Quantum Inspired SIBFA Many-Body Polarizable Force Field: Enabling Condensed Phase Molecular Dynamics Simulations

Balancing Group I Monoatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field