AMBER Archive (2007)Subject: Re: AMBER: adding polar hydrogens
From: j j (fantvamp_at_gmail.com)
Date: Thu May 24 2007 - 06:09:04 CDT
Dear Amber User,
I'm trying to add polar Hydrogens by using protonate by doing:
protonate -d PROTON_INFO.kollua_polH < prot_in.pdb > prot_out_polH.pdb
When I check results I think HG is missing in CYS residues:
charge CYS_prot_out_polH.pdbq
ATOM 101 N CYS B 12 -13.623 -0.254 8.103 0.00 0.00 -0.520
ATOM 102 HN CYS B 12 -14.244 0.084 7.382 0.00 0.00 0.248
ATOM 103 CA CYS B 12 - 12.630 -1.259 7.786 0.00 0.00
0.146
ATOM 104 C CYS B 12 -11.306 -0.844 8.366 0.00 0.00 0.526
ATOM 105 O CYS B 12 -10.816 -1.447 9.314 0.00 0.00 -0.500
ATOM 106 CB CYS B 12 - 12.475 -1.416 6.290 0.00 0.00
0.100
ATOM 107 SG CYS B 12 -11.079 -2.451 5.880 0.00 0.00 -0.135
Total Charge -0.1350 in CYS_prot_out_polH.pdbq
I expected to have:
charge CYS
q[ "CYS", "N" ] = -0.5200
q[ "CYS", "CA" ] = 0.1460
q[ "CYS", "HN" ] = 0.2480
q[ "CYS", "C" ] = 0.5260
q[ "CYS", "O" ] = -0.5000
q[ "CYS", "CB" ] = 0.1000
q[ "CYS", "SG" ] = -0.1350
q[ "CYS", "HG" ] = 0.1350
Total Charge 0.0000 in CYS
as especified in the q.kollua.amber
To do this I though in modifying the PROTON_INFO.kollua_polH by adding:
29,30c29,30
< CYS 6
< CA N C O CB SG
---
> CYS 7
> CA N C O CB SG HG
As specified in the protonate.f:
c*****add all protons to a protein structure according to formulae
c*****specified in a PROTON_INFO file.
I'm missing again the HG so do not seem to work.
I saw this message about it in the reflector,
so I'm writting.
Thank you,
JJ
2006/10/11, Sean Rathlef <sean_at_syncitium.net>:
>
> I'm trying to add polar hydrogens to a PDB file using PROTONATE.exe(oringinally from autodock, I think). I tried the command lines with the
> executable, but could not generate the desired output PDB file. If anyone
> knows the proper command line calls (perhaps an example), please advise.
>
> I tried the following:
> protonate.exe < name.pdb > out.pdb
>
> which generated:
> PROTON_INFO, protonate.exe.stackdump, and out.pdb
>
> The only file with anything in it was the stackdump
> Stack trace:
> Frame Function Args
> End of stack trace
>
> Do I need to use PROTON_INFO.kollua_polH, as in:
> protonate.exe -k PROTON_INFO.kollua_polH < name.pdb > out.pdb
>
> Because I tried that, also with the "-d" flag an that didn't work either.
>
> Help greatly appreciated.
> Sean
>
>
--
fantvamp
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|