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.

Reliable and Accurate Prediction of Single-Residue Pka Values Through Free Energy Perturbation Calculations

Accurately Predicting Protein p Ka Values Using Nonequilibrium Alchemy

Cysteine and histidine pKa prediction: alchemical free energy calculations and force field reparameterization

Very fast empirical prediction and rationalization of protein pKa values

A Statistical Approach To The Prediction Of Pk(A) Values In Proteins

Accurate Prediction of GPCR Ligand Binding Affinity with Free Energy Perturbation

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

Robust prediction of relative binding energies for protein-protein complex mutations using free energy perturbation calculations

Accurate calculation of absolute protein-ligand binding free energies

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

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

Uni-pKa: An Accurate and Physically Consistent pKa Prediction through Protonation Ensemble Modeling

Accurate and Rapid Prediction of Protein pKa: Protein Language Models Reveal the Sequence-pKa Relationship

Accurate Protein Pka Prediction with Physical Organic Chemistry Guided 3D Protein Representation.

Pbsa_e: A Pbsa-Based Free Energy Estimator for Protein-Ligand Binding Affinity

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

Accurate physics-based prediction of binding affinities of RNA and DNA targeting ligands

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

An Efficient Method to Predict Protein Thermostability in Alanine Mutation

PKAD: a database of experimentally measured pKa values of ionizable groups in proteins