Error on running prepare_flexreceptor.py for flexible docking about autodock-vina HOT 8 CLOSED

Hi @eunos-1128 !
I looked at the PDBQT.txt you shared. Quick question: Based on the sequence, is it human dopamine D3 receptor (but unfolded)? How does it compare to a conventional, folded structure, PDB 3PBL, for example? I visualized it and found a gap (long bond) in the backbone between ASP246 and PRO247. It might be helpful to double check the integrity of your starting structure, unless you were trying to perform docking on the unfolded structure

and some other files produce flex PDBQT files but not rigid PDBQT ones due to the same error.
Could you share your command so we can reproduce it?

I want to know about flexibility model used in AutoDock Vina.
In AutoDock Vina, as you might have already noticed, the flexibility model of flex residues is the same as applied to ligands. The sidechain atoms (CB and further from the backbone) are the only mobile portion of the residues. Some residues, however, cannot be made flexible. If I remember correctly these include Pro, Gly and Ala. Alanine's sidechain is like a methyl group (-CH3) in any other ligands which under AutoDock's atom typing will be merged into an aliphatic C, and it doesn't have a torsional degree of freedom.

from autodock-vina.

I'm not sure why rigid files couldn't be generated while flex ones could,though.

I'm not sure, either :< But please let us know if you are still experiencing the same problem with a visually-ok structure of the protein!

from autodock-vina.

Hi @eunos-1128
Do you think that might have something to do with selection of flexible residues (some residues can't be made flexible under AutoDock's flexibility model)? Can you share an example of the input files? We can look at the input files and will let you know if we're able to reproduce this

from autodock-vina.

@rwxayheee

Thanks a lot.

The attached PDBQT file is receptor's one the script produce no rigid/flex PDBQT file due to the error.

PDBQT.txt

from autodock-vina.

@rwxayheee
and some other files produce flex PDBQT files but not rigid PDBQT ones due to the same error.

Do you have any ideas about that?

from autodock-vina.

@rwxayheee
Sorry for continuous posts

AutoDock's flexibility model

I want to know about flexibility model used in AutoDock Vina. Could you tell me the paper describing that mechanism?(I didn't find out that myself)

Hi @eunos-1128 Do you think that might have something to do with selection of flexible residues (some residues can't be made flexible under AutoDock's flexibility model)? Can you share an example of the input files? We can look at the input files and will let you know if we're able to reproduce this

from autodock-vina.

@rwxayheee

I looked at the PDBQT.txt you shared. Quick question: Based on the sequence, is it human dopamine D3 receptor (but unfolded)? How does it compare to a conventional, folded structure, PDB 3PBL, for example? I visualized it and found a gap (long bond) in the backbone between ASP246 and PRO247. It might be helpful to double check the integrity of your starting structure, unless you were trying to perform docking on the unfolded structure

Thanks. As you notice they are the dopamine D3 receptor predicted by alphafold2 with custom settings.

I didn't check the structure using graphical tool. As you say I found it unfolded and so it might not needed for my study.

Could you share your command so we can reproduce it?

Commands are same to the initial one.

Below are those files.
no-rigid_pdbqt1.txt
no-rigid_pdbqt2.txt

The files are generated by alphafold2 and as one I uploaded before, they are also partially unfolded, which could lead to the error. I'm not sure why rigid files couldn't be generated while flex ones could,though.

In AutoDock Vina, as you might have already noticed, the flexibility model of flex residues is the same as applied to ligands. The sidechain atoms (CB and further from the backbone) are the only mobile portion of the residues. Some residues, however, cannot be made flexible. If I remember correctly these include Pro, Gly and Ala. Alanine's sidechain is like a methyl group (-CH3) in any other ligands which under AutoDock's atom typing will be merged into an aliphatic C, and it doesn't have a torsional degree of freedom.

Thank you! I understand.

from autodock-vina.

I'm not sure why rigid files couldn't be generated while flex ones could,though.

I'm not sure, either :< But please let us know if you are still experiencing the same problem with a visually-ok structure of the protein!

OK, Thanks a lot!

from autodock-vina.