Structural Determinants of Buprenorphine Partial Agonism at the μ-Opioid Receptor

The μ-opioid receptor (μOR) is a class A G Protein-Coupled Receptor (GPCR) targeted by natural and synthetic ligands to provide analgesia to patients with pain of various etiologies. Available opioid medications present several unwanted side effects, stressing the need for safer pain therapies. Despite the attractive proposal that biasing μOR signaling toward G protein pathways would lead to fewer side effects, recent studies indicate that low-efficacy opioid drugs, such as buprenorphine, may represent a safer alternative. In the present work, we combine molecular docking, microsecond-time scale molecular dynamics (MD) simulations, and metadynamics to investigate the conformational dynamics of the μOR bound to morphine or buprenorphine. Our objective was to determine structural aspects associated with the unique pharmacological effects caused by the latter, taking morphine as a reference. MD simulations identified a salt bridge with D1493.32 as crucial for stabilizing both ligands into the μOR orthosteric site, with this interaction being weaker in buprenorphine. The morphinan-scaffold of both ligands shared contacts with transmembrane (TM) helix residues of the receptor, including TM3, TM5, TM6, and TM7. Conversely, while morphine showed stronger interactions with a few TM3 residues, additional chemical groups of buprenorphine showed stronger interactions with TM2, extracellular loop 2 (ECL2), and TM7 residues. We also observed distinct TM arrangements induced by these ligands, with buprenorphine causing an extracellular outward movement of TM7 and morphine provoking intracellular inward movements of TM5 and TM7 of the receptor. In addition, we found that buprenorphine tends to explore deeper regions in the μOR orthosteric site, further supported by funnel-metadynamics, resulting in diverse side chain orientations of W2956.48. Metadynamics also unveiled distinct intermediate states for morphine and buprenorphine, with the latter accessing a secondary binding site associated with partial μOR agonists. Our results indicate that the weakened salt bridge of buprenorphine with D1493.32, along with the strong TM7 interaction through its cyclopropyl group, may explain its low efficacy and consequent partial μOR agonism. Furthermore, ECL2 interactions may contribute to explaining the biased agonism of buprenorphine, a common feature shared with other opioid modulators with similar functional effects. Our study sheds light on the complex pharmacology of buprenorphine, identifying structural aspects associated with its partial and biased μOR agonism. These results can provide valuable information for the design of new effective and safer opioid drugs.