Comments (7)
Then it's running. But with that many sidechains, it just takes extremely long.
Molecular viewers show gaps because the protein is split between different files, but all atoms should be there. Each atom should be in either _flex or _rigid, not both.
from autodock-vina.
Hi @farmaceut Forgive me for interjecting again, one more point to consider with your selection of flexible residues is the size of the box. If the box does not cover all the mobile atoms in the flexible residues (for example, Asp100, Glu133, Tyr199), you will get enormously large INTRA energies in the docking output (and those residues will appear to be strained, in the output. So please make sure to check if the size of the box is sufficient to contain residues' sidechains, preferably with some extra room for alternate sidechain conformations.
from autodock-vina.
For people to be able to help you, you need to provide more detailed information:
- exact vina command and all input files to reproduce the problem. You posted only part of what's needed.
- "the job doesn't proceed" can mean a number of things. Does it print something?
- "some residues are missing". What residues are missing?
from autodock-vina.
@diogomart
Command lines:
prepare_flexreceptor.py -r 4lxd.pdbqt -s ALA97_ASP100_PHE101_SER102_TYR105_ASP108_PHE109_MET112_VAL130_GLU133_LEU134_ASN140_TRP141_GLY142_ARG143_VAL145_ALA146_GLU149_PHE150_VAL153_PHE195_TYR199
prepare_gpf.py -l H4Fst.pdbqt -r 4lxd_rigid.pdbqt
autogrid4 -p 4lxd_rigid.gpf -l 4lxd_rigid.glg
vina --flex 4lxd_flex.pdbqt --ligand H4Fst.pdbqt --maps 4lxd_rigid --scoring ad4 --exhaustiveness 6 --out output.pdbqt
The output:
#################################################################
# If you used AutoDock Vina in your work, please cite: #
# #
# O. Trott, A. J. Olson, #
# AutoDock Vina: improving the speed and accuracy of docking #
# with a new scoring function, efficient optimization and #
# multithreading, Journal of Computational Chemistry 31 (2010) #
# 455-461 #
# #
# DOI 10.1002/jcc.21334 #
# #
# Please see http://vina.scripps.edu for more information. #
#################################################################
Scoring function : ad4
Flex receptor: 4lxd_flex.pdbqt
Ligand: H4Fst.pdbqt
Exhaustiveness: 6
CPU: 0
Verbosity: 1
Reading AD4.2 maps ... done.
Performing docking (random seed: 165862369) ...
0% 10 20 30 40 50 60 70 80 90 100%
|----|----|----|----|----|----|----|----|----|----|
*
It's just stuck at *. It doesn't change even though I waited an hour.
The files are in archive:
archive.zip
I use Chimera. Opening simulatenously 4lxd_flex.pdbqt with 4lxd_rigid.pdbqt you can see there are gaps between residues which are filled by AAs in case of original 4lxd.pdbqt file.
from autodock-vina.
Hi @farmaceut
The inputs look correct. I was able to run the calculation, too; The job does proceed very slowly as @diogomart pointed out. If you're not running the job on a cluster, I recommend starting with fewer flexible residues :>
also because your ligand is much smaller than the known ligand in 4lxd. It can only occupy a subsite of this binding site, so having fewer flexible residues might be worth a try. Maybe you could run some jobs with different selection of flex residues and then compare.. or summarize them?
About the missing residues: Gly, Pro and Ala cannot be made flexible, because their sidechains do not have rotatable torsions (for Ala, the nonpolar hydrogens are excluded from PDBQT)
from autodock-vina.
@diogomart @rwxayheee
Oh well. I though, it was me doing something wrong. At least I am aware my old PC is not enough :( Wonder if running on M3 would make it notably faster?
from autodock-vina.
Possibly. But not as fast as having fewer sidechains.
from autodock-vina.
Related Issues (20)
- Running Autodock Vina with batch Ligand occur error basic_string::_M_replace_aux HOT 3
- Clarification on hydrophobic interaction contribution HOT 2
- Charges in Vina HOT 1
- How can I add 2+ to metal atoms using `prepare_receptor`? HOT 9
- Can Autodock Vina perform Constraint Docking? HOT 2
- How to use vina for small molecule and small molecule docking HOT 1
- Release v1.2.5 HOT 13
- Problems with Ruthenium HOT 2
- Flexible docking of Autodock using vina HOT 7
- Can't install 1.2.5 on Win11 HOT 2
- Python bindings vina v.s vina executable HOT 13
- Details on scoring function terms HOT 3
- Autogrid4 Atom Type Error with Large Ligand Database HOT 6
- An internal error occurred in ../../../src/lib/parse_pdbqt.cpp(69). HOT 3
- Autogrid4 raise error "dyld: Library not loaded:" HOT 8
- Enormous free energy of binding HOT 6
- Is AutoDock-Vina applicable to peptides? HOT 3
- Docking to zinc atoms in a non-3,1 tetrahedral conformation HOT 4
- Vina flexible docking HOT 4
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from autodock-vina.