Alan schrieb:
> Dear Vedat,
>
> On Wed, May 19, 2010 at 20:36, <gmx-users-requ...@gromacs.org
> <mailto:gmx-users-requ...@gromacs.org>> wrote:
>
> @rui
> acpypi -i ch3cn_210.pdb
> says: "cannot find template for residue C3N in our library". and
> indeed,
> there's no residue C3N in my ffamber99sb.rtp file
>
> (and i don't know, how to use it in order to generate my topology
> or even
> an rtp file?!)
>
>
> Have a bit of patience and try to read a bit more about ACPYPE.
>
> 1) Read the info there, I am sure you'll find it useful;
> 2) You'll learn that you need to have just one molecule in a pdb and
> not the whole box if you want the topologi of C3N.
> 3) It took me 2s to get the topology with acpype but months to write
> the code, so if you'd take some few minutes to read and use an updated
> version (it's not acpypi anymore BTW)...
>
thanks for your helpful hints. i updated acpype, created a pdb file with
a single molecule and ran
acpype -i ch3cn_210_single.pdb
which generated an .itp and other interesting files. that's
nice. (remember, i want to use gromacs with amber99sb force field and i
downloaded 3 files from the amber site: ch3cn_210.pdb,
frcmod.ch3cn,prep.ch3cn.have you ever seen their content?)
1) the charges do not match the ones listed in the prep.ch3cn file.
shall i just change them by hand accordingly?
2) dummy atoms as listed in the prep.ch3cn are not present in the new
.itp file.
3) the force constants seem totally different. shall i again just adjust
them to the original file obtained from the amber site?
is there another way of using acpype, with a proper args list, that i
should use in this situation?
i know, that's many questions, but is has to be done!
> BTW, how did you get this message "cannot find template for residue
> C3N in our library"?
i got that message *within* the following output when running:
>acpype -i ch3cn_210.pdb
[...]
Warning: cannot find template for residue C3N in our library.
You will not be able to save prmtop for this molecule.
Warning: cannot find template for residue C3N in our library.
You will not be able to save prmtop for this molecule.
[gtkleap]$ #check C3N
[gtkleap]$ saveamberparm C3N ch3cn_210_AC.prmtop ch3cn_210_AC.inpcrd
Error: dparm pchg does not exist!
++++++++++end_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
ERROR: Sleap failed
==> Removing temporary files...
ACPYPE FAILED: [Errno 2] No such file or directory: 'ch3cn_210_AC.inpcrd'
Total time of execution: 7s
>
> acpype.googlecode.com <http://acpype.googlecode.com>
>
> Regards,
> Alan
>
> --
> Alan Wilter S. da Silva, D.Sc. - CCPN Research Associate
> Department of Biochemistry, University of Cambridge.
> 80 Tennis Court Road, Cambridge CB2 1GA, UK.
> >>http://www.bio.cam.ac.uk/~awd28 <http://www.bio.cam.ac.uk/%7Eawd28><<
--
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
