MBG wiki | RecentChanges | Blog | 2021-10-20 | 2021-10-19

Preparing a protein-only DCD file (removing waters)

Difference (from prior major revision)

Changed: 17c17

< <source>

to

> <code>

Changed: 22c22

< </source>

to

> </code>

Changed: 24c24

< <source>

to

> <code>

Changed: 29c29

< </source>

to

> </code>


Using catdcd

Use catdcd with the -i flag to prepare a protein-only DCD. The unix way to prepare the index file is

awk '{for (i=0 ; i < 434 ; i++) print i}'

where 434 is the number of atoms of the solute.


Using carma

If your macromolecular component has segids, say, A and B, then you can in one step remove waters, translations and rotations with carma :

carma -v -fit -atmid "ALLID" -segid " A   " -segid " B   " my.dcd my.psf


Preparation of a protein-only PSF

You will also need a protein-only PSF file. You can do this with an xplor script in the spirit of :

structure @ionized.psf end
delete selection=( resname TIP3 or resname SOD or resname CLA) end
write structure output=protein.psf end
stop

, or,

structure @ionized.psf end
delete selection=( not (segid="A   " or segid="B   ")) end
write structure output=protein.psf end
stop