Substantial evidence in support of the formation of opioid receptor (OR) di-/oligomers suggests previously unknown mechanisms used by these proteins to exert their biological functions. In an attempt to guide experimental assessment of the identity of the minimal signaling unit for ORs, we conducted extensive coarse-grained (CG) molecular dynamics (MD) simulations of different combinations of the three major OR subtypes, i.e., xcexbc-OR, xcexb4-OR, and xcexba-OR, in an explicit lipid bilayer. Specifically, we ran multiple, independent MD simulations of each homomeric xcexbc-OR/xcexbc-OR, xcexb4-OR/xcexb4-OR, and xcexba-OR/xcexba-OR complex, as well as two of the most studied heteromeric complexes, i.e., xcexb4-OR/xcexbc-OR and xcexb4-OR/xcexba-OR, to derive the preferred supramolecular organization and dimer interfaces of ORs in a cell membrane model. These simulations yielded over 250 microseconds of accumulated data, which correspond to approximately 1 millisecond of effective simulated dynamics according to established scaling factors of the CG model we employed. Analysis of these data indicates similar preferred supramolecular organization and dimer interfaces of ORs across the different receptor subtypes, but also important differences in the kinetics of receptor association at specific dimer interfaces. We also investigated the kinetic properties of interfacial lipids, and explored their possible role in modulating the rate of receptor association and in promoting the formation of filiform aggregates, thus supporting a distinctive role of the membrane in OR oligomerization and, possibly, signaling. 