Molecular dynamics simulation of bovine prothrombin fragment 1 in the presence of calcium ions

Early solvation-induced structural reorganization of calcium prothrombin fragment 1 is simulated with molecular dynamics. Initial coordinates are those of the 2.2-A resolution crystal structure [Soriano-Garcia, M., Padmanabhan, K., de Vos, A. M., & Tulinsky, A. (1992) Biochemistry 31, 2554-2556]. The molecular dynamics code AMBER, appropriately modified to include long-range (less than or equal to 22.0 A) ionic forces, was employed. The solution structure appears to equilibrate within 100 ps. Although minor changes are seen in various structural domains, the early solution structure basically maintains an intricate network of nine gamma-carboxyglutamic acid (Gla) residues encapsulating seven calcium ions. However, the Gla domain moves with respect to the kringle domain. This motion is mainly due to the movement of Ser34-Leu35 that appears to be a flexible hinge between the domains. The N-terminus of Ala 1 is in a tightly bound complex with three Gla residues that remains stable in the solution structure when the long-range electrostatic cutoff is employed and the near planar alignment of the seven calcium ions is only slightly distorted. The simulation structure is discussed in terms of experiments that studied calcium ion-induced quenching of the intrinsic fluorescence, protection of the N-terminal amino group from acetylation by calcium ions, chemical modification of the N-terminus to a trinitrophenyl derivative, and the possibility of a calcium-binding site(s) in the kringle domain.