Importance of ligand reorganization free energy in protein-ligand binding-affinity prediction.

Accurate prediction of the binding affinities of small-molecule ligands to their biological targets is fundamental for structure-based drug design but remains a very challenging task. In this paper, we have performed computational studies to predict the binding models of 31 small-molecule Smac (the second mitochondria-derived activator of caspase) mimetics to their target, the XIAP (X-linked inhibitor of apoptosis) protein, and their binding affinities. Our results showed that computational docking was able to reliably predict the binding models, as confirmed by experimentally determined crystal structures of some Smac mimetics complexed with XIAP. However, all the computational methods we have tested, including an empirical scoring function, two knowledge-based scoring functions, and MM-GBSA (molecular mechanics and generalized Born surface area), yield poor to modest prediction for binding affinities. The linear correlation coefficient (r(2)) value between the predicted affinities and the experimentally determined affinities was found to be between 0.21 and 0.36. Inclusion of ensemble protein-ligand conformations obtained from molecular dynamic simulations did not significantly improve the prediction. However, major improvement was achieved when the free-energy change for ligands between their free- and bound-states, or "ligand-reorganization free energy", was included in the MM-GBSA calculation, and the r(2) value increased from 0.36 to 0.66. The prediction was validated using 10 additional Smac mimetics designed and evaluated by an independent group. This study demonstrates that ligand reorganization free energy plays an important role in the overall binding free energy between Smac mimetics and XIAP. This term should be evaluated for other ligand-protein systems and included in the development of new scoring functions. To our best knowledge, this is the first computational study to demonstrate the importance of ligand reorganization free energy for the prediction of protein-ligand binding free energy.

[1]  U. Ryde,et al.  Molecular dynamics simulations of alcohol dehydrogenase with a four‐ or five‐coordinate catalytic zinc ion , 1995, Proteins.

[2]  Oliver F. Lange,et al.  Recognition Dynamics Up to Microseconds Revealed from an RDC-Derived Ubiquitin Ensemble in Solution , 2008, Science.

[3]  M. Gilson,et al.  Calculation of protein-ligand binding affinities. , 2007, Annual review of biophysics and biomolecular structure.

[4]  Sheng Jiang,et al.  Design, synthesis, and characterization of a potent, nonpeptide, cell-permeable, bivalent Smac mimetic that concurrently targets both the BIR2 and BIR3 domains in XIAP. , 2007, Journal of the American Chemical Society.

[5]  Shaomeng Wang,et al.  Structure-based design, synthesis and biochemical testing of novel and potent Smac peptido-mimetics. , 2005, Bioorganic & medicinal chemistry letters.

[6]  Brian K Shoichet,et al.  Prediction of protein-ligand interactions. Docking and scoring: successes and gaps. , 2006, Journal of medicinal chemistry.

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

[8]  Shaomeng Wang,et al.  Structure-based design of potent, conformationally constrained Smac mimetics. , 2004, Journal of the American Chemical Society.

[9]  D. Case,et al.  Generalized born models of macromolecular solvation effects. , 2000, Annual review of physical chemistry.

[10]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[11]  Geng Wu,et al.  Structural basis of IAP recognition by Smac/DIABLO , 2000, Nature.

[12]  W. Fairbrother,et al.  Design, synthesis, and biological activity of a potent Smac mimetic that sensitizes cancer cells to apoptosis by antagonizing IAPs. , 2006, ACS chemical biology.

[13]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

[14]  Stephen F. Betz,et al.  Structural basis for binding of Smac/DIABLO to the XIAP BIR3 domain , 2000, Nature.

[15]  Su Qiu,et al.  Structure-based design, synthesis, evaluation, and crystallographic studies of conformationally constrained Smac mimetics as inhibitors of the X-linked inhibitor of apoptosis protein (XIAP). , 2008, Journal of medicinal chemistry.

[16]  N. Munshi,et al.  Targeting mitochondrial factor Smac/DIABLO as therapy for multiple myeloma (MM). , 2004, Blood.

[17]  Shaomeng Wang,et al.  Design, synthesis, and evaluation of potent, nonpeptidic mimetics of second mitochondria-derived activator of caspases. , 2009, Journal of medicinal chemistry.

[18]  Shaomeng Wang,et al.  M-score: a knowledge-based potential scoring function accounting for protein atom mobility. , 2006, Journal of medicinal chemistry.

[19]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[20]  Mark A Olson,et al.  Calculation of absolute protein-ligand binding affinity using path and endpoint approaches. , 2006, Biophysical journal.

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

[22]  Irwin D Kuntz,et al.  Free energy calculations for theophylline binding to an RNA aptamer: Comparison of MM-PBSA and thermodynamic integration methods. , 2003, Biopolymers.

[23]  Luhua Lai,et al.  Further development and validation of empirical scoring functions for structure-based binding affinity prediction , 2002, J. Comput. Aided Mol. Des..

[24]  J. Snyder,et al.  The Conformations of Taxol in Chloroform , 2000 .

[25]  Shaomeng Wang,et al.  Design, synthesis, and evaluation of tricyclic, conformationally constrained small-molecule mimetics of second mitochondria-derived activator of caspases. , 2008, Journal of medicinal chemistry.

[26]  Vickie Tsui,et al.  Orally bioavailable antagonists of inhibitor of apoptosis proteins based on an azabicyclooctane scaffold. , 2009, Journal of medicinal chemistry.

[27]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[28]  C. E. Peishoff,et al.  A critical assessment of docking programs and scoring functions. , 2006, Journal of medicinal chemistry.

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

[30]  J. Bajorath,et al.  Docking and scoring in virtual screening for drug discovery: methods and applications , 2004, Nature Reviews Drug Discovery.

[31]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[32]  Chaohong Sun,et al.  Non-peptidic small molecule inhibitors of XIAP. , 2005, Bioorganic & medicinal chemistry letters.

[33]  Renxiao Wang,et al.  Comparative evaluation of 11 scoring functions for molecular docking. , 2003, Journal of medicinal chemistry.

[34]  Dajun Yang,et al.  Potent, orally bioavailable diazabicyclic small-molecule mimetics of second mitochondria-derived activator of caspases. , 2008, Journal of medicinal chemistry.