Quantum Mechanics-Based Properties for 3D-QSAR

We have used a set of four local properties based on semiempirical molecular orbital calculations (electron density (ρ), hydrogen bond donor field (HDF), hydrogen bond acceptor field (HAF), and molecular lipophilicity potential (MLP)) for 3D-QSAR studies to overcome the limitations of the current force field-based molecular interaction fields (MIFs). These properties can be calculated rapidly and are thus amenable to high-throughput industrial applications. Their statistical performance was compared with that of conventional 3D-QSAR approaches using nine data sets (angiotensin converting enzyme inhibitors (ACE), acetylcholinesterase inhibitors (AchE), benzodiazepine receptor ligands (BZR), cyclooxygenase-2 inhibitors (COX2), dihydrofolate reductase inhibitors (DHFR), glycogen phosphorylase b inhibitors (GPB), thermolysin inhibitors (THER), thrombin inhibitors (THR), and serine protease factor Xa inhibitors (fXa)). The 3D-QSAR models generated were tested thoroughly for robustness and predictive ability. The average performance of the quantum mechanical molecular interaction field (QM-MIF) models for the nine data sets is better than that of the conventional force field-based MIFs. In the individual data sets, the QM-MIF models always perform better than, or as well as, the conventional approaches. It is particularly encouraging that the relative performance of the QM-MIF models improves in the external validation. In addition, the models generated showed statistical stability with respect to model building procedure variations such as grid spacing size and grid orientation. QM-MIF contour maps reproduce the features important for ligand binding for the example data set (factor Xa inhibitors), demonstrating the intuitive chemical interpretability of QM-MIFs.

[1]  J. Murray,et al.  Relationships of critical constants and boiling points to computed molecular surface properties , 1993 .

[2]  John H. Van Drie,et al.  Pharmacophore Discovery - Lessons Learned , 2003 .

[3]  Kenneth M Merz,et al.  Semiempirical Comparative Binding Energy Analysis (SE-COMBINE) of a Series of Trypsin Inhibitors. , 2006, Journal of chemical theory and computation.

[4]  P. Diaconis,et al.  Computer-Intensive Methods in Statistics , 1983 .

[5]  David R. Lowis.  Molecular Hologram QSAR , 1999 .

[6]  Peter Politzer,et al.  The complementary roles of molecular surface electrostatic potentials and average local ionization energies with respect to electrophilic processes , 2002 .

[7]  Alexander D. MacKerell,et al.  Computational ligand-based rational design: Role of conformational sampling and force fields in model development. , 2011, MedChemComm.

[8]  M. Akamatsu,et al.  Quantum chemical- and 3-D-QSAR (CoMFA) studies of benzalacetones and 1,1,1-trifluoro-4-phenyl-3-buten-2-ones. , 2002, Bioorganic & medicinal chemistry letters.

[9]  Glen Eugene Kellogg,et al.  HINT: A new method of empirical hydrophobic field calculation for CoMFA , 1991, J. Comput. Aided Mol. Des..

[10]  P. Goodford.  Multivariate characterization of molecules for QSAR analysis , 1996 .

[11]  W. G. Richards,et al.  Computers in drug design , 1993 .

[12]  Connecting traditional QSAR and molecular simulations of papain hydrolysis--importance of charge transfer. , 2005, Bioorganic & medicinal chemistry.

[13]  Timothy Clark,et al.  Predicting the Sites and Energies of Noncovalent Intermolecular Interactions Using Local Properties , 2012, J. Chem. Inf. Model..

[14]  Daniel Rinaldi,et al.  Polarisabilites moléculaires et effet diélectrique de milieu à l'état liquide. Étude théorique de la molécule d'eau et de ses diméres , 1973 .

[15]  D. E. Patterson,et al.  Crossvalidation, Bootstrapping, and Partial Least Squares Compared with Multiple Regression in Conventional QSAR Studies , 1988 .

