Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria.
Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria.
Biophys J. 2020 Aug 4;119(3):652-666. doi: 10.1016/j.bpj.2020.06.032. Epub 2020 Jul 10.
Biomolecular recognition between proteins follows complex mechanisms, the understanding of which can substantially advance drug discovery efforts. Here, we track each step of the binding process in atomistic detail with molecular dynamics simulations using trypsin and its inhibitor bovine pancreatic trypsin inhibitor (BPTI) as a model system. We use umbrella sampling to cover a range of unbinding pathways. Starting from these simulations, we subsequently seed classical simulations at different stages of the process and combine them to a Markov state model. We clearly identify three kinetically separated states (an unbound state, an encounter state, and the final complex) and describe the mechanisms that dominate the binding process. From our model, we propose the following sequence of events. The initial formation of the encounter complex is driven by long-range interactions because opposite charges in trypsin and BPTI draw them together. The encounter complex features the prealigned binding partners with binding sites still partially surrounded by solvation shells. Further approaching leads to desolvation and increases the importance of van der Waals interactions. The native binding pose is adopted by maximizing short-range interactions. Thereby side-chain rearrangements ensure optimal shape complementarity. In particular, BPTI's P1 residue adapts to the S1 pocket and prime site residues reorient to optimize interactions. After the paradigm of conformation selection, binding-competent conformations of BPTI and trypsin are already present in the apo ensembles and their probabilities increase during this proposed two-step association process. This detailed characterization of the molecular forces driving the binding process includes numerous aspects that have been discussed as central to the binding of trypsin and BPTI and protein complex formation in general. In this study, we combine all these aspects into one comprehensive model of protein recognition. We thereby contribute to enhance our general understanding of this fundamental mechanism, which is particularly critical as the development of biopharmaceuticals continuously gains significance.
蛋白质之间的生物分子识别遵循复杂的机制,深入理解这些机制可以极大地推动药物发现的进展。在这里,我们使用分子动力学模拟,以胰蛋白酶及其抑制剂牛胰蛋白酶抑制剂(BPTI)为模型系统,以原子细节跟踪结合过程的每一步。我们使用伞形采样来覆盖一系列非结合途径。从这些模拟开始,我们随后在该过程的不同阶段播种经典模拟,并将它们组合成马科夫状态模型。我们清楚地识别出三个动力学上分离的状态(未结合状态、相遇状态和最终复合物),并描述了主导结合过程的机制。从我们的模型中,我们提出了以下事件序列。最初,由于胰蛋白酶和 BPTI 中的相反电荷将它们吸引在一起,因此形成了长程相互作用,从而导致了相遇复合物的形成。该相遇复合物具有预对准的结合伙伴,其结合部位仍部分被溶剂壳包围。进一步接近导致去溶剂化,并增加范德华相互作用的重要性。通过最大化短程相互作用,采用天然结合构象。因此,侧链重排确保了最佳的形状互补性。特别是,BPTI 的 P1 残基适应 S1 口袋,并且前导位点残基重新定向以优化相互作用。在构象选择范例之后,结合竞争构象的 BPTI 和胰蛋白酶已经存在于无配体的混合物中,并且在这个建议的两步结合过程中,它们的概率增加。这种对驱动结合过程的分子力的详细描述包括了许多被认为是胰蛋白酶和 BPTI 结合以及蛋白质复合物形成的核心方面。在这项研究中,我们将所有这些方面结合到一个蛋白质识别的综合模型中。因此,我们有助于增强对这一基本机制的普遍理解,因为生物制药的不断发展具有重要意义。