GROMACS version: 2019
GROMACS modification: Yes/No
Here post your question
Hello Gromacs developers and users,
I’m dealing with a highly negative charged protein-RNA-ligand complex which happens to create a lot of problems with lots of counter ions being added to the system to effectively neutralize it, and during dynamics the ions happen to go and interact with the ligand which is not what i want. Have tried using the Gromacs documentation recommended van der waal and coulombic cut-offs for PME and CHARMM36 force field with no success.
Is there any way to scale down the partial charges of charged groups in Gromacs, which can be done using CHARMM and have been performed in a few papers Why base tautomerization does not cause errors in mRNA decoding on the ribosome. Here in methods section they write, “Phosphate groups beyond 22Å from the simulation center were neutralized by scaling down the partial charges.” Although the authors of this paper have used their own simulation codes, can in anyway this method be performed in Gromacs.?