| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
* Institute of Molecular Pediatric Sciences, Gordon Center for Integrative Science, University of Chicago, Chicago, Illinois; and Bioscience Division, Mathematics and Computer Science Division, Argonne National Laboratory, Argonne, Illinois
Correspondence: Address reprint requests to Prof. Benoît Roux, Tel.: 773-834-3557; E-mail: roux{at}uchicago.edu.
The absolute (standard) binding free energy of eight FK506-related ligands to FKBP12 is calculated using free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent. A number of features are implemented to improve the accuracy and enhance the convergence of the calculations. First, the absolute binding free energy is decomposed into sequential steps during which the ligand-surrounding interactions as well as various biasing potentials restraining the translation, orientation, and conformation of the ligand are turned "on" and "off." Second, sampling of the ligand conformation is enforced by a restraining potential based on the root mean-square deviation relative to the bound state conformation. The effect of all the restraining potentials is rigorously unbiased, and it is shown explicitly that the final results are independent of all artificial restraints. Third, the repulsive and dispersive free energy contribution arising from the Lennard-Jones interactions of the ligand with its surrounding (protein and solvent) is calculated using the Weeks-Chandler-Andersen separation. This separation also improves convergence of the FEP/MD calculations. Fourth, to decrease the computational cost, only a small number of atoms in the vicinity of the binding site are simulated explicitly, while all the influence of the remaining atoms is incorporated implicitly using the generalized solvent boundary potential (GSBP) method. With GSBP, the size of the simulated FKBP12/ligand systems is significantly reduced, from 
25,000 to 2500. The computations are very efficient and the statistical error is small (1 kcal/mol). The calculated binding free energies are generally in good agreement with available experimental data and previous calculations (within 2 kcal/mol). The present results indicate that a strategy based on FEP/MD simulations of a reduced GSBP atomic model sampled with conformational, translational, and orientational restraining potentials can be computationally inexpensive and accurate.
This article has been cited by other articles:
 |
|
 |
|
  O. Khoruzhii, A. G. Donchev, N. Galkin, A. Illarionov, M. Olevanov, V. Ozrin, C. Queen, and V. Tarasov Application of a polarizable force field to calculations of relative protein-ligand binding affinities PNAS, July 29, 2008; 105(30): 10378 - 10383. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 L. Zhou and S. A. Siegelbaum Pathway and Endpoint Free Energy Calculations for Cyclic Nucleotide Binding to HCN Channels Biophys. J., June 15, 2008; 94(12): L90 - L92. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
D. Jiao, P. A. Golubkov, T. A. Darden, and P. Ren Calculation of protein-ligand binding free energy by using a polarizable potential PNAS, April 29, 2008; 105(17): 6290 - 6295. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
M. Ander, V. B. Luzhkov, and J. Aqvist Ligand Binding to the Voltage-Gated Kv1.5 Potassium Channel in the Open State Docking and Computer Simulations of a Homology Model Biophys. J., February 1, 2008; 94(3): 820 - 831. [Abstract] [Full Text] [PDF]
|
|
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
