Abstract:The computational prediction of relative binding free energies is a crucial goal for drug discovery, and G protein-coupled receptors (GPCRs) are arguably the most important drug target class. However, they present increased complexity to model compared to soluble globular proteins. Despite breakthroughs, experimental X-ray crystal and cryo-EM structures are challenging to attain, meaning computational models of the receptor and ligand binding mode are sometimes necessary. This leads to uncertainty in understanding ligand–protein binding induced changes such as, water positioning and displacement, side chain positioning, hydrogen bond networks, and the overall structure of the hydration shell around the ligand and protein. In other words, the very elements that define structure activity relationships (SARs) and are crucial for accurate binding free energy calculations are typically more uncertain for GPCRs. In this work we use free energy perturbation (FEP) to predict the relative binding free energies for ligands of two different GPCRs. We pinpoint the key aspects for success such as the important role of key water molecules, amino acid ionization states, and the benefit of equilibration with specific ligands. Initial calculations following typical FEP setup and execution protocols delivered no correlation with experiment, but we show how results are improved in a logical and systematic way. This approach gave, in the best cases, a coefficient of determination (<i>R</i><sup>2</sup>) compared with experiment in the range of 0.6–0.9 and mean unsigned errors compared to experiment of 0.6–0.7 kcal/mol. We anticipate that our findings will be applicable to other difficult-to-model protein ligand data sets and be of wide interest for the community to continue improving FE binding energy predictions.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jcim.0c00449?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jcim.0c00449</a>.Correlation of calculated and experimental activities for additional studies and the calculated free energy differences for specific perturbation networks (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.0c00449/suppl_file/ci0c00449_si_001.pdf">PDF</a>)Input PDB and ligand structures (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.0c00449/suppl_file/ci0c00449_si_002.zip">ZIP</a>)This article has not yet been cited by other publications.

Accurate determination of protein:ligand standard binding free energies from molecular dynamics simulations

BFEE2: Automated, Streamlined, and Accurate Absolute Binding Free-Energy Calculations

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

Accurate Estimation Of Protein-Ligand Binding Free Energies Based On Geometric Restraints

Equilibrium and Non-equilibrium Ensemble Methods for Accurate, Precise and Reproducible Absolute Binding Free Energy Calculations

Calculation of Relative Binding Free Energy for Mutations in Protein Complexes: The Alchemical Path

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

The maximal and current accuracy of rigorous protein-ligand binding free energy calculations

Accurate and Efficient Calculation of Protein-Protein Binding Free Energy-Interaction Entropy with Residue Type-Specific Dielectric Constants.

An Accurate Free Energy Estimator: Based on MM/PBSA Combined with Interaction Entropy for Protein-Ligand Binding Affinity.

Equilibrium and Nonequilibrium Ensemble Methods for Accurate, Precise and Reproducible Absolute Binding Free Energy Calculations

Standard Binding Free-Energy Calculations: How Far Are We from Automation?

Accurate calculation of absolute protein-ligand binding free energies

Achieving Accurate Standard Protein-Protein Binding Free Energy Calculations through the Geometrical Route and Ergodic Sampling

Calculation of Absolute Protein-Ligand Binding Affinity Using Path and Endpoint Approaches

Accurate Prediction of GPCR Ligand Binding Affinity with Free Energy Perturbation

Computational prediction of protein–protein binding affinities

Meta-Analysis Reveals That Absolute Binding Free-Energy Calculations Approach Chemical Accuracy

Binding Free Energy Estimation for Protein-Ligand Complex Based on MM-PBSA with Various Partial Charge Models.

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

Absolute Binding Free Energy Calculation and Design of a Subnanomolar Inhibitor of Phosphodiesterase-10