<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Stefan Hoorman wrote:<br>
&gt;<br>
&gt;<br>
&gt;     Stefan Hoorman wrote:<br>
&gt;      &gt;<br>
&gt;      &gt;     Stefan Hoorman wrote:<br>
&gt;      &gt;      &gt; I have tried using g_bundle in order to analyse helix axes<br>
&gt;     in my<br>
&gt;      &gt;      &gt; transmembran helices. I created two groups in my ndx file that<br>
&gt;      &gt;     included<br>
&gt;      &gt;      &gt; the alpha carbons of the first half of my helix and a second<br>
&gt;      &gt;     group for<br>
&gt;      &gt;      &gt; the alpha carbons for the second half. When I try using<br>
&gt;     g_bundle the<br>
&gt;      &gt;      &gt; following error comes up:<br>
&gt;      &gt;      &gt; &quot;Fatal error:<br>
&gt;      &gt;      &gt; The size of one of your index groups is not a multiple of n&quot;<br>
&gt;      &gt;      &gt; What am I doing wrong?<br>
&gt;      &gt;<br>
&gt;      &gt;     Not a clue - What is your exact command line?  How many atoms<br>
&gt;     are in<br>
&gt;      &gt;     each index<br>
&gt;      &gt;     group?  The value of n is the number of parts into which the axes<br>
&gt;      &gt;     are divided<br>
&gt;      &gt;     and there are several criteria in the code that must be met, else<br>
&gt;      &gt;     you get this<br>
&gt;      &gt;     error.  So without this information, the best guess is: some<br>
&gt;      &gt;     criterion is not<br>
&gt;      &gt;     being met.<br>
&gt;      &gt;<br>
&gt;      &gt;     -Justin<br>
&gt;      &gt;<br>
&gt;      &gt;      &gt; Thank you<br>
&gt;      &gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;     --<br>
&gt;      &gt;<br>
&gt;      &gt;     g_bundle analyzes bundles of axes - and therefore number of atoms<br>
&gt;      &gt;     each group<br>
&gt;      &gt;     have to be divisible by number of axes (entered by -na). And I<br>
&gt;      &gt;     suppose, you<br>
&gt;      &gt;     are interested in analyzing only one axis, don&#39;t you?<br>
&gt;      &gt;<br>
&gt;      &gt;     Milan<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt; Well, the exact command line was : g_bundle -f md.xtc -s md.tpr -z -n<br>
&gt;      &gt; bundle.ndx<br>
&gt;      &gt; As for analyzing one axis, yes, that is true, I am only interested in<br>
&gt;      &gt; analysing my helix principal axis.<br>
&gt;      &gt; I have tried using different group sizes for the index group. I have<br>
&gt;      &gt; tried using only carbon alpha for the first half and then for the<br>
&gt;     second<br>
&gt;      &gt; half (each had 15 alpha carbons), then I tried with all sorts of<br>
&gt;      &gt; different groups, first five residues and last five residues, first<br>
&gt;      &gt; fifteen residues and so on.<br>
&gt;      &gt; All analysis came up with the same error.<br>
&gt;      &gt;<br>
&gt;<br>
&gt;     You&#39;re not setting the -na flag, which by default is zero.  In the<br>
&gt;     code, if the<br>
&gt;     value from -na (int n) is &lt;= 0, you get the fatal error.  Maybe the<br>
&gt;     error<br>
&gt;     message should be made more specific, so that users actually know<br>
&gt;     what n is.<br>
&gt;     You need to explicitly set a value for -na &gt;= 1.<br>
&gt;<br>
&gt;     -Justin<br>
&gt;<br>
&gt;<br>
&gt; Right, thank you Justin, things worked out nicely now. I was able to<br>
&gt; calculatemy helix tilt with respect to the z axis. Would there be a way<br>
&gt; to calculate the tilt between two helix axes? I cannot simply calculate<br>
&gt; the tilt of each with respect to the z axis and subtract because the<br>
&gt; interhelical tilt does not depend only on the z axis. Can g_bundle give<br>
&gt; me that kind of information?<br>
&gt;<br>
<br>
It seems like it should:<br>
<br>
<a href="http://lists.gromacs.org/pipermail/gmx-users/2007-October/030415.html" target="_blank">http://lists.gromacs.org/pipermail/gmx-users/2007-October/030415.html</a><br>
<br>
There are other programs that calculate angles, like g_sgangle and g_angle, both<br>
of which could potentially be used to calculate what you want.<br>
<br>
-Justin<br>
<font color="#888888"><br></font></blockquote><div>I have tried using g_sgangle, but the problem is that it does not
calculate the principal axis of the helix. Since g_sgangle allows me to
choose at maximum 3 atoms per analysis, it is quite difficult to
perform such an analysis if one is to do this for 4 or five different
systems and dimers.<br>So, in this case, since g_bundle asks me to select a group of top and a group of bottom atoms, I would need to create an index group with my first group (top) being the top residues of both my helices and the second group being the bottom residues of my helices. Or would I create a first group for my first helix and a second group for my second helix?<br>
</div></div><br>