Prediction of binding characteristics between von Willebrand factor and platelet glycoprotein Ibα with various mutations by molecular dynamic simulation

Abstract

IntroductionBinding of platelet glycoprotein (GP)Ibα with von-Willebrand factor (VWF) exclusively mediates the initial platelet adhesion to injured vessel wall. To understand the mechanism of biomedical functions, we calculated the dynamic fluctuating three-dimensional (3D) structures and dissociation energy for GPIbα with various single amino-acid substitution at G233, which location is known to cause significant changes in platelet adhesive characteristics. Material and methodsMolecular dynamics (MD) simulation was utilized to calculate 3D structures and Potential of Mean Force (PMF) for wild-type VWF bound with wild-type, G233A (equal function), G233V (gain of function), and G233D (loss of function) GPIbα. Simulation was done on water-soluble condition with time-step of 2 × 10−15 s using NAnoscale Molecular Dynamics (NAMD) with Chemistry at HARvard Molecular Mechanics (CHARMM) force field. Initial structure for each mutant was obtained by inducing single amino-acid substitution to the stable water-soluble binding structure of wild-type. ResultsThe most stable structures of wild-type VWF bound to GPIbα in wild-type or any mutant did not differ. However, bond dissociation energy defined as difference of PMF between most stable structure and the structure at 65 Å mass center distances in G233D was 4.32 kcal/mol (19.5%) lower than that of wild-type. Approximately, 2.07 kcal/mol energy was required to dissociate VWF from GPIbα with G233V at mass center distance from 48 to 52 Å, which may explain the apparent “gain of function” in G233V. ConclusionThe mechanism of substantially different biochemical characteristics of GPIbα with mutations in G233 location was predicted from physical movement of atoms constructing these proteins.