Before revision:  Our method advances the development of force fields for metal ions and the protein binding sites in dynamical environments.
After revision:  Our method advances the development of force fields for metal ions and protein binding sites in dynamical environments.