<div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Thank you Olivier for your response. <br></blockquote><div><br>SA <br>
</div><div><br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
------------------------------<br>
<br>
Message: 2<br>
Date: Thu, 28 Jul 2011 11:20:44 +0100<br>
From: Oliver Grant &lt;<a href="mailto:olivercgrant@gmail.com">olivercgrant@gmail.com</a>&gt;<br>
Subject: Re: [gmx-users] Conversion of the GLYCAM dihedral parameters<br>
To: Discussion list for GROMACS users &lt;<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>
Message-ID:<br>
        &lt;CAFd05fG5yAzPHr=<a href="mailto:sSfiCaA0FF0Q7UEt3ib_YJrs70qME3N3afg@mail.gmail.com">sSfiCaA0FF0Q7UEt3ib_YJrs70qME3N3afg@mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;windows-1252&quot;<br>
<br>
Hi SA,<br>
<br>
I think this is due to the original <a href="http://amb2gmx.pl" target="_blank">amb2gmx.pl</a> coder only considering AMBER<br>
type dihedrals. He/She didn&#39;t expect any negative values. You&#39;ll notice<br>
certain groups like NAc and COO- on sugars are affected when you run your<br>
simulation. Below is a similar post with a fix to the <a href="http://amb2gmx.pl" target="_blank">amb2gmx.pl</a> code. It&#39;s<br>
short so I pasted it rather than linking.<br>
<br>
You may already know this but note also the default is to correctly fudge<br>
AMBERS 1-4 interactions. So if you run a sugar using the GLYCAM forcefield<br>
they will be incorrectly scaled/fudged. An easy fix is to copy the<br>
ffamber99sb.itp (may be called something else now) to your local directory<br>
where you are running the md and change to:<br>
<br>
[ defaults ]<br>
; nbfunc        comb-rule       gen-pairs       fudgeLJ fudgeQQ<br>
1               2                         yes                     1.0<br>
1.0<br>
<br>
If you have a sugar and a protein you cannot mix the scaling. Mixed scaling<br>
is possible in amber11 if you really need it. I did all this when 4.0.7 was<br>
the latest release so things may have changed and it may now be possible to<br>
mix scaling or it may be that something else has changed :)<br>
<br>
If you find anything when testing the parameters please do post back.<br>
All the best,<br>
Oliver<br>
<br>
<br>
<br>
*Hi all,<br>
<br>
There is a problem that I encountered when I was trying to manually<br>
verify the proper dihedral conversion from AMBER topology to GROMACS<br>
topology using amb2gmx perl script.<br>
Some of the dihedrals were set to zero by amb2gmx even if in the<br>
prmtop file they were not zero. This was happening for all the lines<br>
that had PK with negative values!. Reading the script I came up with<br>
the following lines for V[i] calculation (lines 749  to 755 in the<br>
script file):<br>
<br>
...........................<br>
<br>
 # get all force constants for each line of a dihedral #<br>
 my $lines = $i -1 +$numijkl;<br>
 for(my $j=$i;$j&lt;=$lines;$j++){<br>
   my $period = abs($pn{$j});<br>
   if($pk{$j}&gt;0) {<br>
     $V[$period] = 2*$pk{$j}*$cal/$idivf{$j};<br>
   }<br>
<br>
...........................<br>
<br>
It seems from here that only PK values &gt; 0 are considered when<br>
computing the RB constants.<br>
After I change the &quot;&gt;&quot; sign to &quot;!=&quot; (i.e. not equal to) everything<br>
goes fine and ALL the dihedral are transformed correctly.<br>
<br>
While this is OK with different  AMBER sets and GAFF if one wish to<br>
convert a GLYCAM (which comes also with Amber package) generated<br>
topology, in the gromacs resultant file there will be missing<br>
parameters for dihedrals. That&#39;s because GLYCAM does not use phase<br>
shift and have also negative values for several PK<br>
<br>
Is this a bug or there is a reason for considering only the positive<br>
values of PK or I am missing something (as I am a begginer with<br>
AMBER)?<br>
(or maybe amb2gmx was designed only to work an AMBER ff conversion and<br>
not for example GLYCAM)<br>
<br>
Thanks for any comment,<br>
Andrei*<br>
<br>
<br>
On 21 July 2011 10:48, sa &lt;<a href="mailto:sagmx.mail@gmail.com">sagmx.mail@gmail.com</a>&gt; wrote:<br>
<br>
&gt; Yeah I have also found this page. Finally I did the conversion. I am<br>
&gt;&gt; currently testing the parameters.<br>
&gt;&gt;<br>
&gt;<br>
&gt; Thank you  Austin and Mark for your help<br>
&gt;<br>
&gt; SA-<br>
&gt;<br>
&gt;&gt;<br>
&gt;&gt; ------------------------------<br>
&gt;&gt;<br>
&gt;&gt; Message: 4<br>
&gt;&gt; Date: Wed, 20 Jul 2011 07:17:18 -0700 (PDT)<br>
&gt;&gt; From: &quot;Austin B. Yongye&quot; &lt;<a href="mailto:ybausty@yahoo.com">ybausty@yahoo.com</a>&gt;<br>
&gt;&gt; Subject: Re: [gmx-users] Conversion of the GLYCAM dihedral parameters<br>
&gt;&gt;        in the  GROMACS format.<br>
&gt;&gt; To: Discussion list for GROMACS users &lt;<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>
&gt;&gt; Message-ID:<br>
&gt;&gt;        &lt;<a href="mailto:1311171438.14853.YahooMailClassic@web161422.mail.bf1.yahoo.com">1311171438.14853.YahooMailClassic@web161422.mail.bf1.yahoo.com</a>&gt;<br>
&gt;&gt; Content-Type: text/plain; charset=&quot;utf-8&quot;<br>
&gt;&gt;<br>
&gt;&gt; &gt; So my questions are: how to convert it in the gromacs format ?<br>
&gt;&gt; I haven&#39;t done AMBER to GROMACS conversions. But a google search led me to<br>
&gt;&gt; this page. It&#39;s just a matter of knowing the different formats.<br>
&gt;&gt;<br>
&gt;&gt; <a href="http://ffamber.cnsm.csulb.edu" target="_blank">http://ffamber.cnsm.csulb.edu</a><br>
&gt;&gt;<br>
&gt;&gt; &gt; Is it correct to use the absolute value of the multiplicity in my<br>
&gt;&gt; parameters and use the negative values of barrier heights when they are<br>
&gt;&gt; exist in the AMBER parameters? Since, I have noticed that in the AMBER<br>
&gt;&gt; ff port in GROMACS, the multiplicity values are always set to &gt; 0<br>
&gt;&gt; and   the barrier heights have sometime a negative value.<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; So my questions are: how to convert it in the gromacs format ?  Is it<br>
&gt;&gt; correct to use the absolute value of the multiplicity in my parameters and<br>
&gt;&gt; use the negative values of barrier heights when they are exist in the AMBER<br>
&gt;&gt; parameters? Since, I have noticed that in the AMBER ff port in GROMACS, the<br>
&gt;&gt; multiplicity values are always set to &gt; 0 and   the barrier heights have<br>
&gt;&gt; sometime a negative value.<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; Thank you again for your advice.<br>
&gt;&gt;<br>
&gt;&gt; SA-<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; Message: 1<br>
&gt;&gt;<br>
&gt;&gt; Date: Tue, 19 Jul 2011 11:08:55 -0700 (PDT)<br>
&gt;&gt;<br>
&gt;&gt; From: &quot;Austin B. Yongye&quot; &lt;<a href="mailto:ybausty@yahoo.com">ybausty@yahoo.com</a>&gt;<br>
&gt;&gt;<br>
&gt;&gt; Subject: Re: [gmx-users] Conversion of the GLYCAM dihedral parameters<br>
&gt;&gt;<br>
&gt;&gt;         in the  GROMACS format.<br>
&gt;&gt;<br>
&gt;&gt; To: Discussion list for GROMACS users &lt;<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>
&gt;&gt;<br>
&gt;&gt; Message-ID:<br>
&gt;&gt;<br>
&gt;&gt;         &lt;<a href="mailto:1311098935.93905.YahooMailClassic@web161426.mail.bf1.yahoo.com">1311098935.93905.YahooMailClassic@web161426.mail.bf1.yahoo.com</a>&gt;<br>
&gt;&gt;<br>
&gt;&gt; Content-Type: text/plain; charset=&quot;utf-8&quot;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; It is normal to have combinations of negative and positive values for the<br>
&gt;&gt; barrier heights. Those are just the best coefficients to reproduce some QM<br>
&gt;&gt; rotational energy curve during the parameterization.  The negative<br>
&gt;&gt; periodicities are a convention from AMBER. They simply indicate that the<br>
&gt;&gt; dihedral angle potential has more than one term. For your example below:<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; O2-P<br>
&gt;&gt;<br>
&gt;&gt;           -OS-CP   1    0.10          0.0            -3.         Dimethyl<br>
&gt;&gt; phosphate<br>
&gt;&gt;<br>
&gt;&gt;                              1   -0.50          0.0            -2.<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;           1    0.10<br>
&gt;&gt;<br>
&gt;&gt;           0.0             1<br>
&gt;&gt;<br>
&gt;&gt; -3, -2 and 1 signify V3, V2 and V1 terms, respectively. Once a positive<br>
&gt;&gt; value is reached, terms for the O2-P-OS-CP potential have been completely<br>
&gt;&gt; accounted for.<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; Hope that helps.<br>
&gt;&gt;<br>
&gt;&gt; Austin-<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; --- On Tue, 7/19/11, Mark Abraham &lt;<a href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>&gt; wrote:<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; From: Mark Abraham &lt;<a href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>&gt;<br>
&gt;&gt;<br>
&gt;&gt; Subject: Re: [gmx-users] Conversion of the GLYCAM dihedral parameters in<br>
&gt;&gt; the GROMACS format.<br>
&gt;&gt;<br>
&gt;&gt; To: &quot;Discussion list for GROMACS users&quot; &lt;<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>&gt;<br>
&gt;&gt;<br>
&gt;&gt; Date: Tuesday, July 19, 2011, 7:08 AM<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;     On 19/07/2011 11:56 PM, sa wrote:<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;       Dear<br>
&gt;&gt;<br>
&gt;&gt;           GROMCS users,<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;       I<br>
&gt;&gt;<br>
&gt;&gt;           am trying to convert some GLYCAM parameters in GROMACS format.<br>
&gt;&gt;<br>
&gt;&gt;           For this<br>
&gt;&gt;<br>
&gt;&gt;           purpose, I am using the latest GLYCAM parameters downloaded<br>
&gt;&gt;<br>
&gt;&gt;           from the RJ. Woods’<br>
&gt;&gt;<br>
&gt;&gt;           website and the examples given in the acpype code (here for<br>
&gt;&gt;<br>
&gt;&gt;           the dihedral angles) :<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; <a href="http://code.google.com/p/acpype/source/browse/trunk/acpypi.py?spec=svn9&amp;r=9&amp;redir=1" target="_blank">http://code.google.com/p/acpype/source/browse/trunk/acpypi.py?spec=svn9&amp;r=9&amp;redir=1</a><br>

&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;       -------------------<br>
&gt;&gt;<br>
&gt;&gt; # dihedral    idivf        barrier hight/2 kcal/mol  phase degrees<br>
&gt;&gt; periodicity     comments<br>
&gt;&gt;<br>
&gt;&gt;        X -ca-ca-X    4           14.500*                     180.000<br>
&gt;&gt;          2.000             intrpol.bsd.on C6H6<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; * to convert to gmx: 14.5/4 * 4.184 * 2 (?) (yes in amb2gmx, no in<br>
&gt;&gt; topolbuild, why?) = 30.334 or 15.167 kJ/mol<br>
&gt;&gt;<br>
&gt;&gt; # X -CA-CA-X    4   14.50        180.0             2.<br>
&gt;&gt; intrpol.bsd.on C6H6 (from parm99.dat)<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; # X-CA-CA-X     3    30.33400     0.00000   -30.33400     0.00000<br>
&gt;&gt; 0.00000     0.00000   ; intrpol.bsd.on C6H6<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; -----------<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;       I have no problems with the parameters for proteins. But, in case of<br>
&gt;&gt; the GLYCAM parameters, I am a little confused<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; about the conversion of dihedral force constants (DFC), especially when<br>
&gt;&gt; the DFC and the periodicity values are &lt; 0 for example<br>
&gt;&gt;<br>
&gt;&gt; for this torsion:<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;       O2-P<br>
&gt;&gt;<br>
&gt;&gt;           -OS-CP   1    0.10          0.0            -3.         Dimethyl<br>
&gt;&gt; phosphate<br>
&gt;&gt;<br>
&gt;&gt;                     1   -0.50          0.0            -2.<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;           1    0.10<br>
&gt;&gt;<br>
&gt;&gt;           0.0             1<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;     Where only a positive value makes sense, sometimes people use<br>
&gt;&gt;<br>
&gt;&gt;     negative values to indicate some special functional form. This can<br>
&gt;&gt;<br>
&gt;&gt;     be easier to code. Regardless, you&#39;ll have to check out the GLYCAM<br>
&gt;&gt;<br>
&gt;&gt;     documentation and see what is meant, before you can address how to<br>
&gt;&gt;<br>
&gt;&gt;     convert it into a GROMACS format. Obviously the contents of parts of<br>
&gt;&gt;<br>
&gt;&gt;     chapter 4 and 5 of the manual will be important.<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;     Mark<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt; -----Inline Attachment Follows-----<br>
&gt;&gt;<br>
&gt;&gt; --<br>
&gt;&gt; gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt;&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt;&gt; Please search the archive at<br>
&gt;&gt; <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
&gt;&gt; Please don&#39;t post (un)subscribe requests to the list. Use the<br>
&gt;&gt; www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
&gt;&gt; Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
&gt;&gt; -------------- next part --------------<br>
&gt;&gt; An HTML attachment was scrubbed...<br>
&gt;&gt; URL:<br>
&gt;&gt; <a href="http://lists.gromacs.org/pipermail/gmx-users/attachments/20110720/c963e049/attachment.html" target="_blank">http://lists.gromacs.org/pipermail/gmx-users/attachments/20110720/c963e049/attachment.html</a><br>

&gt;&gt;<br>
&gt;&gt; ------------------------------<br>
&gt;&gt;<br>
&gt;&gt; --<br>
&gt;&gt; gmx-users mailing list<br>
&gt;&gt; <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt;&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt;&gt; Please search the archive at<br>
&gt;&gt; <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
&gt;&gt;<br>
&gt;&gt; End of gmx-users Digest, Vol 87, Issue 109<br>
&gt;&gt; ******************************************<br>
&gt;&gt;<br>
&gt;<br>
&gt;<br>
&gt; --<br>
&gt; gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
&gt; <a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt; Please search the archive at<br>
&gt; <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
&gt; Please don&#39;t post (un)subscribe requests to the list. Use the<br>
&gt; www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>
&gt; Can&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
&gt;<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://lists.gromacs.org/pipermail/gmx-users/attachments/20110728/0d95a512/attachment.html" target="_blank">http://lists.gromacs.org/pipermail/gmx-users/attachments/20110728/0d95a512/attachment.html</a><br>

<br>
------------------------------<br>
<font color="#888888"><br>
--<br>
gmx-users mailing list<br>
<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
<br>
End of gmx-users Digest, Vol 87, Issue 161<br>
******************************************<br>
</font></blockquote></div><br>