Importance of accurate charges in molecular docking: Quantum mechanical/molecular mechanical (QM/MM) approach

The extent to which accuracy of electric charges plays a role in protein‐ligand docking is investigated through development of a docking algorithm, which incorporates quantum mechanical/molecular mechanical (QM/MM) calculations. In this algorithm, fixed charges of ligands obtained from force field parameterization are replaced by QM/MM calculations in the protein environment, treating only the ligands as the quantum region. The algorithm is tested on a set of 40 cocrystallized structures taken from the Protein Data Bank (PDB) and provides strong evidence that use of nonfixed charges is important. An algorithm, dubbed “Survival of the Fittest” (SOF) algorithm, is implemented to incorporate QM/MM charge calculations without any prior knowledge of native structures of the complexes. Using an iterative protocol, this algorithm is able in many cases to converge to a nativelike structure in systems where redocking of the ligand using a standard fixed charge force field exhibits nontrivial errors. The results demonstrate that polarization effects can play a significant role in determining the structures of protein‐ligand complexes, and provide a promising start towards the development of more accurate docking methods for lead optimization applications. © 2005 Wiley Periodicals, Inc. J Comput Chem 26: 915–931, 2005

[1]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[2]  Robert W. Hogg,et al.  Crystallographic data of an l-arabinose-binding protein from Escherichia coli , 1974 .

[3]  Richard A. Friesner,et al.  Mixed ab initio QM/MM modeling using frozen orbitals and tests with alanine dipeptide and tetrapeptide , 1999 .

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

[5]  R Abagyan,et al.  High-throughput docking for lead generation. , 2001, Current opinion in chemical biology.

[6]  Thomas Lengauer,et al.  Computational methods for biomolecular docking. , 1996, Current opinion in structural biology.

[7]  Benny G. Johnson,et al.  The performance of a family of density functional methods , 1993 .

[8]  Hege S. Beard,et al.  Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. , 2004, Journal of medicinal chemistry.

[9]  Charles W. Bauschlicher,et al.  A comparison of the accuracy of different functionals , 1995 .

[10]  Matthew P. Repasky,et al.  Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. , 2004, Journal of medicinal chemistry.

[11]  Florante A. Quiocho,et al.  Novel stereospecificity of the L-arabinose-binding protein , 1984, Nature.

[12]  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..

[13]  Charles Collyer,et al.  Inhibition of bovine β-trypsin, human α-thrombin and porcine pancreatic β-kallikrein-B by 4′,6-diamidino-2-phenylindole, 6-amidinoindole and benzamidine: a comparative thermodynamic and X-ray structural study , 1995 .

[14]  R. Friesner,et al.  Frozen orbital QM/MM methods for density functional theory , 2000 .

[15]  Richard D. Taylor,et al.  Improved protein–ligand docking using GOLD , 2003, Proteins.

[16]  F. A. Quiocho,et al.  Substrate specificity and affinity of a protein modulated by bound water molecules , 1989, Nature.

[17]  P. Strevens.  Iii , 1985 .

[18]  Richard A. Friesner,et al.  A mixed quantum mechanics/molecular mechanics (QM/MM) method for large‐scale modeling of chemistry in protein environments , 2000, J. Comput. Chem..

[19]  R. Harrison,et al.  Prediction of new serine proteinase inhibitors , 1994, Nature Structural Biology.

[20]  L. Curtiss,et al.  Assessment of Gaussian-2 and density functional theories for the computation of enthalpies of formation , 1997 .

[21]  I. Kuntz.  Structure-Based Strategies for Drug Design and Discovery , 1992, Science.

[22]  William A. Goddard,et al.  The MPSim‐Dock hierarchical docking algorithm: Application to the eight trypsin inhibitor cocrystals , 2005, J. Comput. Chem..

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

[24]  A. Becke.  A New Mixing of Hartree-Fock and Local Density-Functional Theories , 1993 .

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

[26]  Todd J. A. Ewing,et al.  DOCK 4.0: Search strategies for automated molecular docking of flexible molecule databases , 2001, J. Comput. Aided Mol. Des..