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.

Determination of the Protonation State of the Asp Dyad: Conventional Molecular Dynamics Versus Thermodynamic Integration.

Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures

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

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

Conformational Activation of a Transmembrane Proton Channel from Constant Ph Molecular Dynamics.

Proton-regulated electron transfers from tyrosine to tryptophan in proteins: through-bond mechanism versus long-range hopping mechanism.

A Comprehensive Study of Factors Affecting the Prediction of the p K a Shift of Asp 26 in Thioredoxin Protein

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

Application of Explicit-Solvent Constant Ph Molecular Dynamics to Proteins

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

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

Proton first: rationalizing a proton transfer in a protein-fragment complex

Selective Protonation of Catalytic Dyad for γ-Secretase-Mediated Hydrolysis Revealed by Multiscale Simulations

Constant-pH MD Simulations of the Protonation-Triggered Conformational Switching in Diphtheria Toxin Translocation Domain

Role of Resultant Dipole Moment in Mechanical Dissociation of Biological Complexes

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

Effects Of System Net Charge And Electrostatic Truncation On All-Atom Constant Ph Molecular Dynamics

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

Joint neutron crystallographic and NMR solution studies of Tyr residue ionization and hydrogen bonding: Implications for enzyme-mediated proton transfer

Computational Study of pKa shift of Aspartate residue in Thioredoxin: Role of Configurational Sampling and Solvent Model

Constant pH Simulation with FMM Electrostatics in GROMACS. (A) Design and Applications