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.

Cryptate Binding Energies Towards High Throughput Chelator Design: Metadynamics Ensembles with Cluster-Continuum Solvation

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

Taming Multiple Binding Poses in Alchemical Binding Free Energy Prediction: the $β$-cyclodextrin Host-Guest SAMPL9 Blinded Challenge

Obtaining Consistent Free Energies for Ion Binding at Surfaces from Solution: Pathways versus Alchemy for Determining Kink Site Stability

Scalable Indirect Free Energy Method Applied to Divalent Cation-Metalloprotein Binding

Automated Fragmentation QM/MM Calculation of NMR Chemical Shifts for Protein-Ligand Complexes.

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

Alchemical Transfer Approach to Absolute Binding Free Energy Estimation

Combining Alchemical Transformation with Physical Pathway to Accurately Compute Absolute Binding Free Energy

A Fragmenting Protocol with Explicit Hydration for Calculation of Binding Enthalpies of Target-Ligand Complexes at a Quantum Mechanical Level

Performance and Analysis of the Alchemical Transfer Method for Binding Free Energy Predictions of Diverse Ligands

Relative Binding Free Energy Calculations for Ligands with Diverse Scaffolds with the Alchemical Transfer Method

Chemical accuracy for ligand-receptor binding Gibbs energies through multi-level SQM/QM calculations

Crystal Polymorph Search in the NPT Ensemble via a Deposition/Sublimation Alchemical Path

Concurrent and Adaptive Extreme Scale Binding Free Energy Calculations

Binding free energy of protein/ligand complexes calculated using dissociation Parallel Cascade Selection Molecular Dynamics and Markov state model

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

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

Physical binding energies using the electron localization function in 4-hydroxyphenylboronic acid co-crystals with aza donors

Multistate Method to Efficiently Account for Tautomerism and Protonation in Alchemical Free-Energy Calculations