CHARMM force field parameters for metal ions Cd(II) and Pb(II)

GROMACS version: 2018
GROMACS modification: Yes/No
I want to simulate an enzyme with toxic metal ion Pb(II) using CHARMM force field. Does any one know for FF parameters for Pb(II) ion? Or how to generate parameters for metal ions?
Thank you.

Both of these ions are already supported in the current version of the CHARMM36 port for GROMACS. The parameters come from but please note none of these have been extensively validated for interactions with biomolecules and may require adjustment (e.g. NBFIX) to preserve appropriate coordination geometries.

I have a similar task - to minimize a protein molecule with zinc and cadmium ions. Is it necessary to generate separate topology for each of this atoms and include their parameters to the topology file as we would do for protein-ligand complex? Or using pdb2gmx command is enough in the case when metal ions are already included in the primary PDB file as HETATMs of the protein chain?