Department of Molecular Biology and Genetics, Democritus University of Thrace, University campus, 68100 Alexandroupolis, Greece.
J Chem Inf Model. 2024 Aug 12;64(15):6081-6091. doi: 10.1021/acs.jcim.4c00575. Epub 2024 Jul 25.
The effect of mutations on protein structures is usually rather localized and minor. Finding a mutation that can single-handedly change the fold and/or topology of a protein structure is a rare exception. The A31P mutant of the homodimeric Repressor of primer (Rop) protein is one such exception: This single mutation ─and as demonstrated by two independent crystal structure determinations─ can convert the canonical (left-handed/all-antiparallel) 4-α-helical bundle of Rop to a new form (right-handed/mixed parallel and antiparallel bundle) displaying a previously unobserved "bisecting U" topology. The main problem with understanding the dramatic effect of this mutation on the folding of Rop is to understand its very existence: Most computational methods appear to agree that the mutation should have had no appreciable effect, with the majority of energy minimization methods and protein structure prediction protocols indicating that this mutation is fully consistent with the native Rop structure, requiring only a local and minor change at the mutation site. Here we use two long (10 μs each) molecular dynamics simulations to compare the stability and dynamics of the native Rop versus a hypothetical structure that is identical with the native Rop but is carrying this single Alanine to Proline mutation. Comparative analysis of the two trajectories convincingly shows that, in contrast to the indications from energy minimization ─but in agreement with the experimental data─, this hypothetical native-like A31P structure is unstable, with its turn regions almost completely unfolding, even under the relatively mild 320 K simulations that we have used for this study. We discuss the implication of these findings for the folding of the A31P mutant, especially with respect to the proposed model of a double-funneled energy landscape.
突变对蛋白质结构的影响通常是局部的和微小的。找到一种能够单独改变蛋白质结构折叠和/或拓扑结构的突变是一个罕见的例外。二聚体引物抑制剂(Rop)蛋白的 A31P 突变体就是这样一个例外:这种单一突变——正如两个独立的晶体结构测定所证明的那样——可以将 Rop 的典型(左手/全反平行)4-α-螺旋束转换为一种新形式(右手/混合平行和反平行束),呈现出以前未观察到的“二分 U”拓扑结构。理解这种突变对 Rop 折叠的巨大影响的主要问题是理解它的存在:大多数计算方法似乎都认为这种突变应该没有明显的影响,大多数能量最小化方法和蛋白质结构预测协议都表明这种突变与天然 Rop 结构完全一致,只需要在突变部位进行局部和微小的改变。在这里,我们使用两个长(每个 10 μs)分子动力学模拟来比较天然 Rop 的稳定性和动力学与一个假设的结构,该结构与天然 Rop 相同,但带有这个单一的丙氨酸到脯氨酸突变。对这两个轨迹的比较分析令人信服地表明,与能量最小化的指示相反——但与实验数据一致——这个假设的天然样 A31P 结构是不稳定的,其转弯区域几乎完全展开,即使在我们用于这项研究的相对温和的 320 K 模拟下也是如此。我们讨论了这些发现对 A31P 突变体折叠的影响,特别是对提出的双漏斗能景观模型的影响。