From: Tristan Croll (tristan.croll_at_qut.edu.au)
Date: Fri Jul 24 2015 - 18:00:32 CDT

You would be much better off using AutoPSF for a task like this. PSFgen alone will treat your protein like a single connected chain (so residues either side of a gap will become connected to each other) and, as you have already discovered, won't automatically handle ligands. AutoPSF is the front-end that handles all of that before feeding to PSFgen, and can be called from a non-interactive script.

Cheers,

Tristan

Tristan Croll
Lecturer
Faculty of Health
School of Biomedical Sciences
Institute of Health and Biomedical Engineering
Queensland University of Technology
60 Musk Ave
Kelvin Grove QLD 4059 Australia
+61 7 3138 6443

This email and its attachments (if any) contain confidential information intended for use by the addressee and may be privileged. We do not waive any confidentiality, privilege or copyright associated with the email or the attachments. If you are not the intended addressee, you must not use, transmit, disclose or copy the email or any attachments. If you receive this email by mistake, please notify the sender immediately and delete the original email.

On 24 Jul 2015, at 3:25 pm, Akshay Bhatnagar <akshaybhatnagar2790_at_gmail.com<mailto:akshaybhatnagar2790_at_gmail.com>> wrote:

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--_000_EF433A3AECDD45F9A0BE1A2783576F4Aquteduau_--