<html>
  <head>

    <meta http-equiv="content-type" content="text/html; charset=UTF-8">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    <blockquote type="cite">
      <pre>Hsin-Lin Chiang wrote:
&gt;<i> Hi,
</i>&gt;<i> 
</i>&gt;<i> I'm trying to heat a protein.
</i>&gt;<i> There are two chains, A-chain and B-chain.
</i>&gt;<i> Two disulfide bonds are between A-chain and B-chain.
</i>&gt;<i> As I know, I should let A-chain and B-chain belong to the same [molecule 
</i>&gt;<i> type] in .top file if I want to have the two inter-bonds.
</i>&gt;<i> 
</i>&gt;<i> So, I delete the TER line between A-chain and B-chain in .pdb file.
</i>&gt;<i> And then I use "pdb2gmx -chainsep interactive" on it .
</i>&gt;<i> 
</i>&gt;<i> I got the two bonds successfully.
</i>&gt;<i> But after heating, the length of A-chain and B-chain in VMD was changed.
</i>&gt;<i> The number of residues in A-chain increase and the ones of B-chain 
</i>&gt;<i> decrease.
</i>&gt;<i> It means some residues move from B-chain to A-chain.
</i>&gt;<i> How should I prevent this kind of error happen.
</i>&gt;<i> 
</i>
The labeling changes?  Is this something VMD is doing or something you can 
actually demonstrate in the coordinate file?
</pre>
    </blockquote>
    Labeling didn't change.<br>
    I don't know what VMD did.<br>
    It seems to some residues disappear, A-chain was elongated, and
    B-chain was shortened.<br>
    <br>
    I also tried Pymol.<br>
    It gave me several small crash peptides.<br>
    <br>
    But I just tried Ramol, everything looks correct.<br>
    <br>
    Something wrong in VMD and Pymol may just because A-chain and
    B-chain curved extremely after heating.<br>
    <br>
    But I'm still afraid<br>
    Is my output structure wrong by GROMACS?<br>
    <br>
    Sincerely yours,<br>
    Hsin-Lin<br>
    <blockquote type="cite">
      <pre>
&gt;<i> Besides, I found even though I only have protein_A which consist both of 
</i>&gt;<i> A-chain and B-chain.
</i>&gt;<i> The number of first residue of B-chain is started from 1 in .gro file.
</i>&gt;<i> I thought it should be equal to "the number of A-chain's residues plus 1".
</i>&gt;<i> 
</i>&gt;<i> And it was also strange that the number of first water molecules was 
</i>&gt;<i> equal to "the number of B-chain plus 1"
</i>&gt;<i> I thought it should be equal to "the number of A-chain's residues plus 
</i>&gt;<i> the number of B-chain's residues plus 1"
</i>&gt;<i> 
</i>&gt;<i> Is it a bug?
</i>&gt;<i> 
</i>
No, it is the expected output.  By default, pdb2gmx does not renumber residues 
consecutively.  This feature was introduced some time ago so that output 
selections could easily be made in multiple chain molecules like this one.  If 
you want residues renumbered consecutively, either post-process the coordinate 
files with genconf -renumber or start over and use pdb2gmx -renum.
</pre>
    </blockquote>
    Thank you for your explanation.<br>
    <blockquote type="cite">
      <pre>
-Justin
</pre>
    </blockquote>
  </body>
</html>