Hi Justin,<br><br>Thanks a lot for the replies.<br><br>I wonder what are the &quot;newer versions&quot; you indicated, but I find one as:<br><h1 style="font-weight: normal;" class="articleTitle"><i><font size="2">A reoptimized GROMOS force field for 
hexopyranose-based carbohydrates accounting for the relative free 
energies of ring conformers, anomers, epimers, hydroxymethyl rotamers, 
and glycosidic linkage conformers</font></i></h1>. This new version 56A<sub>CARBO</sub>  is claimed to be &quot;nearly equivalent to 53A6 for non-carbohydrate systems&quot;, and I find, at least, atomic charges for C1 and C5 look more chemically reasonable.<br>

<br>But changes for a couple of 1-4 van der Waals parameter C12 were made to 53a6 force field, and additional atom types were added. So could you tell me which file contain those LJ parameters with respect to different atom types?<br>

<br>Alternatively, would it be reasonably good enough to use your Lysozyme-JZ4 strategy, i.e., to build initial topology of carbohydrate molecules by PRODRG, and then modify the atomic charges according to the 56A<sub>CARBO </sub>paper?<br>

<br>Thanks,<br><br>Yun Shi <br><br><br><br>Yun Shi wrote:<br>
&gt; Hi all,<br>
&gt;<br>
&gt; I am doing MD simulation of some carbohydrate-protein complex with this<br>
&gt; 53a6 force-field.<br>
&gt;<br>
&gt; I noted that in any oligosaccharide, the charge assigned for anomeric<br>
&gt; carbon is 0.232 while C5 (i.e. in a non-terminal mannopyranose) is<br>
&gt; 0.376. This is kind of counter-intuitive for me since the anomeric<br>
&gt; carbon C1 is linked to two oxygen atoms O5 and &quot;O1&quot; while C5 is linked<br>
&gt; to only one oxygen atom O5, and both C1 and C5 are CH type atom. And in<br>
&gt; other force filed like GLYCAM06, C1 always has a larger positive partial<br>
&gt; charge than C5.<br>
&gt;<br>
&gt; I understand that GROMOS force fields are derived in quite different<br>
&gt; ways, but I still wonder if it makes sense to switch the charges on C1<br>
&gt; and C5?<br>
&gt;<br>
<br>
No, but the charges used may not be the best ones anyway.  There are newer<br>
revisions of 53A6 that are much better.<br>
<br>
&gt; It seems these carbohydrate parameters did not change in the new 54a7<br>
&gt; ff, which was tested against nOe and J-coupling data rather than free<br>
&gt; enthalpy of hydration and solvation, which really puzzles me.<br>
&gt;<br>
<br>
54A7 was principally a revision of protein atom types and one lipid type;<br>
nothing in the sugar parameters was modified.<br>
<br>
-Justin<br>
<br>
--<br>
==============================<div id=":d">==========<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>