Difference between revisions of "Pdbset (ccp4)"
From Christoph's Personal Wiki
(→Example usage) |
|||
Line 17: | Line 17: | ||
Or, | Or, | ||
echo "BFACTOR ALWAYS 0" | pdbset xyzin input.pdb xyzout output.pdb | echo "BFACTOR ALWAYS 0" | pdbset xyzin input.pdb xyzout output.pdb | ||
+ | |||
+ | ===Extended examples=== | ||
+ | ''Note: These were taken directly from the man page.'' | ||
+ | <pre> | ||
+ | ######################## Convert PDB file to mmCIF format | ||
+ | #!/bin/csh -f | ||
+ | # | ||
+ | pdbset xyzin toxd.pdb xyzout toxd.cif << eof-1 | ||
+ | output cif | ||
+ | end | ||
+ | eof-1 | ||
+ | |||
+ | ######################## Take output from O into a form suitable for refinement | ||
+ | #!/bin/csh -f | ||
+ | # | ||
+ | pdbset xyzin bst_113m.pdb xyzout temp1.pdb << eof-1 | ||
+ | cell 132.02 115.21 96.20 90.00 90.00 90.00 | ||
+ | spacegroup P212121 | ||
+ | eof-1 | ||
+ | |||
+ | ################### Take output from Xplor into a form suitable for refinement | ||
+ | #!/bin/csh -f | ||
+ | # | ||
+ | pdbset xyzin bst_113m.pdb xyzout temp1.pdb << eof-1 | ||
+ | cell 132.02 115.21 96.20 90.00 90.00 90.00 | ||
+ | spacegroup P212121 | ||
+ | xplor | ||
+ | eof-1 | ||
+ | |||
+ | |||
+ | |||
+ | ######################## Expand dimer to tetramer, rename chains, transform | ||
+ | #!/bin/csh -f | ||
+ | # | ||
+ | # Make tetramer from dimer | ||
+ | # | ||
+ | pdbset xyzin ecrproducts268.pdb xyzout ecrprodpqrtet.pdb <<eof-1 | ||
+ | remark Tetramer generated from AB dimer | ||
+ | remark rotated to pqr frame | ||
+ | remark | ||
+ | ! Generate other dimer by z-dyad in P21212 | ||
+ | symgen x,y,z | ||
+ | symgen -x,-y,z | ||
+ | ! Rename chains in second dimer: V & W are water chains | ||
+ | chain symmetry 2 A C | ||
+ | chain symmetry 2 B D | ||
+ | chain symmetry 2 V X | ||
+ | chain symmetry 2 W Y | ||
+ | ! transform to molecular frame | ||
+ | transform - | ||
+ | 0.87831 0.47808 0 - | ||
+ | 0 0 -1. - | ||
+ | -0.47808 0.87831 0 - | ||
+ | 0.0 -2.713 0.0 | ||
+ | eof-1 | ||
+ | </pre> | ||
==External links== | ==External links== |
Latest revision as of 01:34, 9 July 2012
pdbset — is part of the ccp4 package and has various useful manipulations on coordinate files (e.g., PDB files).
Syntax
pdbset XYZIN foo_in.pdb XYZOUT foo_out.pdb [Keyworded input]
Keyworded input
- The available keywords are:
BFACTOR, CELL, CHAIN, COM, ELEMENT, EXCLUDE, OCCUPANCY, ORTHOGONALIZATION, OUTPUT, PICK, REMARK, RENUMBER, REORTHOGONALIZE, REPLACE, ROTATE, SELECT, SEQUENCE, SHIFT, SPACEGROUP, SYMGEN, TRANSFORM, UTOB, XPLOR, NOISE, ATRENUMBER
Example usage
#!/bin/sh -f pdbset xyzin input.pdb xyzout output.pdb << eof OCCUPANCY RESET 90 BFACTOR ALWAYS 0 eof
Or,
echo "BFACTOR ALWAYS 0" | pdbset xyzin input.pdb xyzout output.pdb
Extended examples
Note: These were taken directly from the man page.
######################## Convert PDB file to mmCIF format #!/bin/csh -f # pdbset xyzin toxd.pdb xyzout toxd.cif << eof-1 output cif end eof-1 ######################## Take output from O into a form suitable for refinement #!/bin/csh -f # pdbset xyzin bst_113m.pdb xyzout temp1.pdb << eof-1 cell 132.02 115.21 96.20 90.00 90.00 90.00 spacegroup P212121 eof-1 ################### Take output from Xplor into a form suitable for refinement #!/bin/csh -f # pdbset xyzin bst_113m.pdb xyzout temp1.pdb << eof-1 cell 132.02 115.21 96.20 90.00 90.00 90.00 spacegroup P212121 xplor eof-1 ######################## Expand dimer to tetramer, rename chains, transform #!/bin/csh -f # # Make tetramer from dimer # pdbset xyzin ecrproducts268.pdb xyzout ecrprodpqrtet.pdb <<eof-1 remark Tetramer generated from AB dimer remark rotated to pqr frame remark ! Generate other dimer by z-dyad in P21212 symgen x,y,z symgen -x,-y,z ! Rename chains in second dimer: V & W are water chains chain symmetry 2 A C chain symmetry 2 B D chain symmetry 2 V X chain symmetry 2 W Y ! transform to molecular frame transform - 0.87831 0.47808 0 - 0 0 -1. - -0.47808 0.87831 0 - 0.0 -2.713 0.0 eof-1