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.

AlphaFold3 versus experimental structures: assessment of the accuracy in ligand-bound G protein-coupled receptors

AlphaFold2 versus experimental structures: evaluation on G protein-coupled receptors

Multi‐State Modeling of G‐protein Coupled Receptors at Experimental Accuracy

AlphaFold prediction and analysis of Adhesion-family G protein coupled receptor (aGPCR) structures.

A Focus on Selected Class A GPCRs

Evaluation of AlphaFold2 structures as docking targets

How accurately can one predict drug binding modes using AlphaFold models?

Conservation of Hot Spots and Ligand Binding Sites in Protein Models by AlphaFold2

AlphaFold2 structures guide prospective ligand discovery

Accurate structure prediction of biomolecular interactions with AlphaFold 3

Comparative analysis of RNA 3D structure prediction methods: towards enhanced modeling of RNA–ligand interactions

Reliability of Docking-Based Virtual Screening for GPCR Ligands with Homology Modeled Structures: A Case Study of the Angiotensin II Type I Receptor.

Accurate Prediction of GPCR Ligand Binding Affinity with Free Energy Perturbation

Assessment of AI-Based Protein Structure Prediction for the NLRP3 Target

Benchmarking AlphaFold for protein complex modeling reveals accuracy determinants

Evaluation of AlphaFold2 Structures for Hit Identification across Multiple Scenarios

AlphaFold2 structures template ligand discovery

Revolutionizing GPCR–ligand predictions: DeepGPCR with experimental validation for high-precision drug discovery

Revolutionizing GPCR-Ligand Predictions: DeepGPCR with experimental Validation for High-Precision Drug Discovery

Comparative studies of AlphaFold, RoseTTAFold and Modeller: a case study involving the use of G-protein-coupled receptors

A Benchmarking Study on Virtual Ligand Screening Against Homology Models of Human GPCRs