VMD-L Mailing List

From: Akshay Bhatnagar (akshaybhatnagar2790_at_gmail.com)
Date: Thu Jul 23 2015 - 23:33:56 CDT

Thank you very much for your suggestions.

As you said, i have started with the tcl script (shown below), i am able to
upload all the molecules in one go, but for some reason the output psf and
pdb do not show any atoms. can anyone please help me to correct the script.

the script is:

set filelist [glob *.PDB]
  foreach file $filelist {
    mol new $file
package require psfgen
topology top_all27_prot_lipid.inp
pdbalias residue HIS HSE
pdbalias atom ILE CD1 CD
writepdb $file.pdb
writepsf $file.psf
}

Your help is greatly appreciated.

With Regards
Akshay Bhatnagar
PhD Student
BITS Pilani Hyderabad Campus

On Wed, Jul 22, 2015 at 12:14 AM, Josh Vermaas <vermaas2_at_illinois.edu>
wrote:

> Just an FYI, this usually won't work particularly well unless the
> proteins are well behaved. Usually there are missing pieces in pdb files
> that need manual intervention. That being said, autopsf has a tcl interface:
>
> package require autopsf
> autopsf -mol top -top topologyfilehere
>
> Check the other autopsf options to see if they might be useful to you.
>
> -Josh Vermaas
>
>
> On 07/21/2015 02:33 AM, Ashar Malik wrote:
>
> To start you off, make a list of all the pdb. Read them into a tcl script.
> Run a for loop for the number of pdb files you have and in each loop use
> the psfgen package.
>
> Try it. If you are still stuck .. Drop another email.
> On Jul 21, 2015 12:19 PM, "Akshay Bhatnagar" <
> akshaybhatnagar2790_at_gmail.com> wrote:
>
>> Hello everyone
>>
>> I have around 5000 proteins and i want to generate 5000 psf and
>> subsequently generated pdb files for all these 5000 proteins in one shot. I
>> have always use automatic psf generate method to generate psf. Now please
>> anyone help in generating a tcl script through which i can generate all the
>> 5000 psf and pdb files in one go using the script.
>>
>>
>> With Regards
>> Akshay Bhatnagar
>> PhD Student
>> BITS Pilani Hyderabad Campus
>>
>>
>