I acetylated 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, Stephane. I have simulated the heterodimer before. I have also done in silico site directed mutagenesis and simulated the mutants. So I am not worried about the computational resources part. What I am worried about, however, is that I am getting a "fatal error" with all 15 of the GROMACS force fields, as soon as I acetylate the molecule.
You are getting error message because available force fields in the GROMACS do not support acetylated residue. You can add parameters for the modified residue in the parameter files. Please read the given article.