Theoretical Analysis of Activity Cliffs among Benzofuranone-Class Pim1 Inhibitors Using the Fragment Molecular Orbital Method with Molecular Mechanics Poisson-Boltzmann Surface Area (FMO+MM-PBSA) Approach

Significant activity changes due to small structural changes (i.e., activity cliffs) of serine/threonine kinase Pim1 inhibitors were studied theoretically using the fragment molecular orbital method with molecular mechanics Poisson-Boltzmann surface area (FMO+MM-PBSA) approach. This methodology enables quantum-chemical calculations for large biomolecules with solvation. In the course of drug discovery targeting Pim1, six benzofuranone-class inhibitors were found to differ only in the position of the indole-ring nitrogen atom. By comparing the various qualities of complex structures based on X-ray, classical molecular mechanics (MM)-optimized, and quantum/molecular mechanics (QM/MM)-optimized structures, we found that the QM/MM-optimized structures provided the best correlation (R2 = 0.85) between pIC50 and the calculated FMO+MM-PBSA binding energy. Combining the classical solvation energy with the QM binding energy was important to increase the correlation. In addition, decomposition of the interaction energy into various physicochemical components by pair interaction energy decomposition analysis suggested that CH-π and electrostatic interactions mainly caused the activity differences.

[1]  Yuka Koyama,et al.  Affinity of HIV-1 antibody 2G12 with monosaccharides: A theoretical study based on explicit and implicit water models , 2014, Comput. Biol. Chem..

[2]  R. Reeves,et al.  Why target PIM1 for cancer diagnosis and treatment? , 2010, Future oncology.

[3]  Kaori Fukuzawa,et al.  Charge clamps of lysines and hydrogen bonds play key roles in the mechanism to fix helix 12 in the agonist and antagonist positions of estrogen receptor α: intramolecular interactions studied by the ab initio fragment molecular orbital method. , 2014, The journal of physical chemistry. B.

[4]  S. Dhanasekaran,et al.  Delineation of prognostic biomarkers in prostate cancer , 2001, Nature.

[5]  Jaime Pérez-Villanueva,et al.  Activity cliffs and activity cliff generators based on chemotype-related activity landscapes , 2015, Molecular Diversity.

[6]  Kotoko Nakata,et al.  Ab initio quantum mechanical study of the binding energies of human estrogen receptor α with its ligands: An application of fragment molecular orbital method , 2005, J. Comput. Chem..

[7]  K. Kitaura,et al.  Fragment molecular orbital method: an approximate computational method for large molecules , 1999 .

[8]  K. Kitaura,et al.  Energy decomposition analysis in solution based on the fragment molecular orbital method. , 2012, The journal of physical chemistry. A.

[9]  J. Tomasi,et al.  Quantum mechanical continuum solvation models. , 2005, Chemical reviews.

[10]  Wei Wang,et al.  Characterization of domain-peptide interaction interface: prediction of SH3 domain-mediated protein-protein interaction network in yeast by generic structure-based models. , 2012, Journal of proteome research.

[11]  G Klebe,et al.  Impact of ligand and protein desolvation on ligand binding to the S1 pocket of thrombin. , 2012, Journal of molecular biology.

[12]  Youyong Li,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 4. Accuracies of MM/PBSA and MM/GBSA methodologies evaluated by various simulation protocols using PDBbind data set. , 2014, Physical chemistry chemical physics : PCCP.

[13]  S. Yokoyama,et al.  A novel Pim-1 kinase inhibitor targeting residues that bind the substrate peptide. , 2012, Journal of molecular biology.

[14]  Kazuo Kitaura,et al.  Pair interaction energy decomposition analysis , 2007, J. Comput. Chem..

[15]  Kaori Fukuzawa,et al.  Implementation of Pair Interaction Energy DecompositionAnalysis and Its Applications to Protein-Ligand Systems , 2015 .

[16]  M. Nishio,et al.  CH/pi interactions as demonstrated in the crystal structure of guanine-nucleotide binding proteins, Src homology-2 domains and human growth hormone in complex with their specific ligands. , 1998, Bioorganic & medicinal chemistry.

[17]  J. Bajorath,et al.  Follow up: Advancing the activity cliff concept, part II , 2014, F1000Research.

[18]  Yasuhiro Shigemitsu.  Quantum chemical study on molecular‐level affinity of DJ‐1 binding compounds , 2013 .

[19]  Kaori Fukuzawa,et al.  Explicit solvation modulates intra- and inter-molecular interactions within DNA: Electronic aspects revealed by the ab initio fragment molecular orbital (FMO) method , 2015 .

