All docking calculations were performed with the AutoDock Vina program. The crystal structures of human CMKLR1 was used as starting point. The PDB structures of 6os0.1. A was used as template for the predication of CMKLR1. All Molecular dynamics simulations (MD) simulations were performed with the pmemd module of Amber 18 software suite. Molecular mechanics generalized born surface area (MM/GBSA) was performed with mmpbsa.py program, enclosed in Amber 18 software, to estimate the binding free energy. The translational and rotational entropy was calculated with stardard statistical mechanics formulas. The vibrational entropy was computed from the frequencies of global motions by using normal mode analysis (N-mode) in mmpbsa.py program. The coordinates for these calculations are taken from the MD simulations trajectory. All docking calculations were performed with the AutoDock Vina program.