From: Francesco Pietra (
Date: Sat Jun 23 2018 - 11:46:37 CDT

>From a system TIP3-solvated in a periodic box, I am trying to obtain the
.psf for only the protein and one of the ligands.

Loaded the solvated .psf into vmd, from tk console:

% set sel [atomselect top "segname PROB or segname GDPB"]
> $sel writepsf protein-GDP.psf
> frame -1 out of range for molecule 0

I have noticed this type of error on the vmd list but I can't find a
solution here.

thanks for advice

francesco pietra