Molecular Basis of Inhibitory Mechanism of Naltrexone and Its Metabolites through Structural and Energetic Analyses

Naltrexone is a potent opioid antagonist with good blood–brain barrier permeability, targeting different endogenous opioid receptors, particularly the mu-opioid receptor (MOR). Therefore, it represents a promising candidate for drug development against drug addiction. However, the details of the molecular interactions of naltrexone and its derivatives with MOR are not fully understood, hindering ligand-based drug discovery. In the present study, taking advantage of the high-resolution X-ray crystal structure of the murine MOR (mMOR), we constructed a homology model of the human MOR (hMOR). A solvated phospholipid bilayer was built around the hMOR and submitted to microsecond (µs) molecular dynamics (MD) simulations to obtain an optimized hMOR model. Naltrexone and its derivatives were docked into the optimized hMOR model and submitted to µs MD simulations in an aqueous membrane system. The MD simulation results were submitted to the molecular mechanics–generalized Born surface area (MMGBSA) binding free energy calculations and principal component analysis. Our results revealed that naltrexone and its derivatives showed differences in protein–ligand interactions; however, they shared contacts with residues at TM2, TM3, H6, and TM7. The binding free energy and principal component analysis revealed the structural and energetic effects responsible for the higher potency of naltrexone compared to its derivatives.