Hi,
You don't describe exactly what system you're dealing with, so this
may not be helpful, but
there are parameters for N-terminal amino acids (in all_aminont.lib)
So if you change the name of your protonated residue XXX in the pdb
file to NXXX, it might recognise the NH3+ group. The basic problem is
that leap has to match up atoms in residues in your pdb file to atoms
in the residues defined in the .lib file, and if you input two extra
hydrogens it will ignore them (maybe with a warning).
Dave Evans
London School of Pharmacy
---- Message from Daniel Wetzler <danielwetzler.appleinfo.de> at
2004-11-17 12:59:12 ------
>Dear Amber-users,
>
>I have a problem with a single aminoacid which acts like a inhibitor
for
>a enzyme I'm dealing with.
>My problem is, that this aminoacid should be fully protonated
(NH3-group).
>Regardless how I write this structure into my pdb-files (in its own
>chain, as HETATM
>a.s.o) xleap seems to handle it like a part of my structure and
>protonates my N only with one H.
>
>If I protonate my structure with the protonate command I get the
right
>structure with NH3
>but after using xleap to write the .top / .crd files I get the
>deprotonated form again.
>
>Does anyone know how to change that behaviour ?
>
>Hope anyone can help.
>
>Best wishes,
>
>Daniel
>---------------------------------------------------------------------
--
>The AMBER Mail Reflector
>To post, send mail to amber.scripps.edu
>To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Nov 17 2004 - 13:53:00 PST
