Analysis of structures and sequence similarities of the vitamin K-dependent (VKD) proteins in the blood coagulation pathway has shown that these proteins have evolved through a series of gene duplications and diversification, acquiring a degree of functional diversity in the process [1][2][3]. Protein Z (PZ) is a VKD plasma glycoprotein that is highly conserved across different species. It is homologous to the blood coagulation factors VII (FVII), FIX, FX, and protein C (PC). However, PZ differs from other coagulation proteins in that it lacks the critical histidine and serine residues of the catalytic triad, and is therefore not a zymogen of a serine protease (SP) [4,5].Human PZ was first isolated in 1984 [6], and its gene was characterized in 1998 [7]. PZ is relatively abundant in humans, with a wide plasma concentration range, and circulates as a complex with PZ-dependent protease inhibitor (ZPI) [8][9][10]. It has been shown to function as a cofactor in the inhibition of activated FX (FXa) by ZPI, causing a thousandfold increase in inhibition in a Ca 2+ -dependent manner [11]. The relative abundance in plasma and sequence similarity with other VKD proteins serve to pique our interest in the PZ system.In order to gain a better understanding of the function and evolutionary significance of PZ, it is important to start from structural information at an atomic level. Although PZ was isolated over 20 years ago, its physiologic role in relation to its molecular structure is still not clear. In this study, we propose a solvent-equilibrated structural model of human PZ using homology modeling and molecular dynamics (MD) simulation. The possibility of using this approach is supported by recent findings [12].Human PZ is a single-chain molecule of 360 residues with four domains: a Gla domain (residues 1-46) with 13 ccarboxyglutamic acid (Gla) residues (residues 7, 8, 11, 15, 17, 20, 21, 26, 27, 30, 33, 35 and 40), two epidermal growth factor (EGF)-like domains -an EGF1 domain (residues 47-83) with a b-hydroxyaspartate (Bha) residue at 64, and an EGF2 domain (residues 85-126) -and an SP-like domain (residues 135-360) with high homology to the catalytic domain of other VKD SPs [4,5].We have employed an iterative process of homology modeling and structure evaluation followed by MD simulation to construct a solvent-equilibrated structural model of PZ. Two different schemes were adopted to construct the homology model for PZ: one using multiple sequence alignment of human FVIIa [13] 18-23, 51-62, 56-71, 73-82, 89-101, 97-110, 112-125, 163-179 and 287-301). To obtain a refined solvent-equilibrated model and to remove bad contacts in the homology model, we performed an MD simulation using PMEMD9 in the AMBER9 [20] suite with the ff99SB forcefield. The total system was composed of 99 459 atoms, including 12 calcium ions, four sodium counterions for neutralizing the system, and 31 307 TIP3P water molecules. The seven conserved calcium ions in the Gla domain were placed on the basis of FVIIa, and the additional calcium ions w...