Molecular dynamics simulations were performed to simulate Ca2+-dependent conformational change of calmodulin (CaM). Simulations of the fully Ca2+-bound form of CaM (Holo-CaM) and the Ca2+-free form (Apo-CaM) were performed in solution for 4 ns starting from the X-ray crystal structure of Holo-CaM. A striking difference was observed between the trajectories of Holo-CaM and Apo-CaM: the central helix remained straight in the former but became largely bent in the latter. Also, the flexibility of Apo-CaM was higher than that of Holo-CaM. The results indicated that the bound Ca2+ ions harden the structure of CaM.