Great use of comparative genomics to engineer receptors! Stay tuned for our upcoming publication (fingers crossed) investing similar evolutionary hotspots for all human GPCR families.
Article: Molecular mechanisms of native ligand selectivity in catecholamine G protein-coupled receptors
#research supported by Knut and Alice Wallenberg Foundation
#biology#science@NatureComms
https://t.co/NfZx8QDnuX
Improving AlphaFold2 Performance in Virtual Screens Targeting GPCRs by Enhancing Binding-Site Conformational Sampling
1. The paper introduces AFsample2T, a targeted AlphaFold2 sampling strategy that boosts virtual screening performance for GPCRs by generating diverse binding-site conformations rather than relying on a single “best” structure.
2. Core idea: mask selected MSA columns only in the orthosteric-site region (extracellular-facing TM segments + part of EL2) to weaken local coevolution constraints and encourage alternative pocket geometries, while keeping the rest of the receptor well-constrained.
3. AFsample2T contrasts with global masking (AFsample2), which reduced binding-site accuracy for GPCRs here (AUC 0.38). Targeted masking preserved fold quality while improving pocket modeling, showing that “where to perturb” matters as much as “how much to perturb”.
4. The authors benchmarked 10 class A GPCRs using 61 curated experimental binding-site structures (from an initial pool of 119 PDB structures) and generated 1,000 models per receptor to quantify how often predicted pockets match experimental ones within 1–2 Å side-chain RMSD (symmetry-aware).
5. Moderate targeted masking improved binding-site accuracy: default AF2 AUC 0.54; adding dropout without masking AUC 0.57; targeted masking at 10–30% reached AUC ~0.59–0.61. Too much masking (50%) degraded secondary structure and collapsed performance (AUC 0.43).
6. The best-performing ensemble (AFsample2T) mixes 250 models each from 0%, 10%, 20%, and 30% masking (with dropout), yielding AUC 0.63 and capturing 73.8% of experimental binding sites at 1.5 Å RMSD vs 60.7% for AF2 (a 22% relative gain at that threshold).
7. A key mechanistic improvement is realistic binding-site plasticity. Median binding-site side-chain RMSF increased from 0.15 Å (AF2; overly rigid pockets) to 0.45 Å (AFsample2T), approaching experimental variability (median 0.58 Å). Backbone RMSF similarly moved from 0.10 Å (AF2) to 0.28 Å (AFsample2T), close to experiment (0.30 Å).
8. AFsample2T also mitigates a known docking issue: AF2 often predicts narrow/collapsed pockets. Across receptors, mean pocket volume increased (209 → 218 Å3) and the “top 1% most open” pockets expanded substantially (272 → 389 Å3), closer to experimental pockets (mean 256 Å3). This was especially relevant for μ-opioid receptor, where AF2 pockets were too collapsed.
9. Virtual screening evaluation used DOCK3.8 with ChEMBL actives (52–202 per receptor) and property-matched ZINC20 decoys, totaling extremely large-scale docking (reported as >240 trillion complexes scored). Rigid-receptor docking was used, making pocket microstates critical.
10. Ensemble screening + ligand-guided model selection is the practical win: while median enrichment of AF2-based models remained below experimental structures, the top 1% AFsample2T models improved early enrichment (mean aLogAUC top 1%: 10.8 → 12.9; mean EF1% top 1%: 7.5 → 9.6). In some targets (e.g., TAAR1, μ-opioid receptor), best AFsample2T models approached top experimental-structure performance.
11. The paper provides a workflow for prospective use: generate ≥250 AFsample2T models for the relevant receptor state, dock a ligand/decoy control set to compute enrichment, select ~top 1% models, manually inspect key interactions/poses, then proceed to large-library prospective screening.
12. Modeling receptor state is handled explicitly: inactive models use receptor sequence alone; active models are generated by cofolding receptor with heterotrimeric G protein sequences via AF2-Multimer, capturing hallmark TM6 movements and separating “state sampling” from “pocket microstate sampling”.
💻Code: https://t.co/DMW991M9ey
📜Paper: https://t.co/jKVzIdt6bx
#AlphaFold2 #GPCR #VirtualScreening #Docking #StructureBasedDrugDesign #ComputationalChemistry #Bioinformatics #ProteinStructure #DrugDiscovery #MachineLearning
Article: Improving #AlphaFold2 Performance in Virtual Screens Targeting GPCRs by Enhancing Binding-Site Conformational Sampling
#research supported by Knut and Alice Wallenberg Foundation
#science https://t.co/5MbZOapOh6
Read our just published paper on the molecular mechanisms of native ligand selectivity in G protein-coupled receptors: https://t.co/CMN2EsdI3q
@NatureComms
@UU_University
@Scilifelab
@KAW
In this paper (https://t.co/6XKoXIKdpo) we introduce AFsample2T, which increases structural diversity in AlphaFold2 models of binding sites. AFsample2T captures multiple relevant binding-site conformations and exploring diverse models can improve virtual screening performance.
Ultra-large virtual screening discovers agonists of the orphan receptor GPR139, a target for schizophrenia. AlphaFold3 fails to model receptor-ligand complexes for understudied targets! New paper: https://t.co/mmhB0Z2Ssb
@NatureComms
@UU_University
@Scilifelab
@CarlssonLab proves it: GPCR modeling only works with real collaboration—chemists, pharmacologists, modelers in sync.
🎧 Hear how it’s done: https://t.co/LvDSRh6Xu5
#DrGPCR#GPCR#TeamScience#DrugDesign
Control duration, penetration & PK/PD separation. Plus: Jens Carlsson on predictive modeling + this week’s Premium sneak peek.
Read & share ➤ https://t.co/KRQiOselhb
#DrGPCR#GPCR