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

User discussions
1

GROMACS version: 2023
GROMACS modification: Yes
Here post your question :

Recently, I want to make the simulation of a metalloenzyme with ligand. The metalloenzyme has two Mn ions in its active site. When producing the topol file for this metalloenzyme, gromacs cannot find the Mn ion in the rtp file. I have tried many methods. The one is i change the residue name of Mn in the pdb file to Zn, then the topool file of metalloenzyme was produced well using the Gromos 53a6 forcefield. Then i used the ATB server to produce the topol file for the ligand. However, when conducted the gmx genion, it cannot pass.
How can I make a suitable topol file for the metalloenzyme and ligand?
Thanks.

2

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”.

3

Thanks for your help!

4

I also have the same problem. I am new in GROMACS and do not know how to handle Fe in my protein. i try to contact to @alevilla but her account is not available. if you know any other source kindly let me know. Thank you