Tautomerism in chemical information management systems

Transcription

1 Tautomerism in chemical information management systems Dr. Wendy A. Warr

2 Tautomerism in chemical information management systems Author: Wendy A. Warr DOI: /s Perspectives Issue Devoted to Tautomerism in Molecular Design Edited by Yvonne Martin

3 Chemical Information Aspects Registration procedures Storage of tautomers Exact and substructure search Depiction of results

4 Software and Database Vendors Accelrys ACD/Labs Beilstein/Reaxys CambridgeSoft CAS CCDC CCG ChemAxon ChemoSoft ChemSpider CWM Global Search Daylight Dialog IDBS InfoChem InhibOx John Wiley & Sons Molecular Networks NCI/CADD OpenEye Thieme PubChem Questel Schrödinger SciTouch Symyx Thomson Reuters Xemistry (CACTVS)

5 Not Included ABCD (J&J) BioRad (KnowItAll) CDK emolecules SimBioSys Tripos ZINC

6 Chemical Structure Representation

7 Morgan Algorithm Morgan, H. L. The generation of a unique machine description for chemical structures - a technique developed at Chemical Abstracts Service. J. Chem. Doc. 1965, 5(2),

8 CTfile

9 SMILES CC1=CC(Br)CCC1

10 SMILES OpenEye canonical SMILES Daylight canonical SMILES SciTouch canonical SMILES ChemAxon canonical SMILES

11 IUPAC International Chemical Identifier (InChI) InChI=1S/C8H10N4O2/c (10)7(13)12(3)8(14)11(6)2/h4H,1-3H3 InChIKey=RYYVLZVUVIJVGH-UHFFFAOYSA-N

12 NCI/CADD Identifiers (CACTVS Hashcodes) 9850FD9F9E2B4E25-FICTS FD9F9E2B4E25-FICuS FD9F9E2B4E25-uuuuu-01-27

13 Definition of Tautomerism M=Q-ZH HM-Q=Z Q = C, N, S, P, Sb, As, Se, Te, Br, Cl or I M, Z = trivalent N, bivalent O, S, Se or Te [Either M or Z = C] H = H, D, T [or + or -] Extended system, ring/chain, etc.

14 Straightforward 1,3 shift 1,5 shift 1,7 shift

15 More Complex

16 Degree of Unsaturation

17 Ring Opening

18 Fluxional structures

19 Mesomers Different NEMA Keys NEMA Key=6P1SUP7NENNHV4V61WRZP5S2ES8NZF InChI=1S/C16H18N3S.ClH/c1-18(2) (9-11) (19(3)4)6-8-14(16)17-13;/h5-10H,1-4H3;1H/q+1;/p-1 NEMA Key=CKGEHDBX4KZPW3VV6DXTVM5BB689GB InChI=1S/C16H18N3S.ClH/c1-18(2) (9-11) (19(3)4)6-8-14(16)17-13;/h5-10H,1-4H3;1H/q+1;/p-1 CXKWCBBOMKCUKX-UHFFFAOYSA-M Same InChIKey

20 Tautomers Different NEMA Keys NEMA Key=CU3YSHT7DX8KUTKGRNS5GH3B4UQBFA InChI=1S/C10H13N5O3/c (14-10(17)13-8)15(4-12-7) (3-16)18-6/h4-6,16H,1-3H2, (H3,11,13,14,17)/t5-,6+/m0/s1 NEMA Key=CTDBHWQW8CQJHC3S5AH6X4QJWAVMKD InChI=1S/C10H13N5O3/c (14-10(17)13-8)15(4-12-7) (3-16)18-6/h4-6,16H,1-3H2, (H3,11,13,14,17)/t5-,6+/m0/s1 KITPKMKMNZXFDK-NTSWFWBYSA-N Same InChIKey

21 Unreasonable

22 5 Multiple Overlapping 6 7 8

23 Overlapping

24 Registration

25 Registration Objectives Corporate database Stock room database Predicting spectra Reaction mechanisms Ultra-low temperature lab

26 Registration Options Enumerate all tautomers; store all tautomers Calculate canonical tautomer; store canonical tautomer Enumerate all tautomers Rank [as major, minor, or conditions dependent (ACD/Labs)] Allow user to choose which form to store

27 Epik Schrödinger Enumerate all energetically reasonable tautomers Enumerate all energetically reasonable ionization states Store all tautomers and ionization states Canvas identifies duplicates by canonical SMILES

28 Are A and B Tautomers? If A and B are identical, accept If the total number of hydrogen atoms or charges is not identical, reject Examine the heavy-atom skeletons; reject if not identical Enumerate all tautomers for A; if any is the same as B, accept Enumerate all tautomers for B; if any is the same as A, accept Otherwise reject.

29 Enumeration of tautomers Sayle, R. A.; Delany, J. J. Canonicalization and enumeration of tautomers. Paper presented at EuroMUG99, Cambridge, UK, Oct 1999 Oellien, F.; Cramer, J.; Beyer, C.; Ihlenfeldt, W-D.; Selzer, P. M. (2006) The impact of tautomer forms on pharmacophore-based virtual screening. J. Chem. Inf. Model. 2006, 46, Greenwood, J. R.; Calkins, D.; Sullivan, A. P.; Shelley, J. C. Towards the comprehensive, rapid, and accurate prediction of the favorable tautomeric states of drug-like molecules in aqueous solution. J. Comput.-Aided Mol. Des. 2010, published online March 31, 2010

30 Storage of Tautomers

31 Concept A Generate all tautomers Impossible to calculate lowest energy tautomer Use rules for consistent generation [of a low energy form] Store this form [as canonical SMILES]

32 Concept B Generate all tautomers Impossible to calculate lowest energy tautomer Store all tautomers [Store all protomers]

33 Structure Search

34 Structure Search Exact matches done by flexmatch, SMILES, hashcodes etc. Substructure search Hard to perceive all tautomers for a substructure

35 Approaches to Substructure Search Address problem at registration stage store all tautomers Address problem at search stage enumerate database structures on the fly or enumerate query structure or user takes care specifying query Combine methods Ignore the problem

36 Depiction of Results

37 Depicting Results Display input, registered structure Display matched tautomer good approach if substructure is highlighted Display standard form Let user choose Experimental results match displayed tautomer

38 ChemAxon Approaches Normalize the structure ( generic tautomers ) Allow for tautomers at search time Choose a preferred tautomer Customize preferences in Standardizer

39 ChemAxon Tautomerization Plugin Generates all, dominant and canonical tautomers Calculates canonical tautomer by empirical rules Tries to make canonical tautomer the dominant tautomer (includes pk a filter) Handles dearomatization and stereochemistry

40 Customization Choose dominant tautomer set operating ph set maximum distance (# bonds) of a single proton migration protect structural features aromaticity, charge, stereochemistry, stable functional groups exclude unstable antiaromatic compounds

41 ChemAxon software Stores canonical form, or all tautomers Enumerates query tautomers (as far as possible) Usually displays structure originally input Optionally displays standard tautomer

42 Observations Computational chemistry companies does the ligand match the receptor? ligand preparation pk a algorithms, rules, energetics rigorous approaches Informatics companies does the compound match the patent? building registries and inventories graph theory examples (structures) pragmatic approaches Hybrid companies

43 Acknowledgments All 28 vendors including ChemAxon Jonathan Brecher Geoff Skillman Russ Hillard Keith Taylor