A consistent S-Adenosylmethionine force field improved by dynamic Hirshfeld-I atomic charges for biomolecular simulation.

S-Adenosylmethionine (AdoMet) is involved in many biological processes as cofactor in enzymes transferring its sulfonium methyl group to various substrates. Additionally, it is used as drug and nutritional supplement to reduce the pain in osteoarthritis and against depression. Due to the biological relevance of AdoMet it has been part of various computational simulation studies and will also be in the future. However, to our knowledge no rigorous force field parameter development for its simulation in biological systems has been reported. Here, we use electronic structure calculations combined with molecular dynamics simulations in explicit solvent to develop force field parameters compatible with the AMBER99 force field. Additionally, we propose new dynamic Hirshfeld-I atomic charges which are derived from the polarized electron density of AdoMet in aqueous solution to describe its electrostatic interactions in biological systems. The validation of the force field parameters and the atomic charges is performed against experimental interproton NOE distances of AdoMet in aqueous solution and crystal structures of AdoMet in the cavity of three representative proteins.