QM/MM linear response method distinguishes ligand affinities for closely related metalloproteins

Design of selective ligands for closely related targets is becoming one of the most important tasks in the drug development. New tools, more precise than fast scoring functions and less demanding than sophisticated Free Energy Perturbation methods, are necessary to help accomplish this goal. The methods of intermediate complexity, characterizing individual contributions to the binding energy, have been an area of intense research in the past few years. Our recently developed quantum mechanical/molecular mechanical (QM/MM) modification of the Linear Response (LR) method describes the binding free energies as the sum of empirically weighted contributions of the QM/MM interaction energies and solvent‐accessible surface areas for the time‐averaged structures of hydrated complexes, obtained by molecular dynamics (MD) simulations. The method was applied to published data on 27 inhibitors of matrix metalloproteinase‐3 (MMP‐3). The two descriptors explained 90% of variance in the inhibition constants with RMSE of 0.245 log units. The QM/MM treatment is indispensable for characterization of the systems lacking suitable force‐field expressions. In this case, it provided characteristics of H‐bonds of the inhibitors to Glu202, charges of binding site atoms, and accurate coordination geometries of the ligands to catalytic zinc. The geometries were constrained during the MD simulations, which characterized conformational flexibility of the complexes and helped in the elucidation of the binding differences for related compounds. A comparison of the presented QM/MM LR results with those previously published for inhibition of MMP‐9 by the same set of ligands showed that the QM/MM LR approach was able to distinguish subtle differences in binding affinities for MMP‐3 and MMP‐9, which did not exceed one order of magnitude. This precision level makes the approach a useful tool for design of selective ligands to similar targets, because the results can be safely extrapolated to maximize selectivity. Proteins 2007. © 2007 Wiley‐Liss, Inc.

[1]  B. Nordén,et al.  Computational modelling of inhibitor binding to human thrombin. , 2001, European journal of pharmaceutical sciences : official journal of the European Federation for Pharmaceutical Sciences.

[2]  W. L. Jorgensen,et al.  Binding affinities for sulfonamide inhibitors with human thrombin using Monte Carlo simulations with a linear response method. , 1997, Journal of medicinal chemistry.

[3]  J. Aqvist,et al.  A new method for predicting binding affinity in computer-aided drug design. , 1994, Protein engineering.

[4]  Thomas Lengauer,et al.  A fast flexible docking method using an incremental construction algorithm. , 1996, Journal of molecular biology.

[5]  K V Damodaran,et al.  Binding preferences of hydroxamate inhibitors of the matrix metalloproteinase human fibroblast collagenase. , 1999, Journal of medicinal chemistry.

[6]  R. S. Mulliken.  Electronic Population Analysis on LCAO–MO Molecular Wave Functions. I , 1955 .

[7]  Soumyendu Raha,et al.  Similarity of Binding Sites of Human Matrix Metalloproteinases*[boxs] , 2004, Journal of Biological Chemistry.

[8]  J. Åqvist,et al.  Ligand binding affinities from MD simulations. , 2002, Accounts of chemical research.

[9]  Michael K. Gilson,et al.  Tork: Conformational analysis method for molecules and complexes , 2003, J. Comput. Chem..

[10]  Johan Åqvist,et al.  Binding affinity prediction with different force fields: Examination of the linear interaction energy method , 2004, J. Comput. Chem..

[11]  Angelo Vedani,et al.  A new force field for modeling metalloproteins , 1990 .

[12]  Yuan-Ping Pang,et al.  Proton Dissociation Energies of Zinc-Coordinated Hydroxamic Acids and Their Relative Affinities for Zinc: Insight into Design of Inhibitors of Zinc-Containing Proteinases , 2000 .

[13]  P A Kollman,et al.  What determines the van der Waals coefficient β in the LIE (linear interaction energy) method to estimate binding free energies using molecular dynamics simulations? , 1999, Proteins.

[14]  A. Skalka,et al.  The retroviral enzymes. , 1994, Annual review of biochemistry.

