Species-specific protein thermal stability is closely correlated to the living conditions of the organism, especially to its body temperature. In this research, human and zebrafish muscle-type creatine kinases (MMCKs) were taken as model proteins to investigate the molecular adaptation of proteins in poikilothermal and homoiothermal animals. Both the optimal temperature for catalysis and the thermal stability of human MMCK was much higher than those of zebrafish MMCK. Sequence alignment identified 9 amino acid variations conserved in either the teleost MMCKs or the mammal and electric ray MMCKs. Bidirectional mutations were performed to find the residues with beneficial mutations. The results showed that two residues close to the dimer interface of MMCK, the 46th and 146th residue, were crucial for species-specific thermal stability.