Thanks again.<br><br>Since I am running MD simulation on a public cluster, I cannot modify atomtypes.atp and ffnonbonded.itp on that cluster.<br><br>But I could get every parameter ready and prepare the xxx.tpr files locally, and pass them to the cluster. So the practical question is, does the tpr file contain all those C12 and C6 LJ parameters already? Or mdrun need to retrieve these parameters according to the atom types from corresponding files within the gromos53a6.ff folder?<br>

<br>Regards,<br><br>Yun Shi<br><br><br>Yun Shi wrote:<br>
&gt; Hi Justin,<br>
&gt;<br>
&gt; Thanks a lot for the replies.<br>
&gt;<br>
&gt; I wonder what are the &quot;newer versions&quot; you indicated, but I find one as:<br>
&gt;<br>
&gt;<br>
&gt;   /A reoptimized GROMOS force field for hexopyranose-based carbohydrates<br>
&gt;   accounting for the relative free energies of ring conformers, anomers,<br>
&gt;   epimers, hydroxymethyl rotamers, and glycosidic linkage conformers/<br>
&gt;<br>
&gt; . This new version 56A_CARBO   is claimed to be &quot;nearly equivalent to<br>
&gt; 53A6 for non-carbohydrate systems&quot;, and I find, at least, atomic charges<br>
&gt; for C1 and C5 look more chemically reasonable.<br>
&gt;<br>
&gt; But changes for a couple of 1-4 van der Waals parameter C12 were made to<br>
&gt; 53a6 force field, and additional atom types were added. So could you<br>
&gt; tell me which file contain those LJ parameters with respect to different<br>
&gt; atom types?<br>
&gt;<br>
<br>
New atom types need to be added to atomtypes.atp and their parameters defined in<br>
ffnonbonded.itp.<br>
<br>
&gt; Alternatively, would it be reasonably good enough to use your<br>
&gt; Lysozyme-JZ4 strategy, i.e., to build initial topology of carbohydrate<br>
&gt; molecules by PRODRG, and then modify the atomic charges according to the<br>
&gt; 56A_CARBO paper?<br>
&gt;<br>
<br>
Maybe, but check the bonded parameters carefully.  Sugars are particularly<br>
sensitive to incorrect angle and dihedral parameters, especially.<br>
<br>
-Justin<br>
<br>
--<br>
==============================<div id=":1b5">==========<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu/" target="_blank">vt.edu</a> | <a href="tel:%28540%29%20231-9080" value="+15402319080">(540) 231-9080</a><br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================</div>