I have not looked at the source but I believe g_saltbr looks through the topology and lists all possible charge groups. It then analyzes the trajectory and outputs interactions between charge groups within the cutoff. The default cuttoff is 1000 nm so every charge group interaction should be accounted for. <div>
<br></div><div>Now that being said, I think I found the bug. I created a subset trajectory and tpr file that only included the ARG,LYS,APS,GLU residues. When I run g_saltbr on the trajectories containing all the residues this bug goes away. </div>
<div><br></div><div>Ilya</div><div><br><br><div class="gmail_quote">On Tue, Jul 27, 2010 at 1:39 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im"><br>
<br>
Ilya Chorny wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
The code would not know that /a priori/. This list is generated before the frames are analyzed.<br>
<br>
</blockquote>
<br></div>
I'm not implying that the code is going to pre-judge your system. If what you've shown is your output, wouldn't it make sense that g_saltbr only prints out the information for groups actually participating in putative salt bridges? Otherwise, the tool will print all charge groups and render itself useless. Do you for some reason expect the quoted arginines to be equivalent?<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Thanks,<br>
<br>
Ilya<div><div></div><div class="h5"><br>
<br>
On Tue, Jul 27, 2010 at 12:40 PM, Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a> <mailto:<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>>> wrote:<br>
<br>
<br>
<br>
Ilya Chorny wrote:<br>
<br>
I have a question about the g_saltbr output. I looked through<br>
the mailing list and could not find a suitable answer.<br>
The initial output lists all the charge group i.e.<br>
<br>
CG: ARG24 Q: 0.07 Atoms: 58 59 60<br>
CG: ARG24 Q: 0.31 Atoms: 61 62 63<br>
CG: ARG24 Q: 0.38 Atoms: 64 65 66<br>
CG: ARG24 Q: 0.12 Atoms: 67 68 69<br>
CG: ARG24 Q: 0.12 Atoms: 70 71 72<br>
CG: ASP28 Q: -0.1 Atoms: 79 80 81<br>
CG: ASP28 Q: -0.9 Atoms: 82 83 84<br>
CG: ARG36 Q: 0.07 Atoms: 94 95 96<br>
CG: ARG39 Q: 0.31 Atoms: 97 98 99<br>
CG: ARG42 Q: 0.38 Atoms: 100 101 102<br>
CG: ARG45 Q: 0.12 Atoms: 103 104 105<br>
CG: ARG48 Q: 0.12 Atoms: 106 107 108<br>
CG: ARG60 Q: 0.07 Atoms: 118 119 120<br>
CG: ARG63 Q: 0.31 Atoms: 121 122 123<br>
<br>
In some cases it lists 4 charge groups for ARG (i.e. ARG24) and<br>
in others it list 1 charge group (i.e. ARG 36). Why is this?<br>
<br>
<br>
Perhaps because not all the charge groups of Arg36 are participating<br>
in salt bridges, whereas in Arg24 they are.<br>
<br>
-Justin<br>
<br>
<br>
Thanks,<br>
<br>
Ilya<br>
<br>
<br>
<br>
-- Ilya Chorny Ph.D.<br>
<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></div></div>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> <<a href="http://vt.edu" target="_blank">http://vt.edu</a>> | (540) 231-9080<div class="im"><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>
========================================<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 <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before<br>
posting!<br>
Please don't post (un)subscribe requests to the list. Use the www<br>
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/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
<br>
<br>
<br>
<br>
-- <br>
Ilya Chorny Ph.D.<br>
<br>
</div></blockquote><div><div></div><div class="h5">
<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> | (540) 231-9080<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>
========================================<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/search" target="_blank">http://www.gromacs.org/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/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Ilya Chorny Ph.D.<br><br>
</div>