Virtual Libraries and Virtual Screening in Drug Discovery Processes using KNIME

Transcription

1 Virtual Libraries and Virtual Screening in Drug Discovery Processes using KNIME Iván Solt Solutions for Cheminformatics

2 Drug Discovery Strategies for known targets High-Throughput Screening (HTS) Cells or recombinant protein Fluorescent or luminescent readout Automated, miniaturized Thousands of samples / day Number of primer actives: ~1% Virtual Screening (VS) Ligand or structure based Virtual or real libraries Similarity search, 2D or 3D Can lead to thousands of possible actives: further processing needed Measurement: Enrichment ratio, ROC curves for known actives

3 Virtual Library Design Workflow DB DB Databases Reactions Molecules Queries Fragmentation R-group decomposition Fragmentation Reagent clipping Compound selection Similarity searches Substructure searches Enumeration Fuse fragments R-group composition Reaction enumeration Library analysis Clustering 2D similarity screen 3D Shape similarity screen

4 Find or Virtually Create Candidates Virtual screening of existing compounds Pros: Fast Hits are readily available for in vitro experiments Cons: Limitation on available compounds De novo design Pros: No limitation on virtual compound space Structural novelty Cons: Are hits synthetically available?

5 Find or Virtually Create Candidates Virtual screening of existing compounds Pros: Fast Hits are readily available for in vitro experiments Cons: Limitation on available compounds De novo design Pros: No limitation on virtual compound space Structural novelty Cons: Are hits synthetically available?

6 Virtual Screening Workflow DB DB Molecules in-house or commercially available 1. Reactions virtual synthetic path Synthetically Accessible Compounds 2. Filtering in vivo experiment? 5. Clustering 4. 3D alignment 3. Similarity Search

7 Step 1: Reaction Enumeration Reaction schema for accessible syntheses Combinatorial or sequential enumeration Reaction rules: phrase + apply public and in-house chemical knowledge Selectivity with tolerance Reactivity Exclusion rules EXCLUDE: match(reactant(1), "[Cl,Br,I]C(=[O,S])C=C") or match(reactant(0), "[H][O,S]C=[O,S]") or match(reactant(0), "[P][H]") or (max(pka(reactant(0), filter(reactant(0), "match('[o,s;h1]')"), "acidic")) > 14.5) or (max(pka(reactant(0), filter(reactant(0), "match('[#7:1][h]', 1)"), "basic")) > 0)

8 Step 1: Reaction Enumeration

9 Step 1: Reaction Enumeration Reaction rules ON Fewer results than theoretical Unfeasible starting materials eliminated Feasible products only Custom rules can be added to increase selectivity Reaction rules OFF More results Best for debugging purposes Prodcts may be incorrect due to neglecting chemical rules

10 Step 2: Filtering Lead likeness, drug likeness Chemical Terms Could it fit to the active centre? Basic analysis: size, mass... Could it get to the active centre? ADME properties: solubility, pka, polar surface, partition coefficients... Structural filtering e.g. reactive groups Toxicity, environmental concerns, etc... Calculator plugins Elemental Analysis Elemental Analysis IUPAC Name Structure to Name Protonation pk a Microspecies Isoelectric Point Partitioning logp logd Charge Charge Polarizability Orbital Electronegetivity Isomers Tautomerization Stereoisomer Conformation Conformer Flexible 3D Alignment Molecular Dynamics Geometry Topology Analysis Geometry Polar Surface Area (2D) Molecular Surface Area (3D) Markush Markush Enumeration Other Hydrogen Bond Donor- Acceptor Huckel Analysis Refractivity Structural Framework Resonance

11 Step 3: Similarity search Screen 2D + Descriptor package Screen against known bioactives Chemical Fingerprints Topology Pharmacophore Fingerprints: Custom atomic properties + their topological relationship H-bond donors / acceptors Cationic / anionic groups Hydrophobic groups Aromatic groups etc. ECFP/FCFP Similarity searches Tanimoto, Eucledian, Tversky metrics Metrics optimization regular Tanimoto optimized Tanimoto

12 Step 4: Screen 3D Align the candidates to the known active in 3D Treat the candidate flexible! Consider pharmacophore atom types (align cationic to cationic, etc.)! Problem: complicated conformational space

13 Step 4: Screen 3D Simple sampling of the conformational space: Minimum and maximum distance between atom pairs in the full torsion space Select atoms Colors (e.g. pharmacophore types ) Topological features (e.g.:longest chain start/end/center) Ring centers (aromatic, aliphatic) Calculate Min/max internal distance ranges Distance histograms for selected atoms Only once for each molecule

14 Step 4: Screen 3D Hybrid alignment: Separate translation&rotation from torsions Robust and goes fast Needs good guess on atomatom mapping: Same colors Distance ranges must be allowed for all mapped pairs Triangle inequality must be fulfilled for any atom triplet

15 % of the actives retrieved Screen 3D: Test on DUD 30 Average of 1% Enrichments Giganti et al. J. Chem. Inf. Model. 2010, 50, 992

16 % of the actives retrieved Screen 3D: Test on DUD 100 Average of 10% enrichments Giganti et al. J. Chem. Inf. Model. 2010, 50, 992

17 Screen 3D: Test on DUD Average time per compound (without precalculations) ChemAxon Screen3D 0.07 ROCS 0.5 FRED 1.0 ICMsim 2.4 Surflex-sim 6.7 FlexS 6.9 Surflex-dock 14.6 FLEXX 15.6 ICM 17.7 Speed Intel Q GHz Intel Xeon 2.4 GHz Giganti et al. J. Chem. Inf. Model. 2010, 50, 992

18 Step 5: Clustering, library analysis JKlustor Wide range of methods Unsupervised, agglomerative clustering Hierarchical and non-hierarchical methods Similarity based and structure based techniques Flexible search options Tanimoto and Euclidean metrics, weighting Maximum common substructure identification chemical property matching including atom type, bond type, hybridization, charge

19 JChem Extensions in KNIME Worklflow management in KNIME JChem extension nodes developed by InfoCom, Japan Constantly developing palette of available JChem tools

20 JChem Extensions in KNIME IO molecule and reaction import, export, drawing Visualization Manipulators Calculator plugins Reactor Similarity and structure-based search Fingerprint calculation Fragmentation Clustering R-group composition, decompozition Standardization... Database management Molecular format conversion Web search services

21 Step 1: Reaction Enumeration

22 Step 2: Filtering

23 Step 2: Filtering

24 Step 3: Similarity search

25 JChem Extensions in KNIME DB DB 1. Reactions virtual synthetic path Synthetically Accessible Compounds 2. Filtering in vivo experiment? 4. 3D alignment 3. Similarity Search 1. Import reactants 2. Enumerate reaction Carry out topology analysis 3. Calculate properties Filter 4. Screen for similarity against known active 5. Export results

26 Conclusions Virtual libraries and virtual screening are essential tools in modern Drug Discovery No special hardware, short experiment cycles, variety of approaches Database of synthetically accessible compounds can be designed with reaction libraries and custom in-house synthetic knowledge Powerful 3D alignment techniques allow highthroughput conformational screening with great efficiency Straightforward integration into KNIME

27 Contributors Tímea Polgár Attila Tajti

28