I want to study the interaction of cytochrome c with cosolvent by using GROMACS software package ,with OPLS/AA force field. But the forcefield parameter for HEME group is not included in OPLS . So how can I take this parameter?
You can look for atoms of the same chemical environment but in other molecules in the force field file.
Another method maybe is not very accurate.
You can get RESP charge from gaussian and then use Antechamber and acpype to obtain .gro、.itp、.top files for Gromacs. It will also output OPLS file of the molecule.