<br /><br /><span>On 13/05/11, <b class="name">Anna Marabotti </b> &lt;anna.marabotti@isa.cnr.it&gt; wrote:</span><blockquote cite="mid:709E9A6BF36A4FC2B09AA5F9D51AF6D2@AnnaNetbook" class="iwcQuote" style="border-left: 1px solid rgb(0, 0, 255); padding-left: 13px; margin-left: 0pt;" type="cite"><div class="mimepart text html">
<span><p>
<meta content="text/html; charset=us-ascii" http-equiv="Content-Type" />
<meta content="MSHTML 8.00.6001.19046" name="GENERATOR" />
<table><tbody><tr><td><p>
</p><div><font face="Arial" size="2"><span class="218405310-13052011">Dear 
gmx-users,</span></font></div>
<div><font face="Arial" size="2"><span class="218405310-13052011">I'm simulating a 
homodimeric protein obtained by homology modelling. In the .pdb file, I have 2 
chain identifiers (A and B). I used Gromacs 4.0.7 to prepare the .gro files for 
simulation (up to nvt); then I was forced to switch to Gromacs 4.5.4 on another 
machine to continue with npt and full MD. I used the files from Gromacs 4.0.7 to 
continue, and it seems to me that all went OK. Now I have the files from the 
full MD, and I started to analyze them.</span></font></div>
<div><font face="Arial" size="2"><span class="218405310-13052011">I used g_rmsf to 
analyze the fluctuations of the residues, and I found that the results for the 
two subunits of the protein are superimposed. In the file .gro that comes from 
the MD as well as in the file .gro that comes from the npt MD made with Gromacs 
4.5.4, the number of residues of the second subunit does not continue after the 
first subunit (i.e. after the last 359th residue of the first subunit, I see 
that the first residue of the second subunit is numbered 1 instead of 360), but 
in the nvt file produced by Gromacs 4.0.7 that I provided to Gromacs 
4.5.4, the numbering of the two subunits is 
sequential. </span></font></div>
<div><font face="Arial" size="2"><span class="218405310-13052011">My questions 
are:</span></font></div>
<div><font face="Arial" size="2"><span class="218405310-13052011">1) why Gromacs 4.5.4 
changed the numbering of the residues in the two subunits, even if it started 
from a .gro file?</span></font></div></td></tr></tbody></table></p></span></div></blockquote><br />The trajectory files do not number the atoms, so the only way you can generate such numbering is by matching them with a structure file that also has such numbering. It's hard to be any more specific in the absence of command lines. If you just match the trajectory file with a more suitably numbered structure file (with atom ordering preserved) then everything is fine.<br /><br /><blockquote cite="mid:709E9A6BF36A4FC2B09AA5F9D51AF6D2@AnnaNetbook" class="iwcQuote" style="border-left: 1px solid rgb(0, 0, 255); padding-left: 13px; margin-left: 0pt;" type="cite"><div class="mimepart text html"><span><p><table><tbody><tr><td><div><font face="Arial" size="2"><span class="218405310-13052011"></span></font></div>
<div><font face="Arial" size="2"><span class="218405310-13052011">2) Does this affect 
the result of the MD?</span></font></div></td></tr></tbody></table></p></span></div></blockquote><br />I highly doubt it.<br /><br /><blockquote cite="mid:709E9A6BF36A4FC2B09AA5F9D51AF6D2@AnnaNetbook" class="iwcQuote" style="border-left: 1px solid rgb(0, 0, 255); padding-left: 13px; margin-left: 0pt;" type="cite"><div class="mimepart text html"><span><p><table><tbody><tr><td><div><font face="Arial" size="2"><span class="218405310-13052011"></span></font></div>
<div><font face="Arial" size="2"><span class="218405310-13052011">3) how can I 
separate, or renumber, the resulting .gro files in order to have a graph in 
which the fluctuations are not superimposed?</span></font></div></td></tr></tbody></table></p></span></div></blockquote><br />I rather expect that the input file to your NPT grompp had this &quot;feature&quot;. In any case, you need to identify where it first arose in your workflow in order to work out why it happened and what to do about it.<br /><br />Mark<br /><blockquote cite="mid:709E9A6BF36A4FC2B09AA5F9D51AF6D2@AnnaNetbook" class="iwcQuote" style="border-left: 1px solid rgb(0, 0, 255); padding-left: 13px; margin-left: 0pt;" type="cite"><div class="mimepart text html"><span><p>
</p></span></div></blockquote>