Hi! I acetylated K40 my heterodimer's PDB file using PyTMs in PyMOL. Now I need to simulate it using GROMACS. Which force field would be the best for this purpose? Thank you in advance! Regards, Raag