<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><p class="yiv665091845MsoNormal">Hi everyone, <br>
<br>
I am using g53a6 force filed for the simulation of cellulose. However it seems
that parametrization works only for the sugar blocked linked together and no
terminal atoms for terminal residues are defined there (I mean the start and
end atom, H and O atom at the beginning linked to c4 and H atom at the last
Sugar residue linked to O1). I tried to solve the problem by making the
topology file removing the start OH and end H group from the initial PDB file to&nbsp;
have the topology file made by g53a6 and later add the parameters of OH and H
to the topology file. The problem is at the time of minimization almost no
change in energy (not converged even after 7325 step) with the maximum force
(about 9e+03) on the H of new added OH group (the start OH group). and when I
checked the output of minimization step I found the HO position and C2 and C3,
O3 atoms are deformed and the defined angles/bonds/dihedral of this new OH are
not kept (the other H group the the other side, end, is fine). I know the
problem is due to the way I defined the parameters but I do not know what the
problem exactly is. Since the original topology file without the terminal atoms
start from 1, the terminal atoms are added at the end of gro file and topology
file (although OH group is part of first residue): <br>
<br>
the molecule is a cellotetetraose and here is the added part to top file (Atom
number: H at end: 57 ,&nbsp; O at start:58, &nbsp; H at start:59):<br>
<br>
<span style="font-family:courier, monaco, monospace, sans-serif;">; residue 1001 GLCB rtp GLCB q&nbsp;
0.0<br>
&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
CH1&nbsp;&nbsp; 1001&nbsp;&nbsp; GLCB&nbsp;&nbsp;&nbsp;&nbsp;
C4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.332&nbsp;&nbsp;&nbsp;&nbsp; 13.019&nbsp;&nbsp; ; qtot 0.232<br>
&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
CH1&nbsp;&nbsp; 1001&nbsp;&nbsp; GLCB&nbsp;&nbsp;&nbsp;&nbsp;
C3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.232&nbsp;&nbsp;&nbsp;&nbsp; 13.019&nbsp;&nbsp; ; qtot 0.464<br>
&nbsp;&nbsp;&nbsp;&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
OA&nbsp;&nbsp; 1001&nbsp;&nbsp; GLCB&nbsp;&nbsp;&nbsp;&nbsp;
O3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;
-0.642&nbsp;&nbsp;&nbsp; 15.9994&nbsp;&nbsp; ; qtot -0.178<br>
&nbsp;&nbsp;&nbsp;&nbsp;
4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;
1001&nbsp;&nbsp; GLCB&nbsp;&nbsp;&nbsp; HO3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.41&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
1.008&nbsp;&nbsp; ; qtot 0.232</span><br style="font-family:courier, monaco, monospace, sans-serif;"><span style="font-family:courier, monaco, monospace, sans-serif;">
&nbsp; &nbsp; &nbsp; &nbsp; .</span><br style="font-family:courier, monaco, monospace, sans-serif;"><span style="font-family:courier, monaco, monospace, sans-serif;">
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; .</span><br style="font-family:courier, monaco, monospace, sans-serif;"><span style="font-family:courier, monaco, monospace, sans-serif;">
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; .</span><br style="font-family:courier, monaco, monospace, sans-serif;">
<span style="font-family:courier, monaco, monospace, sans-serif;">&nbsp;&nbsp;&nbsp;
56&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; OA&nbsp;&nbsp;
1004&nbsp;&nbsp; GLCB&nbsp;&nbsp;&nbsp;&nbsp; O1&nbsp;&nbsp;&nbsp;&nbsp;
20&nbsp;&nbsp;&nbsp;&nbsp; -0.642&nbsp;&nbsp;&nbsp; 15.9994&nbsp;&nbsp; ; qtot
-0.20<br>
&nbsp;&nbsp;&nbsp; 57&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
H&nbsp;&nbsp; 1004&nbsp;&nbsp; GLCB&nbsp;&nbsp;&nbsp;
HO1&nbsp;&nbsp;&nbsp;&nbsp; 20&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.282&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.008&nbsp;&nbsp; ; qtot 0<br>
;added later to 1000 GLCB<br>
&nbsp;&nbsp;&nbsp; 58&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
OA&nbsp;&nbsp; 1001&nbsp;&nbsp; GLCB&nbsp;&nbsp;&nbsp;&nbsp;
O4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;
-0.542&nbsp;&nbsp;&nbsp; 15.9994&nbsp;&nbsp; ;<br>
&nbsp;&nbsp;&nbsp; 59&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
H&nbsp;&nbsp; 1001&nbsp;&nbsp; GLCB&nbsp;&nbsp;&nbsp;
HO4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.442&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.008&nbsp;&nbsp; ;<br>
<br>
[ bonds ]<br>
;&nbsp; ai&nbsp;&nbsp;&nbsp; aj
funct&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; c2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c3<br>
&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;
2&nbsp;&nbsp;&nbsp; gb_26<br>
&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; 11&nbsp;&nbsp;&nbsp;&nbsp;
2&nbsp;&nbsp;&nbsp; gb_26<br>
.<br>
.<br>
.<br>
&nbsp; 55&nbsp;&nbsp;&nbsp; 56&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;
gb_20<br>
&nbsp;&nbsp; 56&nbsp;&nbsp;&nbsp; 57&nbsp;&nbsp;&nbsp;&nbsp;
2&nbsp;&nbsp;&nbsp; gb_1<br>
&nbsp;&nbsp; 58&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;
2&nbsp;&nbsp;&nbsp; gb_20<br>
&nbsp;&nbsp; 58&nbsp;&nbsp;&nbsp; 59&nbsp;&nbsp;&nbsp;&nbsp;
2&nbsp;&nbsp;&nbsp; gb_1<br>
<br>
[ angles ]<br>
;&nbsp; ai&nbsp;&nbsp;&nbsp; aj&nbsp;&nbsp;&nbsp; ak
funct&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; c3<br>
&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;
11&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; ga_8<br>
&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;
3&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; ga_9<br>
&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;
5&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; ga_8<br>
.<br>
.<br>
.<br>
&nbsp;&nbsp; 55&nbsp;&nbsp;&nbsp; 56&nbsp;&nbsp;&nbsp;
57&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; ga_12<br>
&nbsp;&nbsp; 11&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;
58&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; ga_9<br>
&nbsp;&nbsp; 58&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;
2&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; ga_9<br>
&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; 58&nbsp;&nbsp;&nbsp;
59&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; ga_12<br>
<br>
[ dihedrals ]<br>
;&nbsp; ai&nbsp;&nbsp;&nbsp; aj&nbsp;&nbsp;&nbsp; ak&nbsp;&nbsp;&nbsp; al
funct&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; c5<br>
&nbsp;&nbsp; 11&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;
2&nbsp;&nbsp;&nbsp;&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; gd_17<br>
&nbsp;&nbsp; 11&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;
2&nbsp;&nbsp;&nbsp;&nbsp; 5&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; gd_34<br>
&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;
11&nbsp;&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; gd_34<br>
.<br>
.<br>
.<br>
&nbsp;&nbsp; 59&nbsp;&nbsp;&nbsp; 58&nbsp;&nbsp;&nbsp;&nbsp;
1&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; gd_30<br>
&nbsp;&nbsp; 58&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;
2&nbsp;&nbsp;&nbsp;&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; gd_18<br>
&nbsp;&nbsp; 58&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;
2&nbsp;&nbsp;&nbsp;&nbsp; 5&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; gd_17<br>
&nbsp;&nbsp; 58&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;
11&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; gd_17<br>
<br>
[ dihedrals ]<br>
;&nbsp; ai&nbsp;&nbsp;&nbsp; aj&nbsp;&nbsp;&nbsp; ak&nbsp;&nbsp;&nbsp; al
funct&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
c2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; c3<br>
.<br>
.<br>
.<br>
&nbsp;&nbsp; 11&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;
58&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; gi_2</span><br style="font-family:courier, monaco, monospace, sans-serif;">
<br>
<br>
I prepare this parameters based on the other OH group in the sugars and
parameters of link oxygen (the shared oxygen at the position of glycosidic bond
that is expected to be O4 group of the next residue). <a rel="nofollow" name="_MailAutoSig"></a></p>