[16]  G Klebe,et al.  Three-dimensional quantitative structure-activity relationship analyses using comparative molecular field analysis and comparative molecular similarity indices analysis to elucidate selectivity differences of inhibitors binding to trypsin, thrombin, and factor Xa. , 1999, Journal of medicinal chemistry.

[17]  Brendan Twamley,et al.  The nature of halogen...halogen synthons: crystallographic and theoretical studies. , 2006, Chemistry.

[18]  M. E. Beck,et al.  Quantitative structure–activity relations based on quantum theory and wavelet transformations , 2009 .

[19]  Johann Gasteiger,et al.  Multivariate structure‐activity relationships between data from a battery of biological tests and an ensemble of structure descriptors: The PLS method , 1984 .

[20]  Peter Willett,et al.  Evaluation of a novel infrared range vibration-based descriptor (EVA) for QSAR studies. 1. General application , 1997, J. Comput. Aided Mol. Des..

[21]  G E Kellogg,et al.  Allosteric modifiers of hemoglobin. 2. Crystallographically determined binding sites and hydrophobic binding/interaction analysis of novel hemoglobin oxygen effectors. , 1991, Journal of medicinal chemistry.

[22]  J. Murray,et al.  Computed electrostatic potentials and average local ionization energies on the molecular surfaces of some tetracyclines , 2001 .

[23]  Mahmoud A A Ibrahim,et al.  AMBER empirical potential describes the geometry and energy of noncovalent halogen interactions better than advanced semiempirical quantum mechanical method PM6-DH2X. , 2012, The journal of physical chemistry. B.

[24]  Thomas A. Halgren,et al.  Maximally diagonal force constants in dependent angle-bending coordinates. II. Implications for the design of empirical force fields , 1990 .

[25]  Timothy Clark,et al.  Local molecular properties and their use in predicting reactivity , 2003, Journal of molecular modeling.

[26]  Li Zhang,et al.  Quantitative structure–activity relationships for phenyl triazolinones of protoporphyrinogen oxidase inhibitors: A density functional theory study , 2004, J. Comput. Chem..

[27]  F. Pasha,et al.  DFT-based QSAR study of testosterone and its derivatives. , 2004, Bioorganic & medicinal chemistry.

[28]  Maurizio Recanatini,et al.  Target-related applications of first principles quantum chemical methods in drug design. , 2006, Chemical reviews.

[29]  Stefan Güssregen,et al.  Evidence for C-Cl/C-Br...pi interactions as an important contribution to protein-ligand binding affinity. , 2009, Angewandte Chemie.

[30]  B. Berne,et al.  Role of the active-site solvent in the thermodynamics of factor Xa ligand binding. , 2008, Journal of the American Chemical Society.

[31]  Bing Wang,et al.  The role of quantum mechanics in structure-based drug design. , 2007, Drug discovery today.

[32]  Pavel Hobza,et al.  On Extension of the Current Biomolecular Empirical Force Field for the Description of Halogen Bonds. , 2012, Journal of chemical theory and computation.

[33]  Timothy Clark,et al.  σ-Holes: σ-Holes , 2013 .

[34]  Timothy Clark,et al.  3D-QSAR Based on Quantum-Chemical Molecular Fields: Toward an Improved Description of Halogen Interactions , 2012, J. Chem. Inf. Model..

[35]  Patrick Bultinck,et al.  Conceptual DFT properties‐based 3D QSAR: Analysis of inhibitors of the nicotine metabolizing CYP2A6 enzyme , 2009, J. Comput. Chem..

[36]  Maurizio Sironi,et al.  Halogen bonding in ligand-receptor systems in the framework of classical force fields. , 2011, Physical chemistry chemical physics : PCCP.

[37]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[38]  P. Geladi.  Notes on the history and nature of partial least squares (PLS) modelling , 1988 .

[39]  Hans Matter,et al.  Structural requirements for factor Xa inhibition by 3-oxybenzamides with neutral P1 substituents: combining X-ray crystallography, 3D-QSAR, and tailored scoring functions. , 2005, Journal of medicinal chemistry.

