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.

Effect of pH and ligand charge state on BACE-1 fragment docking performance

Prediction of the Enantioselectivity of Lipases and Esterases by Molecular Docking Method with Modified Force Field Parameters.

Binding of Inhibitors to BACE1 Affected by pH-Dependent Protonation: An Exploration from Multiple Replica Gaussian Accelerated Molecular Dynamics and MM-GBSA Calculations

Molecular mechanism underlying effect of D93 and D289 protonation states on inhibitor-BACE1 binding: exploration from multiple independent Gaussian accelerated molecular dynamics and deep learning

Protonation effects in protein-ligand complexes - a case study of endothiapepsin and pepstatin A with computational and experimental methods

D3R Grand Challenge 4: ligand similarity and MM-GBSA-based pose prediction and affinity ranking for BACE-1 inhibitors

Docking guidance with experimental ligand structural density improves docking pose prediction and virtual screening performance

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

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

Coupling enhanced sampling of the apo-receptor with template-based ligand conformers selection: Performance in pose prediction in the D3R-GC4

Docking Challenge: Protein Sampling and Molecular Docking Performance

Effects of Disulfide Bonds on Bindings of Inhibitors to BACE1 Decoded by Multiple Replica Accelerated Molecular Dynamics Simulations.

Activity prediction of substrates in NADH-dependent carbonyl reductase by docking requires catalytic constraints and charge parameterization of catalytic zinc environment

Molecular Mechanism of Binding Selectivity of Inhibitors Toward BACE1 and BACE2 Revealed by Multiple Short Molecular Dynamics Simulations and Free-Energy Predictions.

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

Predicting binding events in very flexible, allosteric, multi-domain proteins

Determining Structural Changes for Ligand Recognition between Human and Rat Phosphorylated BACE1 in Silico and Its Phosphorylation by GSK3β at Thr252 by in Vitro Studies

Effects of ligand binding on the energy landscape of acyl-CoA-binding protein

Computational Modeling of β-Secretase 1 (BACE-1) Inhibitors Using Ligand Based Approaches

Identifying Artifacts from Large Library Docking

A Combination of Computational and Experimental Approaches to Investigate the Binding Behavior of B.sub Lipase A Mutants with Substrate Pnpp