[gmx-users] selecting res CYS from G53a5 using pdb2gmx

Payman Pirzadeh ppirzade at ucalgary.ca
Sun Jan 6 22:45:09 CET 2013

I am using gromacs 4.5.4, and I am generating a topology file using G53a5
force field. My protein has a Cys residue which is supposed to be in
thiolate state (negatively charged side chain). During issuing the pdb2gmx
command I use the switch -inter to be prompted about the charge state of the
ionisable side chains, but it asks nothing about Cys side chains except
about they being involved in disulfide bonds. My questions are:
1. how can I get pdb2gmx ask me about the charge status of Cys residues? Or
do I need to manually change their status in the topology file and output
gro file (removing the HG)?
2. I checked the aminoacid.rtp of G53a5, and it seems that for CYS, the
total charge of the residue does not add to -1 but -0.5. Is anyone aware of
any modifications to this force field parameters?


More information about the gromacs.org_gmx-users mailing list