Abstract:Accurate prediction of binding free energies is critical to streamlining the drug development and protein design process. With the advent of GPU acceleration, absolute alchemical methods, which simulate the removal of ligand electrostatics and van der Waals interactions with the protein, have become routinely accessible and provide a physically rigorous approach that enables full consideration of flexibility and solvent interaction. However, standard explicit solvent simulations are unable to model protonation or electronic polarization changes upon ligand transfer from water to the protein interior, leading to inaccurate prediction of binding affinities for charged molecules. Here, we perform extensive simulation totaling ∼540 μs to benchmark the impact of modeling conditions on predictive accuracy for absolute alchemical simulations. Binding to urokinase plasminogen activator (UPA), a protein frequently overexpressed in metastatic tumors, is evaluated for a set of 10 inhibitors with extended flexibility, highly charged character, and titratable properties. We demonstrate that the alchemical simulations can be adapted to utilize the MBAR/PBSA method to improve the accuracy upon incorporating electronic polarization, highlighting the importance of polarization in alchemical simulations of binding affinities. Comparison of binding energy prediction at various protonation states indicates that proper electrostatic setup is also crucial in binding affinity prediction of charged systems, prompting us to propose an alternative binding mode with protonated ligand phenol and Hid-46 at the binding site, a testable hypothesis for future experimental validation.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jctc.0c01305?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jctc.0c01305</a>.Discussion of apparent binding free energies with one titratable group, MBAR/PBSA accuracy, binding free energy predictions for all conditions, illustration of restraint definitions, decharging and VDW free energies from alchemical simulations, illustration of the UPA active site with ligand, B-Factor, GAFF and GAFF2 structure analysis, and correlation plots of binding affinity for ligand restraints, differing protonation states, with and without outlier and with MBAR/PBSA postprocessing (<a class="ext-link" href="/doi/suppl/10.1021/acs.jctc.0c01305/suppl_file/ct0c01305_si_001.pdf">PDF</a>)This article has not yet been cited by other publications.

Accurate Prediction of Relative Binding Affinities of a Series of HIV-1 Protease Inhibitors Using Semi-Empirical Quantum Mechanical Charge

Prediction of Potency of Protease Inhibitors Using Free Energy Simulations with Polarizable Quantum Mechanics-Based Ligand Charges and a Hybrid Water Model

Improving the accuracy of predicting protein-ligand binding-free energy with semiempirical quantum chemistry charge.

Combining Molecular Dynamic Information and an Aspherical-Atom Data Bank in the Evaluation of the Electrostatic Interaction Energy in Multimeric Protein-Ligand Complex: A Case Study for HIV-1 Protease

On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities

Estimating the Roles of Protonation and Electronic Polarization in Absolute Binding Affinity Simulations

The impact of interior dielectric constant and entropic change on HIV-1 complex binding free energy prediction

Rigorous Free Energy Perturbation Approach to Estimating Relative Binding Affinities Between Ligands with Multiple Protonation and Tautomeric States.

Formal Estimation of Errors in Computed Absolute Interaction Energies of Protein-Ligand Complexes

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

Proteochemometric Modeling Of The Bioactivity Spectra Of Hiv-1 Protease Inhibitors By Introducing Protein-Ligand Interaction Fingerprint

Elucidating the Inhibiting Mode of AHPBA Derivatives Against HIV-1 Protease and Building Predictive 3D-QSAR Models

Fast and accurate predictions of binding free energies using MM-PBSA and MM-GBSA

Binding Free Energy Calculation Based on the Fragment Molecular Orbital Method and Its Application in Designing Novel SHP-2 Allosteric Inhibitors

Binding Affinity and Dissociation Pathway Predictions for a Series of USP7 Inhibitors with Pyrimidinone Scaffold by Multiple Computational Methods.

Fast and Effective Prediction of the Absolute Binding Free Energies of Covalent Inhibitors of SARS-CoV-2 Main Protease and 20S Proteasome

An Extensive Test of 14 Scoring Functions Using the Pdbbind Refined Set of 800 Protein-Ligand Complexes

Development of QSAR-Improved Statistical Potential for the Structure-Based Analysis of ProteinPeptide Binding Affinities

Rapid and Reliable Binding Affinity Prediction of Bromodomain Inhibitors: A Computational Study

Computing Relative Binding Affinity of Ligands to Receptor: An Effective Hybrid Single-Dual-Topology Free-Energy Perturbation Approach in NAMD

Relative Binding Affinity Prediction of Charge-Changing Sequence Mutations with FEP in Protein–Protein Interfaces