DUDE: Difference between revisions
m (Beginning documentation) |
mNo edit summary |
||
Line 19: | Line 19: | ||
* subset_decoys.py in the target directory can covert the full dudgen_ecfp4 set into a dudgen_clustered type set given a list of molregno ids. | * subset_decoys.py in the target directory can covert the full dudgen_ecfp4 set into a dudgen_clustered type set given a list of molregno ids. | ||
= FAQ = | |||
== Q1. How duplicates are removed == | |||
In the paper, you said "We then remove duplicate decoys from the ligand set by sorting decoys from | |||
least to most duplicated and assigned each decoy to the protonated ligand | |||
which has the least number of decoys already assigned." I don't know which | |||
molecules you considered as duplicates, and how to define and calculate the | |||
similarity. Please explain! | |||
== A1. == | |||
Decoys are uniquely identified by their ZINC protonation ids. We | |||
ensure that a particular protonation id (prot_id) is only assigned to | |||
one ligand of a given target. | |||
* 0) After filtering the for the 25% most dissimilar decoys | |||
* 1) map each prot_id to the number of different ligands it could be assigned to | |||
* 2) sort from the prot_ids that hit the fewest ligands to those that hit the most | |||
* 3) loop over that sorted list, starting with the most constrained decoys | |||
* 4) assign each prot_id to the the ligand with the fewest other prot_ids so far | |||
The effect is to spread the decoys as evenly as possible among the ligands they could belong to. | |||
END OF FAQ | |||
[[Category:DUDE]] | [[Category:DUDE]] | ||
[[Category: | [[Category:FAQ]] |
Revision as of 15:32, 12 July 2012
This is the Wiki Page for DUD-E, a directory of useful decoys - enhanced. DUD-E is on the web at http://dude.docking.org.
This page may be used for posting errors and ommissions in DUD-E, and also for commenting on the database's design and usefulness.
File explanation
The clustered sets live in dudgen_clustered while the full raw sets are in dudgen_ecfp4. Inside, the ligands.charge file contains the mapping from chembl ids to unique property sets, and the search/decoys.*.picked contain the actual decoys for that protonation set.
The file formats are as follows:
- ligands.charge - gives unique protonation states of input ligands.
Format: one ligand protonation form per line SMILES input_id protonation_id mwt logp rb hba hbd charge
- search/
- decoys.<protonation_id>.picked - contains matched decoys for each unique ligand protonation state
Format: ligand protomer and then 50 matched decoys first line: ligand SMILES input_id protonation_id SMILES ZINC_ID ZINC_Protonation_ID
- subset_decoys.py in the target directory can covert the full dudgen_ecfp4 set into a dudgen_clustered type set given a list of molregno ids.
FAQ
Q1. How duplicates are removed
In the paper, you said "We then remove duplicate decoys from the ligand set by sorting decoys from least to most duplicated and assigned each decoy to the protonated ligand which has the least number of decoys already assigned." I don't know which molecules you considered as duplicates, and how to define and calculate the similarity. Please explain!
A1.
Decoys are uniquely identified by their ZINC protonation ids. We ensure that a particular protonation id (prot_id) is only assigned to one ligand of a given target.
- 0) After filtering the for the 25% most dissimilar decoys
- 1) map each prot_id to the number of different ligands it could be assigned to
- 2) sort from the prot_ids that hit the fewest ligands to those that hit the most
- 3) loop over that sorted list, starting with the most constrained decoys
- 4) assign each prot_id to the the ligand with the fewest other prot_ids so far
The effect is to spread the decoys as evenly as possible among the ligands they could belong to.
END OF FAQ