Deep learning in GPCR drug discovery: benchmarking the path to accurate peptide binding

Deep learning (DL) methods have drastically advanced structure-based drug discovery by directly predicting protein structures from sequences. Recently, these methods have become increasingly accurate in predicting complexes formed by multiple protein chains. We evaluated these advancements to predict and accurately model the largest receptor family and its cognate peptide hormones. We benchmarked DL tools, including AlphaFold 2.3 (AF2), AlphaFold 3 (AF3), Chai-1, NeuralPLexer, RoseTTAFold-AllAtom, Peptriever, ESMFold, and D-SCRIPT, to predict interactions between G protein-coupled receptors (GPCRs) and their endogenous peptide ligands. Our results showed that structure-aware models outperformed language models in peptide binding classification, with the top-performing model achieving an area under the curve of 0.86 on a benchmark set of 124 ligands and 1240 decoys. Rescoring predicted structures on local interactions further improved the principal ligand discovery among decoy peptides, whereas DL-based approaches did not. We explored a competitive tournament approach for modeling multiple peptides simultaneously on a single GPCR, which accelerates the performance but reduces true-positive recovery. When evaluating the binding poses of 67 recent complexes, AF2 reproduced the correct binding modes in nearly all cases (94%), surpassing those of both AF3 and Chai-1. Confidence scores correlate with structural binding mode accuracy, which provides a guide for interpreting interface predictions. These results demonstrated that DL models can reliably rediscover peptide binders, aid peptide drug discovery, and guide the selection of optimal tools for GPCR-targeted therapies. To this end, we provided a practical guide for selecting the best models for specific applications and an independent benchmarking set for future model evaluation.