From: Brian Radak (bradak_at_anl.gov)
Date: Thu Dec 15 2016 - 14:28:56 CST
I'm not sure that "pdbalias residue" is doing what is expected - can
someone enlighten me?
I am reading a pre-made psf/pdb/coor/vel set from files, but want to
change some of the topologies that are used. Specifically I want every
instance of the standard C36 residues HSD, HSE, and HSP to become HIS (a
new topology that I made and loaded an RTF for). This is essentially the
same thing as opening the PSF and PDB and running sed to replace
HSD/HSD/HSP with HIS. I thought the following would do what I want:
pdbalias residue HSD HIS
pdbalias residue HSE HIS
pdbalias residue HSP HIS
readpsf <blah>
writepsf <some new psf>
writepdb <some new pdb>
but this does not seem work as I wanted it to (the output pdb still
contains HSD residues in my test). Does pdbalias only work when using
the segment command? If so, is that a bug or a feature?
Brian
-- Brian Radak Postdoctoral Appointee Leadership Computing Facility Argonne National Laboratory 9700 South Cass Avenue, Bldg. 240 Argonne, IL 60439-4854 (630) 252-8643 brian.radak_at_anl.gov
This archive was generated by hypermail 2.1.6 : Sun Dec 31 2017 - 23:20:53 CST