Prediction of SAMPL3 host-guest affinities with the binding energy distribution analysis method (BEDAM)

BEDAM calculations are described to predict the free energies of binding of a series of anaesthetic drugs to a recently characterized acyclic cucurbituril host. The modeling predictions, conducted as part of the SAMPL3 host-guest affinity blind challenge, are generally in good quantitative agreement with the experimental measurements. The correlation coefficient between computed and measured binding free energies is 70% with high statistical significance. Multiple conformational stereoisomers and protonation states of the guests have been considered. Better agreement is obtained with high statistical confidence under acidic modeling conditions. It is shown that this level of quantitative agreement could have not been reached without taking into account reorganization energy and configurational entropy effects. Extensive conformational variability of the host, the guests and their complexes is observed in the simulations, affecting binding free energy estimates and structural predictions. A conformational reservoir technique is introduced as part of the parallel Hamiltonian replica exchange molecular dynamics BEDAM protocol to fully capture conformational variability. It is shown that these advanced computational strategies lead to converged free energy estimates for these systems, offering the prospect of utilizing host-guest binding free energy data for force field validation and development.

[1]  B. Roux,et al.  Computations of standard binding free energies with molecular dynamics simulations. , 2009, The journal of physical chemistry. B.

[2]  Ronald M. Levy,et al.  AGBNP: An analytic implicit solvent model suitable for molecular dynamics simulations and high‐resolution modeling , 2004, J. Comput. Chem..

[3]  Emilio Gallicchio,et al.  The Binding Energy Distribution Analysis Method (BEDAM) for the Estimation of Protein-Ligand Binding Affinities. , 2010, Journal of chemical theory and computation.

[4]  Shoshana J Wodak.  From the Mediterranean coast to the shores of Lake Ontario: CAPRI's premiere on the American continent , 2007, Proteins.

[5]  J. Guthrie,et al.  A blind challenge for computational solvation free energies: introduction and overview. , 2009, The journal of physical chemistry. B.

[6]  A. Tramontano,et al.  Critical assessment of methods of protein structure prediction (CASP)—round IX , 2011, Proteins.

[7]  Ian W. Wyman,et al.  Host-guest complexations of local anaesthetics by cucurbit[7]uril in aqueous solution. , 2010, Organic & biomolecular chemistry.

[8]  Emilio Gallicchio,et al.  The AGBNP2 Implicit Solvation Model. , 2009, Journal of chemical theory and computation.

[9]  David L Mobley,et al.  Alchemical free energy methods for drug discovery: progress and challenges. , 2011, Current opinion in structural biology.

[10]  M. Gilson,et al.  Ligand configurational entropy and protein binding , 2007, Proceedings of the National Academy of Sciences.

[11]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

[12]  Ronald M. Levy,et al.  Conformational populations of ligand‐sized molecules by replica exchange molecular dynamics and temperature reweighting , 2009, J. Comput. Chem..

[13]  Richard A Friesner,et al.  Serial replica exchange. , 2007, The journal of physical chemistry. B.

[14]  Emilio Gallicchio,et al.  Conformational Transitions and Convergence of Absolute Binding Free Energy Calculations. , 2012, Journal of chemical theory and computation.

[15]  David L. Mobley,et al.  Drug Design: Free-energy calculations in structure-based drug design , 2010 .

[16]  Lyle Isaacs,et al.  Acyclic cucurbit[n]uril congeners are high affinity hosts. , 2010, The Journal of organic chemistry.

[17]  Manish Parashar,et al.  Asynchronous replica exchange for molecular simulations , 2008, J. Comput. Chem..

[18]  David L. Mobley,et al.  Chapter 4 Alchemical Free Energy Calculations: Ready for Prime Time? , 2007 .

[19]  Michael R. Shirts,et al.  Statistically optimal analysis of samples from multiple equilibrium states. , 2008, The Journal of chemical physics.

[20]  Anthony K. Felts,et al.  Temperature weighted histogram analysis method, replica exchange, and transition paths. , 2005, The journal of physical chemistry. B.

[21]  A. Pohorille,et al.  Free energy calculations : theory and applications in chemistry and biology , 2007 .

[22]  K. Dill,et al.  Binding of small-molecule ligands to proteins: "what you see" is not always "what you get". , 2009, Structure.

[23]  A. Roitberg,et al.  Coupling of replica exchange simulations to a non-Boltzmann structure reservoir. , 2007, The journal of physical chemistry. B.

[24]  Christophe Chipot,et al.  Good practices in free-energy calculations. , 2010, The journal of physical chemistry. B.

[25]  Michael K Gilson,et al.  New ultrahigh affinity host-guest complexes of cucurbit[7]uril with bicyclo[2.2.2]octane and adamantane guests: thermodynamic analysis and evaluation of M2 affinity calculations. , 2011, Journal of the American Chemical Society.

[26]  Michael K. Gilson,et al.  Structure and Thermodynamics of Molecular Hydration via Grid Inhomogeneous Solvation Theory , 2011, 1108.4876.

[27]  David L Mobley,et al.  Predicting ligand binding affinity with alchemical free energy methods in a polar model binding site. , 2009, Journal of molecular biology.

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

[29]  Jonathan C. Horton,et al.  What you see ... , 2001, Nature.

[30]  Emilio Gallicchio,et al.  Advances in all atom sampling methods for modeling protein-ligand binding affinities. , 2011, Current opinion in structural biology.

[31]  Richard A. Friesner,et al.  Integrated Modeling Program, Applied Chemical Theory (IMPACT) , 2005, J. Comput. Chem..

[32]  B. Rost,et al.  Critical assessment of methods of protein structure prediction—Round VIII , 2009, Proteins.

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

[34]  J. Nielsen,et al.  The pKa Cooperative: A collaborative effort to advance structure‐based calculations of pKa values and electrostatic effects in proteins , 2011, Proteins.

[35]  Emilio Gallicchio,et al.  Recent theoretical and computational advances for modeling protein-ligand binding affinities. , 2011, Advances in protein chemistry and structural biology.

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

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

[38]  Michael K. Gilson,et al.  Blind prediction of host–guest binding affinities: a new SAMPL3 challenge , 2012, Journal of Computer-Aided Molecular Design.

[39]  Christophe Chipot,et al.  Comprar Free Energy Calculations · Theory and Applications in Chemistry and Biology | Chipot, Christophe | 9783540736172 | Springer , 2007 .

[40]  Christophe Chipot,et al.  Free Energy Calculations , 2008 .

[41]  W. L. Jorgensen.  The Many Roles of Computation in Drug Discovery , 2004, Science.

[42]  R. Levy,et al.  Simple continuous and discrete models for simulating replica exchange simulations of protein folding. , 2008, The journal of physical chemistry. B.