Prothrombin is a key zymogen of the coagulation process and can be converted to thrombin by the prothrombinase complex, which consists of factor Xa (FXa), cofactor Va (FVa), and phospholipids. Prothrombin can be activated at two cleavage sites, R271 and R320, which generates two intermediates: prethrombin-2 via the initial cleavage at R271, and meizothrombin via the first cleavage at R320. Several mechanisms have been proposed to explain this activation preference, but the role of cleavage site sequences in prothrombin activation has not been thoroughly investigated. Here, we used an advanced sampling technique, parallel tempering metadynamics with a well-tempered ensemble (PTMetaDWTE), to study the binding modes of prothrombin cleavage site sequences R266AIEGRTATSEY277 (denoted as Pep271) and S315YIDGRIVEGSD326 (denoted as Pep320) to the FXa catalytic triad. Our study indicates that there exist three binding modes for Pep271 to the FXa catalytic triad but only one binding mode for Pep320 to the FXa catalytic triad. Further molecular dynamics simulations revealed that due to the strong electrostatic interactions, especially the H-bond interactions and salt bridges formed between Pep320 and FXa, the binding mode in the Pep320-FXa system is more stable than the binding modes in the Pep271-FXa system. In view of experimental observations and our results that there exists only one binding mode for Pep320 to the FXa catalytic triad and especially R320 in Pep320 can stably bind to the FXa catalytic triad, we believe that the first cleavage at R320 is favored.