Dear all,
In recent, I download a pdb file from RCSB (PDB ID: 3HBT) and want to generate a psf file. According some texts , I used the topology file (top_all27_prot_lipid_na.inp and top_all27_prot_lipid.rtf) to creat the psf file by a tcl script. The tcl script as follow:
resetpsf
package require psfgen
topology ./top_all27_prot_lipid_na.inp
topology ./top_all27_prot_lipid.rtf
pdbalias atom ILE CD1 CD
pdbalias atom HOH O OH2
pdbalias residue HOH TIP3
pdbalias residue HIC HSE
segment A {
auto none
pdb 3HBT.pdb
}
coordpdb 3HBT.pdb A
guesscoord
writepdb actin.pdb
writepsf actin.psf
exit
However, VMD had error in running the script like this :
psfgen) reading residues from pdb file 3HBT.pdb
psfgen) unknown residue type CA
psfgen) extracted 424 residues from pdb file
psfgen) Info: generating structure...psfgen) unknown atom type OG2P1
psfgen) add atom failed in residue SO4:500
failed!
ERROR: failed on end of segment
MOLECULE DESTROYED BY FATAL ERROR! Use resetpsf to start over.
In another one, I used the charmm-gui to creat a topology file and psf file. But the psf file could not be used in NAMD2.9.
I do not know how to creat a psf file correctly.
Looking forward for your help and suggestion.