Computational study of coagulation factor VIIa’s affinity for phospholipid membranes

Abstract

The interaction between the gamma-carboxyglutamic acid-rich domain of coagulation factor VIIa (FVIIa), a vitamin-K-dependent enzyme, and phospholipid membranes plays a major role in initiation of blood coagulation. However, despite a high sequence and structural similarity to the Gla domain of other vitamin-K-dependent enzymes with a high membrane affinity, its affinity for negatively charged phospholipids is poor. A few amino acid differences are responsible for this observation. Based on the X-ray structure of lysophosphatidylserine (lysoPS) bound to the Gla domain of bovine prothrombin (Prth), models of the Gla domain of wildtype FVIIa and mutated FVIIa Gla domains in complex with lysoPS were built. Molecular dynamics (MD) and steered molecular dynamics (SMD) simulations on the complexes were applied to investigate the significant difference in the binding affinity. The MD simulation approach provides a structural and dynamic support to the role of P10Q and K32E mutations in the improvement of the membrane contact. Hence, rotation of the Gly11 main chain generated during the MD simulation results in a hydrogen bond with Q10 side chain as well as the appearance of a hydrogen bond between E32 and Q10 forcing the loop harbouring Arg9 and Arg15 to shrink and thereby enhances the accessibility of the phospholipids to the calcium ions. Furthermore, the application of the SMD simulation method to dissociate C6-lysoPS from a series of Gla domain models exhibits a ranking of the rupture force that can be useful in the interpretation of the PS interaction with Gla domains. Finally, adiabatic mapping of Gla6 residue in FVIIa with or without insertion of Tyr4 confirms the critical role of the insertion on the conformation of the side chain Gla6 in FVIIa and the corresponding Gla7 in Prth.