DUDE: Difference between revisions

From DISI
Jump to navigation Jump to search
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:DUD]]
[[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