Dear Vedat,<div><br></div><div>Thank you very much indeed. Things only get better with constructive feedbacks, either positive or negative.</div><div><br></div><div>Glad you&#39;re managing well your problem. However, as I see, your problem now is in the tleap level.</div>

<div><br></div><div>So, are you in amber mailing list? If so I would suggest to post your question there. If not, let me know and I try to post the question there (you may want to consider joining it, very useful, even for those who uses only ambertools).</div>

<div><br></div><div>However, at first, I would go ahead and neglect the error. In the end, if you not happy, you can always adjust this particular parameter (in the itp file), since in principle a triple bond would mean a stronger bond, so implying in higher K-string constant. For how much to adjust? Well, you cannot even ever say that you got is &quot;correct&quot;; all you can get is a better or worse model and in the end it&#39;s up to you to set that.</div>

<div><br></div><div>Nevertheless, be sure that you have the right protonation state and try to re-run your case putting CONECT info in the pdb just to check it.</div><div><br></div><div>I hope it helps.</div><div><br></div>

<div>Best,</div><div><br></div><div>Alan</div><div><br><div class="gmail_quote">On Fri, Jun 4, 2010 at 11:00,  <span dir="ltr">&lt;<a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>&gt;</span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div id=":vm" class="ii gt">hi alan and all others,<br>
<br>
i love it. i&#39;ve probably never seen such a helpful and constuctive<br>
answer to a question asked over the gmx user mailing list. following<br>
your instructions was very informative on the one hand and provided the<br>
results that i needed on the other.<br>
<br>
slightly modifying the the ch3cn.frcmod and -.prep files from the<br>
manchester site in accordance with a newer solvent model for<br>
acetonitrile of the year 2007 and feeding tleap with that data resulted<br>
in a ch3cn.lib file and afterwards in ch3cn.prmtop and ch3cn.inpcrd<br>
which i could use as input for acpype in order to generate .gro and .itp<br>
files. perfect.<br>
<br>
almost perfect, because in tleap i was not able to add the YC-YN triple<br>
bond (with the &quot;bond ...&quot; command) without generating an error message<br>
when saving my parameter files (&quot;saveAmberParm ...&quot; command). see the<br>
following output:<br>
<br>
<br>
&gt; bond C3N.1.C2 C3N.1.N1 &quot;T&quot;<br>
&gt; saveAmberParm C3N ch3cn.prmtop ch3cn.inpcrd<br>
Checking Unit.<br>
Building topology.<br>
Building atom parameters.<br>
Building bond parameters.<br>
Building angle parameters.<br>
Building proper torsion parameters.<br>
!FATAL ERROR----------------------------------------<br>
!FATAL:    In file [unitio.c], line 1778<br>
!FATAL:    Message: 1-4: cannot add bond 1 2<br>
This may be caused by duplicate bond specifications;<br>
for example, explicit bond commands in addition to PDB conect records.<br>
!<br>
!ABORTING.<br>
<br>
where could the &quot;second&quot; bond specification be, that avoids the<br>
specification of a triple bond &quot;T&quot;? there&#39;s no bonding information in<br>
the pdb file.<br>
<br>
a step earlier in tleap, i had to addAtomTypes YN and YC, however, there<br>
was no way to specify them with &quot;sp&quot; hybridization, but &quot;sp2&quot; resulting<br>
in two lines each containing a single bond between these two atoms<br>
within the &quot;!entry.C3N.unit.connectivity table&quot; section of the ch3cn.lib<br>
file (seems to be a double bond).<br>
<br>
is this issue negligible, since the force constants and equilibrium<br>
angles (180 degrees) are already describing the right molecule?<br>
<br>
vedat<br>
<br>
<br>
<br>
<br>
Alan schrieb:<br>
&gt; Nice, glad you did progress. See below.<br>
&gt;<br>
&gt; On Thu, May 20, 2010 at 12:38, &lt;<a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a><br>
&gt; &lt;mailto:<a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>&gt;&gt; wrote:<br>
&gt;<br>
&gt;     thanks for your helpful hints. i updated acpype, created a pdb<br>
&gt;     file with<br>
&gt;     a single molecule and ran<br>
&gt;<br>
&gt;     acpype -i ch3cn_210_single.pdb<br>
&gt;<br>
&gt;     which generated an .itp and other interesting files. that&#39;s<br>
&gt;     nice. (remember, i want to use gromacs with amber99sb force field<br>
&gt;     and i<br>
&gt;     downloaded 3 files from the amber site: ch3cn_210.pdb,<br>
&gt;     frcmod.ch3cn,prep.ch3cn.have you ever seen their content?)<br>
&gt;<br>
&gt;     1) the charges do not match the ones listed in the prep.ch3cn file.<br>
&gt;     shall i just change them by hand accordingly?<br>
&gt;<br>
&gt;<br>
&gt; It doesn&#39;t match because it&#39;s using am1bcc, which was parametrised to<br>
&gt; reproduce the RESP charges, but obviously (sqm is semi-empirical<br>
&gt; method, not like gaussian)  it won&#39;t be accurate.<br>
&gt;<br>
&gt; However, you&#39;re right, if you have the RESP charges in prep.ch3cn just<br>
&gt; copy them by hand accordingly.<br>
&gt;<br>
&gt; Or even better, if you want to learn more about the whole stuff,<br>
&gt; double check if the parameters you got from the Manchester site are<br>
&gt; OK, why not trying <a href="http://q4md-forcefieldtools.org/RED/" target="_blank">q4md-forcefieldtools.org/RED/</a><br>
&gt; &lt;<a href="http://q4md-forcefieldtools.org/RED/" target="_blank">http://q4md-forcefieldtools.org/RED/</a>&gt;? Once you got the charges (they<br>
&gt; should be very close if not the same from  prep.ch3cn), you can use<br>
&gt; acpype just to generate the topology by providing a c3n.MOL2 file with<br>
&gt; the charges calculated by RED and then using &quot;acpype -di c3n.mol2 -c<br>
&gt; user&quot;.<br>
&gt;<br>
&gt;<br>
&gt;     2) dummy atoms as listed in the prep.ch3cn are not present in the new<br>
&gt;     .itp file.<br>
&gt;<br>
&gt;<br>
&gt; I guess you don&#39;t know how a prep file works, so<br>
&gt; see <a href="http://ambermd.org/doc/prep.html" target="_blank">http://ambermd.org/doc/prep.html</a>.<br>
&gt;<br>
&gt;<br>
&gt;     3) the force constants seem totally different. shall i again just<br>
&gt;     adjust<br>
&gt;     them to the original file obtained from the amber site?<br>
&gt;<br>
&gt;<br>
&gt; If using acpype with default mode, so you&#39;d get GAFF parameters. You<br>
&gt; may want to try:<br>
&gt;<br>
&gt; acpype -di c3n.mol2 -c user -a amber<br>
&gt;<br>
&gt; However, it still may diff. If you read Jaime&#39;s paper and you agree<br>
&gt; with what he did, so you can &quot;copy&amp;paste&quot; his parameters as well.<br>
&gt;<br>
&gt;<br>
&gt;     is there another way of using acpype, with a proper args list, that i<br>
&gt;     should use in this situation?<br>
&gt;<br>
&gt;<br>
&gt; Read the Wikis in the acpype site and &#39;acpype -h&#39;. I am always keen<br>
&gt; for suggestions.<br>
&gt;<br>
&gt; Another possible way, would be using tleap from AmberTools, generate<br>
&gt; just one molecule, save parameters and use acpype to convert from<br>
&gt; amber to gromacs, something like<br>
&gt;<br>
&gt; acpype -p c3n.prmtop -x c3n.inpcrd<br>
&gt;<br>
&gt; If doing so, you&#39;d get the exactly Jaime&#39;s topology but in gromacs<br>
&gt; format (gro and top file, not itp, so you may need to adjust things in<br>
&gt; the top file in order to create a itp, but should be a simple task).<br>
&gt;<br>
&gt;     &gt; BTW, how did you get this message &quot;cannot find template for residue<br>
&gt;     &gt; C3N in our library&quot;?<br>
&gt;     i got that message *within* the following output when running:<br>
&gt;     &gt;acpype -i ch3cn_210.pdb<br>
&gt;     [...]<br>
&gt;     Warning: cannot find template for residue C3N in our library.<br>
&gt;             You will not be able to save prmtop for this molecule.<br>
&gt;     Warning: cannot find template for residue C3N in our library.<br>
&gt;             You will not be able to save prmtop for this molecule.<br>
&gt;     [gtkleap]$ #check C3N<br>
&gt;     [gtkleap]$ saveamberparm C3N ch3cn_210_AC.prmtop ch3cn_210_AC.inpcrd<br>
&gt;     Error: dparm pchg does not exist!<br>
&gt;<br>
&gt;     ++++++++++end_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++<br>
&gt;     ERROR: Sleap failed<br>
&gt;     ==&gt; Removing temporary files...<br>
&gt;     ACPYPE FAILED: [Errno 2] No such file or directory:<br>
&gt;     &#39;ch3cn_210_AC.inpcrd&#39;<br>
&gt;     Total time of execution: 7s<br>
&gt;<br>
&gt;<br>
&gt; Ah, ok, I should&#39;ve know this... It&#39;s a fall back routine to try to<br>
&gt; use &#39;sleap&#39;, but sleap is broken in AmberTools 1.3 and 1.4, unfortunately.<br>
&gt;<br>
&gt; Thanks for trying acpype.<br>
&gt;<br>
&gt; Cheers,<br>
&gt;<br>
&gt; Alan<br></div></blockquote></div><br><br clear="all"><br>-- <br>Alan Wilter S. da Silva, D.Sc. - CCPN Research Associate<br>Department of Biochemistry, University of Cambridge. <br>80 Tennis Court Road, Cambridge CB2 1GA, UK.<br>

&gt;&gt;<a href="http://www.bio.cam.ac.uk/~awd28">http://www.bio.cam.ac.uk/~awd28</a>&lt;&lt;<br>
</div>