Comparison of Carbohydrate Force Fields Using Gaussian Accelerated Molecular Dynamics Simulations and Development of Force Field Parameters for Heparin-Analogue Pentasaccharides

Computational description of conformational and dynamic properties of anticoagulant heparin analogue pentasaccharides is of crucial importance in understanding their biological activities. We designed and synthesized idraparinux derivatives modified with sulfonatomethyl moieties at the D, F, and H glucose units that display varied potencies depending on the exact nature of the substitution. In this report we examined the capability of molecular dynamics (MD) simulations to describe the conformational behavior of these novel idraparinux derivatives. We used Gaussian accelerated MD (GAMD) simulations on the parent compound, idraparinux, to choose the most suitable carbohydrate force field for these type of compounds. GAMD provided significant acceleration of conformational transitions compared to classical MD. We compared descriptors obtained from GAMD with NMR spectroscopic parameters related to geometrical descriptors such as scalar couplings and nuclear Overhauser effects (NOE) measured on idraparinux. We found that the experimental data of idraparinux is best reproduced by the CHARMM carbohydrate force field. Furthermore, we propose a torsion angle parameter for the sulfonatomethyl group, which was developed for the chosen CHARMM force field using quantum chemical calculations and validated by comparison with NMR data. The work lays down the foundation of using MD simulations to gain insight into the conformational properties of sulfonato-methyl group modified idraparinux derivatives and to understand their structure-activity relationship thus enabling rational design of further modifications.