Iron sulfur cluster force field parameters

Hi! I need to simulate a protein with iron-sulfur clusters. I manually added the bonded parameters for iron-sulfur clusters. How should I deal with non-bonded parameters, and how are co-ordinate covalent bonds taken into could?
In general, could anyone please help with generating topology for these het. atoms?
Any suggestions would be appreciated.
Thanks!