SAMPL6 host–guest blind predictions using a non equilibrium alchemical approach

In this paper, we compute, by means of a non equilibrium alchemical technique, called fast switching double annihilation methods (FSDAM), the absolute standard dissociation free energies of the the octa acids host–guest systems in the SAMPL6 challenge initiative. FSDAM is based on the production of canonical configurations of the bound and unbound states via enhanced sampling and on the subsequent generation of hundreds of fast non-equilibrium ligand annihilation trajectories. The annihilation free energies of the ligand when bound to the receptor and in bulk solvent are obtained from the collection of work values using an estimate based on the Crooks theorem for driven non equilibrium processes. The FSDAM blind prediction, relying on the normality assumption for the annihilation work distributions, ranked fairly well among the submitted blind predictions that were not adjusted with a linear corrections obtained from retrospective data on similar host guest systems. Improved results for FSDAM can be obtained by post-processing the work data assuming mixtures of normal components.

[1]  G. Crooks Nonequilibrium Measurements of Free Energy Differences for Microscopically Reversible Markovian Systems , 1998 .

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

[3]  B. Gibb,et al.  Molecular containers assembled through the hydrophobic effect. , 2015, Chemical Society reviews.

[4]  J. Kirkwood Statistical Mechanics of Fluid Mixtures , 1935 .

[5]  D. Darling,et al.  A Test of Goodness of Fit , 1954 .

[6]  Piero Procacci,et al.  Hybrid MPI/OpenMP Implementation of the ORAC Molecular Dynamics Program for Generalized Ensemble and Fast Switching Alchemical Simulations , 2016, J. Chem. Inf. Model..

[7]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[8]  Alessio Rocchi,et al.  The role of medium size facilities in the HPC ecosystem: the case of the new CRESCO4 cluster integrated in the ENEAGRID infrastructure , 2014, 2014 International Conference on High Performance Computing & Simulation (HPCS).

[9]  Tanweer Ul Islam,et al.  Stringency-based ranking of normality tests , 2017, Commun. Stat. Simul. Comput..

[10]  B. Gibb,et al.  Guest-mediated switching of the assembly state of a water-soluble deep-cavity cavitand. , 2013, Chemical communications.

[11]  Michael K. Gilson,et al.  Theory of free energy and entropy in noncovalent binding. , 2009, Chemical reviews.

[12]  D. Lauffenburger,et al.  Receptors: Models for Binding, Trafficking, and Signaling , 1993 .

[13]  Lee G. Pedersen,et al.  Ionic charging free energies: Spherical versus periodic boundary conditions , 1998 .

[14]  Piero Procacci,et al.  I. Dissociation free energies of drug-receptor systems via non-equilibrium alchemical simulations: a theoretical framework. , 2016, Physical chemistry chemical physics : PCCP.

[15]  P. Procacci Unbiased free energy estimates in fast nonequilibrium transformations using Gaussian mixtures. , 2015, The Journal of chemical physics.

[16]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[17]  K. Sharp,et al.  On the calculation of absolute macromolecular binding free energies , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[18]  Y. B. Wah,et al.  Power comparisons of Shapiro-Wilk , Kolmogorov-Smirnov , Lilliefors and Anderson-Darling tests , 2011 .

[19]  K. Schulten,et al.  Calculating potentials of mean force from steered molecular dynamics simulations. , 2004, The Journal of chemical physics.

[20]  Michael R Shirts,et al.  Linear basis function approach to efficient alchemical free energy calculations. 2. Inserting and deleting particles with coulombic interactions. , 2015, Journal of chemical theory and computation.

[21]  Michael P Eastwood,et al.  Minimizing thermodynamic length to select intermediate states for free-energy calculations and replica-exchange simulations. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[23]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[24]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[25]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[26]  P. Procacci,et al.  Binding Free Energies of Host-Guest Systems by Nonequilibrium Alchemical Simulations with Constrained Dynamics: Theoretical Framework. , 2017, Journal of chemical theory and computation.

[27]  P. Procacci,et al.  Energy-Driven Undocking (EDU-HREM) in Solute Tempering Replica Exchange Simulations. , 2014, Journal of chemical theory and computation.

[28]  R GuptaMaya,et al.  Theory and Use of the EM Algorithm , 2011 .

