Department of Chemistry, Pusan National University, Busan 609-735, Korea.
J Chem Phys. 2009 Nov 21;131(19):195102. doi: 10.1063/1.3266510.
We performed an all-atom ab initio native structure prediction of 1FME, which is one of the computationally challenging mixed fold beta beta alpha miniproteins, by combining a novel conformational search algorithm (multiplexed Q-replica exchange molecular dynamics scheme) with a well-balanced all-atom force field employing a generalized Born implicit solvation model (param99MOD5/GBSA). The nativelike structure of 1FME was identified from the lowest free energy minimum state and in excellent agreement with the NMR structure. Based on the interpretation of the free energy landscape, the structural properties as well as the folding behaviors of 1FME were compared with other beta beta alpha miniproteins (1FSD, 1PSV, and BBA5) that we have previously studied with the same force field. Our simulation showed that the 28-residue beta beta alpha miniproteins (1FME, 1FSD, and 1PSV) share a common feature of the free energy topography and exhibit the three local minimum states on each computed free energy map, but the 23-residue miniprotein (BBA5) follows a downhill folding with a single minimum state. Also, the structure and stability changes resulting from the two point mutation (Gln1-->Glu1 and Ile7-->Tyr7) of 1FSD were investigated in details for direct comparison with the experiment. The comparison shows that upon mutation, the experimentally observed turn type switch from an irregular turn (1FSD) to type I(') turn (1FME) was well reproduced with the present simulation.
我们通过将一种新颖的构象搜索算法(多路 Q 复制交换分子动力学方案)与一种平衡的全原子力场(采用广义 Born 隐式溶剂模型 (param99MOD5/GBSA))相结合,对 1FME 进行了全原子从头算天然结构预测。1FME 的天然结构从最低自由能最小状态中确定,与 NMR 结构吻合良好。基于自由能景观的解释,比较了 1FME 的结构特性和折叠行为与我们之前使用相同力场研究的其他ββα 小蛋白(1FSD、1PSV 和 BBA5)。我们的模拟表明,28 残基ββα 小蛋白(1FME、1FSD 和 1PSV)具有相似的自由能地形特征,在每个计算的自由能图谱上都表现出三个局部最小状态,但 23 残基小蛋白(BBA5)则呈单一最小状态的下坡折叠。此外,还详细研究了 1FSD 两点突变(Gln1-->Glu1 和 Ile7-->Tyr7)引起的结构和稳定性变化,以便与实验进行直接比较。比较表明,在突变后,实验观察到的从不规则环(1FSD)到 I 型(')环(1FME)的转角类型转变在本模拟中得到了很好的再现。