Hyaluronic acid (HA) has a wide range of biomedical applications including the formation of hydrogels, microspheres, sponges, and films. The modeling of HA to understand its behavior and interaction with other biomolecules at the atomic level is of considerable interest. The atomistic representation of long HA polymers for the study of the macroscopic structural formation and its interactions with other polyelectrolytes is computationally demanding. To overcome this limitation, we developed a coarse grained (CG) model for HA adapting the Martini scheme. A very good agreement was observed between the CG model and all-atom simulations for both local (bonded interactions) and global properties (end-to-end distance, a radius of gyration, RMSD). Our CG model successfully demonstrated the formation of HA gel and its structural changes at high salt concentrations. We found that the main role of CaCl2 is screening the electrostatic repulsion between chains. HA gel did not collapse even at high CaCl2 concentrations, and the osmotic pressure decreased, which agrees well with the experimental results. This is a distinct property of HA from other proteins or polynucleic acids which ensures the validity of our CG model. Our HA CG model is compatible with other CG biomolecular models developed under the Martini scheme, which allows for large-scale simulations of various HA-based complex systems.
Keywords: Martini; coarse-grained simulation; gel formation; hyaluronic acid.