[40]  G. Klebe,et al.  Molecular similarity indices in a comparative analysis (CoMSIA) of drug molecules to correlate and predict their biological activity. , 1994, Journal of medicinal chemistry.

[41]  Woody Sherman,et al.  Contribution of Explicit Solvent Effects to the Binding Affinity of Small‐Molecule Inhibitors in Blood Coagulation Factor Serine Proteases , 2011, ChemMedChem.

[42]  Eamonn F. Healy,et al.  Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model , 1985 .

[43]  P. Furet,et al.  3D molecular lipophilicity potential profiles: a new tool in molecular modeling , 1988 .

[44]  F. Diederich,et al.  Interactions with aromatic rings in chemical and biological recognition. , 2003, Angewandte Chemie.

[45]  Kenneth M. Merz,et al.  QMQSAR: Utilization of a semiempirical probe potential in a field‐based QSAR method , 2005, J. Comput. Chem..

[46]  Jonathan M. Goodman,et al.  Hydrogen Bonding and π-Stacking: How Reliable are Force Fields? A Critical Evaluation of Force Field Descriptions of Nonbonded Interactions , 2009, J. Chem. Inf. Model..

[47]  QM and Pharmacophore based 3D-QSAR of MK886 Analogues against mPGES-1 , 2008 .

[48]  Andreas Klamt,et al.  COSMOsar3D: Molecular Field Analysis Based on Local COSMO σ-Profiles , 2012, J. Chem. Inf. Model..

[49]  J. Sutherland,et al.  A comparison of methods for modeling quantitative structure-activity relationships. , 2004, Journal of medicinal chemistry.

[50]  Peter Politzer,et al.  An overview of halogen bonding , 2007, Journal of molecular modeling.

[51]  D. A. Dougherty,et al.  The Cationminus signpi Interaction. , 1997, Chemical reviews.

[52]  Bernd Wendt,et al.  Pushing the boundaries of 3D-QSAR , 2007, J. Comput. Aided Mol. Des..

[53]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[54]  Bernard Testa,et al.  Lipophilicity in Molecular Modeling , 1996, Pharmaceutical Research.

[55]  Amedeo Caflisch,et al.  Quantum mechanical methods for drug design. , 2010, Current topics in medicinal chemistry.

[56]  R. Cramer,et al.  Comparative molecular field analysis (CoMFA). 1. Effect of shape on binding of steroids to carrier proteins. , 1988, Journal of the American Chemical Society.

