Hi Hiromichi,<br><br>It is likely that you have some groups which are not covered by &quot;Protein&quot; and &quot;Other&quot;. You can use &quot;Protein&quot; and &quot;Non-Protein&quot; in stead. But make sure that the two groups correspond to two distinct groups in your system. You may run into trouble if you have some non-standard residue in your protein which will end up in the &quot;Non-Protein&quot; group.
<br><br>Cheers,<br><br>Tsjerk<br><br><div><span class="gmail_quote">On 6/14/06, <b class="gmail_sendername">Hiromichi Tsurui</b> &lt;<a href="mailto:tsurui@med.juntendo.ac.jp">tsurui@med.juntendo.ac.jp</a>&gt; wrote:</span>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div>





<div bgcolor="#ffffff">
<div><font face="MS UI Gothic" size="2">
</font><p style="margin: 0mm 0mm 0pt;"><span lang="EN-US"><font face="MS UI Gothic" size="2"><font size="3">Dear all;</font></font></span></p>
<p style="margin: 0mm 0mm 0pt;"><span lang="EN-US"><font face="MS UI Gothic" size="2"><font size="3">&nbsp;</font></font></span></p>
<p style="margin: 0mm 0mm 0pt;"><span lang="EN-US"><font face="MS UI Gothic" size="2"><font size="3">I 
am trying to calculate the binding free energy of an MHC-molecule and a peptide 
by thermodynamic integration method. The system contains binding domain of MHC 
molecule, a peptide presented on it, water-molecules (about 13,000), and ions 
(Na+ and Cl- at PBS-concentration).</font></font></span></p>
<p style="margin: 0mm 0mm 0pt;"><span lang="EN-US"><font face="MS UI Gothic" size="2"><font size="3">After energy minimization, I tried to equilibrate the system with the 
mdp-file as follows, and an error-message appeared as 
such;</font></font></span></p>
<p style="margin: 0mm 0mm 0pt;"><span lang="EN-US"><font face="MS UI Gothic" size="2"><font size="3">"Fatal error.</font></font></span></p>
<p style="margin: 0mm 0mm 0pt; text-indent: 10.5pt;"><span lang="EN-US"><font face="MS UI Gothic" size="2"><font size="3">Not enough ref_t and tau_t values."</font></font></span></p>
<p style="margin: 0mm 0mm 0pt; text-indent: 10.5pt;"><span lang="EN-US"><font face="MS UI Gothic" size="2"><font size="3">How should I change the conditions ?</font></font></span></p>
<p style="margin: 0mm 0mm 0pt;"><span lang="EN-US"><font face="MS UI Gothic" size="2"><font size="3"></font></font></span><font face="MS UI Gothic" size="2">&nbsp;</font></p>
<p style="margin: 0mm 0mm 0pt;"><span lang="EN-US"><font face="MS UI Gothic" size="2"><font size="3"><font face="MS UI Gothic">Many thanks for your help.</font>&nbsp;<font face="MS UI Gothic"> </font></font></font></span></p>

<p style="margin: 0mm 0mm 0pt;"><span lang="EN-US"><font face="MS UI Gothic" size="2"><font face="MS UI Gothic"></font></font></span><font face="MS UI Gothic" size="2">&nbsp;</font></p>
<p style="margin: 0mm 0mm 0pt;"><span lang="EN-US"><font face="MS UI Gothic" size="2"><font size="3">&nbsp;</font></font></span></p>
<p style="margin: 0mm 0mm 0pt;"><span lang="EN-US"><font face="MS UI Gothic" size="2"><font size="3">cpp<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>/lib/cpp<br>define<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span><br>integrator<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>md<br>tinit<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>0.0<br>dt<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>0.002<span>&nbsp;&nbsp;&nbsp; </span><br>nsteps<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>100000<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span><br>nstcomm<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>1<br>nstxout<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>50000<br>nstvout<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>50000<br>nstfout<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>0<br>nstlog<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>10000<br>nstenergy<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>250<br>nstxtcout<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; 
</span>10000<br>xtc_precision<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>=<span>&nbsp; </span>1000<br>xtc_grps<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>Protein<br>energygrps<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>Protein 
other<br>nstlist<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>5<br>ns_type<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>grid<br>pbc<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>xyz<br>rlist<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>0.9<br>coulombtype<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; 
</span>Reaction-Field<br>epsilon_r<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>78.0<br>rcoulomb<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>1.4<br>vdw-type<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>Cut-off<br>rvdw<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>1.4<br>tcoupl<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>Berendsen<br>tc-grps<span>&nbsp;&nbsp; </span><span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</span>=<span>&nbsp; </span>Protein other<br>ref_t<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>300 300<br>tau_t<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>0.1 0.1<br>pcoupl<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>Berendsen<br>tau_p<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; 
</span>1.0<br>compressibility<span>&nbsp;&nbsp;&nbsp;&nbsp; </span>=<span>&nbsp; </span>4.6e-5<br>ref_p<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>1.0<br>gen-vel<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>yes<br>gen-temp<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; </span>300<br>gen-seed<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; 
</span>$RANDOM<br>constraints<span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
</span>=<span>&nbsp; 
</span>all-bonds<br>&nbsp;</font></font></span></p></div>
<div><font face="MS UI Gothic" size="2">Hiromichi Tsurui MD, PhD<br>Dept. 
Pathology, <br>Juntendo 
University<br>Phone:81-3-5802-1039<br>Fax:81-3-3813-3164<br><a href="mailto:tsurui@med.juntendo.ac.jp" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">tsurui@med.juntendo.ac.jp</a></font></div></div>


</div><br>_______________________________________________<br>gmx-users mailing list &nbsp; &nbsp;<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">
http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please don't post (un)subscribe requests to the list. Use the<br>www interface or send it to <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users-request@gromacs.org">
gmx-users-request@gromacs.org</a>.<br>Can't post? Read <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php
</a><br><br></blockquote></div><br><br clear="all"><br>-- <br><br>Tsjerk A. Wassenaar, M.Sc.<br>Groningen Biomolecular Sciences and Biotechnology Institute (GBB)<br>Dept. of Biophysical Chemistry<br>University of Groningen
<br>Nijenborgh 4<br>9747AG Groningen, The Netherlands<br>+31 50 363 4336<br>