Docking of flexible ligands to flexible receptors in solution by molecular dynamics simulation

In this paper, a method of simulating the docking of small flexible ligands to flexible receptors in water is reported. The method is based on molecular dynamics simulations and is an extension of an algorithm previously reported by Di Nola et al. (Di Nola et al., Proteins 1994;19:174–182). The method allows a fast exploration of the receptor surface, using a high temperature of the center of mass translational motion, while the ligand internal motions, the solvent, and the receptor are simulated at room temperature. In addition, the method allows a fast center of mass motion of the ligand, even in solution. The dampening effect of the solvent can be overcome by applying different weights to the interactions between system subsets (solvent, receptor, and ligand). Specific ligand–receptor distances have been used to compare the results of the simulations with the crystal structure. The method is applied, as a test system, to the docking of the phosphocholine to the immunoglobulin McPC603. The results show the similarity of structure between the complex in solution and in the crystal. Proteins 1999;35:153–162. © 1999 Wiley‐Liss, Inc.

[1]  Huajun Wang.  Grid‐search molecular accessible surface algorithm for solving the protein docking problem , 1991 .

[2]  C. Brooks.  Computer simulation of liquids , 1989 .

[3]  A. Leach,et al.  Ligand docking to proteins with discrete side-chain flexibility. , 1994, Journal of molecular biology.

[4]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[5]  I. Kuntz,et al.  Molecular docking to ensembles of protein structures. , 1997, Journal of molecular biology.

[6]  Kevin P. Clark,et al.  Flexible ligand docking without parameter adjustment across four ligand–receptor complexes , 1995, J. Comput. Chem..

[7]  E. Katchalski‐Katzir,et al.  Molecular surface recognition: determination of geometric fit between proteins and their ligands by correlation techniques. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[8]  Charles DeLisi,et al.  Peptide docking using dynamic programming , 1996, J. Comput. Chem..

[9]  The joys and vicissitudes of protein science , 1993, Protein science : a publication of the Protein Society.

[10]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[11]  I. Kuntz,et al.  Using shape complementarity as an initial screen in designing ligands for a receptor binding site of known three-dimensional structure. , 1988, Journal of medicinal chemistry.

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

[13]  M. Mizutani,et al.  Rational automatic search method for stable docking models of protein and ligand. , 1994, Journal of molecular biology.

[14]  H J Berendsen,et al.  Molecular dynamics simulation of the docking of substrates to proteins , 1994, Proteins.

[15]  Ron Elber,et al.  MOIL-View - A Program for Visualization of Structure and Dynamics of Biomolecules and STO - A Program for Computing Stochastic Paths , 1995 .

[16]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[17]  T. Stouch,et al.  Affinity and specificity of serine endopeptidase-protein inhibitor interactions. Empirical free energy calculations based on X-ray crystallographic structures. , 1993, Journal of molecular biology.

[18]  D. Goodsell,et al.  Automated docking of substrates to proteins by simulated annealing , 1990, Proteins.

[19]  M. Lewis,et al.  Calculation of the free energy of association for protein complexes , 1992, Protein science : a publication of the Protein Society.

[20]  Pieter F. W. Stouten,et al.  A molecular mechanics/grid method for evaluation of ligand–receptor interactions , 1995, J. Comput. Chem..

[21]  G. Cohen,et al.  The three-dimensional structure of a phosphorylcholine-binding mouse immunoglobulin Fab and the nature of the antigen binding site. , 1974, Proceedings of the National Academy of Sciences of the United States of America.

[22]  P. Goodford.  A computational procedure for determining energetically favorable binding sites on biologically important macromolecules. , 1985, Journal of medicinal chemistry.

[23]  A. N. Jain,et al.  Hammerhead: fast, fully automated docking of flexible ligands to protein binding sites. , 1996, Chemistry & biology.

[24]  I. Kuntz,et al.  Protein docking and complementarity. , 1991, Journal of molecular biology.

[25]  R. Bruccoleri,et al.  On the attribution of binding energy in antigen-antibody complexes McPC 603, D1.3, and HyHEL-5. , 1989, Biochemistry.

[26]  I. Kuntz,et al.  Conformational analysis of flexible ligands in macromolecular receptor sites , 1992 .

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

[28]  A Aitken,et al.  Expression and structural analysis of 14-3-3 proteins. , 1995, Journal of molecular biology.

[29]  D S Goodsell,et al.  Automated docking of flexible ligands: Applications of autodock , 1996, Journal of molecular recognition : JMR.

[30]  Amedeo Caflisch,et al.  Docking small ligands in flexible binding sites , 1998, J. Comput. Chem..

[31]  M J Sternberg,et al.  New algorithm to model protein-protein recognition based on surface complementarity. Applications to antibody-antigen docking. , 1992, Journal of molecular biology.

[32]  R. Glen,et al.  Molecular recognition of receptor sites using a genetic algorithm with a description of desolvation. , 1995, Journal of molecular biology.

[33]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[34]  J. Scott Dixon,et al.  Flexible ligand docking using a genetic algorithm , 1995, J. Comput. Aided Mol. Des..