Abstract

Although molecular dynamics simulations allow for the study of interactions among virtually all biomolecular entities, metal ions still pose significant challenges in achieving an accurate structural and dynamical description of many biological assemblies, particularly to coarse-grained (CG) models. Although the reduced computational cost of CG methods often makes them the technique of choice for the study of large biomolecular systems, the parameterization of metal ions is still very crude or not available for the vast majority of CG force fields. Here, we show that incorporating statistical data retrieved from the Protein Data Bank (PDB) to set specific Lennard-Jones interactions can produce structurally accurate CG molecular dynamics simulations using the SIRAH force field. We provide a set of interaction parameters for calcium, magnesium, and zinc ions, which cover more than 80% of the metal-bound structures reported in the PDB. Simulations performed on several proteins and DNA systems show that it is possible to preclude the use of topological constraints by modifying specific Lennard-Jones interactions.