Binding affinity prediction with different force fields: Examination of the linear interaction energy method

A systematic study of the linear interaction energy (LIE) method and the possible dependence of its parameterization on the force field and system (receptor binding site) is reported. We have calculated the binding free energy for nine different ligands in complex with P450cam using three different force fields (Amber95, Gromos87, and OPLS‐AA). The results from these LIE calculations using our earlier parameterization give relative free energies of binding that agree remarkably well with the experimental data. However, the absolute energies are too positive for all three force fields, and it is clear that an additional constant term (γ) is required in this case. Out of five examined LIE models, the same one emerges as the best for all three force fields, and this, in fact, corresponds to our earlier one apart from the addition of the constant γ, which is almost identical for the three force fields. Thus, the present free energy calculations clearly indicate that the coefficients of the LIE method are independent of the force field used. Their relation to solvation free energies is also demonstrated. The only free parameter of the best model is γ, which is found to depend on the hydrophobicity of the binding site. We also attempt to quantify the binding site hydrophobicity of four different proteins which shows that the ordering of γ's for these sites reflects the fraction of hydrophobic surface area. © 2004 Wiley Periodicals, Inc. J Comput Chem 10: 1242–1254, 2004

[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]  M L Lamb,et al.  Estimation of the binding affinities of FKBP12 inhibitors using a linear response method. , 1999, Bioorganic & medicinal chemistry.

[3]  O. Alvarez,et al.  Ion-selective properties of a small ionophore in methanol studied by free energy perturbation simulations , 1992 .

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

[5]  J Andrew McCammon,et al.  Protein flexibility and computer-aided drug design. , 2003, Annual review of pharmacology and toxicology.

[6]  S. Sligar,et al.  The roles of active site hydrogen bonding in cytochrome P-450cam as revealed by site-directed mutagenesis. , 1988, The Journal of biological chemistry.

[7]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

[8]  J. Åqvist,et al.  C2-symmetric inhibitors of Plasmodium falciparum plasmepsin II: synthesis and theoretical predictions. , 2003, Bioorganic & medicinal chemistry.

[9]  J. Åqvist,et al.  Computational analysis of binding of P1 variants to trypsin , 2001, Protein science : a publication of the Protein Society.

[10]  Arieh Warshel,et al.  A local reaction field method for fast evaluation of long‐range electrostatic interactions in molecular simulations , 1992 .

[11]  J. Åqvist,et al.  Q: a molecular dynamics program for free energy calculations and empirical valence bond simulations in biomolecular systems. , 1998, Journal of molecular graphics & modelling.

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

[13]  M G Ford,et al.  Binding constants of neuraminidase inhibitors: An investigation of the linear interaction energy method. , 1999, Journal of medicinal chemistry.

[14]  Heather A. Carlson,et al.  Free energies of solvation in chloroform and water from a linear response approach , 1997 .

[15]  K N Houk,et al.  Quantitative calculations of antibody--antigen binding: steroid--DB3 binding energies by the linear interaction energy method. , 2001, The Journal of organic chemistry.

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

[17]  R C Wade,et al.  Improved binding of cytochrome P450cam substrate analogues designed to fill extra space in the substrate binding pocket. , 1996, Biochemistry.

[18]  J. Åqvist,et al.  Cyclic HIV-1 protease inhibitors derived from mannitol: synthesis, inhibitory potencies, and computational predictions of binding affinities. , 1997, Journal of medicinal chemistry.

[19]  Arieh Warshel,et al.  A surface constrained all‐atom solvent model for effective simulations of polar solutions , 1989 .

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

[21]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[22]  T. Hansson,et al.  On the Validity of Electrostatic Linear Response in Polar Solvents , 1996 .

[23]  S. Sligar,et al.  Temperature jump relaxation kinetics of the P-450cam spin equilibrium. , 1987, Biochemistry.

[24]  R. Raag,et al.  The structural basis for substrate-induced changes in redox potential and spin equilibrium in cytochrome P-450CAM. , 1991, Biochemistry.

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

[26]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[27]  R. Ornstein,et al.  Binding free energy calculations for P450cam-substrate complexes. , 1996, Protein engineering.

[28]  Johan Åqvist,et al.  Computation of affinity and selectivity: Binding of 2,4-diaminopteridine and 2,4-diaminoquinazoline inhibitors to dihydrofolate reductases , 1998, J. Comput. Aided Mol. Des..

[29]  T. Poulos,et al.  High-resolution crystal structure of cytochrome P450cam. , 1987, Journal of molecular biology.

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

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

[32]  R. Ornstein,et al.  Stereoselective Hydroxylation of Norcamphor by Cytochrome P450cam , 1995, The Journal of Biological Chemistry.

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

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

[35]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

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

[37]  György M. Keserü,et al.  A virtual high throughput screen for high affinity cytochrome P450cam substrates. Implications for in silico prediction of drug metabolism , 2001, J. Comput. Aided Mol. Des..

[38]  G. Klebe,et al.  Approaches to the description and prediction of the binding affinity of small-molecule ligands to macromolecular receptors. , 2002, Angewandte Chemie.

[39]  J. Åqvist,et al.  Sensitivity of an empirical affinity scoring function to changes in receptor-ligand complex conformations. , 2001, European journal of pharmaceutical sciences : official journal of the European Federation for Pharmaceutical Sciences.

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

[41]  P A Kollman,et al.  Ranking ligand binding affinities with avidin: a molecular dynamics‐based interaction energy study , 1999, Proteins.

[42]  H. Nar,et al.  Structural basis for inhibition promiscuity of dual specific thrombin and factor Xa blood coagulation inhibitors. , 2001, Structure.

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

[44]  Sherry L. Mowbray,et al.  Sugar Recognition by a Glucose/Galactose Receptor , 1995, The Journal of Biological Chemistry.

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

[46]  R. Raag,et al.  Crystal structures of cytochrome P-450CAM complexed with camphane, thiocamphor, and adamantane: factors controlling P-450 substrate hydroxylation. , 1991, Biochemistry.

[47]  Martin Almlöf,et al.  Free energy calculations and ligand binding. , 2003, Advances in protein chemistry.

[48]  M L Lamb,et al.  Computational approaches to molecular recognition. , 1997, Current opinion in chemical biology.

[49]  D. Osguthorpe,et al.  Structure and energetics of ligand binding to proteins: Escherichia coli dihydrofolate reductase‐trimethoprim, a drug‐receptor system , 1988, Proteins.

[50]  R. Huber,et al.  The Geometry of the Reactive Site and of the Peptide Groups in Trypsin, Trypsinogen and its Complexes with Inhibitors , 1983 .

[51]  James Andrew McCammon,et al.  Ligand-receptor interactions , 1984, Comput. Chem..

[52]  G. Malpeli,et al.  Structure of pig plasma retinol-binding protein at 1.65 A resolution. , 1998, Acta crystallographica. Section D, Biological crystallography.