GROMACS version: 2018.4
GROMACS modification: No
Hi everyone,
I’m trying to use pdb2gmx to set up my protein while controlling what residues are protonated. I of course use the -inter flag for this. However, since manually entering each residue charge state isn’t feasible for my application, I’m using a txt file with all the protonation states that pdb2gmx can read from. The full command I use is:
gmx_mpi pdb2gmx -f protein.pdb -o protein.gro -water none -ff oplsaa -inter < Charges.txt
Charges.txt is a txt file that only contains the appropriate number on each line, like so:
1
1
0
0
1
etc.
This was working perfectly on a system with gromacs version 5.0.4 installed. That system is being retired though, so I’ve moved to the new cluster with the same OS (a Unix build) but with version 2018.4 installed. When I try and run this command with the version gromacs will not read the file for the inter flag. Specifically the error I get is:
Fatal error:
Answer for me res LYS 16!
Where LYS 16 is the first protonatable residue. I’ve also tried using a here file, terminating the command with:
<<EOF
1
1
0
etc
EOF
But I still get the same error. Trying to use echo and expect terminology also doesn’t work. I’m out of ideas, I’ve tried all the combinations of the above I can think of. I find it really odd that this works just fine on one version, but not at all on a newer one. Any help or ideas you guys can provide is greatly appreciated.