[20]  Kaori Fukuzawa,et al.  Hydration of ligands of influenza virus neuraminidase studied by the fragment molecular orbital method. , 2016, Journal of molecular graphics & modelling.

[21]  Toshio Fujita,et al.  Novel Quantitative Structure-Activity Studies of HIV-1 Protease Inhibitors of the Cyclic Urea Type Using Descriptors Derived from Molecular Dynamics and Molecular Orbital Calculations , 2009 .

[22]  Kaori Fukuzawa,et al.  Fragment molecular orbital method: use of approximate electrostatic potential , 2002 .

[23]  T. Kigawa,et al.  Improving cell-free protein synthesis for stable-isotope labeling , 2007, Journal of biomolecular NMR.

[24]  Kaori Fukuzawa,et al.  Molecular interactions between estrogen receptor and its ligand studied by the ab initio fragment molecular orbital method. , 2006, The journal of physical chemistry. B.

[25]  Chang-Guo Zhan,et al.  Binding free energies for nicotine analogs inhibiting cytochrome P450 2A6 by a combined use of molecular dynamics simulations and QM/MM-PBSA calculations. , 2014, Bioorganic & medicinal chemistry.

[26]  S. Yokoyama,et al.  Rational evolution of a novel type of potent and selective proviral integration site in Moloney murine leukemia virus kinase 1 (PIM1) inhibitor from a screening-hit compound. , 2012, Journal of medicinal chemistry.

[27]  Araz Jakalian,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method , 2000 .

[28]  Alexander Rohe,et al.  Application of Docking and QM/MM-GBSA Rescoring to Screen for Novel Myt1 Kinase Inhibitors , 2014, J. Chem. Inf. Model..

[29]  Randy J. Read,et al.  Acta Crystallographica Section D Biological , 2003 .

[30]  Hiroshi Chuman,et al.  QSAR Study of Cyclic Urea Type HIV‐1 PR Inhibitors Using Ab Initio MO Calculation of Their Complex Structures with HIV‐1 PR , 2008 .

[31]  博文 渡邉,et al.  Comparison of binding affinity evaluations for FKBP ligands with state-of-the-art computational methods: FMO, QM/MM, MM-PB/SA and MP-CAFEE approaches , 2010 .

[32]  Kevin Cowtan,et al.  research papers Acta Crystallographica Section D Biological , 2005 .

[33]  G. Basarab,et al.  Structure-based design of potent inhibitors of scytalone dehydratase: displacement of a water molecule from the active site. , 1998, Biochemistry.

[34]  S. Knapp,et al.  Structural analysis identifies imidazo[1,2-b]pyridazines as PIM kinase inhibitors with in vitro antileukemic activity. , 2007, Cancer research.

[35]  Kazuo Kitaura,et al.  Molecular recognition mechanism of FK506 binding protein: An all‐electron fragment molecular orbital study , 2007, Proteins.

[36]  Hui Li,et al.  The polarizable continuum model (PCM) interfaced with the fragment molecular orbital method (FMO) , 2006, J. Comput. Chem..

[37]  Kaori Fukuzawa,et al.  Partial geometry optimization with FMO-MP2 gradient: Application to TrpCage , 2012 .

[38]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[39]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[40]  J. Bajorath,et al.  Advancing the activity cliff concept , 2013 .

[41]  Kazuo Kitaura,et al.  Structure-based design of novel potent protein kinase CK2 (CK2) inhibitors with phenyl-azole scaffolds. , 2012, Journal of medicinal chemistry.

[42]  Jürgen Bajorath,et al.  Recent progress in understanding activity cliffs and their utility in medicinal chemistry. , 2014, Journal of medicinal chemistry.

[43]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[44]  Kenneth M. Merz,et al.  Fully Quantum Mechanical Description of Proteins in Solution. Combining Linear Scaling Quantum Mechanical Methodologies with the Poisson−Boltzmann Equation , 1999 .

[45]  Randy J. Read,et al.  Overview of the CCP4 suite and current developments , 2011, Acta crystallographica. Section D, Biological crystallography.

[46]  Jürgen Bajorath,et al.  Exploring activity cliffs in medicinal chemistry. , 2012, Journal of medicinal chemistry.

[47]  Takanori Kigawa,et al.  Preparation of Escherichia coli cell extract for highly productive cell-free protein expression , 2004, Journal of Structural and Functional Genomics.

[48]  Yuji Mochizuki,et al.  Acceleration of fragment molecular orbital calculations with Cholesky decomposition approach , 2010 .

