On the calculation of absolute macromolecular binding free energies

The standard framework for calculating the absolute binding free energy of a macromolecular association reaction A + B → AB with an association constant KAB is to equate chemical potentials of the species on the left- and right-hand sides of this reaction and evaluate the chemical potentials from theory. This theory involves (usually hidden) assumptions about what constitutes the bound species, AB, and where the contribution of the solvent appears. We present here an alternative derivation that can be traced back to Bjerrum, in which the expectation value of KAB is obtained directly through the statistical mechanical method of evaluating its ensemble (Boltzmann-weighted) average. The generalized Bjerrum approach more clearly delineates: (i) the different contributions to binding; (ii) the origin of the much-discussed and somewhat controversial association entropy term; and (iii) where the solvent contribution appears. This approach also allows approximations required for practical evaluation of the binding constant in complex macromolecular systems, to be introduced in a well defined way. We provide an example, with application to test cases that illustrate a range of binding behavior.

[1]  B Honig,et al.  Association entropy in adsorption processes. , 2000, Biophysical journal.

[2]  M. Karplus,et al.  Evaluation of the configurational entropy for proteins: application to molecular dynamics simulations of an α-helix , 1984 .

[3]  Yizhak Marcus,et al.  A simple empirical model describing the thermodynamics of hydration of ions of widely varying charges, sizes, and shapes , 1994 .

[4]  Klotz Im,et al.  The thermodynamics of transfer of amides from an apolar to an aqueous solution. , 1969 .

[5]  William H. Press,et al.  Numerical recipes , 1990 .

[6]  K. Sharp,et al.  Energetics of cyclic dipeptide crystal packing and solvation. , 1997, Biophysical journal.

[7]  M. Karplus,et al.  Method for estimating the configurational entropy of macromolecules , 1981 .

[8]  Effects of liquid structures on chemical reactions and conformational changes of non-rigid molecules in condensed phases , 1978 .

[9]  B. Roux,et al.  Statistical mechanical equilibrium theory of selective ion channels. , 1999, Biophysical journal.

[10]  R. Bruccoleri,et al.  Heavy chain position 50 is a determinant of affinity and specificity for the anti-digoxin antibody 26-10. , 1993, The Journal of biological chemistry.

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

[12]  Barry Honig,et al.  A local dielectric constant model for solvation free energies which accounts for solute polarizability , 1992 .

[13]  M. Gilson,et al.  Coordinate Systems and the Calculation of Molecular Properties , 2002 .

[14]  K. Sharp,et al.  Accurate Calculation of Hydration Free Energies Using Macroscopic Solvent Models , 1994 .

[15]  R L Campbell,et al.  26-10 Fab-digoxin complex: affinity and specificity due to surface complementarity. , 1993, Proceedings of the National Academy of Sciences of the United States of America.

[16]  Barry Honig,et al.  Reevaluation of the Born model of ion hydration , 1985 .

[17]  S. Sridharan,et al.  A rapid method for calculating derivatives of solvent accessible surface areas of molecules , 1995, J. Comput. Chem..

[18]  J. Wendoloski,et al.  Structural origins of high-affinity biotin binding to streptavidin. , 1989, Science.

[19]  K. Sharp,et al.  Calculation of HyHel10‐lysozyme binding free energy changes: Effect of ten point mutations , 1998, Proteins.

[20]  M Karplus,et al.  The contribution of vibrational entropy to molecular association. The dimerization of insulin. , 1994, Journal of molecular biology.

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

[22]  P A Kollman,et al.  Absolute and relative binding free energy calculations of the interaction of biotin and its analogs with streptavidin using molecular dynamics/free energy perturbation approaches , 1993, Proteins.

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

[24]  K. Sharp,et al.  Electrostatic contributions to heat capacity changes of DNA-ligand binding. , 1998, Biophysical journal.

[25]  C. Tanford.  Macromolecules , 1994, Nature.

[26]  Faraday Discuss , 1985 .

[27]  T. L. Hill,et al.  Cooperativity Theory in Biochemistry , 1985, Springer Series in Molecular and Cell Biology.

[28]  L M Amzel,et al.  Loss of translational entropy in binding, folding, and catalysis , 1997, Proteins.

[29]  Ray Luo,et al.  Ligand-receptor docking with the Mining Minima optimizer , 2001, J. Comput. Aided Mol. Des..

[30]  R. Hodges,et al.  Contribution of translational and rotational motions to molecular association in aqueous solution. , 2001, Biophysical journal.

[31]  M. Hamermesh.  Group theory and its application to physical problems , 1962 .

[32]  William L. Jorgensen,et al.  Efficient computation of absolute free energies of binding by computer simulations. Application to the methane dimer in water , 1988 .

[33]  M. Gilson,et al.  The statistical-thermodynamic basis for computation of binding affinities: a critical review. , 1997, Biophysical journal.

[34]  A. Tamura,et al.  The entropy cost of protein association. , 1997, Journal of molecular biology.

[35]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

[36]  Michael K. Gilson,et al.  ''Mining minima'': Direct computation of conformational free energy , 1997 .

[37]  Jeremy C. Smith,et al.  Thermodynamic stability of water molecules in the bacteriorhodopsin proton channel: a molecular dynamics free energy perturbation study. , 1996, Biophysical journal.

[38]  Thomas Simonson,et al.  Solvation Free Energies Estimated from Macroscopic Continuum Theory: An Accuracy Assessment , 1994 .

[39]  L. Amzel,et al.  Calculation of entropy changes in biological processes: folding, binding, and oligomerization. , 2000, Methods in enzymology.

[40]  George M. Whitesides,et al.  Estimating the Entropic Cost of Self-Assembly of Multiparticle Hydrogen-Bonded Aggregates Based on the Cyanuric Acid·Melamine Lattice , 1998 .