Institute of Biomechanics/School of Bioscience and Bioengineering, South China University of Technology, Guangzhou, China.
PLoS One. 2012;7(7):e42263. doi: 10.1371/journal.pone.0042263. Epub 2012 Jul 30.
Binding of platelet receptor glycoprotein Ibα (GPIbα) to the A1 domain of von Willebrand factor (vWF) is a critical step in both physiologic hemostasis and pathologic thrombosis, for initiating platelet adhesion to subendothelium of blood vessels at sites of vascular injury. Gain-of-function mutations in GPIbα contribute to an abnormally high-affinity binding of platelets to vWF and can lead to thrombosis, an accurate complication causing heart attack and stroke. Of various antithrombotic monoclonal antibodies (mAbs) targeting human GPIbα, 6B4 is a potent one to inhibit the interaction between GPIbα and vWF-A1 under static and flow conditions. Mapping paratope to epitope with mutagenesis experiments, a traditional route in researches of these antithrombotic mAbs, is usually expensive and time-consuming. Here, we suggested a novel computational procedure, which combines with homology modeling, rigid body docking, free and steered molecular dynamics (MD) simulations, to identify key paratope residues on 6B4 and their partners on GPIbα, with hypothesis that the stable hydrogen bonds and salt bridges are the important linkers between paratope and epitope residues. Based on a best constructed model of 6B4 bound with GPIbα, the survival ratios and rupture times of all detected hydrogen bonds and salt bridges in binding site were examined via free and steered MD simulations and regarded as indices of thermal and mechanical stabilizations of the bonds, respectively. Five principal paratope residues with their partners were predicted with their high survival ratios and/or long rupture times of involved hydrogen bonds, or with their hydrogen bond stabilization indices ranked in top 5. Exciting, the present results were in good agreement with previous mutagenesis experiment data, meaning a wide application prospect of our novel computational procedure on researches of molecular of basis of ligand-receptor interactions, various antithrombotic mAbs and other antibodies as well as theoretically design of biomolecular drugs.
血小板受体糖蛋白 Ibα(GPIbα)与血管性血友病因子(vWF)的 A1 结构域的结合是生理性止血和病理性血栓形成的关键步骤,它启动血小板在血管损伤部位与血管内皮下的黏附。GPIbα 的功能获得性突变导致血小板与 vWF 的异常高亲和力结合,并可导致血栓形成,这是导致心脏病发作和中风的一种准确并发症。在针对人 GPIbα 的各种抗血栓形成单克隆抗体(mAbs)中,6B4 是一种有效的抗体,可以在静态和流动条件下抑制 GPIbα 与 vWF-A1 的相互作用。通过突变实验进行表位作图(研究这些抗血栓形成 mAbs 的传统途径)通常既昂贵又耗时。在这里,我们提出了一种新的计算程序,该程序结合了同源建模、刚体对接、自由和导向分子动力学(MD)模拟,以确定 6B4 上的关键表位残基及其在 GPIbα 上的伙伴,假设稳定的氢键和盐桥是表位和抗原决定簇残基之间的重要连接。基于 6B4 与 GPIbα 结合的最佳构建模型,通过自由和导向 MD 模拟检查结合位点中所有检测到的氢键和盐桥的存活比例和断裂时间,并分别将其视为氢键和盐桥热和机械稳定性的指标。通过自由和导向 MD 模拟检测到的氢键和盐桥的存活比例和/或涉及的氢键的断裂时间较高的 5 个主要表位残基及其伙伴,或者它们的氢键稳定指数排名前 5,被预测为主要表位残基及其伙伴。令人兴奋的是,目前的结果与以前的突变实验数据吻合较好,这意味着我们的新计算程序在配体-受体相互作用、各种抗血栓形成 mAbs 以及其他抗体的分子基础研究以及生物分子药物的理论设计方面具有广阔的应用前景。