eplop, or Easy PLOP, is a wrapper script for PLOP. It prepares the ligand and the receptor for use with PLOP, and can also invokes PLOP to performing re-scoring. There is option to turn on and off the ligand flexibility, and it could also save all working files so that you could manually tune the PLOP control file to perform more sophisticated analysis. The scoring output format is consistent with that of MUD - Michael's Utilities for Docking and SCORE.
The Simplest Approach
eplop with paths to the receptor and the ligand:
~ycao/eplop/eplop receptor.pdb ligand.mol2
Note that the receptor structure must be in PDB format and that of the ligand must be in mol2 format. You may use
convert.py to convert between different formats. Please be cautious when using
This will prepare the receptor using the first ligand in the
ligand.mol2 file, and then use the prepared receptor to score all ligands in the
ligand.mol2 file. The output will be written to the standard output as well as an output file named
The Two-Step Approach
Typically, you should prepare the receptor using the crystal ligand, and then use the prepared receptor to score the other ligands.
To prepare, run:
~ycao/eplop/eplop prepare receptor.pdb ligand.mol2
It will use only the first ligand in the
ligand.mol2 no matter how many ligands are present in that file. This will generated a PDB file whose name has a suffix
-prepared.pdb. Now use that to score:
~ycao/eplop/eplop score receptor-prepared.pdb other-ligands.mol2
Configuring the Scoring Methods
If no methods are specified, the preparing step by default will add
opt yes in the PLOP control file to prepare the protein, and the scoring step by default will assume ligand to be flexible using
minimize residue in the PLOP control file. You could use the
--method option to change the methods.
An alternative preparing method is
pka, which uses
pka clust in the PLOP control file. This is (much much) more time consuming, but may give you better results. For example:
~ycao/eplop/eplop prepare --method pka receptor.pdb ligand.mol2
Alternative scoring methods are
rigid, which treats ligand as rigid and only changes its orientation, and
no_optim which only performs energy calculation and no optimization using the given coordinates of ligands. For example:
~ycao/eplop/eplop score --method no_optim receptor.pdb ligand.mol2
If you are using the one-step simple approach to do both prepare and score, you could specify both the method of preparation and that of scoring:
~ycao/eplop/eplop --method pka:no_optim receptor.pdb ligand.mol2
or just the preparation method:
~ycao/eplop/eplop --method pka: receptor.pdb ligand.mol2
or just the scoring method:
~ycao/eplop/eplop --method :no_optim receptor.pdb ligand.mol2
All the options are very useful. Run
to see a list of supported options. Here are some important ones.
-ctells eplop to write the output in format compatible with Michael's script.
-stells eplop to save the pose of the ligands.
-ptells eplop to keep the working directory. This is useful if you want to take a took at the log or the mol2 file containing the poses of ligands, or if you want to reuse the working directory to save time when used with
-otells eplop to write the score to an alternative file instead of standard
-atells eplop to use an existing
plop.scoresfile in the current folder and do not recalculate scores that already exist in that file. Please note that if you also set output file using
-o, it will read from the designated output file instead of
-qchanges the verbosity of the script.