Naltrexone (NTX) 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 NTX 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. NTX 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 Molecular Mechanics Generalized-Born surface area (MMGBSA) binding free energy calculations and principal component analysis. Our results revealed that NTX and its derivatives showed differences in protein-ligand interactions; however, they shared contact 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 NTX compared to its derivatives.