Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John's, Newfoundland, Canada.
Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John's, Newfoundland, Canada.
Biophys J. 2020 Mar 24;118(6):1370-1380. doi: 10.1016/j.bpj.2020.01.020. Epub 2020 Jan 28.
Experiments have compared the folding of proteins with different amino acid sequences but the same basic structure, or fold. Results indicate that folding is robust to sequence variations for proteins with some nonlocal folds, such as all-β, whereas the folding of more local, all-α proteins typically exhibits a stronger sequence dependence. Here, we use a coarse-grained model to systematically study how variations in sequence perturb the folding energy landscapes of three model sequences with 3α, 4β + α, and β-barrel folds, respectively. These three proteins exhibit folding features in line with experiments, including expected rank order in the cooperativity of the folding transition and stability-dependent shifts in the location of the free-energy barrier to folding. Using a generalized-ensemble simulation approach, we determine the thermodynamics of around 2000 sequence variants representing all possible hydrophobic or polar single- and double-point mutations. From an analysis of the subset of stability-neutral mutations, we find that folding is perturbed in a topology-dependent manner, with the β-barrel protein being the most robust. Our analysis shows, in particular, that the magnitude of mutational perturbations of the transition state is controlled in part by the size or "width" of the underlying conformational ensemble. This result suggests that the mutational robustness of the folding of the β-barrel protein is underpinned by its conformationally restricted transition state ensemble, revealing a link between sequence and topological effects in protein folding.
实验比较了具有不同氨基酸序列但具有相同基本结构(或折叠)的蛋白质的折叠情况。结果表明,对于具有某些非局部折叠(如全-β)的蛋白质,折叠对序列变化具有鲁棒性,而更局部的全-α蛋白质的折叠通常表现出更强的序列依赖性。在这里,我们使用粗粒度模型系统地研究了序列变化如何扰动具有 3α、4β+α 和β桶折叠的三个模型序列的折叠能景观。这三种蛋白质的折叠特征与实验一致,包括折叠跃迁协同性的预期排序以及稳定性相关的自由能势垒折叠位置的变化。使用广义系综模拟方法,我们确定了约 2000 种序列变体的热力学性质,这些变体代表了所有可能的疏水或极性单点和双点突变。从对稳定性中性突变的子集的分析中,我们发现折叠以拓扑依赖性的方式受到干扰,β桶蛋白是最稳健的。我们的分析特别表明,过渡态的突变扰动程度部分受到基础构象集合的大小或“宽度”的控制。这一结果表明,β桶蛋白折叠的突变稳健性是由其构象受限的过渡态集合支撑的,揭示了序列和拓扑效应对蛋白质折叠的影响之间的联系。