Force field parameters that use a combination of Lennard-Jones and electrostatic interactions are developed for divalent zinc and tested in solution and protein simulations. It is shown that the parameter set gives free energies of solution in good agreement with experiment. Molecular dynamics simulations of carboxypeptidase A and carbonic anhydrase are performed with these zinc parameters and the CHARMM 22 beta all-atom parameter set. The structural results are as accurate as those obtained in published simulations that use specifically bonded models for the zinc ion and the AMBER force field. The inclusion of longer-range electrostatic interactions by use of the Extended Electrostatics model is found to improve the equilibrium conformation of the active site It is concluded that the present parameter set, which permits different coordination geometries and ligand exchange for the zinc ion, can be employed effectively for both solution and protein simulations of zinc-containing systems.