Understanding and Predicting Ligand Efficacy in the μ-Opioid Receptor through Quantitative Dynamical Analysis of Complex Structures

The μ-opioid receptor (MOR) is a G-protein coupled receptor involved in nociception and the primary target of opioid drugs. Understanding the relationships among the ligand structure, receptor dynamics, and efficacy in activating MOR is crucial for drug discovery and development. Here, we use coarse-grained normal-mode analysis to predict ligand-induced changes in receptor dynamics with the Quantitative Dynamics Activity Relationship (QDAR) DynaSig-ML methodology, training a LASSO regression model on the entropic signatures (ESs) computed from ligand-receptor complexes. We train and validate the methodology using a data set of 179 MOR ligands with experimentally measured efficacies split into strictly chemically different cross-validation sets. By analyzing the coefficients of the ES LASSO model, we identified key residues involved in MOR activation, several of which have mutational data supporting their role in MOR activation. Additionally, we explored a contact-only LASSO model based on ligand-protein interactions. While the model showed predictive power, it failed at predicting efficacy for ligands with low structural similarity to the training set, emphasizing the importance of receptor dynamics for predicting ligand-induced receptor activation. Moreover, the low computational cost of our approach, at 3 CPU s per ligand-receptor complex, opens the door to its application in large-scale virtual screening contexts. Our work contributes to a better understanding of dynamics-function relationships in the μ-opioid receptor and provides a framework for predicting ligand efficacy based on ligand-induced changes in receptor dynamics.