Therapeutics targeting cytokines such as the oncostatin M (OSM)-mediated inflammation represent a potential strategy for the treatment of inflammatory bowel disease (IBD). Despite the investigation of the specific role of the interactions between OSM and the receptor (OSMR) in IBD pathogenesis, the 3D structure of the OSM-OSMR complex remains elusive. In this work, the interaction mode between OSM and OSMR at atomic level was predicted by computational simulation approach. The interaction domain of the OSMR was built with the homology modeling method. The near-native structure of the OSM-OSMR complex was obtained by docking, and long-time scale molecular dynamics (MD) simulation in an explicit solvent was further performed to sample the conformations when OSM binds to the OSMR. After getting the equilibrated states of the simulation system, per-residue energy contribution was calculated to characterize the important residues for the OSM-OSMR complex formation. Based on these important residues, eight residues (OSM: Arg100, Leu103, Phe160, and Gln161; OSMR: Tyr214, Ser223, Asp262, and Trp267) were identified as the "hot spots" through computational alanine mutagenesis analysis and verified by additional MD simulation of R100A (one of the identified "hotspots") mutant. Moreover, six cavities were detected at the OSM-OSMR interface through the FTMap analysis, and they were suggested as important binding sites. The predicted 3D structure of the OSM-OSMR complex and the identified "hot spots" constituting the core of the binding interface provide helpful information in understanding the OSM-OSMR interactions, and the detected sites serve as promising targets in designing small molecules to block the interactions.