Mamonova Tatyana, Kurnikova Maria
Chemistry Department, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA.
J Phys Chem B. 2006 Dec 14;110(49):25091-100. doi: 10.1021/jp065009n.
The ion channel protein alpha-hemolysin (alphaHL) forms supramolecular complexes with the polysaccharide beta-cyclodextrin (betaCD). This system has potential uses in nanoscale device engineering. It has been found recently that betaCD formed longer- or shorter-lived complexes with some engineered alphaHL mutants then with a wild type protein (Gu et al. J. Gen. Physiol. 2001, 118, 481-493). However, how changes in the protein sequence affect complex lifetime was not completely understood in part due to the lack of knowledge of structures of these metastable complexes. In this paper, we present an extensive molecular modeling study of the betaCD-alphaHL and selected mutant complexes to gain insights into the betaCD-alphaHL interaction mechanisms and to predict possible structures and energetics of the complexes. Thermodynamic integration (TI) and umbrella sampling (US) techniques (with the weighted histogram analysis method (WHAM)) were used to calculate the relative binding affinities of the complexes formed with the wild type alphaHL and the M113N, M113E, M113A, and M113V mutants. Our results are in excellent agreement with experiment. While betaCD-M113N and betaCD-M113A complexes were stable in the configuration of the wild type complex, the equilibrium configuration of the betaCD-M113V and betaCD-M113E complexes was significantly different. In these cases, TI alone was insufficient to accurately calculate the corresponding free energy differences. By utilizing a TI/US combination in a novel manner, we were able to accurately calculate free energy changes in these flexible systems. The betaCD-M113A and betaCD-M113E complexes, which exhibited shorter lifetimes than other complexes in an experiment, in simulations exhibited greater flexibility and higher water solvation of the betaCD adapter. MD simulations of the betaCD-M113N complex with betaCD in a downward orientation were also performed.
离子通道蛋白α-溶血素(αHL)可与多糖β-环糊精(βCD)形成超分子复合物。该系统在纳米级器件工程中具有潜在应用价值。最近发现,βCD与一些工程化的αHL突变体形成的复合物,其寿命比与野生型蛋白形成的复合物更长或更短(Gu等人,《普通生理学杂志》,2001年,第118卷,481 - 493页)。然而,由于缺乏对这些亚稳复合物结构的了解,蛋白质序列的变化如何影响复合物寿命尚不完全清楚。在本文中,我们对βCD - αHL及选定的突变体复合物进行了广泛的分子模拟研究,以深入了解βCD - αHL的相互作用机制,并预测复合物可能的结构和能量学。采用热力学积分(TI)和伞形抽样(US)技术(结合加权直方图分析方法(WHAM))来计算与野生型αHL以及M113N、M113E、M113A和M113V突变体形成的复合物的相对结合亲和力。我们得到的结果与实验结果高度吻合。虽然βCD - M113N和βCD - M113A复合物在野生型复合物的构型中是稳定的,但βCD - M113V和βCD - M113E复合物的平衡构型却显著不同。在这些情况下,仅用TI不足以准确计算相应的自由能差。通过以一种新颖的方式结合使用TI/US,我们能够准确计算这些柔性系统中的自由能变化。在实验中,βCD - M113A和βCD - M113E复合物的寿命比其他复合物短,而在模拟中,βCD适配体表现出更大的柔性和更高的水合作用。还对βCD - M113N复合物与向下取向的βCD进行了分子动力学(MD)模拟。