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.

Similar questions and discussions