Thrombotic thrombocytopenic purpura (TTP) is a rare and life-threatening disease. One hallmark is severe ADAMTS13 deficiency, causing ultra-large von Willebrand factor (VWF) multimers to accumulate, leading to microthrombi and lastly to microangiopathic hemolytic anemia and severe thrombocytopenia. Despite great success in recent decades, the molecular picture of the interaction between VWF and ADAMTS13 remains vague. Here, we utilized modern replica-exchange molecular dynamics simulations with the TIGER2h method to sample a vast configurational space of the isolated ADAMTS13-MDTCS domains and the exposure to its substrate and activating cofactor - the unraveled VWF-A2 domain. The sampling of binding sites and conformations was guided and filtered in agreement with available experimental evidence. We provide comprehensive information on exosites for each domain and direct pairs of interacting amino acids, for the first time. The major binding cluster for the active site of the MP domain contrasts the previous mapping of VWF-A2 residues and reciprocal binding pockets. Two major binding modes are revealed and provide access to conformational changes of an extended gatekeeper tetrad upon overcoming local latency during substrate binding and to a dedicated recruitment mechanism. Our work adds the first molecular interaction model that places previous experimental results in perspective to better understand disease-related mutations towards improved therapies. Numerous empirical targets are proposed to verify the given binding modes, to refine the overall picture of MP binding pockets, the role of Dis binding in MP activation and the passage of the Cys-rich domain through VWF-A2, thus deepening the understanding of a highly dynamic interplay.Communicated by Ramaswamy H. Sarma.
Keywords: REMD; allosteric change; molecular dynamics simulations; protein-protein interaction; structural sampling.