Accurate Calculation of Relative Binding Free Energies between Ligands with Different Net Charges.

In drug discovery programs, modifications that change the net charge of the ligands are often considered to improve the binding potency and solubility, or to address other ADME/Tox problems. Accurate calculation of the binding free-energy changes associated with charge-changing perturbations remains a great challenge of central importance in computational drug discovery. The finite size effects associated with periodic boundary condition and lattice summation employed in common molecular dynamics simulations introduce artifacts in the electrostatic potential energy calculations, which need to be carefully handled for accurate free-energy calculations between systems with different net charges. The salts in the buffer solution of experimental binding affinity assays also have a strong effect on the binding free energies between charged species, which further complicates the modeling of the charge-changing perturbations. Here, we extend our free-energy perturbation (FEP) algorithm, which has been extensively applied to many drug discovery programs for relative binding free-energy calculations between ligands with the same net charge (charge-conserving perturbation), to enable charge-changing perturbations. We have investigated three different approaches to correct the finite size effects and tested them on 10 protein targets and 31 charge-changing perturbations. We have found that all three methods are able to successfully eliminate the box-size dependence of calculated binding free energies associated with brute force FEP. Moreover, inclusion of salts matching the ionic strength of experimental buffer solution significantly improves the calculated binding free energies. For ligands with multiple possible protonation states, we applied the p Ka correction to account for the ionization equilibrium of the ligands and the results are significantly improved. Finally, the calculated binding free energies from these methods agree with each other, and also agree well with the experimental results. The root-mean-square error between the calculated binding free energies and experimental data is 1.1 kcal/mol, which is on par with the accuracy of charge-conserving perturbations. We anticipate that the outstanding accuracy demonstrated here across a broad range of target classes may have significant implications for drug discovery projects, where charge-changing modifications must be considered.

[1]  Robert Abel,et al.  Advancing Drug Discovery through Enhanced Free Energy Calculations. , 2017, Accounts of chemical research.

[2]  Jacob D. Durrant,et al.  Molecular dynamics simulations and drug discovery , 2011, BMC Biology.

[3]  B. Brooks,et al.  Hydronium Ions Accompanying Buried Acidic Residues Lead to High Apparent Dielectric Constants in the Interior of Proteins. , 2018, The journal of physical chemistry. B.

[4]  David L Mobley,et al.  Calculating the binding free energies of charged species based on explicit-solvent simulations employing lattice-sum methods: an accurate correction scheme for electrostatic finite-size effects. , 2013, The Journal of chemical physics.

[5]  Yuzhu Chen,et al.  N2-substituted O6-cyclohexylmethylguanine derivatives: potent inhibitors of cyclin-dependent kinases 1 and 2. , 2004, Journal of medicinal chemistry.

[6]  N. Farrow,et al.  Optimization of piperidyl-ureas as inhibitors of soluble epoxide hydrolase. , 2010, Bioorganic & medicinal chemistry letters.

[7]  William L Jorgensen,et al.  Efficient drug lead discovery and optimization. , 2009, Accounts of chemical research.

[8]  Adam R. Johnson,et al.  Lead identification of novel and selective TYK2 inhibitors. , 2013, European journal of medicinal chemistry.

[9]  Christophe Chipot,et al.  Can Absolute Free Energies of Association Be Estimated from Molecular Mechanical Simulations? The Biotin−Streptavidin System Revisited , 2001 .

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

[11]  Patrick G. Blachly,et al.  Protocols Utilizing Constant pH Molecular Dynamics to Compute pH-Dependent Binding Free Energies , 2014, The journal of physical chemistry. B.

[12]  M. Edwards,et al.  Lipophilic Efficiency as an Important Metric in Drug Design. , 2018, Journal of medicinal chemistry.

[13]  Proton-Coupled Conformational Allostery Modulates the Inhibitor Selectivity for β-Secretase. , 2017, The journal of physical chemistry letters.

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

[15]  Lingle Wang,et al.  On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities , 2012, Proceedings of the National Academy of Sciences.

[16]  John P. Overington,et al.  Probing the links between in vitro potency, ADMET and physicochemical parameters , 2011, Nature Reviews Drug Discovery.

[17]  Wei Chen,et al.  Introducing titratable water to all-atom molecular dynamics at constant pH. , 2013, Biophysical journal.

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

[19]  D. Tulshian,et al.  Discovery and Structure Enabled Synthesis of 2,6-Diaminopyrimidin-4-one IRAK4 Inhibitors. , 2015, ACS medicinal chemistry letters.

[20]  Scott A. Erickson,et al.  Imidazo[2,1-b]thiazoles: multitargeted inhibitors of both the insulin-like growth factor receptor and members of the epidermal growth factor family of receptor tyrosine kinases. , 2010, Bioorganic & medicinal chemistry letters.

[21]  P. McEwan,et al.  Tetrahydroindazoles as Interleukin-2 Inducible T-Cell Kinase Inhibitors. Part II. Second-Generation Analogues with Enhanced Potency, Selectivity, and Pharmacodynamic Modulation in Vivo. , 2015, Journal of medicinal chemistry.

[22]  Jana K. Shen,et al.  Charge-leveling and proper treatment of long-range electrostatics in all-atom molecular dynamics at constant pH. , 2012, The Journal of chemical physics.

