<br><br>----- Original Message -----<br>From: onetwo&nbsp; &lt;twoone21@rediffmail.com&gt;<br>Date: Monday, August 16, 2010 18:55<br>Subject: [gmx-users] Running MD on a dimeric protein<br>To: gmx-users@gromacs.org<br><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font>Hello Sir,<br> <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> I am simulating a protein which is a homodimer, when I did pdb2gmx, it generated a topolgy file  in which there were two chain topologies for each chain as :              <br> <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> ; Include chain topologies<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> #include "topol_A.itp"<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> #include "topol_B.itp"       <br> <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> and in Compound section gave ;<br> <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> Compound        #mols<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> Protein_A           1<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> Protein_B           1<br> <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> I want to ask that if it will consider the system as one or two different proteins.<br> <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> Also in production MD phase, in md.mdp file I mentioned<br> <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> Berendsen temperature coupling is on in two groups<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> Tcoupl              = V-rescale<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> tc-grps             = Protein Non-Protein<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> tau_t               = 0.1     0.1<br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> ref_t               = 300     300<br> <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font><br><font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;">&gt; </font> I have this doubt that as I have given tc-grps as Protein and Non-Protein, so if "Protein" will consider both the chains of the protein or not.<br><br> <font style="font-style: normal; font-weight: normal; background-color: rgb(245, 248, 240); font-size: 14px;"></font>The GROMACS index groups (of which tc-grps are perhaps the most commonly used, see manual for more info) are either provided to grompp via a file supplied with the -n parameter, or (in its absence) automatically generated by grompp. Either way, its output should provide enough clues to work out what has been regarded as Protein (e.g. residue or atom count). In extremis, you can use gmxdump on the .tpr to see the contents of the groups.<br><br>Unless you've gone out of your way to break things (e.g. changed the database of amino acid atom names), the default index groups generated by grompp will regard both your dimer parts as Protein.<br><br>If the dimer parts are adjacent, then tc-grps of Protein and Non-Protein are probably best. If they're apart, then you may want to couple the two parts separately. See http://www.gromacs.org/Documentation/Terminology/Thermostats in the first instance.<br><br>Mark