[15]  Akash Khandelwal,et al.  A combination of docking, QM/MM methods, and MD simulation for binding affinity estimation of metalloprotein ligands. , 2005, Journal of medicinal chemistry.

[16]  A. D. McLachlan,et al.  Solvation energy in protein folding and binding , 1986, Nature.

[17]  Hans-Joachim Böhm,et al.  The development of a simple empirical scoring function to estimate the binding constant for a protein-ligand complex of known three-dimensional structure , 1994, J. Comput. Aided Mol. Des..

[18]  W. R. Wadt,et al.  Ab initio effective core potentials for molecular calculations. Potentials for main group elements Na to Bi , 1985 .

[19]  Xiaojie Xu,et al.  Binding Affinities for a Series of Selective Inhibitors of Gelatinase-A Using Molecular Dynamics with a Linear Interaction Energy Approach , 2001 .

[20]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[21]  V. Lukacova,et al.  A comparison of the binding sites of matrix metalloproteinases and tumor necrosis factor-alpha converting enzyme: implications for selectivity. , 2005, Journal of medicinal chemistry.

[22]  W. Bode,et al.  Structure determination and analysis of human neutrophil collagenase complexed with a hydroxamate inhibitor. , 1995, Biochemistry.

[23]  Y. Martin,et al.  A general and fast scoring function for protein-ligand interactions: a simplified potential approach. , 1999, Journal of medicinal chemistry.

[24]  Robert J. Deeth,et al.  The ligand field molecular mechanics model and the stereoelectronic effects of d and s electrons , 2001 .

[25]  M. Browner,et al.  Matrilysin-inhibitor complexes: common themes among metalloproteases. , 1996, Biochemistry.

[26]  K. Merz,et al.  A quantum mechanics-based scoring function: study of zinc ion-mediated ligand binding. , 2004, Journal of the American Chemical Society.

[27]  R. Huber,et al.  X-ray structures of human neutrophil collagenase complexed with peptide hydroxamate and peptide thiol inhibitors. Implications for substrate binding and rational drug design. , 1995, European journal of biochemistry.

[28]  Johan Åqvist,et al.  Calculation of absolute binding free energies for charged ligands and effects of long‐range electrostatic interactions , 1996, J. Comput. Chem..

[29]  Rebecca C. Wade,et al.  COMPUTATIONAL ALCHEMY TO CALCULATE ABSOLUTE PROTEIN-LIGAND BINDING FREE ENERGY , 1998 .

[30]  Michael Bräuer,et al.  Molecular mechanics for zinc complexes with fluctuating atomic charges , 2001 .

[31]  Peter A. Kollman,et al.  Free energy calculation methods: A theoretical and empirical comparison of numerical errors and a new method qualitative estimates of free energy changes , 1997, J. Comput. Chem..

[32]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[33]  T. Hansson,et al.  Estimation of binding free energies for HIV proteinase inhibitors by molecular dynamics simulations. , 1995, Protein engineering.

[34]  Peter A. Kollman,et al.  Free energy calculation methods: A theoretical and empirical comparison of numerical errors and a new method qualitative estimates of free energy changes , 1997, J. Comput. Chem..

[35]  Hirosato Kondo,et al.  New type of metalloproteinase inhibitor: design and synthesis of new phosphonamide-based hydroxamic acids. , 2002, Journal of medicinal chemistry.

[36]  M L Lamb,et al.  Estimation of the binding affinities of FKBP12 inhibitors using a linear response method. , 1999, Bioorganic & medicinal chemistry.

[37]  Ajay N. Jain.  Scoring noncovalent protein-ligand interactions: A continuous differentiable function tuned to compute binding affinities , 1996, J. Comput. Aided Mol. Des..

[38]  M Karplus,et al.  Improving the accuracy of protein pKa calculations: Conformational averaging versus the average structure , 1998, Proteins.

[39]  A. Warshel,et al.  Examining methods for calculations of binding free energies: LRA, LIE, PDLD‐LRA, and PDLD/S‐LRA calculations of ligands binding to an HIV protease , 2000, Proteins.

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

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