<p class="yiv665091845MsoNormal" style=""><span style="">&nbsp;</span></p>

<p class="yiv665091845MsoNormal" style="margin-bottom:12.0pt;"><span style="">I
did not use the PRODRG since i need 53a6 and I found for carbohydrate there is
too much different between the parameters in 53a6 and 43a1. Somewhere it ahs
been mentioned to change the terminal database file (.tdb), however I think it
may make the problem more complicated. <br>
<br>
Your help is appreciated.<br>
-Somaye</span></p>

<p class="yiv665091845MsoNormal" style=""><span style=""><i style=""><span style="">&nbsp;</span></i></span></p>

<p class="yiv665091845MsoNormal" style=""><span style=""><i style=""><span style="">...................................................................</span></i></span></p>

<p class="yiv665091845MsoNormal" style=""><span style=""><b style=""><i style=""><span style="color:#943634;">Somayesadat Badieyan</span></i></b></span></p>

<p class="yiv665091845MsoNormal" style=""><span style=""><b style=""><i style=""><span style="color:#E36C0A;">PhD Candidate and Research Assistant </span></i></b></span></p>

<p class="yiv665091845MsoNormal" style=""><span style=""><b style=""><i style=""><span style="color:#943634;">Biological Syatems Engineering</span></i></b></span></p>

<p class="yiv665091845MsoNormal" style=""><span style=""><b style=""><i style=""><span style="color:#E36C0A;">201 Seitz Hall, Virginia Tech</span></i></b></span></p>

<p class="yiv665091845MsoNormal" style=""><span style=""><b><i style=""><span style="color:#943634;">Blacksburg,
VA 24060</span></i></b></span></p>

<p class="yiv665091845MsoNormal" style=""><span style=""><i style=""><span style="">...................................................................</span></i></span></p></td></tr></table>