<P>
&nbsp; <BR>
Thanks Justin for your reply <BR>
I am interested to get file which contain how many hydrogen bonds are there in my protein and the bonds present between which atoms?<BR>
Is there anyway to get this file <BR>
can you please tell me,<BR>
Thanks in advance for your kind suggestion&nbsp; <BR>
<BR>
On Thu, 28 Aug 2008 Justin A.Lemkul wrote :<BR>
&gt;<BR>
&gt;<BR>
&gt;minnale wrote:<BR>
&gt;&gt;&nbsp; <BR>
&gt;&gt;Thanks Justin for your quick reply<BR>
&gt;&gt;the g_hbond command ran succesfully and gave me .xvg file and it contain 3 columns<BR>
&gt;&gt;that is like<BR>
&gt;&gt;<BR>
&gt;&gt;@&nbsp; &nbsp; title &quot;Hydrogen Bonds&quot;<BR>
&gt;&gt;@&nbsp; &nbsp; xaxis&nbsp; label &quot;Time&quot;<BR>
&gt;&gt;@&nbsp; &nbsp; yaxis&nbsp; label &quot;Number&quot;<BR>
&gt;&gt;@TYPE xy<BR>
&gt;&gt;@ view 0.15, 0.15, 0.75, 0.85<BR>
&gt;&gt;@ legend on<BR>
&gt;&gt;@ legend box on<BR>
&gt;&gt;@ legend loctype view<BR>
&gt;&gt;@ legend 0.78, 0.8<BR>
&gt;&gt;@ legend length 2<BR>
&gt;&gt;@ s0 legend &quot;Hydrogen bonds&quot;<BR>
&gt;&gt;@ s1 legend &quot;Pairs within 0.35 nm&quot;<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp; &nbsp;  0&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 79&nbsp; &nbsp; &nbsp; &nbsp; 674<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp;  0.2&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 87&nbsp; &nbsp; &nbsp; &nbsp; 687<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp;  0.4&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 80&nbsp; &nbsp; &nbsp; &nbsp; 693<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp;  0.6&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 69&nbsp; &nbsp; &nbsp; &nbsp; 690<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp;  0.8&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 79&nbsp; &nbsp; &nbsp; &nbsp; 691<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp; &nbsp;  1&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 74&nbsp; &nbsp; &nbsp; &nbsp; 690<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp;  1.2&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 81&nbsp; &nbsp; &nbsp; &nbsp; 700<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp;  1.4&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 75&nbsp; &nbsp; &nbsp; &nbsp; 681<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp;  1.6&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 83&nbsp; &nbsp; &nbsp; &nbsp; 687<BR>
&gt;&gt;<BR>
&gt;&gt;the first column is time then could please tell about 2nd and 3rd columns<BR>
&gt;&gt;Thanks in advance.<BR>
&gt;<BR>
&gt;Read the header of the file, specifically:<BR>
&gt;<BR>
&gt;@ s0 legend &quot;Hydrogen bonds&quot;<BR>
&gt;@ s1 legend &quot;Pairs within 0.35 nm&quot;<BR>
&gt;<BR>
&gt;Gromacs always prints legends for what the data mean.<BR>
&gt;<BR>
&gt;-Justin<BR>
&gt;<BR>
&gt;&gt;<BR>
&gt;&gt;On Thu, 28 Aug 2008 Justin A.Lemkul wrote :<BR>
&gt;&gt;&nbsp; &gt;<BR>
&gt;&gt;&nbsp; &gt;<BR>
&gt;&gt;&nbsp; &gt;minnale wrote:<BR>
&gt;&gt;&nbsp; &gt;&gt;&nbsp; Hi all,<BR>
&gt;&gt;&nbsp; &gt;&gt;&nbsp; I am new to gromacs, I am interested in calculate Hydrogen bonds of my protein. So I have issued the *g_hbond command like this<BR>
&gt;&gt;&nbsp; &gt;&gt;g_hbond -f .xtc -s .tpr -n .ndx(Index file contain all residues mainchain+H atoms) -num .xvg -g .log<BR>
&gt;&gt;&nbsp; &gt;&gt;its showed<BR>
&gt;&gt;&nbsp; &gt;&gt;<BR>
&gt;&gt;&nbsp; &gt;&gt;Program g_hbond, VERSION 3.3.1<BR>
&gt;&gt;&nbsp; &gt;&gt;Source code file: statutil.c, line: 799<BR>
&gt;&gt;&nbsp; &gt;&gt;<BR>
&gt;&gt;&nbsp; &gt;&gt;Invalid command line argument:<BR>
&gt;&gt;&nbsp; &gt;&gt;-g<BR>
&gt;&gt;&nbsp; &gt;&gt;<BR>
&gt;&gt;&nbsp; &gt;&gt;Can you please tell where I have done mistake?<BR>
&gt;&gt;&nbsp; &gt;&gt;Thanks in advance for your valuable suggestions.<BR>
&gt;&gt;&nbsp; &gt;&gt;<BR>
&gt;&gt;&nbsp; &gt;&gt;<BR>
&gt;&gt;&nbsp; &gt;<BR>
&gt;&gt;&nbsp; &gt;The error tells you exactly what the problem is.&nbsp; There is no such option.&nbsp; The -g flag appears in newer versions of Gromacs.<BR>
&gt;&gt;&nbsp; &gt;<BR>
&gt;&gt;&nbsp; &gt;-Justin<BR>
&gt;&gt;&nbsp; &gt;<BR>
&gt;&gt;&nbsp; &gt;&gt;<BR>
&gt;&gt;&nbsp; &gt;&gt;<BR>
&gt;&gt;&nbsp; &gt;&gt;------------------------------------------------------------------------<BR>
&gt;&gt;&nbsp; &gt;&gt;<BR>
&gt;&gt;&nbsp; &gt;&gt;_______________________________________________<BR>
&gt;&gt;&nbsp; &gt;&gt;gmx-users mailing list&nbsp; &nbsp; gmx-users@gromacs.org<BR>
&gt;&gt;&nbsp; &gt;&gt;http://www.gromacs.org/mailman/listinfo/gmx-users<BR>
&gt;&gt;&nbsp; &gt;&gt;Please search the archive at http://www.gromacs.org/search before posting!<BR>
&gt;&gt;&nbsp; &gt;&gt;Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-request@gromacs.org.<BR>
&gt;&gt;&nbsp; &gt;&gt;Can't post? Read http://www.gromacs.org/mailing_lists/users.php<BR>
&gt;&gt;&nbsp; &gt;<BR>
&gt;&gt;&nbsp; &gt;-- ========================================<BR>
&gt;&gt;&nbsp; &gt;<BR>
&gt;&gt;&nbsp; &gt;Justin A. Lemkul<BR>
&gt;&gt;&nbsp; &gt;Graduate Research Assistant<BR>
&gt;&gt;&nbsp; &gt;Department of Biochemistry<BR>
&gt;&gt;&nbsp; &gt;Virginia Tech<BR>
&gt;&gt;&nbsp; &gt;Blacksburg, VA<BR>
&gt;&gt;&nbsp; &gt;jalemkul[at]vt.edu | (540) 231-9080<BR>
&gt;&gt;&nbsp; &gt;http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin<BR>
&gt;&gt;&nbsp; &gt;<BR>
&gt;&gt;&nbsp; &gt;========================================<BR>
&gt;&gt;<BR>
&gt;&gt;<BR>
&gt;&gt;<BR>
&gt;&gt;<BR>
&gt;&gt;------------------------------------------------------------------------<BR>
&gt;&gt;<BR>
&gt;&gt;_______________________________________________<BR>
&gt;&gt;gmx-users mailing list&nbsp; &nbsp; gmx-users@gromacs.org<BR>
&gt;&gt;http://www.gromacs.org/mailman/listinfo/gmx-users<BR>
&gt;&gt;Please search the archive at http://www.gromacs.org/search before posting!<BR>
&gt;&gt;Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-request@gromacs.org.<BR>
&gt;&gt;Can't post? Read http://www.gromacs.org/mailing_lists/users.php<BR>
&gt;<BR>
&gt;-- ========================================<BR>
&gt;<BR>
&gt;Justin A. Lemkul<BR>
&gt;Graduate Research Assistant<BR>
&gt;Department of Biochemistry<BR>
&gt;Virginia Tech<BR>
&gt;Blacksburg, VA<BR>
&gt;jalemkul[at]vt.edu | (540) 231-9080<BR>
&gt;http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin<BR>
&gt;<BR>
&gt;========================================<BR>

</P>
<br><br>
<Table border=0 Width=644 Height=57 cellspacing=0 cellpadding=0 style='font-family:Verdana;font-size:11px;line-height:15px;'><TR><td><a href='http://adworks.rediff.com/cgi-bin/AdWorks/click.cgi/www.rediff.com/signature-home.htm/1050715198@Middle5/2206641_2199021/2201649/1?PARTNER=3&OAS_QUERY=null' target=new ><img src ='http://imadworks.rediff.com/cgi-bin/AdWorks/adimage.cgi/2206641_2199021/creative_2201649.gif'  alt='Rediff Shopping'  border=0></a></td></TR></Table>