We have developed new force field and parameters for copper(I) and mercury(II) to be used in molecular dynamics simulations of metalloproteins. Parameters have been derived from fitting of ab initio interaction potentials calculated at the MP2 level of theory, and results compared to experimental data when available. Nonbonded parameters for the metals have been calculated from ab initio interaction potentials with TIP3P water. Due to high charge transfer between Cu(I) or Hg(II) and their ligands, the model is restricted to a linear coordination of the metal bonded to two sulfur atoms. The experimentally observed asymmetric distribution of metal ligand bond lengths (r) is accounted for by the addition of an anharmonic (r3) term in the potential. Finally, the new parameters and potential, introduced into the CHARMM force field, are tested in short molecular dynamics simulations of two metal thiolates fragments in water. (Brooks BR et al. J Comput Chem 1983, 4, 1987.[1]) © 2006 Wiley Periodicals, Inc. J Comput Chem 27: 837-856, 2006