[29]  Stefan Boresch,et al.  Absolute Binding Free Energies: A Quantitative Approach for Their Calculation , 2003 .

[30]  S. Menichetti,et al.  Efficient nonequilibrium method for binding free energy calculations in molecular dynamics simulations. , 2015, Journal of chemical theory and computation.

[31]  Piero Procacci,et al.  Fast Switching Alchemical Transformations in Molecular Dynamics Simulations. , 2014, Journal of chemical theory and computation.

[32]  Wei Yao,et al.  Overview of the SAMPL6 host-guest binding affinity prediction challenge , 2018, bioRxiv.

[33]  M. Stephens Tests of fit for the logistic distribution based on the empirical distribution function , 1979 .

[34]  M. Parrinello,et al.  Crystal structure and pair potentials: A molecular-dynamics study , 1980 .

[35]  D. Chandler,et al.  Introduction To Modern Statistical Mechanics , 1987 .

[36]  Corinne L. D. Gibb,et al.  Binding of cyclic carboxylates to octa-acid deep-cavity cavitand , 2014, Journal of Computer-Aided Molecular Design.

[37]  P. Procacci,et al.  Binding Free Energies of Host-Guest Systems by Nonequilibrium Alchemical Simulations with Constrained Dynamics: Illustrative Calculations and Numerical Validation. , 2017, Journal of chemical theory and computation.

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

[39]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[40]  Riccardo Chelli,et al.  II. Dissociation free energies in drug-receptor systems via nonequilibrium alchemical simulations: application to the FK506-related immunophilin ligands. , 2016, Physical chemistry chemical physics : PCCP.

[41]  Simone Marsili,et al.  Crooks equation for steered molecular dynamics using a Nosé-Hoover thermostat. , 2006, The Journal of chemical physics.

[42]  Piero Procacci,et al.  Coordinates scaling and multiple time step algorithms for simulation of solvated proteins in the NPT ensemble , 1998 .

[43]  Piero Procacci,et al.  ORAC: A Molecular dynamics program to simulate complex molecular systems with realistic electrostatic interactions , 1997 .

[44]  Jian Yin,et al.  Overview of the SAMPL5 host–guest challenge: Are we doing better? , 2016, Journal of Computer-Aided Molecular Design.

[45]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[46]  R. Chelli,et al.  Combining path-breaking with bidirectional nonequilibrium simulations to improve efficiency in free energy calculations. , 2014, The Journal of chemical physics.

[47]  Lu Wang,et al.  Inclusion of Loss of Translational and Rotational Freedom in Theoretical Estimates of Free Energies of Binding. Application to a Complex of Benzene and Mutant T4 Lysozyme , 1997 .

[48]  S. E. Ahmed,et al.  Handbook of Statistical Distributions with Applications , 2007, Technometrics.

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

[50]  Eamonn F. Healy,et al.  Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model , 1985 .

[51]  P. Procacci PrimaDORAC: A Free Web Interface for the Assignment of Partial Charges, Chemical Topology, and Bonded Parameters in Organic or Drug Molecules , 2017, J. Chem. Inf. Model..

[52]  David L. Mobley,et al.  The SAMPL4 host–guest blind prediction challenge: an overview , 2014, Journal of Computer-Aided Molecular Design.

[53]  R. Zwanzig High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[54]  P. Procacci,et al.  Statistical Mechanics of Ligand-Receptor Noncovalent Association, Revisited: Binding Site and Standard State Volumes in Modern Alchemical Theories. , 2017, Journal of chemical theory and computation.

[55]  Jennifer L. Knight,et al.  Accurate and reliable prediction of relative ligand binding potency in prospective drug discovery by way of a modern free-energy calculation protocol and force field. , 2015, Journal of the American Chemical Society.

[56]  Simone Marsili,et al.  ORAC: A molecular dynamics simulation program to explore free energy surfaces in biomolecular systems at the atomistic level , 2009, J. Comput. Chem..

[57]  Christopher I. Bayly,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation , 2002, J. Comput. Chem..

[58]  A. Mark,et al.  Avoiding singularities and numerical instabilities in free energy calculations based on molecular simulations , 1994 .

[59]  Maya R. Gupta,et al.  Theory and Use of the EM Algorithm , 2011, Found. Trends Signal Process..