[42]  Thomas Lengauer,et al.  Evaluation of the FLEXX incremental construction algorithm for protein–ligand docking , 1999, Proteins.

[43]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[44]  Irene Luque,et al.  Structural parameterization of the binding enthalpy of small ligands , 2002, Proteins.

[45]  P. Kollman,et al.  Binding of a diverse set of ligands to avidin and streptavidin: an accurate quantitative prediction of their relative affinities by a combination of molecular mechanics and continuum solvent models. , 2000, Journal of medicinal chemistry.

[46]  G. V. Paolini,et al.  Empirical scoring functions: I. The development of a fast empirical scoring function to estimate the binding affinity of ligands in receptor complexes , 1997, J. Comput. Aided Mol. Des..

[47]  Angelo Vedani,et al.  YETI: An interactive molecular mechanics program for small‐molecule protein complexes , 1988 .

[48]  Olivier Michielin,et al.  Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors , 2003, J. Comput. Aided Mol. Des..

[49]  Robert J. Deeth,et al.  Molecular Mechanics for Coordination Complexes: The Impact of Adding d-Electron Stabilization Energies , 1995 .

[50]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

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

[52]  Kenneth M. Merz,et al.  Force Field Design for Metalloproteins , 1991 .

[53]  M. Gilson,et al.  Free energy, entropy, and induced fit in host-guest recognition: calculations with the second-generation mining minima algorithm. , 2004, Journal of the American Chemical Society.

[54]  U. Ryde,et al.  Carboxylate binding modes in zinc proteins: a theoretical study. , 1999, Biophysical journal.

[55]  Pedro Alexandrino Fernandes,et al.  Protein–ligand docking: Current status and future challenges , 2006, Proteins.

[56]  E. Shakhnovich,et al.  SMall Molecule Growth 2001 (SMoG2001): an improved knowledge-based scoring function for protein-ligand interactions. , 2002, Journal of medicinal chemistry.

[57]  U. Ryde,et al.  Molecular dynamics simulations of alcohol dehydrogenase with a four‐ or five‐coordinate catalytic zinc ion , 1995, Proteins.

[58]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[59]  S. Balaz,et al.  A practical approach to docking of zinc metalloproteinase inhibitors. , 2004, Journal of molecular graphics & modelling.

[60]  B. Barnett,et al.  Crystal structure of the stromelysin catalytic domain at 2.0 A resolution: inhibitor-induced conformational changes. , 1999, Journal of molecular biology.

[61]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

[62]  M. Browner,et al.  Crystal structures of matrilysin-inhibitor complexes , 1995 .

[63]  R. Cramer,et al.  Validation of the general purpose tripos 5.2 force field , 1989 .

[64]  W. R. Wadt,et al.  Ab initio effective core potentials for molecular calculations , 1984 .

[65]  Stephen R. Johnson,et al.  A Novel Method to Simulate the Hydrophobic Effect and Its Application to the Ranking of Host/Guest Complexes , 2006, J. Chem. Inf. Model..

[66]  Y. Pang,et al.  Novel Stable Configurations and Tautomers of the Neutral and Deprotonated Hydroxamic Acids Predicted from High-Level ab Initio Calculations , 1999 .

[67]  Johan Åqvist,et al.  Ligand binding affinity prediction by linear interaction energy methods , 1998, J. Comput. Aided Mol. Des..

[68]  J. Åqvist,et al.  The linear interaction energy method for predicting ligand binding free energies. , 2001, Combinatorial chemistry & high throughput screening.

[69]  W. L. Jorgensen,et al.  AN EXTENDED LINEAR RESPONSE METHOD FOR DETERMINING FREE ENERGIES OF HYDRATION , 1995 .

[70]  J M Blaney,et al.  A geometric approach to macromolecule-ligand interactions. , 1982, Journal of molecular biology.

[71]  Nohad Gresh,et al.  Inclusion of the ligand field contribution in a polarizable molecular mechanics: SIBFA‐LF , 2003, J. Comput. Chem..

[72]  G. Klebe,et al.  Knowledge-based scoring function to predict protein-ligand interactions. , 2000, Journal of molecular biology.