[49]  Tingjun Hou,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 5. Improved docking performance using high solute dielectric constant MM/GBSA and MM/PBSA rescoring. , 2014, Physical chemistry chemical physics : PCCP.

[50]  Kaori Fukuzawa,et al.  Three- and four-body corrected fragment molecular orbital calculations with a novel subdividing fragmentation method applicable to structure-based drug design. , 2013, Journal of molecular graphics & modelling.

[51]  Stefan Knapp,et al.  PIM serine/threonine kinases in the pathogenesis and therapy of hematologic malignancies and solid cancers , 2010, Haematologica.

[52]  Yuto Komeiji,et al.  Electron-correlated fragment-molecular-orbital calculations for biomolecular and nano systems. , 2014, Physical chemistry chemical physics : PCCP.

[53]  Hiroshi Chuman,et al.  Correlation Analyses on Binding Affinity of Sialic Acid Analogues with Influenza Virus Neuraminidase-1 Using ab Initio MO Calculations on Their Complex Structures , 2010, J. Chem. Inf. Model..

[54]  M. Nishio,et al.  The CH/π hydrogen bond in chemistry. Conformation, supramolecules, optical resolution and interactions involving carbohydrates. , 2011, Physical chemistry chemical physics : PCCP.

[55]  Christopher I. Bayly,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation , 2002, J. Comput. Chem..

[56]  M. Fleming,et al.  Pim-1 Ligand-bound Structures Reveal the Mechanism of Serine/Threonine Kinase Inhibition by LY294002* , 2005, Journal of Biological Chemistry.

[57]  D. Boschelli,et al.  4-Anilino-6,7-dialkoxyquinoline-3-carbonitrile inhibitors of epidermal growth factor receptor kinase and their bioisosteric relationship to the 4-anilino-6,7-dialkoxyquinazoline inhibitors. , 2000, Journal of medicinal chemistry.

[58]  T. Nakano,et al.  Fragment interaction analysis based on local MP2 , 2007 .

[59]  Yuto Komeiji,et al.  Change in a protein's electronic structure induced by an explicit solvent: An ab initio fragment molecular orbital study of ubiquitin , 2007, J. Comput. Chem..

[60]  Kaori Fukuzawa,et al.  Development of the four-body corrected fragment molecular orbital (FMO4) method , 2012 .

[61]  Mark Whittaker,et al.  Prediction of cyclin-dependent kinase 2 inhibitor potency using the fragment molecular orbital method , 2011, J. Cheminformatics.

[62]  Kohji Itoh,et al.  Correlation Analyses on Binding Affinity of Sialic Acid Analogues and Anti-Influenza Drugs with Human Neuraminidase Using ab Initio MO Calculations on Their Complex Structures - LERE-QSAR Analysis (IV) , 2011, J. Chem. Inf. Model..

[63]  Shigenori Tanaka,et al.  Incorporation of solvation effects into the fragment molecular orbital calculations with the Poisson–Boltzmann equation , 2010 .

[64]  Travis V. Harris,et al.  The ONIOM Method and Its Applications. , 2015, Chemical reviews.

[65]  S. Yokoyama,et al.  Flexibility of the P-loop of Pim-1 kinase: observation of a novel conformation induced by interaction with an inhibitor. , 2012, Acta crystallographica. Section F, Structural biology and crystallization communications.

[66]  P. Kollman,et al.  Continuum Solvent Studies of the Stability of DNA, RNA, and Phosphoramidate−DNA Helices , 1998 .

[67]  Yutaka Akiyama,et al.  Fragment molecular orbital method: application to polypeptides , 2000 .

[68]  Youyong Li,et al.  Assessing the performance of MM/PBSA and MM/GBSA methods. 3. The impact of force fields and ligand charge models. , 2013, The journal of physical chemistry. B.

[69]  D. Moustakas,et al.  The Necessary Nitrogen Atom: A Versatile High-Impact Design Element for Multiparameter Optimization. , 2017, Journal of medicinal chemistry.

[70]  D. Leonidas,et al.  3′‐Axial CH2OH Substitution on Glucopyranose does not Increase Glycogen Phosphorylase Inhibitory Potency. QM/MM‐PBSA Calculations Suggest Why , 2012, Chemical biology & drug design.

[71]  F. Sigaux,et al.  The human protooncogene product p33pim is expressed during fetal hematopoiesis and in diverse leukemias. , 1989, Proceedings of the National Academy of Sciences of the United States of America.