[57]  S. C. Nyburg.  `Polar flattening': non-spherical effective shapes of atoms in crystals , 1979 .

[58]  P. Bultinck,et al.  3D QSAR based on conceptual DFT molecular fields : antituberculotic activity , 2010 .

[59]  Rajarshi Guha,et al.  Structure-Activity Landscape Index: Identifying and Quantifying Activity Cliffs , 2008, J. Chem. Inf. Model..

[60]  Timothy Clark,et al.  Conformation-Dependent QSPR Models: logPOW , 2011, J. Chem. Inf. Model..

[61]  Timothy Clark,et al.  QSAR and QSPR based solely on surface properties? , 2004, Journal of molecular graphics & modelling.

[62]  Jitender Verma,et al.  3D-QSAR in drug design--a review. , 2010, Current topics in medicinal chemistry.

[63]  Gerhard Klebe,et al.  Development of new hydrogen-bond descriptors and their application to comparative molecular field analyses. , 2002, Journal of medicinal chemistry.

[64]  G. Schürer,et al.  Accurate parametrized variational calculations of the molecular electronic polarizability by NDDO‐based methods , 1999 .

[65]  Bodo Martin,et al.  Additive NDDO-based atomic polarizability model , 2000 .

[66]  S. Wold,et al.  Multivariate Data Analysis in Chemistry , 1984 .

[67]  J. Murray,et al.  Statistical analysis of the molecular surface electrostatic potential: an approach to describing noncovalent interactions in condensed phases , 1998 .

[68]  Timothy Clark,et al.  Halogen bonding: an electrostatically-driven highly directional noncovalent interaction. , 2010, Physical chemistry chemical physics : PCCP.

[69]  Kevin E. Riley,et al.  σ-Holes, π-holes and electrostatically-driven interactions , 2012, Journal of Molecular Modeling.

[70]  S. Wold.  Cross-Validatory Estimation of the Number of Components in Factor and Principal Components Models , 1978 .

[71]  Allan M. Ferguson,et al.  EVA: A new theoretically based molecular descriptor for use in QSAR/QSPR analysis , 1997, J. Comput. Aided Mol. Des..

[72]  G. Klebe,et al.  DrugScore meets CoMFA: adaptation of fields for molecular comparison (AFMoC) or how to tailor knowledge-based pair-potentials to a particular protein. , 2002, Journal of medicinal chemistry.

[73]  E. Audry,et al.  Le potentiel de lipophilie molculaire: application une srie d'amines -adrnolytiques , 1989 .

[74]  Colin McMartin,et al.  QXP: Powerful, rapid computer algorithms for structure-based drug design , 1997, J. Comput. Aided Mol. Des..

[75]  J. Murray,et al.  Characterization of aromatic hydrocarbons by means of average local ionization energies on their molecular surfaces , 2000 .

[76]  Stephen R. Johnson,et al.  The Trouble with QSAR (or How I Learned To Stop Worrying and Embrace Fallacy) , 2008, J. Chem. Inf. Model..

[77]  Gerhard Klebe,et al.  Comparative Molecular Similarity Index Analysis (CoMSIA) to study hydrogen-bonding properties and to score combinatorial libraries , 1999, J. Comput. Aided Mol. Des..

[78]  Peter Politzer,et al.  Average local ionization energies on the molecular surfaces of aromatic systems as guides to chemical reactivity , 1990 .

[79]  QSAR analysis of salicylamide isosteres with the use of quantum chemical molecular descriptors. , 2009, European journal of medicinal chemistry.

[80]  Glen Eugene Kellogg,et al.  The effect of physical organic properties on hydrophobic fields , 1994, J. Comput. Aided Mol. Des..

[81]  J. Fauchère,et al.  Estimating and representing hydrophobicity potential , 1988 .

[82]  Patrick Gaillard,et al.  Molecular Lipophilicity Potential, a tool in 3D QSAR: Method and applications , 1994, J. Comput. Aided Mol. Des..

[83]  Timothy Clark,et al.  Directional Weak Intermolecular Interactions: σ-Hole Bonding , 2010 .

[84]  Patric Schyman,et al.  Treatment of Halogen Bonding in the OPLS-AA Force Field; Application to Potent Anti-HIV Agents. , 2012, Journal of chemical theory and computation.

[85]  Garland R. Marshall,et al.  3D-QSAR of angiotensin-converting enzyme and thermolysin inhibitors: A comparison of CoMFA models based on deduced and experimentally determined active site geometries , 1993 .

[86]  A Tropsha,et al.  Cross-validated R2-guided region selection for comparative molecular field analysis: a simple method to achieve consistent results. , 1995, Journal of medicinal chemistry.

[87]  Li Zhang,et al.  Quantitative Structure-Activity Relationship for Cyclic Imide Derivatives of Protoporphyrinogen Oxidase Inhibitors: A Study of Quantum Chemical Descriptors from Density Functional Theory , 2004, J. Chem. Inf. Model..

[88]  Peter Willett,et al.  Similarity Searching in Files of Three-Dimensional Chemical Structures: Evaluation of the EVA Descriptor and Combination of Rankings Using Data Fusion , 1997, J. Chem. Inf. Comput. Sci..

[89]  Timothy Clark,et al.  The local electron affinity for non-minimal basis sets , 2010, Journal of molecular modeling.

[90]  George D. Purvis,et al.  On the use of isovalued surfaces to determine molecule shape and reaction pathways , 1991, J. Comput. Aided Mol. Des..