<SPAN>On 03/04/12, <B class=name>balaji nagarajan </B>&lt;balaji_sethu@hotmail.com&gt; wrote:</SPAN>
<BLOCKQUOTE style="BORDER-LEFT: #00f 1px solid; PADDING-LEFT: 13px; MARGIN-LEFT: 0px" class=iwcQuote cite=mid:SNT123-W17ADA135276CB06B15FE3CE24D0@phx.gbl type="cite">
<TABLE>
<TBODY>
<TR>
<TD class=hmmessage>
<P>
<DIV dir=ltr>Dear Users ! <br /><br />I have done a minimization for a tetrapeptide in a water box<br /><br />with conjugate gradient minimization of 1000 steps<br /><br />when I convert it to pdb as below <br />------------ <br /><br />trjconv -f em.trr -o min.pdb -s em.tpr &lt;&lt; EOF<br />0<br />EOF<br />--------------------<br /><br />the out put pdb has many structures <br />as follows <br /><br />MODEL&nbsp;&nbsp;&nbsp;&nbsp; 0 <br /><br />MODEL&nbsp;&nbsp;&nbsp;&nbsp; 1<br /><br />MODEL&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp; <br /><br />but the minimized pdb is single structre , <br /><br />what does these structures mean. <br /><br />how to obtain a single structre. <br /><br />my mimization input file is as follows<br />--------------------------------------<br />define&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -DFLEXIBLE<br />constraints&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = none<br />integrator&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = cg<br />dt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.001&nbsp;&nbsp;&nbsp;&nbsp; ; ps !<br />nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br />nstlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10<br />ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = grid<br />rlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<br />coulombtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = PME<br />rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<br />vdwtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = cut-off<br />rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<br />optimize_fft&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = yes<br />;<br />;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Energy minimizing stuff<br />;<br />emtol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10.0<br />emstep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.01<br /><br />energygrps = protein SOL<br />----------------------------------------<br /><br />thanks in advance Dear Users ! <br /><br />I have done a minimization for a tetrapeptide in a water box<br /><br />with conjugate gradient minimization of 1000 steps<br /><br />when I convert it to pdb as below <br />------------ <br /><br />trjconv -f em.trr -o min.pdb -s em.tpr &lt;&lt; EOF<br />0<br />EOF<br />--------------------</DIV></TD></TR></TBODY></TABLE></BLOCKQUOTE>
<DIV>&nbsp;</DIV><DIV>That will convert the contents of the trajectory file into a pseudo-trajectory in .pdb format. Hence you see multiple models. So you can simply look only at the last model, or direct trjconv to only convert the last trajectory frame, or perhaps use the file written by mdrun -c. You could have pre-solved the problem by writing the -c file in&nbsp;.pdb format by choosing mdrun -s yourem.tpr -c yourfinalstructure.pdb.&nbsp;</DIV><DIV>&nbsp;</DIV><DIV>Mark&nbsp;</DIV><DIV>&nbsp;</DIV><BLOCKQUOTE style="BORDER-LEFT: #00f 1px solid; PADDING-LEFT: 13px; MARGIN-LEFT: 0px" class=iwcQuote cite=mid:SNT123-W17ADA135276CB06B15FE3CE24D0@phx.gbl type="cite">
<DIV class="mimepart text html"><SPAN>
<P>
<STYLE><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></STYLE>

<TABLE>
<TBODY>
<TR>
<TD class=hmmessage>
<P>
<DIV dir=ltr>Dear Users ! <br /><br />I have done a minimization for a tetrapeptide in a water box<br /><br />with conjugate gradient minimization of 1000 steps<br /><br />when I convert it to pdb as below <br />------------ <br /><br />trjconv -f em.trr -o min.pdb -s em.tpr &lt;&lt; EOF<br />0<br />EOF<br />--------------------<br /><br />the out put pdb has many structures <br />as follows <br /><br />MODEL&nbsp;&nbsp;&nbsp;&nbsp; 0 <br /><br />MODEL&nbsp;&nbsp;&nbsp;&nbsp; 1<br /><br />MODEL&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp; <br /><br />but the minimized pdb is single structre , <br /><br />what does these structures mean. <br /><br />how to obtain a single structre. <br /><br />my mimization input file is as follows<br />--------------------------------------<br />define&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -DFLEXIBLE<br />constraints&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = none<br />integrator&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = cg<br />dt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.001&nbsp;&nbsp;&nbsp;&nbsp; ; ps !<br />nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br />nstlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10<br />ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = grid<br />rlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<br />coulombtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = PME<br />rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<br />vdwtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = cut-off<br />rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<br />optimize_fft&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = yes<br />;<br />;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Energy minimizing stuff<br />;<br />emtol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10.0<br />emstep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.01<br /><br />energygrps = protein SOL<br />----------------------------------------<br /><br />thanks in advance Dear Users ! <br /><br />I have done a minimization for a tetrapeptide in a water box<br /><br />with conjugate gradient minimization of 1000 steps<br /><br />when I convert it to pdb as below <br />------------ <br /><br />trjconv -f em.trr -o min.pdb -s em.tpr &lt;&lt; EOF<br />0<br />EOF<br />--------------------<br /><br />the out put pdb has many structures <br />as follows <br /><br />MODEL&nbsp;&nbsp;&nbsp;&nbsp; 0 <br /><br />MODEL&nbsp;&nbsp;&nbsp;&nbsp; 1<br /><br />MODEL&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp; <br /><br />but the minimized pdb is single structre , <br /><br />what does these structures mean. <br /><br />how to obtain a single structre. <br /><br />my mimization input file is as follows<br />--------------------------------------<br />define&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -DFLEXIBLE<br />constraints&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = none<br />integrator&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = cg<br />dt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.001&nbsp;&nbsp;&nbsp;&nbsp; ; ps !<br />nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br />nstlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10<br />ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = grid<br />rlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<br />coulombtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = PME<br />rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<br />vdwtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = cut-off<br />rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<br />optimize_fft&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = yes<br />;<br />;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Energy minimizing stuff<br />;<br />emtol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10.0<br />emstep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.01<br /><br />energygrps = protein SOL<br />----------------------------------------<br /><br />thanks in advance </DIV></TD></TR></TBODY></TABLE></SPAN></P></DIV></BLOCKQUOTE>
<DIV>&nbsp;</DIV>