Molecular Dynamics Simulations to Investigate How PZM21 Affects the Conformational State of the μ-Opioid Receptor Upon Activation

Opioid analgesics such as morphine have indispensable roles in analgesia. However, morphine use can elicit side effects such as respiratory depression and constipation. It has been reported that G protein-biased agonists as substitutes for classic opioid agonists can alleviate (or even eliminate) these side effects. The compounds PZM21 and TRV130 could be such alternatives. Nevertheless, there are controversies regarding the efficacy and G protein-biased ability of PZM21. To demonstrate a rationale for the reduced biasing agonism of PZM21 compared with that of TRV130 at the molecular level, we undertook a long-term molecular dynamics simulation of the μ-opioid receptor (MOR) upon the binding of three ligands: morphine, TRV130, and PZM21. We found that the delayed movement of the W293(6.48) (Ballesteros–Weinstein numbering) side chain was a factor determining the dose-dependent agonism of PZM21. Differences in conformational changes of W318(7.35), Y326(7.43), and Y336(7.53) in PZM21 and TRV130 explained the observed differences in bias between these ligands. The extent of water movements across the receptor channel was correlated with analgesic effects. Taken together, these data suggest that the observed differences in conformational changes of the studied MOR–ligand complexes point to the low-potency and lower bias effects of PZM21 compared with the other two ligands, and they lay the foundation for the development of G protein-biased agonists.