|
|
Joined: Oct 2003
Posts: 39
Forum Member
|
OP
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,882 Likes: 12
Forum Member
|
Forum Member
Joined: Sep 2003
Posts: 4,882 Likes: 12 |
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
Forum Member
Joined: Oct 2003
Posts: 39 |
Thanks.
Last edited by zhangxd68; 02/15/07 05:15 AM.
|
|
|
|
Joined: Sep 2003
Posts: 4,882 Likes: 12
Forum Member
|
Forum Member
Joined: Sep 2003
Posts: 4,882 Likes: 12 |
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
Forum Member
Joined: Oct 2003
Posts: 39 |
|
|
|
|
Joined: Sep 2003
Posts: 8,650 Likes: 26
Forum Member
|
Forum Member
Joined: Sep 2003
Posts: 8,650 Likes: 26 |
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
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,650 Likes: 26
Forum Member
|
Forum Member
Joined: Sep 2003
Posts: 8,650 Likes: 26 |
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
Forum Member
Joined: Oct 2003
Posts: 39 |
Thanks Prof Lennart and Rick.
|
|
|
|
Joined: Sep 2003
Posts: 8,650 Likes: 26
Forum Member
|
Forum Member
Joined: Sep 2003
Posts: 8,650 Likes: 26 |
COOR HBOND and COOR ANAL do work, and your effort may be better spent learning to use them, rather than "doing it by hand".
|
|
|
|
|
|