Independent-Trajectory Thermodynamic Integration: a practical guide to protein-drug binding free energy calculations using distributed computing.

The Independent-Trajectory Thermodynamic Integration (IT-TI) approach for free energy calculation with distributed computing is described. IT-TI utilizes diverse conformational sampling obtained from multiple, independent simulations to obtain more reliable free energy estimates compared to single TI predictions. The latter may significantly under- or over-estimate the binding free energy due to finite sampling. We exemplify the advantages of the IT-TI approach using two distinct cases of protein-ligand binding. In both cases, IT-TI yields distributions of absolute binding free energy estimates that are remarkably centered on the target experimental values. Alternative protocols for the practical and general application of IT-TI calculations are investigated. We highlight a protocol that maximizes predictive power and computational efficiency.

[1]  B. Roux,et al.  Absolute binding free energy calculations using molecular dynamics simulations with restraining potentials. , 2006, Biophysical journal.

[2]  K. Dill,et al.  On the use of orientational restraints and symmetry corrections in alchemical free energy calculations. , 2006, The Journal of chemical physics.

[3]  Lingle Wang.  Towards more robust and efficient methods for the calculation of Protein-Ligand binding affinities , 2012 .

[4]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[5]  Rommie E. Amaro,et al.  Impact of calcium on N1 influenza neuraminidase dynamics and binding free energy , 2010, Proteins.

[6]  James Andrew McCammon,et al.  Ligand-receptor interactions , 1984, Comput. Chem..

[7]  Ignacio J General.  A Note on the Standard State's Binding Free Energy. , 2010, Journal of chemical theory and computation.

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

[9]  Ken A Dill,et al.  Use of the Weighted Histogram Analysis Method for the Analysis of Simulated and Parallel Tempering Simulations. , 2007, Journal of chemical theory and computation.

[10]  F. Hayden,et al.  In Vitro Characterization of A-315675, a Highly Potent Inhibitor of A and B Strain Influenza Virus Neuraminidases and Influenza Virus Replication , 2002, Antimicrobial Agents and Chemotherapy.

[11]  Y. Cheng,et al.  Relationship between the inhibition constant (K1) and the concentration of inhibitor which causes 50 per cent inhibition (I50) of an enzymatic reaction. , 1973, Biochemical pharmacology.

[12]  J. Åqvist,et al.  Ion-water interaction potentials derived from free energy perturbation simulations , 1990 .

[13]  David J. Stevens,et al.  The structure of H5N1 avian influenza neuraminidase suggests new opportunities for drug design , 2006, Nature.

[14]  T. Hansson,et al.  On the Validity of Electrostatic Linear Response in Polar Solvents , 1996 .

[15]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

[16]  Wilfred F. van Gunsteren,et al.  Computation of free energy , 2002 .

[17]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[18]  Michael R. Shirts,et al.  Direct calculation of the binding free energies of FKBP ligands. , 2005, The Journal of chemical physics.

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

[20]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[21]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

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

[23]  S. Diamond,et al.  Identification of triazinoindol-benzimidazolones as nanomolar inhibitors of the Mycobacterium tuberculosis enzyme TDP-6-deoxy-d-xylo-4-hexopyranosid-4-ulose 3,5-epimerase (RmlC). , 2010, Bioorganic & medicinal chemistry.

[24]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[25]  Donald Hamelberg,et al.  Standard free energy of releasing a localized water molecule from the binding pockets of proteins: double-decoupling method. , 2004, Journal of the American Chemical Society.

[26]  Wei Yang,et al.  Free energy simulations: use of reverse cumulative averaging to determine the equilibrated region and the time required for convergence. , 2004, The Journal of chemical physics.

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

[28]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[29]  H. C. Andersen.  Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

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

[31]  E. Carlstein.  The Use of Subseries Values for Estimating the Variance of a General Statistic from a Stationary Sequence , 1986 .

[32]  Wilfred F. van Gunsteren,et al.  Basic ingredients of free energy calculations: A review , 2009, J. Comput. Chem..

[33]  Ioannis G. Economou,et al.  Effect of the integration method on the accuracy and computational efficiency of free energy calculations using thermodynamic integration , 2010 .

[34]  Yi Wang,et al.  Effects of Biomolecular Flexibility on Alchemical Calculations of Absolute Binding Free Energies. , 2011, Journal of chemical theory and computation.

[35]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[36]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[37]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[38]  Wilfred F van Gunsteren,et al.  Computational Analysis of the Mechanism and Thermodynamics of Inhibition of Phosphodiesterase 5A by Synthetic Ligands. , 2007, Journal of chemical theory and computation.

[39]  T. Straatsma,et al.  Separation‐shifted scaling, a new scaling method for Lennard‐Jones interactions in thermodynamic integration , 1994 .

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

[41]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[42]  García,et al.  Large-amplitude nonlinear motions in proteins. , 1992, Physical review letters.

[43]  Jung-Hsin Lin,et al.  Remarkable loop flexibility in avian influenza N1 and its implications for antiviral drug design. , 2007, Journal of the American Chemical Society.

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

[45]  J. Andrew McCammon,et al.  Independent-Trajectories Thermodynamic-Integration Free-Energy Changes for Biomolecular Systems: Determinants of H5N1 Avian Influenza Virus Neuraminidase Inhibition by Peramivir , 2009, Journal of chemical theory and computation.

[46]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .