Abstract: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 pK(a) 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.

PEARLS: Program for Energetic Analysis of Receptor-Ligand System

CLiBE: a Database of Computed Ligand Binding Energy for Ligand-Receptor Complexes.

LARMD: Integration of Bioinformatic Resources to Profile Ligand-Driven Protein Dynamics with a Case on the Activation of Estrogen Receptor.

Amber Free Energy Tools: Interoperable Software for Free Energy Simulations Using Generalized Quantum Mechanical/molecular Mechanical and Machine Learning Potentials.

Pbsa_e: A Pbsa-Based Free Energy Estimator for Protein-Ligand Binding Affinity

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

Understanding Aldose Reductase-Inhibitors interactions with free energy simulation

CHARMM-GUI Free Energy Calculator for Absolute and Relative Ligand Solvation and Binding Free Energy Simulations

CHARMM-GUI Ligand Binder for Absolute Binding Free Energy Calculations and Its Application.

Implementing Protein-Polarized Ligand Charges into Relative Protein Ligand Binding Affinities Calculations

Effect of Protein-Polarized Ligand Charges on Relative Protein Ligand Binding Affinities

BFEE: A User-Friendly Graphical Interface Facilitating Absolute Binding Free-Energy Calculations.

Hotlig: A Molecular Surface-Directed Approach to Scoring Protein-Ligand Interactions

Alchemical Binding Free Energy Calculations in AMBER20: Advances and Best Practices for Drug Discovery

Insights into the Molecular Mechanisms of Protein-Ligand Interactions by Molecular Docking and Molecular Dynamics Simulation: A Case of Oligopeptide Binding Protein

Ligand-protein database: linking protein-ligand complex structures to binding data.

ProteinLens: a web-based application for the analysis of allosteric signalling on atomistic graphs of biomolecules

Automated relative binding free energy calculations from SMILES to ΔΔG

Accurate calculation of absolute protein-ligand binding free energies

farPPI: a webserver for accurate prediction of protein-ligand binding structures for small-molecule PPI inhibitors by MM/PB(GB)SA methods

Estimating protein-ligand binding free energy: atomic solvation parameters for partition coefficient and solvation free energy calculation.