How to do the MD simulation of a metalloenzyme with ligand?

This topic is relevant: Fe2+ and Mn2+ force field parameters
if it doesn’t help you fix the problem, then @alevilla is the best person to answer since she is expert at force fields and one of the authors of the “GROMOS force-field parameter sets 53A5 and 53A6”.