I want to perform Molecular dynamic simulations of metalloprotein containing cadmium ion. To treat metal ions, there are different force field available, how can we decide which force field is appropriate for a particular metal ion?
An excellent paper has been published in the nature regarding your question by Kevin J. Waldron1, Julian C. Rutherford, Dianne Ford & Nigel J. Robinson as follows:
Nature 460, 823-830 (13 August 2009) | doi:10.1038/nature08300; Published online 12 August 2009
Simulate metals (in general) by MD simulations is a real challenge! Consider yourself that a metal in classical force fields is a point integer fixed charge and a pair of Lennard-Jones parameters...saying this, and far as I know, we can not simulate at the same time the water-exchange (that is an energetic term), the coordination geometry and ion-oxygen distance that define the chemistry of a metal! you must choose to sacrifice one parameter depending of your needs.
Please see our paper about this, in our case we performed classical TI simulations to find the best L-J parameters....but it was a hard work!!!