QM/MM linear response method distinguishes ligand affinities for closely related metalloproteins
- 2 July 2007
- journal article
- research article
- Published by Wiley in Proteins-Structure Function and Bioinformatics
- Vol. 69 (2), 326-339
- https://doi.org/10.1002/prot.21500
Abstract
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.Keywords
This publication has 74 references indexed in Scilit:
- Free Energy, Entropy, and Induced Fit in Host−Guest Recognition: Calculations with the Second-Generation Mining Minima AlgorithmJournal of the American Chemical Society, 2004
- UCSF Chimera—A visualization system for exploratory research and analysisJournal of Computational Chemistry, 2004
- Inclusion of the ligand field contribution in a polarizable molecular mechanics: SIBFA-LFJournal of Computational Chemistry, 2003
- Knowledge-based scoring function to predict protein-ligand interactionsJournal of Molecular Biology, 2000
- The Protein Data BankNucleic Acids Research, 2000
- Crystal structure of the stromelysin catalytic domain at 2.0 Å resolution: inhibitor-induced conformational changesJournal of Molecular Biology, 1999
- Development and validation of a genetic algorithm for flexible docking 1 1Edited by F. E. CohenJournal of Molecular Biology, 1997
- A Fast Flexible Docking Method using an Incremental Construction AlgorithmJournal of Molecular Biology, 1996
- Solvation energy in protein folding and bindingNature, 1986
- A geometric approach to macromolecule-ligand interactionsJournal of Molecular Biology, 1982