Molecular dynamics optimizations of five metalloproteins as an example show that their structures obtained by ABEEMσπ fluctuating charge force field with the calibrated ABEEMσπ parameters is very close to the structure optimized by the ab initio MP2/6-311G method.