1 / 26

Scientific & technical presentation Fragmenter Nóra Máté

Scientific & technical presentation Fragmenter Nóra Máté. Sept 2005. Contents. Fragmenter creating fragments by cleavage rules Fragment statistics sorting fragments by activity values R-group decomposition finding a scaffold with attached ligands

nhung
Télécharger la présentation

Scientific & technical presentation Fragmenter Nóra Máté

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Scientific & technical presentationFragmenter Nóra Máté Sept2005

  2. Contents • Fragmenter • creating fragments by cleavage rules • Fragment statistics • sorting fragments by activityvalues • R-group decomposition • finding a scaffold with attached ligands • web:http://www.chemaxon.com/jchem/doc/user/LibraryProfiling.html

  3. Fragmenter Basics Fragmenter cleaves single bonds to generate molecular fragments. The cleavage rules correspond to chemical reactions in order to enhance synthetic accessibility. Fragmenter fragments molecules based on predefined cleavage rules. The cleavage rules are given in form of reaction molecules in the configuration XML. By default, all non-ring bonds matching the cleavage bonds in the rules are cleaved. However, it is possible to provide a revision algorithm that forbids certain cuts depending on predefined criteria (e.g. the resulting fragment size, the structural environment of the bond, the number of cleaved bonds in the resulting fragments, etc.). Currently one such algorithm is implemented: the RECAP method.

  4. The RECAP Method All non-ring bonds matching the cleavage bonds in the rules are allowed to be cleaved by default.The RECAP algorithmforbids the cleavage of some bonds according to the following rules: • Never cut a hydrogen-connecting bond. • Never cut a bond connecting a ring-carbon and a hetero atom (optional). • Never cut ring bonds. (Fragmenter always keeps this rule, we add it here for completeness.) • Refuse a cut if any of the resulting fragments is on the specified Notlist. • Refuse a cut if the number of open bonds in any of the resulting fragments exceeds the specified limit. • Refuse a cut if the number of atoms in any of the resulting fragments is less than the predefined minimal atom count.

  5. Fragmenter Cleavage Rules The following rules are typically used with the RECAP algorithm, but Fragmenter accepts any custom cleavage rules described by reaction equations. The cleavage points on the fragments are labeled with the cleavage rules:

  6. Configuration Standardization Cleavage reactions Fragmenter parameters RECAP parameters

  7. Fragmentation Example I. An example fragmentation of tamoxifen (left), an oestrogen antagonist and atenolol (right), an anti-hypertension drug:

  8. Fragmentation Example II. An example fragmentation with amine type cleavage bonds: input molecule amine cleavage amine cleavage fragments

  9. Fragmentation Example III. All fragments of the same input molecule (extensive fragmentation):

  10. Fragment Statistics Basics FragmentStatistics creates statistical results from the output of Fragmenter. The simplest usage is to remove duplicate fragments and sort fragments by occurrence, but FragmentStatistics can also sort fragments by a scoring function based on molecule activity or other data read from the input molecules and stored together with the generated fragments.

  11. Fragment Statistics Input / Output • The input of FragmentStatistics is the output of Fragmenter in cxsmiles format with the following fields: • SMILES string • atom labels storing fragment cleavage data • unique ID (used for fragment duplicate check) • input molecule data read from SDFile tag (optional, e.g. molecule activity) • The output of FragmentStatistics is a sorted cxsmiles table with the following data: • SMILES string • atom labels storing fragment cleavage data • atom count • fragment counts per activity categories (number of identical fragments in each activity category, one field for each) • score

  12. The Scoring Function • Fragments are sorted by activity which is calculated in form of a scoring function: • acx (w1*c1 + w2*c2 + ... + wN*cN) • ac is the heavy atom count • w1, w2, ..., wN are the category weights in descending order (default: from +1 to -1, equidistant) • c1, c2, ..., cN are the fragment counts in each category, in descending activity order • x is the exponent of the heavy atom count (default: 1) • If there is no activity data then FragmentStatistics simply removes fragment duplicates and sorts fragments byacx c1wherec1is the fragment count. By default the exponent is1and the score is thusac*c1. • If there are two activity categories then the default scoring function is • ac(c1 - c2),if there are three categories, then it isac(c1 - c3).

  13. Scoring Function Example – single cutoff value Two activity ranges with cutoff value 0.5:

  14. Scoring Function Example – discrete activity range Discrete activity values:

  15. Generating Fragment Statistics I. We start with a large set of input molecules with activity data: Activity = 5 Activity = 4 Activity = 50 Activity = 0.05

  16. Generating Fragment Statistics II. For the purpose of fragment statistics, start with generating a broad set of fragments without the RECAP (or any other) restrictions: Standardization Cleavage reactions Extensive fragmentation

  17. Generating Fragment Statistics III. The generated fragments inherit the activity values from the parent molecule: Fragments (with cleavage and activity data) are stored in cxsmiles format. The activity data is stored in field_1. field_1 = 50

  18. Generating Fragment Statistics IV. First make statistics with duplicate filtering and sorting. These are the 4 most active fragments (byscore = atom count * occurrence): • field_0: • atom count • field_1: • fragment occurrence • field_2: • score • (field_1 * field_2)

  19. Generating Fragment Statistics V. Next include activity data in the scoring, with cutoff value 1. This means that molecules with activity value at least 1 are considered active, while all others are inactive. These are the 4 most active fragments (byscore = atom count * (active occurrence - inactive occurrence)): • field_0: atom count • field_1: fragment occurrence in the active set (score >= 1) • field_2: fragment occurrence in the inactive set (score < 1) • field_3: score • (atom count * (active occurrence - inactive occurrence))

  20. R-group Decomposition – the query R-group decomposition is a special kind of substructure search that aims at finding a central structure - scaffold - and identify its ligands at certain attachment positions. The query molecule consists of the scaffold and ligand attachment points represented by R-groups: The two R1 nodes should match identical structures by default – but this behaviour can be changed by setting the –p (--skip-same-structure-check) parameter

  21. R-group Decomposition – the targets Our sample targets all contain the query (scaffold with R-group attachment points) but not all of them can satisfy the condition of identical R1-ligands: single hit: identical R1 ligands single hit: the same R1 ligand (R-bridge) more hits: all with different R1 ligands more hits: one with identical R1 ligands

  22. R-group Decomposition – decomposition I. Attachment points can be denoted by different symbols, depending on the –a (--attachment-symbol) option: N: none P: attachment point A: any-atom (default) M: atom map L: atom label

  23. R-group Decomposition – decomposition II. SMILES table: the output is written in a SMILES table if the –f (--format) option is omitted: Otherwise the output is written as molecule series in the specified output form, with atom color codes (separated by ‘;’ characters) stored in the DMAP property (SDF/MRV tag). The code is: 0: scaffold atom n: Rn ligand atom (n > 0) -: non-hit atom Example: 0;0;1;1;1;1;2;2;1;0;-;-

  24. R-group Decomposition – decomposition III. The DMAP property can be used in mview to color the atoms according to a color-map file that maps the color codes to colors. Set the –p (--skip-same-structure-check) option to allow the two R1 nodes match different ligands. Finally, use the –A (--allhits) option to see all possible decompositions. In this way our last target will also have two decompositions:

  25. Visit othertechnical presentations MarvinSketch/Viewhttp://www.chemaxon.com/MarvinSketch_View.ppt MarvinSpacehttp://www.chemaxon.com/MarvinSpace.ppt Calculator Pluginshttp://www.chemaxon.com/Calculator_Plugins.ppt JChem Basehttp://www.chemaxon.com/JChem_Base.ppt JChem Cartridgehttp://www.chemaxon.com/JChem_Cartridge.ppt Standardizerhttp://www.chemaxon.com/Standardizer.ppt Screenhttp://www.chemaxon.com/Screen.ppt JKlustorhttp://www.chemaxon.com/JKlustor.ppt Fragmenterhttp://www.chemaxon.com/Fragmenter.ppt Reactor http://www.chemaxon.com/Reactor.ppt

  26. References Fragmenter, fragment statistics: http://www.chemaxon.com/jchem/doc/user/Fragmenter.html http://www.chemaxon.com/jchem/doc/user/FragmentStatistics.html http://www.chemaxon.com/jchem/examples/fragmenter • RECAP - Retrosynthetic Combinatorial Analysis Procedure: A Powerful New Technique for Identifying Privileged Molecular Fragments with Useful Applications in Combinatorical Chemistry In: J. Chem. Inf. Comput. Sci. 1998, 38. 511-522 • Schneider, G. et al.; De novo design of molecular architectures by evolutionary assembly of drug-derived building blocks. J. Comput.-Aided Mol. Des. 2000, 14, 487-494 R-group decomposition: http://www.chemaxon.com/jchem/doc/user/RGroupDecomposition.html

More Related