Vina generates one output pdbqt file with many conformations. How can I generate a single pdb file of my protein complex with one of the conformations that I choose from the output pdbqt file?
I think you are looking for protein-ligand conformation after docking has been performed right? You can use Chimera to create one pdb file i.e all of the docked posed in your protein could be seen in Chimera, or you can also use AutoDock GUI, it also provides you with the option to do so plus in autodock GUI the residue interactions in the active site can also be seen. Discovery studio visualizer is also an option.
You can use the clustering option in Autodock, where you rank the binding sites using histograms by energy. Then you click the histogram and the "average" of the conformations with that energy range will show, you should be able to save that as a pdb. Only one macromolecule and ligand should show on the output.
Thank you but I'm not familiar with Chimera because I'm using PyMol. Do you know how to generate one pdb file of a complex protein-ligand using pymol? For AutoDock GUI, I have been using it to view the ligand conformation and interactions with amino acid but I do not see option that we could build a complex of protein-ligand in one single file. Do you know how to do it?
Thank you. But that is for analyzing Autodock result. We won't be able to see cluster if we are analyzing Vina result in Autodock Tools. I know how to build protein-ligand complex for Autodock4.2 but I don't see the option for Vina.
I used vina_split program to split the multi conformation docking output *.pdbqt file, then wrote a simple Linux bash shell scipt to build the complex pdb file. Writing simple shell scripts you can make your molecular modelling more effective!
Put your xxx.pdbqt output and protein.pdbqt in one folder.
Run the cygwin and go to this folder and follow only two command as -
1. cut -c-66 xxx.pdbqt > xxx.pdb
(this creates pdb file from your output pdbqt. Edit this pdb so as to include only required conformer. e.g. if you want 1st conformer only keep model 1 in pdb and delete other models. Thus your pdb is having only one desired conformer now)
Please use vina_split to split your file. now you have each conformation in complex with receptor. Then use mgltools or openbables to convert complex with pdbqt format to a file with pdb format.
I want to share my personal experience. There are several ways for combining ligand and receptor after docking. It can be done simply by using Open-BABEL. After docking was performed and .pdbqt files were generated one is able to select both ligand.pdbqt and receptor.pdbqt files simultaneously and convert them into one complex.pdb file with preserve coordinates. If one use PaDEL-ADV for screening all different pose of the ligand will be combined with receptor, then the preferred pose can be chosed in complex with receptor by using Chimera. The the complex file can be visualized by any software like Chimera and etc.... It also also can be used for MDS.
i m facing problem in converting docked file into .pdb format. autodock vina save ligand.pdbqt and receptor.pdbqt different and it do not save both molecule in .pdb format single file