Hi, 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? Best,
Payman -- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! * Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
