> ptraj finds the correct number of donors but it can not find any
> acceptors (and there should be at least 1 acceptor based on the ptraj
> scripts that work).
> So, the problem seems to be with the mask in the acceptor line and it
> does not seem to like the wild card or range specifications. (But, a
> wild card/range works with the mask for the donor command so it would
> seem that I've got the format of the mask command correct.) So,
> assuming I do have the format correct, why, with the acceptor command,
> can't I specify/use a wild card or range with mask?
....
> acceptor mask :*.N :*.H
There has to be a 1-1 correspondence between the acceptors, i.e. the first
N matched is paired with the first H matched, the second N with the second
H and so on... as the program is not smart enough (while remaining
completely general) to auto-detect what N should be paired with what H
atom.
The problem is that in this case, with your matching
:1-159.N or .N matches residues 1-159.N
:1-159.H or .H matches residues 2-159.H
(since residue 1, if it is a normal N-terminal replaces the H with H1, H2 and H3).
Therefore, your acceptor list would be
:1.N :2.H
:2.N :3.H
:3.N :4.H
...
and these will never be satisfied; moreover, ptraj will likely chuck out
the acceptors since there are more "N" atoms chosen than "H" atoms and
therefore it doesn't have a complete 1-1 match and will complain about
this.
To get it to work, do
acceptor mask :2-159.N :2-159.H
acceptor mask :1.N :1.H1
acceptor mask :1.N :1.H2
acceptor mask :1.N :1.H3
To get more information, before specifying the donors/acceptors, add
prnlev 4
and then set it back to zero after this (before doing the hbond command
unless you want a lot of debugging info)
prnlev 0
--tom
p.s. also note that :* is redundant; you just need .N, however, as you saw
this will not work since the N-terminal does not have atom H but atoms H1,
H2 and H3 bound to the N and there is no 1-1 match.
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Tue Jul 26 2005 - 21:53:01 PDT
