Thanks Justin, that was helpful.<br><br>I have a following question.<br><br>Since in my system I have both methanol and water and I want to order both of them ( my eventual aim to make a sphere), is there is way to override -na option ( for water na 4, methanol 3).<br>
Is there is way that all the components of the system get ordered with respect to the protein?<br><br>or I have have to play around ordering them one by one?<br><br>Thanks,<br>SN<br><br><br><div class="gmail_quote">On Thu, May 5, 2011 at 4:37 PM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im"><br>
<br>
shivangi nangia wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Hello,<br>
<br>
Came back to the set of calculations I was doing few days back.<br>
<br>
I have a box of water (TIP4P), methanol, protein and DHB anions.<br>
<br>
I minimized this system.<br>
<br>
Trying to use trjorder for water molecules with respect to the protein (eventually I want to make a sphere by ordering and getting rid of the molecules I do not want)<br>
<br>
In the index.ndx file i specified atom "144" (belongs to protein)<br>
<br>
so when I run<br>
<br>
trjorder -f em.gro -s em.tpr -n index.ndx -na 4 -o ordered.gro<br>
<br>
I pick number associated with atom "144" and SOL ( tip4p water) and I generate ordered.gro<br>
<br>
The sequence for SOL changes ( random), but when I use g_dist as I go down the ordered.gro, there is still no trend ( ascending with respect to protein)<br>
<br>
<br>
</blockquote>
<br></div>
I suspect your index file is wrong. Presumably, you have measured distances to individual water molecules using g_dist, which would require a custom index group for each water molecule analyzed, correct? If you use ordered.gro, the residue numbers will indeed be scrambled based on ordering, but the distance will correspond to the original distance, not the ordered one.<br>
<br>
The solution is to use your unordered coordinate file or .tpr file to create the index group, then perform the distance measurements. Alternatively, use genconf to renumber ordered.gro to get sequential residue numbers, then create whatever index file you might need from the renumbered coordinate file.<div class="im">
<br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
<br>
<br>
Also, -da 0 refers to an atom OR COM of the molecule.<br>
using trjorder -f em.gro -s em.tpr -n index.ndx -da 0 -na 4 -o ordered.gro and running it<br>
<br>
and choosing "1" for protein and "14" for SOL means that gromacs automatically understands that protein's COM is to be used to order SOL (tip4p, -na 4) ??<br>
<br>
</blockquote>
<br></div>
As the documentation is written, yes.<div class="im"><br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
If so, then tha talso generates a .gro file which produces random water molecules but there is again no trend.<br>
<br>
</blockquote>
<br></div>
Same problem as above, I suspect.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im">
I am very confused about using trjorder, it will be really helpful if someone or Mark can help me understand.<br>
<br>
Thanks a lot<br>
SN<br>
<br>
<br></div><div class="im">
On Thu, Apr 28, 2011 at 9:43 PM, Mark Abraham <<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a> <mailto:<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a>>> wrote:<br>
<br>
On 4/29/2011 11:29 AM, shivangi nangia wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Hello,<br>
<br>
The manual explaining trjorder says:<br>
<br>
trjorder orders molecules according to the smallest distance to<br>
atoms in a reference group or on z-coordinate (with option -z).<br>
With distance ordering, it will ask for a group of reference atoms<br>
and a group of molecules. For each frame of the trajectory the<br>
selected molecules will be reordered according to the shortest<br>
distance between atom number -da in the molecule and all the atoms<br>
in the reference group. *The center of mass of the molecules can<br>
be used instead of a reference atom by setting -da to 0*<br>
<br>
In order to arrange water molecules in accordance with the COM of<br>
the polypeptide, I chose -da 0.<br>
</blockquote>
<br>
As it says above, -da refers to an atom or COM of the molecule, not<br>
the reference group. This could be worded better in the documentation.<br>
<br>
Be sure you're choosing the groups you think you are choosing - you<br>
not copying relevant parts of your terminal output into emails is<br>
making things difficult.<br>
<br>
Mark<br>
<br>
</div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div class="im">
<br>
Am I wrong?<br>
<br>
Thanks,<br>
SN<br>
<br>
<br>
<br>
<br>
On Thu, Apr 28, 2011 at 8:29 PM, Mark Abraham<br></div><div class="im">
<<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a> <mailto:<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a>>> wrote:<br>
<br>
On 4/29/2011 4:08 AM, shivangi nangia wrote:<br>
<br>
Hello all,<br>
<br>
I am trying to order the TIP4P water molecules in my<br>
system with respect to the polypeptide in my system.<br>
<br>
The command I am using is:<br>
<br>
trjorder -f shape.gro -s shape.tpr -da 0 -na 4 -o<br>
ordered.gro<br>
<br>
This runs without any error and ordered.gro is generated<br>
with random sequence of water molecules.<br>
<br>
Just to cross check I calculated the distances between one<br>
of atoms of the polypeptide and oxyegn atom of different<br>
"ordered" water molecules.<br>
I found, there is no ascendig trend in the distances with<br>
respect to the polypeptide as a go down in the<br>
"ordered.gro" file.<br>
<br>
What could be going wrong?<br>
<br>
<br>
-da 0 has a particular effect - is it appropriate? Did you<br>
choose the right groups? You could use the -nshell option to<br>
probe what trjorder thinks is going on.<br>
<br>
Mark<br>
<br>
<br>
<br>
-- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><div class="im"><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at<br>
<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before<br>
posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.<div class="im"><br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
<br>
<br>
</div></blockquote><div class="im">
<br>
<br>
--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><div class="im"><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at<br>
<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.<div class="im"><br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
<br>
<br>
</div></blockquote>
<br>
-- <br>
========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | <a href="tel:%28540%29%20231-9080" value="+15402319080" target="_blank">(540) 231-9080</a><br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<div><div></div><div class="h5"><br>
-- <br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br>