Previous Thread
Next Thread
Print Thread
Page 1 of 2 1 2
Joined: Oct 2003
Posts: 39
Forum Member
OP Offline
Forum Member
Joined: Oct 2003
Posts: 39
Hi,

Would you like to suggest how to analysis the water molecule's occpuancy
at a speical position during the MD simulations?

Thanks.

Joined: Sep 2003
Posts: 4,840
Likes: 3
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 4,840
Likes: 3
The CORRel module (correl.doc) can be used for an individual water molecule - eg to get its position or distance from a point vs time. To obtain a density map you can use COOR ANALysis (corman.doc).


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden
Joined: Oct 2003
Posts: 39
Forum Member
OP Offline
Forum Member
Joined: Oct 2003
Posts: 39
Thanks.

Last edited by zhangxd68; 02/15/07 05:15 AM.
Joined: Sep 2003
Posts: 4,840
Likes: 3
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 4,840
Likes: 3
Did you plot columns 2 and 3? Did any of them look like a g(r)?
Column 2 is g(r) and 3 is basically its integral, from which you can get the "solvation number". Coluimn 1 is not the number of frames, but it is quite easy to figure out what it is.
Since you are not doing a solvent-solvent g(r), but rather a site (your reference point) to solvent g(r) it may be that you should specify iSdist instead of iGdist - I don't remember exactly how the code works for this combination.


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden
Joined: Oct 2003
Posts: 39
Forum Member
OP Offline
Forum Member
Joined: Oct 2003
Posts: 39
Thanks Prof. Lennart.

Joined: Sep 2003
Posts: 8,596
Likes: 12
rmv Online Content
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 8,596
Likes: 12
To examine water around a peptide or protein site with COOR ANAL, I think you want to use the SITE keyword, and omit the WATER keyword, as in the 3rd g(r) example for water O around backbone HN in corman.doc

If the site is somewhat buried, you may wish to consider exploring the excluded volume options.

COOR HBOND might also be useful to compute how often the site is occupied, and the average lifetime of H-bonds.

Joined: Oct 2003
Posts: 39
Forum Member
OP Offline
Forum Member
Joined: Oct 2003
Posts: 39
Thanks Rick.

I figure out the error, I will do this by hand.

Last edited by zhangxd68; 02/15/07 05:20 AM.
Joined: Sep 2003
Posts: 8,596
Likes: 12
rmv Online Content
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 8,596
Likes: 12
Analysing hydrogen bond patterns for:
#donors #acceptors
1st sel: 767 1744
2nd sel: 0 0


For this command to work, both atom selections need to have either donors or acceptors (or both); the 2nd selection is clearly a problem (0 for both), hence the error message.

Your habit of posting cryptic excerpts out of context makes it difficult to say much more.

Joined: Oct 2003
Posts: 39
Forum Member
OP Offline
Forum Member
Joined: Oct 2003
Posts: 39
Thanks Prof Lennart and Rick.

Joined: Sep 2003
Posts: 8,596
Likes: 12
rmv Online Content
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 8,596
Likes: 12
COOR HBOND and COOR ANAL do work, and your effort may be better spent learning to use them, rather than "doing it by hand".

Page 1 of 2 1 2

Moderated by  BRBrooks, lennart, rmv 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.5
(Release build 20201027)
Responsive Width:

PHP: 7.3.31-1~deb10u1 Page Time: 0.005s Queries: 35 (0.002s) Memory: 0.7792 MB (Peak: 0.8705 MB) Data Comp: Off Server Time: 2022-05-25 12:50:45 UTC
Valid HTML 5 and Valid CSS