Calculation of Absolute Protein-Ligand Binding Affinity Using Path and Endpoint Approaches

Abstract

A comparative analysis is provided of rigorous and approximate methods for calculating absolute binding affinities of two protein-ligand complexes: the FKBP protein bound with small molecules 4-hydroxy-2-butanone and FK506. Our rigorous approach is an umbrella sampling technique where a potential of mean force is determined by pulling the ligand out of the protein active site over several simulation windows. The results of this approach agree well with experimentally observed binding affinities. Also assessed is a commonly used approximate endpoint approach, which separately estimates enthalpy, solvation free energy, and entropy. We show that this endpoint approach has numerous variations, all of which are prone to critical shortcomings. For example, conventional harmonic and quasiharmonic entropy estimation procedures produce disparate results for the relatively simple protein-ligand systems studied in this work.

Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration

“…(1) Small molecule conformational sampling [73][74][75][76]; (2) Protein-folding [77][78][79] and large-scale protein conformational sampling [80][81][82][83]; (3) Protein-protein/peptide-peptide interactions [84][85][86][87][88][89][90][91][92]; (4) DNA conformational changes [93] and DNA-DNA interactions [94][95][96]; (5) Binding and association free-energies [97][98][99][100][101][102][103][104][105][106][107]; (6) Adsorption on and permeation through lipid bilayers [108][109][110][111][112][113][114][115][116][117] …”

Electrostatics in Proteins and Protein–Ligand Complexes

“…Several variations of MM/GBSA models have been tested in recent years and compared with experimental data and microscopic models [75][76][77][78][79][80]. For the majority of applications, MM/GBSA models showed the strong correlation between experimental and computational data.…”

Free energy surfaces for the interaction of D-glucose with planar aromatic groups in aqueous solution

“…In an isotropic, infinitely dilute system, A is assumed to rotate freely in the reference frame of B, and the integration over the Euler angles can be carried out directly. 17 For the case of spherical symmetry, the simplest possible case, the work function is a function of only one variable, the distance r. The free energy then becomes…”