Department of Physics, Applied Physics and Astronomy, Rensselaer Polytechnic Institute, Troy, NY 12180, USA.
Phys Chem Chem Phys. 2014 Apr 14;16(14):6422-9. doi: 10.1039/c3cp54336f. Epub 2014 Jan 29.
Equilibrium thermodynamics of a short beta-hairpin are studied using unbiased all-atom replica exchange molecular dynamics simulations in explicit solvent. An exploratory analysis of the free energy landscape of the system is provided in terms of various structural characteristics, for both the folded and unfolded ensembles. We find that the favorable interactions between the ends introduced by the tryptophan cap, along with the flexibility of the turn region, explain the remarkable stability of the folded state. Charging of the N termini results in effective roughening of the free energy landscape and stabilization of non-native contacts. Folding-unfolding dynamics are further discussed using a set of 2413 independent molecular dynamics simulations, 2 ns to 20 ns long, at the melting temperature of the beta-hairpin. A novel method for the construction of Markov models consisting of an iterative refinement of the discretization in reduced dimensionality is presented and used to generate a detailed kinetic network of the system. The hairpin is found to fold heterogeneously on sub-microsecond timescales, with the relative position of the tryptophan side chains driving the selection of the specific pathway.
使用无偏全原子复制交换分子动力学模拟在显式溶剂中研究了短 β-发夹的平衡热力学。根据各种结构特征,对系统的自由能景观进行了探索性分析,包括折叠和展开两种集合。我们发现,色氨酸帽引入的末端之间的有利相互作用,以及转弯区域的灵活性,解释了折叠状态的显著稳定性。N 末端的带电导致自由能景观的有效粗糙化和非天然接触的稳定化。进一步使用一组 2413 个独立的分子动力学模拟(2 ns 至 20 ns),在 β-发夹的熔化温度下,讨论了折叠-展开动力学。提出了一种新的构建马尔可夫模型的方法,该方法由在降低维数的离散化中的迭代细化组成,并用于生成系统的详细动力学网络。发夹在亚微秒时间尺度上异质折叠,色氨酸侧链的相对位置驱动特定途径的选择。