DeepBAR: A Fast and Exact Method for Binding Free Energy Computation.

The fast and accurate calculation of standard binding free energy has many important applications. Existing methodologies struggle at balancing accuracy and efficiency. We introduce a new method to compute binding free energy using deep generative models and the Bennett acceptance ratio method (DeepBAR). Compared to the rigorous potential of mean force (PMF) approach that requires sampling from intermediate states, DeepBAR is an order-of-magnitude more efficient as demonstrated in a series of host-guest systems. Notably, DeepBAR is exact and does not suffer from approximations for entropic contributions used in methods such as the molecular mechanics energy combined with the generalized Born and surface area continuum solvation (MM/GBSA). We anticipate DeepBAR to be a valuable tool for computing standard binding free energy used in drug design.

[1]  D. Case,et al.  Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. , 2003, Journal of molecular biology.

[2]  Jonah Z. Vilseck,et al.  Fast Solver for Large Scale Multistate Bennett Acceptance Ratio Equations. , 2019, Journal of chemical theory and computation.

[3]  Holger Gohlke,et al.  MMPBSA.py: An Efficient Program for End-State Free Energy Calculations. , 2012, Journal of chemical theory and computation.

[4]  Andrew T. Fenley,et al.  Toward Improved Force-Field Accuracy through Sensitivity Analysis of Host-Guest Binding Thermodynamics. , 2015, The journal of physical chemistry. B.

[5]  Michail Papadourakis,et al.  The SAMPL6 SAMPLing challenge: assessing the reliability and efficiency of binding free energy calculations , 2020, Journal of Computer-Aided Molecular Design.

[6]  David L Mobley,et al.  Predicting Binding Free Energies: Frontiers and Benchmarks. , 2017, Annual review of biophysics.

[7]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[8]  S. Genheden,et al.  The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities , 2015, Expert opinion on drug discovery.

[9]  Jimmy Ba,et al.  Adam: A Method for Stochastic Optimization , 2014, ICLR.

[10]  G. Whitesides,et al.  Self-Assembly of Mesoscale Objects into Ordered Two-Dimensional Arrays , 1997, Science.

[11]  Frank Noé,et al.  Equivariant Flows: sampling configurations for multi-body systems with symmetric energies , 2019, ArXiv.

[12]  P. Kollman,et al.  Combined molecular mechanical and continuum solvent approach (MM-PBSA/GBSA) to predict ligand binding , 2000 .

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

[14]  Hao Wu,et al.  Boltzmann generators: Sampling equilibrium states of many-body systems with deep learning , 2018, Science.

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

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

[17]  B. Roux,et al.  Calculation of absolute protein-ligand binding free energy from computer simulations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[18]  Shakir Mohamed,et al.  Variational Inference with Normalizing Flows , 2015, ICML.

[19]  Samy Bengio,et al.  Density estimation using Real NVP , 2016, ICLR.

[20]  Charles Blundell,et al.  Targeted free energy estimation via learned mappings. , 2020, The Journal of chemical physics.

[21]  Liping Cao,et al.  Cucurbit[7]uril⋅guest pair with an attomolar dissociation constant. , 2014, Angewandte Chemie.

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

[23]  Lyle Isaacs,et al.  The cucurbit[n]uril family: prime components for self-sorting systems. , 2005, Journal of the American Chemical Society.

[24]  Bin Zhang,et al.  Computing Absolute Free Energy with Deep Generative Models , 2020, The journal of physical chemistry. B.

[25]  Andrew T. Fenley,et al.  Computational Calorimetry: High-Precision Calculation of Host–Guest Binding Thermodynamics , 2015, Journal of chemical theory and computation.

[26]  Ajay,et al.  Computational methods to predict binding free energy in ligand-receptor complexes. , 1995, Journal of medicinal chemistry.

[27]  Michael K. Gilson,et al.  A synthetic host-guest system achieves avidin-biotin affinity by overcoming enthalpy–entropy compensation , 2007, Proceedings of the National Academy of Sciences.

[28]  Michael K Gilson,et al.  Host-guest complexes with protein-ligand-like affinities: computational analysis and design. , 2009, Journal of the American Chemical Society.

[29]  Kevin O'Connor,et al.  Permutation Invariant Likelihoods and Equivariant Transformations , 2019, ArXiv.