Abstract

Methadone is a synthetic opioid agonist with notoriously unique properties, such as lower abuse liability and induced relief of withdrawal symptoms and drug cravings, despite acting on the same opioid receptors triggered by classical opioids, in particular the µ-opioid receptor (MOR). Its distinct pharmacologic properties, which have recently been attributed to the preferential activation of β-arrestin over G proteins, make methadone a standard-of-care maintenance medication for opioid addiction. Although a recent biophysical study suggests that methadone stabilizes different MOR active conformations from those stabilized by classical opioid drugs or G protein-biased agonists, how this drug modulates the conformational equilibrium of MOR and what specific active conformation of the receptor it stabilizes is unknown. Here, we report the results of sub-millisecond adaptive sampling molecular dynamics (MD) simulations of a predicted methadone-bound MOR complex, and compare them with analogous data obtained for the classical opioid morphine and the G protein-biased ligand TRV130. The model, which is supported by existing experimental data, is analyzed using Markov State Models (MSMs) and transfer entropy analysis to provide testable hypotheses of methadone-specific conformational dynamics and activation kinetics of MOR.