- Original Message -
From: Jack Shultz
Date: Friday, January 29, 2010 4:15
Subject: [gmx-users] including a custom itp file in topology
To: Discussion list for GROMACS users
Cc: Alan , Andrey Voronkov
> Hi,
>
> I was trying to figure out if there is a short-cut for what I'm doing. I have
> complexes that I'm trying to prep using pdb2gmx. The ligand does not have a
> standard residue name. The way I know this can work is seperating out the
> ligand and protein into seperate files and preping the ligand using acpypi
> and the protein using pdb2gmx. Then incorporating them into a single pdb
> complex and including a reference to the ligand.itp (generated by acpypi)
> into a complex topology file. Is there any shortcut to doing this? any way to
> reference the ligand's itp file when running pdb2gmx?
The procedure looks good, but there's no way to get pdb2gmx to do this step for
you. It will be straightforward to use sed or perl to add the #include and the
[molecules] changes into the .top produced by pdb2gmx - that's what these tools
are for!
Even
head -n -2 topol.top > new.top
echo '#include "ligand.itp"' >> new.top
tail -n 2 topol.top > footer >> new.top
echo 'ligand 1' >> new.top
will do the job if the last two lines of topol.top are of the form
[molecules]
Protein 1
This kind of UNIX-tool scripted problem-solving using head, tail, cut, paste,
grep, sed, uniq, sort, echo, perl, cat and others often requires some
consultation of the man pages and awareness that the tool exists... but most
problems you want to solve have been solved before and the tool is out there
already... Google away!
Mark
--
gmx-users mailing listgmx-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