[23]  J. Schmit,et al.  SLTCAP: A Simple Method for Calculating the Number of Ions Needed for MD Simulation. , 2018, Journal of chemical theory and computation.

[24]  Thomas Simonson,et al.  An Overview of Electrostatic Free Energy Computations for Solutions and Proteins. , 2014, Journal of chemical theory and computation.

[25]  M. Hann.  Molecular obesity, potency and other addictions in drug discovery , 2011 .

[26]  William L. Jorgensen,et al.  A Quantum Mechanical and Molecular Mechanical Method Based on CM1A Charges: Applications to Solvent Effects on Organic Equilibria and Reactions , 1998 .

[27]  Gregory A Ross,et al.  Biomolecular Simulations under Realistic Macroscopic Salt Conditions , 2017, bioRxiv.

[28]  K. Scearce-Levie,et al.  Scaffold-Hopping and Structure-Based Discovery of Potent, Selective, And Brain Penetrant N-(1H-Pyrazol-3-yl)pyridin-2-amine Inhibitors of Dual Leucine Zipper Kinase (DLK, MAP3K12). , 2015, Journal of medicinal chemistry.

[29]  Robert Abel,et al.  Accurate Modeling of Scaffold Hopping Transformations in Drug Discovery. , 2017, Journal of chemical theory and computation.

[30]  B. Roux,et al.  Constant-pH Hybrid Nonequilibrium Molecular Dynamics–Monte Carlo Simulation Method , 2015, Journal of chemical theory and computation.

[31]  Jennifer L. Knight,et al.  OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins. , 2016, Journal of chemical theory and computation.

[32]  J. Mccammon,et al.  Ewald artifacts in computer simulations of ionic solvation and ion–ion interaction: A continuum electrostatics study , 1999 .

[33]  A. Moretto,et al.  Structure-based optimization of protein tyrosine phosphatase 1B inhibitors: from the active site to the second phosphotyrosine binding site. , 2007, Journal of medicinal chemistry.

[34]  B. García-Moreno E.,et al.  The pK(a) values of acidic and basic residues buried at the same internal location in a protein are governed by different factors. , 2009, Journal of molecular biology.

[35]  B. Berne,et al.  Replica exchange with solute scaling: a more efficient version of replica exchange with solute tempering (REST2). , 2011, The journal of physical chemistry. B.

[36]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[37]  Chaohong Sun,et al.  Aminopyridine-based c-Jun N-terminal kinase inhibitors with cellular activity and minimal cross-kinase activity. , 2006, Journal of medicinal chemistry.

[38]  Karl T. Debiec,et al.  Evaluating the Strength of Salt Bridges: A Comparison of Current Biomolecular Force Fields , 2014, The journal of physical chemistry. B.

[39]  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.

[40]  James Andrew McCammon,et al.  Conformational Dynamics and Binding Free Energies of Inhibitors of BACE-1: From the Perspective of Protonation Equilibria , 2015, PLoS Comput. Biol..

[41]  B. García-Moreno E.,et al.  Stabilization of internal charges in a protein: water penetration or conformational change? , 2004, Biophysical journal.

[42]  D. van der Spoel,et al.  A temperature predictor for parallel tempering simulations. , 2008, Physical chemistry chemical physics : PCCP.

[43]  Jing Zhang,et al.  Jaguar: A high-performance quantum chemistry software program with strengths in life and materials sciences , 2013 .

[44]  Robert Abel,et al.  A Critical Review of Validation, Blind Testing, and Real- World Use of Alchemical Protein-Ligand Binding Free Energy Calculations. , 2017, Current topics in medicinal chemistry.

[45]  G. Kontopidis,et al.  Quantification of the effects of ionic strength, viscosity, and hydrophobicity on protein-ligand binding affinity. , 2014, ACS medicinal chemistry letters.

[46]  Adam R. Johnson,et al.  Lead optimization of a 4-aminopyridine benzamide scaffold to identify potent, selective, and orally bioavailable TYK2 inhibitors. , 2013, Journal of medicinal chemistry.

[47]  Robert Abel,et al.  Accurate and Reliable Prediction of the Binding Affinities of Macrocycles to Their Protein Targets. , 2017, Journal of chemical theory and computation.

[48]  Adam R. Johnson,et al.  Novel triazolo-pyrrolopyridines as inhibitors of Janus kinase 1. , 2013, Bioorganic & medicinal chemistry letters.

[49]  Nathan A. Baker,et al.  Electrostatics of nanosystems: Application to microtubules and the ribosome , 2001, Proceedings of the National Academy of Sciences of the United States of America.

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

[51]  Jane R Kenny,et al.  Identification of C-2 hydroxyethyl imidazopyrrolopyridines as potent JAK1 inhibitors with favorable physicochemical properties and high selectivity over JAK2. , 2013, Journal of medicinal chemistry.

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

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

[54]  Jennifer L. Knight,et al.  Modeling Local Structural Rearrangements Using FEP/REST: Application to Relative Binding Affinity Predictions of CDK2 Inhibitors. , 2013, Journal of chemical theory and computation.

[55]  Nathan A. Baker,et al.  APBSmem: A Graphical Interface for Electrostatic Calculations at the Membrane , 2010, PloS one.

[56]  B. García-Moreno E.,et al.  Large shifts in pKa values of lysine residues buried inside a protein , 2011, Proceedings of the National Academy of Sciences.