Assessment of tautomer distribution using the condensed reaction graph approach

We report the first direct QSPR modeling of equilibrium constants of tautomeric transformations (logKT) in different solvents and at different temperatures, which do not require intermediate assessment of acidity (basicity) constants for all tautomeric forms. The key step of the modeling consisted in the merging of two tautomers in one sole molecular graph (“condensed reaction graph”) which enables to compute molecular descriptors characterizing entire equilibrium. The support vector regression method was used to build the models. The training set consisted of 785 transformations belonging to 11 types of tautomeric reactions with equilibrium constants measured in different solvents and at different temperatures. The models obtained perform well both in cross-validation (Q2 = 0.81 RMSE = 0.7 logKT units) and on two external test sets. Benchmarking studies demonstrate that our models outperform results obtained with DFT B3LYP/6-311 ++ G(d,p) and ChemAxon Tautomerizer applicable only in water at room temperature.

[1]  K. Baumann,et al.  Chemoinformatic Classification Methods and their Applicability Domain , 2016, Molecular informatics.

[2]  Igor I. Baskin,et al.  Development of “structure-property” models in nucleophilic substitution reactions involving azides , 2014, Journal of Structural Chemistry.

[3]  D. Horvath,et al.  ISIDA Property‐Labelled Fragment Descriptors , 2010, Molecular informatics.

[4]  Wendy A. Warr,et al.  Tautomerism in chemical information management systems , 2010, J. Comput. Aided Mol. Des..

[5]  Javier Catalán,et al.  A Generalized Solvent Acidity Scale: The Solvatochromism of o-tert-Butylstilbazolium Betaine Dye and Its Homomorph o,o′-Di-tert-butylstilbazolium Betaine Dye† , 1997 .

[6]  J. Andrew Grant,et al.  SAMPL2 and continuum modeling , 2010, J. Comput. Aided Mol. Des..

[7]  David Calkins,et al.  Towards the comprehensive, rapid, and accurate prediction of the favorable tautomeric states of drug-like molecules in aqueous solution , 2010, J. Comput. Aided Mol. Des..

[8]  Gerd Folkers,et al.  Tautomerism in Computer‐Aided Drug Design , 2003, Journal of receptor and signal transduction research.

[9]  Javier Catalán,et al.  Progress towards a generalized solvent polarity scale: The solvatochromism of 2‐(dimethylamino)‐7‐nitrofluorene and its homomorph 2‐fluoro‐7‐nitrofluorene , 1995 .

[10]  Sergei V. Trepalin,et al.  Advanced Exact Structure Searching in Large Databases of Chemical Compounds , 2003, J. Chem. Inf. Comput. Sci..

[11]  J. Tomasi,et al.  The IEF version of the PCM solvation method: an overview of a new method addressed to study molecular solutes at the QM ab initio level , 1999 .

[12]  R. Taft,et al.  The solvatochromic comparison method. 2. The .alpha.-scale of solvent hydrogen-bond donor (HBD) acidities , 1976 .

[13]  Igor I. Baskin,et al.  Structure-reactivity relationships in terms of the condensed graphs of reactions , 2014, Russian Journal of Organic Chemistry.

[14]  R. Taft,et al.  The solvatochromic comparison method. I. The .beta.-scale of solvent hydrogen-bond acceptor (HBA) basicities , 1976 .

[15]  Майре koostaja Тамме.  Таблицы констант скорости и равновесия гетеролитических органических реакций. Доп. том 4, выпуск 1 = Tables of rate and equilibrium constants of heterolytic organic reactions , 1989 .

[16]  J. Riveros,et al.  The Cluster−Continuum Model for the Calculation of the Solvation Free Energy of Ionic Species , 2001 .

[17]  I. Tetko,et al.  ISIDA - Platform for Virtual Screening Based on Fragment and Pharmacophoric Descriptors , 2008 .

[18]  A. Albert,et al.  479. Ionization constants of heterocyclic substances. Part III. Mercapto-derivatives of pyridine, quinoline, and isoquinoline , 1959 .

[19]  Timothy Clark,et al.  Tautomers and reference 3D-structures: the orphans of in silico drug design , 2010, J. Comput. Aided Mol. Des..

[20]  Gilles Marcou,et al.  An Evolutionary Optimizer of libsvm Models , 2014 .

[21]  Maciej Haranczyk,et al.  Quantum Mechanical Energy-Based Screening of Combinatorially Generated Library of Tautomers. TauTGen: A Tautomer Generator Program , 2007, J. Chem. Inf. Model..

[22]  Dragos Horvath,et al.  Models for Identification of Erroneous Atom-to-Atom Mapping of Reactions Performed by Automated Algorithms , 2012, J. Chem. Inf. Model..

[23]  S. Mason.  1002. The tautomerism of N-heteroaromatic hydroxy-compounds. Part II. Ultraviolet spectra , 1957 .

[24]  Marc C. Nicklaus,et al.  Experimental and Chemoinformatics Study of Tautomerism in a Database of Commercially Available Screening Samples , 2016, J. Chem. Inf. Model..

[25]  C. Cramer,et al.  Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. , 2009, The journal of physical chemistry. B.

[26]  Donald G. Truhlar,et al.  Hydride transfer catalyzed by xylose isomerase: Mechanism and quantum effects , 2003, J. Comput. Chem..

[27]  Nicolas Lachiche,et al.  A Representation to Apply Usual Data Mining Techniques to Chemical reactions - Illustration on the Rate Constant of SN2 reactions in water , 2011, Int. J. Artif. Intell. Tools.

[28]  A. Balaban,et al.  Tautomeric Cyclic Thiones. Part IV. The Thione-thiol Equilibrium in some Azoline-2-thiones. , 1969 .

