AMBER: How to calculate hydrogen bond involved with water molecules quickly?

From: haixiao jin <jinhx952.gmail.com>
Date: Mon, 8 Aug 2005 13:07:48 +0800

Here is my input file for calculate hydrogen bond involved with water molecules.

Hbond1.in (the oxygen of water molecules is donor, and the hydrogen of
residues is acceptor)
trajin 1npt.mdcrd
center :1-347 origin mass
image origin center familiar
donor mask :345.O
donor mask :346.O
donor mask :347.O
acceptor ARG NE HE
acceptor ARG NH1 HH11
acceptor ARG NH1 HH12
acceptor ARG NH2 HH21
acceptor ARG NH2 HH22
acceptor ASN ND2 HD21
acceptor ASN ND2 HD22
acceptor CYS SG HG
acceptor GLN NE2 HE21
acceptor GLN NE2 HE22
acceptor HID ND1 HD1
acceptor SER OG HG
acceptor THR OG1 HG1
acceptor TYR OH HH
acceptor ALA N H
acceptor ARG N H
acceptor ASN N H
acceptor ASP N H
acceptor CYS N H
acceptor GLN N H
acceptor GLU N H
acceptor GLY N H
acceptor HIS N H
acceptor ILE N H
acceptor LEU N H
acceptor LYS N H
acceptor MET N H
acceptor PHE N H
acceptor PRO N H
acceptor SER N H
acceptor THR N H
acceptor TRP N H
acceptor TYR N H
acceptor VAL N H
hbond distance 3.5 angle 120.0 donor acceptor neighbor 2 series hbond

this file has been run without any problem. However, running hbond2.in
(the hydrogen of water molecules is acceptor and the heteroatom of
residues is donor) has obtained a empty result. Here is the hbond2.in
file

trajin 1npt.mdcrd
center :1-347 origin mass
image origin center familiar
donor mask :ALA.N,O
donor mask :ARG.N,O,NE,NH1,NH2
donor mask :ASN.N,O,OD1,ND2
donor mask :ASP.N,O,OD1,OD2
donor mask :CYS.N,O,SG
donor mask :GLN.N,O,OE1,NE2
donor mask :GLU.N,O,OE1,OE2
donor mask :GLY.N,O
donor mask :HIS.N,O,ND1
donor mask :ILE.N,O
donor mask :LEU.N,O
donor mask :LYS.N,O
donor mask :MET.N,O,SD
donor mask :PHE.N,O
donor mask :PRO.N,O
donor mask :SER.N,O,OG
donor mask :THR.N,O,OG1
donor mask :TRP.N,O
donor mask :TYR.N,O,OH
donor mask :VAL.N,O
acceptor mask :345.O H1
acceptor mask :345.O H2
acceptor mask :346.O H1
acceptor mask :346.O H2
acceptor mask :347.O H1
acceptor mask :347.O H2
hbond distance 3.5 angle 120.0 donor acceptor neighbor 2 series hbond


is the expression of acceptor in hbond2.in right? I have been try
"acceptor WAT O H1" and "acceptor WAT O H2" in hbond2.in, It fails due
to low internal memory. Any suggestion will be appreciated !
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Mon Aug 08 2005 - 06:53:00 PDT
Custom Search