Liu Shiyong wrote:
Hi,
Chain name is very important for some application. But, when I use
the following script to get a minimized structure from a multi-chain
complex.
But the chain name was changed. I checked the man trjconv , but I did
not find any information on chain name for trjconv.
The .gro file format does not permit chain identifiers. So they were
lost when you wrote ${file}.gro, and re-invented for output.pdb. If I
remember correctly they're also not stored in a .tpr file.
pdb2gmx -ff ${forcefield} -f starting.pdb -p ${file}.top -i
${file}.posre.itp -o ${file}.gro > ${file}.output.pdb2gmx 2>&1
mdrun -nice 0 -v -s ${file}.input.tpr -o ${file}.minim_traj.trr -c
${file}.minimized.gro -e ${file}.minim_ener.edr -g ${file}.emlog.log
> ${file}.output.mdrun 2>&1
trjconv -sep -f ${file}.minim_traj.trr -s ${file}.input.tpr -o
output.pdb << HHH
2
2
HHH
The starting.pdb has four chain: A B D E
The output file output.pdb has a different chain name: A B C D.
There's no way to fix this one with GROMACS (at least, as far as I
know). However if you edit your starting.pdb file to have the same kind
of format as output.pdb, then you can use the standard unix utility
programs head, tail, grep, cut and paste to solve this problem. Use
head -n x output.pdb > head
tail -n y output.pdb > tail
to extract the header and footer to separate files, and then grep out
all of the lines that begin with atom records
grep -e '^ATOM' output.pdb > outputatom
grep -e '^ATOM' starting.pdb > startingatom
Then
cut -c=z startingatom > chainidentifiers
cut -c-1-a outputatom > left
cut -c=b- outputatom > right
paste -d='' left chainidentifiers right > middle
cat head middle tail > fixed.pdb
and to clean up
rm head tail outputatom startingatom chainidentifiers middle
You will need to identify the correct numbers for x, y, and z, and a=z-1
and b=z+1. See the man pages for each of these commands to learn more.
Mark
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.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 [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php