[29]  Roger A. Sayle,et al.  So you think you understand tautomerism? , 2010, J. Comput. Aided Mol. Des..

[30]  R. Bell,et al.  The enol content and acidity of cyclopentanone, cyclohexanone, and acetone in aqueous solution , 1966 .

[31]  Nikolay P. Todorov,et al.  The Influence of Variations of Ligand Protonation and Tautomerism on Protein—Ligand Recognition and Binding Energy Landscape. , 2006 .

[32]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[33]  Loriano Storchi,et al.  Tautomer Enumeration and Stability Prediction for Virtual Screening on Large Chemical Databases. , 2009 .

[34]  S. J. Angyl,et al.  268. The tautomerism of N-hetero-aromatic amines. Part I , 1952 .

[35]  R. Taft,et al.  The solvatochromic comparison method. 6. The .pi.* scale of solvent polarities , 1977 .

[36]  S. Mason.  980. The tautomerism of N-heteroaromatic hydroxy-compounds. Part I. Infrared spectra , 1957 .

[37]  Fiorella Ruggiu,et al.  QSPR modelling -a valuable support in HTS quality control , 2014 .

[38]  Nina Jeliazkova,et al.  Ambit‐Tautomer: An Open Source Tool for Tautomer Generation , 2013, Molecular informatics.

[39]  Donald G. Truhlar,et al.  Implicit Solvation Models: Equilibria, Structure, Spectra, and Dynamics , 1999 .

[40]  D. N. Laikov.  Fast evaluation of density functional exchange-correlation terms using the expansion of the electron density in auxiliary basis sets , 1997 .

[41]  Wolf-Dietrich Ihlenfeldt,et al.  Tautomerism in large databases , 2010, J. Comput. Aided Mol. Des..

[42]  Marc C. Nicklaus,et al.  Enumeration of Ring–Chain Tautomers Based on SMIRKS Rules , 2014, J. Chem. Inf. Model..

[43]  Alexandre Varnek,et al.  Automatized Assessment of Protective Group Reactivity: A Step Toward Big Reaction Data Analysis , 2016, J. Chem. Inf. Model..

[44]  Dimitri N. Laikov,et al.  PRIRODA-04: a quantum-chemical program suite. New possibilities in the study of molecular systems with the application of parallel computing , 2005 .

[45]  Gilles Marcou,et al.  Do Not Hesitate to Use Tversky - and Other Hints for Successful Active Analogue Searches with Feature Count Descriptors , 2013, J. Chem. Inf. Model..

[46]  Chih-Jen Lin,et al.  LIBSVM: A library for support vector machines , 2011, TIST.

[47]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[48]  Igor I. Baskin,et al.  Structure–reactivity relationship in Diels–Alder reactions obtained using the condensed reaction graph approach , 2017, Journal of Structural Chemistry.

[49]  Alexandre Varnek,et al.  Substructural fragments: an universal language to encode reactions, molecular and supramolecular structures , 2005, J. Comput. Aided Mol. Des..

[50]  Fiorella Ruggiu,et al.  Quantitative structure-property relationship modeling: a valuable support in high-throughput screening quality control. , 2014, Analytical chemistry.

[51]  F. Javier Luque,et al.  Performance of the IEF-MST solvation continuum model in the SAMPL2 blind test prediction of hydration and tautomerization free energies , 2010, J. Comput. Aided Mol. Des..

[52]  Donald G. Truhlar,et al.  Prediction of SAMPL2 aqueous solvation free energies and tautomeric ratios using the SM8, SM8AD, and SMD solvation models , 2010, J. Comput. Aided Mol. Des..

[53]  Benjamin Parent,et al.  Fuzzy Tricentric Pharmacophore Fingerprints, 1. Topological Fuzzy Pharmacophore Triplets and Adapted Molecular Similarity Scoring Schemes , 2006, J. Chem. Inf. Model..

[54]  Yvonne C. Martin,et al.  Let’s not forget tautomers , 2009, J. Comput. Aided Mol. Des..

[55]  P. Arnaud,et al.  Binding of the tautomeric forms of isoniazid-NAD adducts to the active site of the Mycobacterium tuberculosis enoyl-ACP reductase (InhA): a theoretical approach. , 2008, Journal of molecular graphics & modelling.

[56]  Arthur Dalby,et al.  Description of several chemical structure file formats used by computer programs developed at Molecular Design Limited , 1992, J. Chem. Inf. Comput. Sci..

[57]  J. Stewart.  Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements , 2007, Journal of molecular modeling.

[58]  Gilles Marcou,et al.  Computational chemogenomics: Is it more than inductive transfer? , 2014, Journal of Computer-Aided Molecular Design.

[59]  Peter A. Kollman,et al.  INSIGHT INTO THE SPECIFICITY OF THYMIDYLATE SYNTHASE FROM MOLECULAR DYNAMICS AND FREE ENERGY PERTURBATION CALCULATIONS , 1995 .

[60]  Jacopo Tomasi,et al.  Quantum Mechanical Continuum Solvation Models , 2005 .

[61]  Paul M. Selzer,et al.  The Impact of Tautomer Forms on Pharmacophore-Based Virtual Screening , 2006, J. Chem. Inf. Model..

[62]  Timur I. Madzhidov,et al.  Structure–reactivity relationship in bimolecular elimination reactions based on the condensed graph of a reaction , 2015, Journal of Structural Chemistry.

[63]  A. Albert,et al.  264. Ionization constants of heterocyclic substances. Part II. Hydroxy-derivatives of nitrogenous six-membered ring-compounds , 1956 .

[64]  S. Mason.  131. The tautomerism of N-heteroaromatic hydroxy-compounds. Part III. Ionisation constants , 1958 .