Prediction of solution structures of the Ca2+-bound gamma-carboxyglutamic acid domains of protein S and homolog growth arrest specific protein 6: use of the particle mesh Ewald method

The solution structures of the N-terminal domains of protein S, a plasma vitamin K-dependent glycoprotein, and its homolog growth arrest specific protein 6 (Gas6) were predicted by molecular dynamics computer simulations. The initial structures were based on the x-ray crystallographic structure of the corresponding region of bovine prothrombin fragment 1. The subsequent molecular dynamics trajectories were calculated using the second-generation AMBER force field. The long-range electrostatic forces were evaluated by the particle mesh Ewald method. The structures that stabilized over a 400-ps time interval were compared with the corresponding region of the simulated solution structure of bovine prothrombin fragment 1. Structural properties of the gamma-carboxyglutamic acid (Gla) domains obtained from simulations and calcium binding were found to be conserved for all three proteins. Analysis of the predicted solution structure of the Gla domain of Gas6 suggests that this domain should bind with negatively charged phospholipid surfaces analogous to bovine prothrombin fragment 1 and protein S.