Could someone tell how can i modify the specific force field, for example gromacs43a1, to build topology file for protein has iron that interact with some ligands (cytochrome, hemoglobin...etc)? which file i should modify?