Nice, glad you did progress. See below.<br><br><div class="gmail_quote">On Thu, May 20, 2010 at 12:38, <span dir="ltr"><<a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div id=":45f" class="ii gt">thanks for your helpful hints. i updated acpype, created a pdb file with<br>
a single molecule and ran<br>
<br>
acpype -i ch3cn_210_single.pdb<br>
<br>
which generated an .itp and other interesting files. that's<br>
nice. (remember, i want to use gromacs with amber99sb force field and i<br>
downloaded 3 files from the amber site: ch3cn_210.pdb,<br>
frcmod.ch3cn,prep.ch3cn.have you ever seen their content?)<br>
<br>
1) the charges do not match the ones listed in the prep.ch3cn file.<br>
shall i just change them by hand accordingly?<br></div></blockquote><div><br></div><div>It doesn't match because it's using am1bcc, which was parametrised to reproduce the RESP charges, but obviously (sqm is semi-empirical method, not like gaussian) it won't be accurate.</div>
<div><br></div><div>However, you're right, if you have the RESP charges in prep.ch3cn just copy them by hand accordingly.</div><div><br></div><div>Or even better, if you want to learn more about the whole stuff, double check if the parameters you got from the Manchester site are OK, why not trying <a href="http://q4md-forcefieldtools.org/RED/">q4md-forcefieldtools.org/RED/</a>? Once you got the charges (they should be very close if not the same from prep.ch3cn), you can use acpype just to generate the topology by providing a c3n.MOL2 file with the charges calculated by RED and then using "acpype -di c3n.mol2 -c user".</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div id=":45f" class="ii gt">
2) dummy atoms as listed in the prep.ch3cn are not present in the new<br>
.itp file.<br></div></blockquote><div><br></div><div>I guess you don't know how a prep file works, so see <a href="http://ambermd.org/doc/prep.html">http://ambermd.org/doc/prep.html</a>.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":45f" class="ii gt">
3) the force constants seem totally different. shall i again just adjust<br>
them to the original file obtained from the amber site?<br></div></blockquote><div><br></div><div>If using acpype with default mode, so you'd get GAFF parameters. You may want to try:</div><div><br></div><div>acpype -di c3n.mol2 -c user -a amber</div>
<div><br></div><div>However, it still may diff. If you read Jaime's paper and you agree with what he did, so you can "copy&paste" his parameters as well.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":45f" class="ii gt">
is there another way of using acpype, with a proper args list, that i<br>
should use in this situation?<br></div></blockquote><div><br></div><div>Read the Wikis in the acpype site and 'acpype -h'. I am always keen for suggestions.</div><div><br></div><div>Another possible way, would be using tleap from AmberTools, generate just one molecule, save parameters and use acpype to convert from amber to gromacs, something like</div>
<div><br></div><div>acpype -p c3n.prmtop -x c3n.inpcrd</div><div><br></div><div>If doing so, you'd get the exactly Jaime's topology but in gromacs format (gro and top file, not itp, so you may need to adjust things in the top file in order to create a itp, but should be a simple task).</div>
<div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div id=":45f" class="ii gt">> BTW, how did you get this message "cannot find template for residue<br>
> C3N in our library"?<br>
i got that message *within* the following output when running:<br>
>acpype -i ch3cn_210.pdb<br>
[...]<br>
Warning: cannot find template for residue C3N in our library.<br>
You will not be able to save prmtop for this molecule.<br>
Warning: cannot find template for residue C3N in our library.<br>
You will not be able to save prmtop for this molecule.<br>
[gtkleap]$ #check C3N<br>
[gtkleap]$ saveamberparm C3N ch3cn_210_AC.prmtop ch3cn_210_AC.inpcrd<br>
Error: dparm pchg does not exist!<br>
<br>
++++++++++end_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++<br>
ERROR: Sleap failed<br>
==> Removing temporary files...<br>
ACPYPE FAILED: [Errno 2] No such file or directory: 'ch3cn_210_AC.inpcrd'<br>
Total time of execution: 7s<br></div></blockquote></div><div><br></div>Ah, ok, I should've know this... It's a fall back routine to try to use 'sleap', but sleap is broken in AmberTools 1.3 and 1.4, unfortunately.<br>
<br><div>Thanks for trying acpype.</div><div><br></div><div>Cheers,</div><div><br></div><div>Alan<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>>><a href="http://www.bio.cam.ac.uk/~awd28">http://www.bio.cam.ac.uk/~awd28</a><<<br>
</div>