Dear Prof. Cheatham:
I used ptraj to perform some water analysis, especially the "hbond" utility.
Very interesting tool; but I have some doubts:
I have a protein with a Glutamic Acid near (Oxygen of an) Asparagine
residue. During a 1.2ns simulation, I imagined that perhaps water molecule
might bridge in a H-bond. So I specified OE1 oxygen of the GLU residue as a
DONOR; and H1 and H2 of WAT as potential ACCEPTORS:
***part of my script***:
trajin boxed6b_alltraj_centrdimg.coord
donor mask :63.OE2
acceptor WAT O H1
acceptor WAT O H2
hbond series HBOND6B time 2 \
out wat_hbond_63OE2.dat solventneighbor 5 \
solventacceptor WAT O H1 solventacceptor WAT O H2 \
This gave me a wat_hbond_63OE2.dat file WHICH I QUITE DON'T UNDERSTAND.
Please find it attached. My questions are:
(1) The distances look OK, but ALL the angles are too low. Since I used a
cutoff of 120 deg (default), why these 20, 30, 40 etc. degrees are showing
up?
(2) What does that last column mean :
|.-xxo**o |
means 8 symbols in these case, and in other cases only a few or none.....
(3) What's the difference between "%occupied", "lifetime" and
"maxocc"....?
(4) Reading the 1st and 5th rows:
975 :63.OE2 | 33347 :9563.H2 33345 :9563.O | 45.83 2.699 ( 0.11)
17.10 ( 8.78) 17.2 ( 13.0) 68 |.-xxo**o |
and
975 :63.OE2 | 33346 :9563.H1 33345 :9563.O | 18.83 2.718 ( 0.12)
18.60 ( 9.00) 9.0 ( 3.5) 19 |x--.- . |
IS IT SAFE TO SAY that WAT9563 showed 68+19=87% maximum occupancy and
45.83+18.83=64.66% occupancy.....if its correct, then what's the difference
between two options.
(5) Reading the ultimate last row:
solvent donor | 5746 :363.H1 5745 :363.O | 156.83 2.695 ( 0.12)
16.56 ( 9.18) 10.6 ( 11.0) 68 |..........|
Seems that WAT363 has so many .=>100% occupancy, why so? What then about
WAT9563.....
May be it is rapidly exchanged....
I Thank you a Zillion in advance for your time and help!
Kind regards,
Vineet
_________________________________________________________________
Express yourself instantly with MSN Messenger! Download today it's FREE!
http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Mar 27 2005 - 15:53:00 PST
