I want to model a metalloprotein with a Fe-containing cofactor coordinated by several cysteine residues. I have obtained FF parameters for the cluster. What is the proper procedure for building the correct input for GROMACS simulation?