[gmx-users] pdb2gmx: can we use a file to specify the protonation instead of using "-inter"?

[Christopher Neale]

Dear Cheng:

You can define protonation states like this:

echo "0 1 2" | pdb2gmx -ff amber99sb-ildn -water tip3p -his

That will give you, in order for three HIS residues:
0. H on ND1 only (HID)
1. H on NE2 only (HIE)
2. H on ND1 and NE2 (HIP)

You can also use -inter, but you'll need to do a dry-run to see what order gmx asks you for these things in. From my tests, it follow an order that is first residue type (in order: LYS, ARG, GLN, ASP, GLU, HIS, disulfide bonds) and for each of these, order by residue number.

Chris.

________________________________________
From: gromacs.org_gmx-users-bounces at maillist.sys.kth.se <gromacs.org_gmx-users-bounces at maillist.sys.kth.se> on behalf of Zhang, Cheng <c.zhang.11 at ucl.ac.uk>
Sent: 21 May 2015 10:38
To: gromacs.org_gmx-users at maillist.sys.kth.se
Cc: Zhang, Cheng
Subject: [gmx-users] pdb2gmx: can we use a file to specify the protonation instead of using "-inter"?

Dear GROMACS experts,

Can I ask is there a more efficient way to deal with the "-inter" option in the pdb2gmx command?


Now, I have to manually assign individual protonation status one by one, which takes a very long time. Sometimes I make a mistake and I have to re-do all of them again.


Would it be possible if the command line can link to a file, which specifies the protonation for individual residues?


Thank you very much.


Yours sincerely

Cheng
--
Gromacs Users mailing list

* Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a mail to gmx-users-request at gromacs.org.