Based on the up-regulation of the proviral integration site of the Moloney murine leukemia virus (Pim) kinase family (Pim1, 2, and 3) observed in several types of leukemias and lymphomas, the development of pan-Pim inhibitors is an attractive therapeutic strategy. While only PIM447 and AZD1208 have entered the clinical stages. To elucidate the interaction mechanisms of three Pim kinases with PIM447 and AZD1208, six Pim/ligand systems were studied by homology modeling, molecular docking, molecular dynamics (MD) simulation and molecular mechanics/generalized Born surface area (MM/GBSA) binding free energy calculation. The residues of the top group (Leu44, Val52, Ala65, Lys67, and Leu120 in Pim1) dominated the pan-Pim inhibitors binding to Pim kinases. The residues of the bottom group (Gln127, Asp128, and Leu174 in Pim1) were crucial for Pims/PIM447 systems, while the contributions of these residues were decreased sharply for Pims/AZD1208 systems. It is likely that the more potent pan-Pim inhibitors should be bound strongly to the top and bottom groups. The residues of the left, right and loop groups were located in the loop regions of the binding pocket, however, the flexibility of these regions triggered the protein interacting with diverse pan-Pim inhibitors efficiently. We hope this work can provide valuable information for the design of novel pan-Pim inhibitors in the future.